?? baseline_200702.m
字號:
% 本程序用于計算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');
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -