%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