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');