SAR雷达回波信号仿真

clear;
clc;
close all;

%% 基本参数
c=3e8;                                      %光速
fstart=4300e6;                              %起始频率
fstop=4600e6;                               %终止频率
fc=(fstart+fstop)/2;                        %载频
lambda=c/fc;                                %波长
Br=fstop-fstart;                            %发射信号带宽
Tr=20e-3;                                   %信号持续时间
K=Br/Tr;                                    %调频率
s=16.8;                                       %雷达平台到场景平面的垂直距离
theta_L=0/180*pi;                          %波束俯仰角%%%%%%%
Rs=s/cos(theta_L);                          %景中心斜距
Yc=sqrt(Rs^2-s^2);                          %场景中心位置（距离方向）
Y0=5;                                       %场景宽度一半
Ymax=Yc+Y0;                                 %距离方向上的最大值
Ymin=Yc-Y0;                                 %距离方向上的最小值%%%%
Xmin=-5;                                    %方位向成像最小距离
Xmax=5;                                     %方位向成像最大距离
D=0.1;                                      %天线尺寸
Lsar=0.886*Rs*lambda/D;                     %合成孔径长度
l=2.4;                                      %%%%%
Rmin=sqrt(Ymin^2+s^2);                      %最短斜距
Rmax=sqrt(Ymax^2+s^2+(Lsar)^2);             %最长斜距
DY=c/2/Br;                                  %距离向理论分辨率
DX=D/2*Lsar/l;                                     %方位向理论分辨率

%% 点目标信息
PY=[3 3 3 0 0 0 -3 -3 -3]+Yc;               %以场景中心为中心的多目标
PX=[-3 0 3 -3 0 3 -3 0 3];
Ntarget=size(PX,2);

%% 计算相关数据
dl=0.04;                                    %雷达步进位移
N=(Lsar+Xmax-Xmin)/dl;                      %雷达步进点数
N=ceil(N);                                  %步进点数取整
Xa=linspace((Xmin-Lsar/2),(Xmax+Lsar/2),N); %雷达位置序列
tao_max=2*Rmax/c;                           %回波信号最大延时
B=K*tao_max;                                %差频信号带宽
Fsr=1.2*B;                                  %差频信号采样率
M=ceil(Fsr*Tr);                             %一个快时间内的采样点数
M=2^nextpow2(M);                            %在正无穷方向上取最大的2次幂数字
dt=Tr/M;                                    %采样间隔
Fsr=1/dt;                                   %更新采样率
t=ones(N,1)*linspace(-Tr/2,Tr/2,M);         %快时间序列
Kr=linspace(4*pi/c*fstart,4*pi/c*fstop,M);  %Kr域序列

%% 回波信号
sig=zeros(N,M);
for k=1:Ntarget
  R=sqrt((PY(k)+s)^2+(Xa'-PX(k)).^2)*ones(1,M);   %距离矩阵
    tao=2*R/c;              
     phase=fc*tao+K*tao.*t;
      Xlimt=(Xa')*ones(1,M)
   sig=sig+exp(-1j*2*pi*phase).*(abs(Xlimt)<(l/2));
end
figure('Name','回波信号幅度'); imagesc(abs(sig));
title('回波信号幅度'); xlabel('距离向（采样点）'); ylabel('方位向（采样点）');
figure('Name','回波信号相位'); imagesc(abs(angle(sig)));
title('回波信号相位'); xlabel('距离向（采样点）'); ylabel('方位向（采样点）');