資源簡介
根據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
- 上一篇:UML食堂飯卡系統報告設計
- 下一篇:解決了關閉死鎖的CSerialPort類項目
評論
共有 條評論