isar成像程序——适合初学者

根据ISAR概念写的，没有考虑包络对齐、初相补偿等。 ISARImage为主程序，rotate为角度旋转函数

%ISAR 成像仿真
%分辨率为0.3m×0.3m

clc;clear;
%%%%%%%%参数设置%%%%%%%%%%%%%%%%%%%%%%%%%%
c=3e8;    %光速
fc=10e9;  %载频10GHz
lamda=c/fc; 
ruo_r=0.3;  %距离分辨率
ruo_a=0.3;  %方位分辨率
B=c/2/ruo_r;   %信号带宽
theta=lamda/2/ruo_a;  %累计旋转角度
pw=1e-6;    %信号脉冲宽度
gama=B/pw;  %调制斜率
ft=2.5*B;   %单个脉冲采样频率
Ns=pw*ft;   %单个脉冲采样点数
t=linspace（-pw/2pw/2Ns）;  %快时间
dt=pw/Ns;   %时间采样间隔
Chirp=exp（j*pi*gama*t.*t）;   %线性调频脉冲

%%%%%%%%目标坐标信息%%%%%%%%%%%%%%%%%%%%%%%%
%目标大小范围5m×5m
Tar=[0 5; -1 4; 1 4; -1 3; 1 3; -1 2; 1 2; 1 1; -1 1; -1 0; 1 0; 2 0; -2 0; 3 0; -3 0; 4 0; -4 0; 4 -1; -4 -1; -3 -2; 3 -2; 2 -3; -2 -3; 1 -4; -1 -4; 0 -4 ];
figure;plot（Tar（:1）Tar（:2）‘*‘）;
axis（[-20 20 -20 20]）;


%%%%%%%针对目标生成回波%%%%%%%%%%%%%%%%%%%%%
%逆时针旋转取128次照射回波
[t1 t2]=size（Tar）;
PR=[0 -1000];   %雷达位置
RSee=[990 1200];   %观察记录范围，距离雷达
Nr=（RSee（2）-RSee（1））/（c*dt/2）;  %距离向点数
Na=128;    %方位向点数
Sr=zeros（NaNr）;
deltat=linspace（0thetaNa）;  %每个重频周期回波时旋转的角度
for ii=1:Na
    ttrar=rotate（Tardeltat（ii））;
    for jj=1:t1
        R=sqrt（（ttrar（jj1）-PR（1））^2+（ttrar（jj2）-PR（2））^2）;
        n=floor（（R-RSee（1））/（c*dt/2））;
        tChirp=Chirp*exp（-j*4*pi*R/lamda）;
        Sr（iin:n+Ns-1）=Sr（iin:n+Ns-1）+tChirp;
    end
end


%%%%%%%%%%%%%%匹配滤波成像%%%%%%%%%%%%%%%%%%%%%%%%%
Hr=exp（-j*pi*gama*t.*t）;
Hr=Hr.*hamming（Ns）‘;
hr=[Hr zeros（1Nr-Ns）];
Sout=Sr;
for ii=1:Na
    Sout（ii:）=fftshift（ifft（（fft（Sr（ii:）））.*（fft（hr））））;
end
tt=Sout‘;
out=zeros（NrNr）;
for ii=1:Nr
    out（ii:）=[tt（ii:） zeros（1Nr-Na）];
    out（ii:）=fftshift（fft（out（ii:）））;
end
Sout=out;
kk=abs（Sout）;
figure;imagesc（kk）;
axis on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1776  2010-10-26 21:18  ISARImage.m

     文件        759  2010-10-24 16:09  rotate.m

     文件         89  2011-04-11 19:22  新建文本文档.txt

----------- ---------  ---------- -----  ----

                 2624                    3