SAR单散射点条带式回波数据仿真程序（带注释）

       场景设置：

航空器方位移动范围：-1500m-1500m

散射点坐标（0,11180m,0）

线性调频波形采样率：100MHz

载频：10GHz

中频中心频率：6MHz

带宽：10MHz

时宽：5us

（1）SAR单散射点条带式回波数据仿真（时间-方位）

（2）SAR单散射点条带式回波数据仿真（距离-方位） 

clear all
close all
clc
fs=100e6;%线性调频波形中频采样率
f0=6e6;%线性调频波形中频波形中心频率
B=10e6;%线性调频波形带宽
t=5e-6;%线性调频波形时宽
miu=B/(2*t);%线性调频波形调频斜率
Nt=round(fs*t);
tt=0:1/fs:t-1/fs;

c=300e6;%电磁波传播速度
fc=10e9;%线性调频波形载频
lamda=c/fc;%线性调频波形波长
Rp=11.18e3;%散射点在xoy轴上的y坐标位置，x坐标位置为0
u=-1500:1:1500;%航空器延x轴方向的运动轨迹
umax=1500;%航空器在xoy轴上的投影与原点的最大距离
umin=0;%航空器在xoy轴上的投影与原点的最小距离
R=(Rp+u.^2/(2*Rp));%航空器在每个位置处与散射点的近似直线距离
td=2*R./c;%航空器在每个位置处对应的回波时延
tdmax=2*(Rp+umax.^2/(2*Rp))/c+t;%最大回波时延对应的回波信号最大边界
tdmin=2*(Rp+umin.^2/(2*Rp))/c-1e-6;%最小回波时延对应的回波信号起始最小边界（-1e-6是扩宽最小边界，提高可视性）
Nsmax=round(fs*tdmax);
Nsmin=round(fs*tdmin);
Ns=round(fs.*td);

tt_d=tdmin:1/fs:tdmax-1/fs;

lfm_td_ind=[];
lfm_td=[];
data=zeros(length(u),Nsmax-Nsmin);%回波数据矩阵
for k=1:length(u);
    lfm_td_ind(k,:)=exp(1i*(2*pi*(tt-td(k)).*(f0-B/2)+pi*miu.*(tt-td(k)).^2));%每一个方位上发射一个线性调频波形
    lamda_x(k)=exp(-1i*4*pi*R(k)/lamda);%每个方位上的回波对应的相位调制
    lfm_td(k,:)=lfm_td_ind(k,:).*lamda_x(k);
    data(k,(Ns(k)-Nsmin):(Ns(k)+Nt-1-Nsmin))=lfm_td(k,:);%每个方位上的回波对应的时延调制
end

%% 单散射点回波数据（时间-方位）
figure
imagesc(tt_d.*1e6,u,real(data))
xlabel('时延（us）')
ylabel('方位（m）')

%% 单散射点回波数据（距离-方位）
figure
imagesc(tt_d.*1e6.*c/2/1e3,u,real(data))
xlabel('距离（km）')
ylabel('方位（m）')