www.gusucode.com > 波前重构SAR图像程序算法源码程序 > 波前重构SAR图像程序算法源码程序/105230337wavefront-algorithm/Untitled4.m

    clear; 
clc;
close all;

r=3e-2;B=60e6;fs=90e6;prf=60;ls=300;R=50e3;c=3e8;T=2e-6;v=100;ts=ls/(2*v); 
pa=r*R/(2*ls);  pr=c/(2*B);  k=B/T; 
qt=-T/2:(1/fs):T/2-1/fs; 
t=-ts:(1/prf):ts-1/prf; 
rt=(R.^2+(v*t).^2).^2; 
f=-fs/2:(fs/256):fs/2-fs/256; 
hf=exp((i*2*pi*f).^2/(4*pi*k)); 
lfm=exp(i*pi*k*qt.^2)'*exp(-i*4*pi*rt/r); 
[l,L]=size(lfm); 
hw1=exp(i*(2*pi*f).^2/(4*pi*k)); 
ka=-2*v.^2/(r*k); 
w=-prf/2:(prf/256):prf/2-prf/256; 
hw2=exp(j*w.^2/(4*pi*ka)); 
s=1:256; 
lfm=lfm'; 
for m=1:L 
    s1(m,s)=fft(lfm(m,1:180),256); 
    s1(m,s)=fftshift(s1(m,s)); 
    s1(m,s)=ifft(s1(m,s).*hw1); 
    s1(m,s)=ifftshift(s1(m,s)); 
end 
s1=s1'; 
for k=1:256 
          s2(k,s)=fft(s1(k,1:180),256); 
          s2(k,s)=fftshift(s2(k,s)); 
          s2(k,s)=ifft(s2(k,s).*hw2); 
          s2(k,s)=ifftshift(s2(k,s)); 
end 
[m,n]=size(s2) 
subplot(2,1,1) 
contour(abs(s2)) 
subplot(2,1,2) 
mesh(abs(s2)) 
clear all