单站模式InSAR成像流程仿真

%% *************************************
% 2013.11.01 编写__单站模式
% 生成目标回波
% 2015.10.09 完成单条方位线上目标的高程提取， 在没有高程模糊的情况下，简化解缠绕流程
% 简单结论有：采用斜平面成像，需要去地平面效应；采用地平面成像，直接得到场景的高程起伏
% 2015.10.24 验证了干涉成像所有流程；然而在相位解缠绕方面稳健性不佳
%            采用更快的相位补偿算法+线性插值 提高成像质量 
% 2015.11.05 增加了对遮挡因素影响的仿真
%% *************************************
%% 校正目标位置坐标
clear all;
close all;
clc;
TrgtNum=1;              %目标数量
TrgtPos=[0*ones(TrgtNum,1),15+(1:TrgtNum)'*0.2,sin((0:TrgtNum-1)/TrgtNum*pi*2)'*0.3]; 
TrgtRCS=[];             %目标后向散射强度
FLAG=0;
%% XXXXXXXXXXXXXXXXX系统参数设置XXXXXXXXXXXXXXXXXXX
AntennaHeight=2;                  %发射天线高度
HeightDiff=0.25;                     %发射天线与接收天线高度差
AntIntvl=0.01*4;                     %天线间距
f0=13e9;                            %发射信号起始频率
B=2.00e9;                            %系统带宽
Delt_f=2e6;                          %频率步进增量
ReAntCord=-0.4:AntIntvl:0.4;        %接收天线水平位置坐标
RxNum=length(ReAntCord);             %接收天线数81
ReAntHght1=2.2-0.35/2;
ReAntHght2=2.2+0.35/2;
BaselineLen=ReAntHght2-ReAntHght1;
%% XXXXXXXXXXXXXXXX脉冲压缩相关参数XXXXXXXXXXXXXXXXX
RC_Filter=[zeros(50,1);hamming(925);zeros(25,1)];    % 距离脉压（RC）窗函数
%% ***************成像参数****************
c=3e8;                               %光速
FreqNum=B/Delt_f;                    %频率步进数
delt_r=c/Delt_f/2;
Rou_r=delt_r/FreqNum;                 %系统距离分辨率
fs=8192*1;                             %采样频率
ImgRngRes=0.0375;%2*Rou_r/3;                  %距离向成像分辨率
ImgAzmRes=0.0375;%2*Rou_r/3;                  %方位向成像分辨率
ImageWide=5;                          %成像区域宽
ImageRangeN=10;                       %成像最近距离 (距离大小，对干涉性能也有影响)
ImageRangeF=15.95;                    %成像最远距离
f_vector=(0:FreqNum-1)'*Delt_f+f0;    %频率矢量
Ix=-(ImageWide-ImgAzmRes)/2:ImgAzmRes:(ImageWide-ImgAzmRes)/2;%图像方位向坐标
Iy=(ImageRangeN:ImgRngRes:ImageRangeF)';                         %图像距离向坐标
IxMesh=repmat(Ix,length(Iy),1);                            IyMesh=repmat(Iy,1,length(Ix));
fc=(f0+B/2);Lamda=3e8/fc; 
NewPhaseCompVec=exp(1i*4*pi*(f0-fc)*delt_r/fs/c*[0:fs-1].');%效率更高
%% **********生成DEM和RCS******
%如果网格单元小于分辨单元，根据相干斑形成原理，图像中将存在严重的相干斑
%如果目标相位相同，根据散射原理，干涉之后，点线阵退化为两个点
%20111208看北理工论文，三种解决途径：加随机高度、加随机相位和加随机位置扰动
%看来三种扰动效果都不好。现采用的解决途径为：降低场景采样率，增大网格间距
Tx=Ix(1):Rou_r*0.8:Ix(end);%大于0.0375
Ty=(Iy(1):Rou_r*1.5:Iy(end))';
TxMesh=repmat(Tx,[length(Ty),1]);
TyMesh=repmat(Ty,[1,length(Tx)]);
DEM=zeros([length(Ty),length(Tx)]);
if(length(Tx)>=length(Ty))
    Z = peaks(length(Tx));
    %Z = cone(length(Tx),length(Tx));


    %%Z_y=Ty(1):(Ty(end)-Ty(1))/(length(Tx)-1):Ty(end);
    Z_y=linspace(Ty(1),Ty(end),length(Tx));
    for kk=1:length(Tx)
        DEM(:,kk)=interp1(Z_y,Z(:,kk),Ty);
    end
    %%在83*83中通过插值取53*83%%
else
    Z = peaks(length(Ty));
    %Z = cone(length(Ty),length(Ty));
    Z_x=Tx(1):(Tx(end)-Tx(1))/(length(Ty)-1):Tx(end);
    for kk=1:length(Ty)
        DEM(kk,:)=interp1(Z_x,Z(kk,:),Tx);
    end
end
DEM=DEM/max(DEM(:))*0.3;
%figure;mesh(DEM);
%关键步骤，增加位置扰动
TxMesh=TxMesh+(rand(size(TxMesh))-0.5)*Rou_r*0.2;