点目标飞机模型iSAR成像

强力推荐的ISAR程序代码，自己调试过，可以运行，对学习雷达成像的朋友很有帮助！

clc;
clear all;
%=============基本参数设置=======================
c=3.0e+8;                   %光速
r0=100;                       %雷达到目标旋转中心的距离为100m 
w=0.125;                     %目标旋转角速度，单位（rad）                      
Mr=128;                     %距离向采样点数
Ma=128;                     %方位向采样点数
%=============基本参数设置=======================

%=============调频信号参数设置====================
f0=10e+9;                    %载频10GHz
lamda=c/f0;                  %波长
pulsewide=25.6e-6;     %发射信号的脉冲宽度为25.6微秒
B=400e+6;                   %发射信号的带宽为400MHz
PRF=400;                      % 脉冲重复频率 400Hz
PRT=1/PRF;                  % 脉冲重复周期为2.5毫秒
K=B/pulsewide ;           % 调频斜率
%=============调频信号参数设置====================

%=====================Load Target ===============
load Target_Plane_SeventyFour_Point.mat posxy
x=posxy（:1）;
y=posxy（:2）;
[M_Target_Plane_SeventyFour_PointN_Target_Plane_SeventyFour_Point]=size（posxy）;
%===================== Load Target ===============
%=====================Load Target ===============
% % load Target_Plane_SeventyFour_Point.mat posxy
% x=[-9 -8 -7 -7 -6 -6 -5 -5 -4 -4 -4 -3 -3 -3 -2 -2 -2 -1 -1 -1 ... 
%     0 0 0 0 0 0 0 0 0 0 1 1 2 2 3 3 3 3 3];
% y=[0 0 0.5 0.5 1 1 1.5 0 -1.5 2 0 -2 2.5 0 -2.5 3 0 -3 3.5 ...
%     0-3.5 4 3.6 2.8 2 1 -1 -2 -2.8 -3.6 -4 0.5 -0.5 1 -1 2 1 0 -1 -2];
% [M_Target_Plane_SeventyFour_PointN_Target_Plane_SeventyFour_Point]=size（x）;
% %===================== Load Target ===============

sigma=5*ones（1M_Target_Plane_SeventyFour_Point）;              %?为什么全是5？为什么不是目标存在点的σ为1，其他为0？
j=sqrt（-1）; 
t=-pulsewide/2:pulsewide/（Ma-1）:pulsewide/2;
Hamming_Windows=hamming（Mr）;
F=4.0*pi*（f0+K*t）/c;

%=====================开辟矩阵便于保存数据 ==========
echos=[MaMr]; 
echos_real=[MaMr];
echos_final=[MaMr];
%=====================开辟矩阵便于保存数据 ==========

%=====================仿真回波数据 =================成像区域化为Ma*Mr
for m=1:Ma
    time=（m-Ma/2）*PRT;     %慢时间
    for n=1:Mr
        s1=0.0;
        for k=1:M_Target_Plane_SeventyFour_Point
            s2=sigma（k）*exp（-j*（（F（n）*（x（k）*cos（w*time）+y（k）*sin（w*time）））））;  %此处为回波
            s1=s1+s2;
        end
        s1=exp（j*F（n））*s1;  %匹配滤波
        s1=s1*Hamming_Windows（n）;  % 对行（距离向）加海明窗 将距离向的波形控制在n个采样点内，
        echos（mn）=s1;
    end
end
%=====================仿真回波数据 =================

%=====================输出显示回波数据 ==================
figure（1）imagesc（abs（echos））;colormap（1-gray）;
xlabel（‘距离向采样点（个）‘）ylabel（‘回波帧数（个）‘）;title（‘仿真回波数据‘）;
set（get（gca‘xlabel‘）‘FontSize‘16）;
set（get（gca‘ylabel‘）‘FontSize‘16）;
set（get（gca‘title‘）‘FontSize‘16）;
set（gca‘FontName‘‘Times New Roman‘‘FontSize‘16）;
save echos   echos  
%=====================输出显示回波数据 ==================

%=====================将回波数据处理为实测数据 =============
for m=1:Ma
    for n=1:Mr
        echos_real（n）=echos（mn）;
    end
    echos_real=ifftshift（ifft（echos_real））;             % 对行（距离向）傅氏逆变换
    for n=1:Mr
        echos_final（mn）=echos_real（n）;
    end
end
%=

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4809  2017-09-30 11:41  点目标的飞机模型ISAR仿真程序\ISAR_Imaging.m
     文件         438  2017-09-27 17:47  点目标的飞机模型ISAR仿真程序\Target_Plane_SeventyFour_Point.mat
     文件         673  2016-01-15 20:43  点目标的飞机模型ISAR仿真程序\Target_Plane_SeventyFour_Point1.m
     文件      253804  2017-09-30 11:41  点目标的飞机模型ISAR仿真程序\echos.mat
     文件        2192  2017-09-30 11:41  点目标的飞机模型ISAR仿真程序\echos_real.mat
     目录           0  2017-09-30 16:32  点目标的飞机模型ISAR仿真程序\