www.pudn.com > isar.zip > isar_image.m, change:2013-04-09,size:1140b


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=[5,5];    %散射点坐标 
 
%目标回波 
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(w,m)=s0(m).*conj(s1(m)); 
end 
 
end 
   
%  strech 脉冲压缩处理 
for q=1:128 
  sr(q,:)=fft(s2(q,:));          %对每个转角的回波数进行距离维的压缩 
  n1(q)=q; 
   
 % plot(p,abs(sr(q,:))) 
end 
 
for f=1:1024 
    sr2(:,f)=fftshift(fft(sr(:,f))); 
    hold on 
    n2(f)=f; 
   %plot(n1,abs(sr2(:,f))) 
end 
  %contour(abs(sr)); 
figure(1) 
contour(abs(sr2))