%tic;
clear;
load?‘F:\實驗室課題\雷達成像\雷達成像\ISAR\data‘;

s1=fft（s）;%fft對矩陣按列進行傅里葉變換
for?i=1:256;
???s1（:i）=fftshift（s1（:i））;%把直流分量移到圖像中
end

figure（1）;
subplot（211）;
contour（abs（s1））;
title（‘雅克42對準前的數據‘）;

%2.?ISAR運動補償距離對準空域法
s2=abs（s1）;?????????????????%從前向后
s2=fft（s2）;
tmp=abs（s1（1:256129））;
TMP=fft（abs（tmp））;
s3=s1;
for?i=1:256
???for?k=1:256
??????out（ki）=TMP（k）*conj（s2（ki））;
???end
??????out1=ifft（out（:i））;
??????out1=abs（out1）;
??????[mj]=max（out1）;
??????%for?i=1:256
?????????if（j>128）
????????????j=j-256;
?????????end??
?????????for?k=1:256
????????????J=k-j;?
	????????if（（J>0）&（J<=256））?%右移或左移
	???????????s3（ki）=s1（Ji）;
	???	????else?
??????	???????s3（ki）=0;
	????????end???
	?????end
?????????%end
??
???if（i~=129）
?????????tmp=tmp+abs（s3（:i））;
?????????if?i<129
?????????????t=i;
?????????else
?????????????t=i-1;?
?????????end?????????
?????????tmpt=tmp/（t+1）;
?????????TMP=fft（tmpt）;
???end
end

figure（5）
subplot（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-05-13?21:05??ISAR\
?????文件?????1048768??2013-02-25?20:14??ISAR\data.mat
?????文件??????????51??2013-02-25?20:14??ISAR\data說明.txt
?????文件????????2802??2013-02-25?20:14??ISAR\isar.m