SAR点目标仿真

clear all;
close all;

%const
c = 3e8;
pi = 3.1415926;
j = sqrt（-1）;

Rn = 20e3;  %点据
Vr = 150;     %速度
Tp = 10e-6;   %脉冲宽度
fc = 5.3e9;      %载频
Kr = 20e12;    %距离向调频率
Bw = Kr*Tp;    %距离向带宽
Fs = 1.2*Bw;   %采样率
Ba = 200;         %多普勒带宽
PRF = 1.2*Ba;  %脉冲重复频率
lamda = c/fc;      %波长

sin_theta = Ba/2*lamda/2/Vr;                %最大斜视角
cos_theta = sqrt（1-sin_theta^2）;
tan_theta = sin_theta/cos_theta;
La = Rn*tan_theta*2;              %合成孔径
Na = fix（La/Vr*PRF）+1;      %方位向采样数
Nr = fix（Tp*Fs）;                   %距离向采样数

%回波数据
data = zeros（Na2*Nr）;
for k = 1:Na;
    L= abs（k-Na/2）/PRF*Vr;
    slantRange = sqrt（L^2+Rn^2）;
    timeOffset = 2*（slantRange-Rn）/c;
    offsetNum = round（timeOffset*Fs）;
    Ni = （0+fix（Nr/2））:（Nr+fix（Nr/2））;
    tr = （Ni-Nr+offsetNum）/Fs;                   %发射信号时间向量
    sig = exp（-j*pi*Kr*（tr-timeOffset）.^2）*exp（-j*4*pi/lamda*slantRange）;
    data（kNi+offsetNum）=sig;
end
figureimagesc（abs（data））;        title（‘点目标回波幅度‘）
figureimagesc（（angle（data）））;   title（‘点目标回波相位‘）
figureplot（real（data（Na/2:）））;

%距离向压缩
f_axis = [-Nr:Nr-1]/2/Nr*Fs;
HfRange = exp（-j*pi/Kr*f_axis.^2）;
% HfRange = fftshift（HfRange）;
figureplot（f_axisabs（HfRange））;
figureplot（abs（myfft（data（Na/2:））））;
for i = 1:Na
    data（i:） = myifft（myfft（data（i:））.*HfRange）;
end
figureimagesc（abs（data））;         title（‘距离压缩数据‘）
figureplot（abs（