ISAR雷达成像

提供ISAR雷达成像matlab代码，有原始数据，可直接运行。文件中还有关于ISAR的ppt,很详细的资料

%tic;
clear all
close all
clc
load ‘D:\文献\ISAR\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对准前的数据‘）;
 
s2=fft（abs（s1））;
for j=1:256
    for i=1:256
        phs2（ij）=conj（s2（ij））*（s2（i129））;%利用互相关求取相位
    end
end
PHs2=ifft（phs2）;%反傅里叶变化确定相位
MaPHs2=abs（PHs2）;%取g模值
[mj]=max（MaPHs2）;% 将 MagPHs2 每列中的最大值赋给m，行数赋给j
for i=1:256
    if（j（i）>128）
        j（i）=j（i）-256;%保证j中的数值在-128——128中
    end
    for k=1:256
        J=k-j（i）;
        if（（J>0）&&（J<=256）） %右移或左移
            s3（ki）=s1（Ji）;
        else
            s3（ki）=0;
        end
    end
end
p=abs（s3）;
subplot（212）;
contour（abs（s3））;
title（‘频域补偿后‘）;
%3. ISAR运动补偿--- 相位补偿
e（1）=1;
for i=2:256
    sum1=0;
    for j=1:256
        sum1=sum1+conj（s3（ji-1））*s3（ji）;
    end
    sum1=sum1/abs（sum1）;
    e（i）=-e（i-1）*sum1;
end
for i=1:256
    w（i）=0.54-0.46*cos（2*pi*i/256）;
%加窗
end
for i=1:256
    for j=1:256
        fi（ji）=w（j）*s3（ij）*conj（e（j））;
    end
end
fi=fft（fi）;
ss=abs（fi）;
ss=rot90（ss）;
figure（2）;
contour（ss8）;
xlabel（‘二维成像‘）;
title（‘运动补偿结果（频域法）‘）;
 
ee=0;
for n=1:256
    for m=1:256
        ee=ee+ss（mn）^2*log（ss（mn）^2）;
    end
end
 
 
%3. ISAR运动补偿--- 相位补偿
e（1）=1;
for i=2:256
   sum1=0;
   for j=1:256
      sum1=sum1+conj（s3（ji-1））*s3（ji）;
   end   
   sum1=sum1/abs（sum1）;
   e（i）=e（i-1）*sum1;  
end
for i=1:256
     w（i）=0.54-0.46*cos（2*pi*i/256）;%加窗
end
 
for i=1:256
   for j=1:256
    fi（ji）=w（j）*s3（ij）*conj（e（j））;
   end
end
 
fi=fft（fi）;             
for k=1:256
    fi（:k）=fftshift（fi（:k））;
end
ss=abs（fi）;
ss=rot90（ss）;
figure（3）;
contour（ss10）;
xlabel（‘雅克42 ISAR 目标二维像‘）;
%title（‘距离对准相位补偿研究‘）;
 
[xy]=max（ss）;
mm=max（x）;
for i=1:256
  for j=1:256
     ss（ij）=20*log10（ss（ij）/mm）+21;
      if（ss（ij）<0）
       ss（ij）=0;
      else
       ss（ij）=ss（ij）+120;
      end
   end
end
figure（）;
contour（ss10）;
%xlabel（‘雅克42 ISAR 目标二维像‘）;
title（‘（整体相关法）经图像处理后的结果‘）;
%toc

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

     文件    1048768  2007-04-17 11:40  ISAR\data.mat

     文件         51  2008-04-08 18:30  ISAR\data说明.txt

     文件   16208896  2007-04-12 09:43  ISAR\ISAR&SAR.ppt

     文件       1788  2009-05-05 17:43  ISAR\kongyu.asv

     文件       2263  2016-05-23 20:26  ISAR\zhengti.m

     文件     606720  2007-04-12 10:58  ISAR\成像过程.doc

     文件        916  2007-07-04 14:29  ISAR\新建 文本文档.txt

     文件     152064  2007-04-12 16:05  ISAR\逆合成孔径雷达成象.ppt

     目录          0  2016-05-23 20:30  ISAR

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

             18021466                    9