% 本程序用于計算InSAR中各種誤差元對于高度誤差的影響曲線
% 請先輸入已知數據

% 平臺相對高度（推薦500，單位km）
H_zero=499.226;

% 投射角（單位degree）
theta_temp=35;
theta=theta_temp*pi/180;

% 基線長度（單位m）
baseline_length=3000;

% 基線角（單位degree）
alpha_temp=30;
alpha=alpha_temp*pi/180;

% 工作頻率（單位GHz）
fre=3.2;
lamda=0.3/fre; % unit:m

r=H_zero/cos(theta);

%====================定義變量（誤差元、橫坐標）====================
del_B_length=zeros(1001,1); % unit m
del_h_B_length=zeros(1001,1);

del_alpha=zeros(1001,1);
del_h_alpha=zeros(1001,1);

del_r=zeros(1001,1); %(unit m)
del_h_r=zeros(1001,1);

del_phi=zeros(1001,1);
del_h_phi=zeros(1001,1);

for i=1:1001
    del_B_length(i)=(i-1)/10000;
    del_alpha(i)=(i-1)/500000;
    del_r(i)=(i-1)/50;
    del_phi(i)=(i-1)/5;
end

%====================基線長度測量誤差對于測高誤差曲線(unit m)====================
for i=1:1001
    del_h_B_length(i)=1000*r*tan(theta-alpha)*sin(theta)*del_B_length(i)/baseline_length;
end
figure(1)
plot(del_B_length,del_h_B_length);
xlabel('m'),ylabel('m');

%====================基線角測量誤差對于測高誤差曲線(unit m)====================
for i=1:1001
    del_h_alpha(i)=1000*r*sin(theta)*del_alpha(i)*pi/180;
end
figure(2)
plot(del_alpha,del_h_alpha);
xlabel('degree'),ylabel('m');

%====================距離測量誤差對于測高誤差曲線（unit m）====================
for i=1:1001
    del_h_r(i)=cos(theta)*del_r(i);
end
figure(3)
plot(del_r,del_h_r);
xlabel('m'),ylabel('m');

%====================相位測量誤差對于測高誤差曲線(unit m)====================
for i=1:1001
    del_h_phi(i)=1000*lamda*r*sin(theta)*del_phi(i)*pi/180/4/pi/baseline_length/cos(theta-alpha);
end
figure(4)
plot(del_phi,del_h_phi);
xlabel('degree'),ylabel('m');