clear

f0=10e9;???%載波頻率
B=500e6;???%?線性調頻信號帶寬
T=100e-6;%102.4e-6;??%脈沖寬度
k=B/T;?????%?調頻斜率
c=3e8;?????%光速
fs=1024/T;%20e6;???%采樣率
R0=15e3-100;???%?目標距離
Rf=15e3;?????%參考距離

R=[55];????%散射點坐標

%目標回波
num=128;???%方位向采樣數
ang1=2;????%總的轉角??度
ang2=ang1/num;??%單位轉角??度
for?w=1:num
????ang（w）=（w-1）*ang2*pi/180;
????R1（w）=R0+R*[cos（ang（w））sin（ang（w））]‘;

delay1=2*Rf/c;??????%?參考延時
delay（w）=2*R1（w）/c;???????%?目標延時
for?m=1:1024%1024
????p（m）=m;
????t（m）=-T/2+（m-1）*1/fs;??%?目標回波信號采樣
????s0（m）=exp（j*2*pi*f0*（t（m）-delay（w））+j*pi*k*（t（m）-delay（w））^2）;??????????%?目標回波信號
????s1（m）=exp（j*2*pi*f0*（t（m）-delay1）+j*pi*k*（t（m）-delay1）^2）;???????????????%?參考本振信號??
????s2（wm）=s0（m）.*conj（s1（m））;
end

end
??
%??strech?脈沖壓縮處理
for?q=1:128
??sr（q:）=fft（s2（q:））;??????????%對每個轉角的回波數進行距離維的壓縮
??n1（q）=q;
??
?%?plot（pabs（sr（q:）））
end

for?f=1:1024
????sr2（:f）=fftshift（fft（sr（:f）））;
????hold?on
????n2（f）=f;
???%plot（n1abs（sr2（:f）））
end
??%contour（abs（sr））;

contour（abs（sr2））


??
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1131??2008-03-18?13:10??isar_image.m

-----------?---------??----------?-----??----

?????????????????1349????????????????????2