条带sar成像实验报告.pdfVIP

clearall

clc

Tp=1.e-6;%脉冲宽度

B=100.e6;%带宽

kr=B/Tp;

f0=10.e9;%载频

c=3.e8;

lamd=c/f0;

V=150.;%雷达飞行速度

PRF=1800.;

Na=2048*4;%方位采样点数

Ta=Na/PRF;%方位采样总时间

ta=linspace(-Ta/2,Ta/2,Na);%方位采样时间定义

beamwidth=9./180*pi;

xr=V*ta;%雷达参数

yr=3000.;

xm1=0.;%目标参数

ym1=0.;

xm2=50.;

ym2=0.;

xm3=0.;

ym3=50.;

Nr=256;%距离采样点数

fs=120.e6;

R0=3000.;

tr=linspace(2*R0/c-1/fs*Nr/2,2*R0/c+1/fs*Nr/2,Nr);%距离采样时间定义

data=zeros(Nr,Na);%二维数组

fork=1:Na

record=zeros(1,Nr);%回波信号，关于每个方位时间数组

Rt=sqrt((xr(k)-xm1)^2+(yr-ym1)^2);%目标1距离雷达的距离

cita=asin((xm1-xr(k))/Rt);%反正弦

%方位限制

if(abs(cita)beamwidth/2)

%计算目标回波信号

record=exp(j*pi*kr*(tr-2*Rt/c).^2).*(abs(tr-2*Rt/c)Tp/2).*exp(-j*2*pi*f0*2*

Rt/c);

end

Rt=sqrt((xr(k)-xm2)^2+(yr-ym2)^2);

cita=asin((xm2-xr(k))/Rt);

if(abs(cita)beamwidth/2)

record=record+exp(j*pi*kr*(tr-2*Rt/c).^2).*(abs(tr-2*Rt/c)Tp/2).*exp(-j*2*p

i*f0*2*Rt/c);

end

Rt=sqrt((xr(k)-xm3)^2+(yr-ym3)^2);

cita=asin((xm3-xr(k))/Rt);

if(abs(cita)beamwidth/2)

record=record+exp(j*pi*kr*(tr-2*Rt/c).^2).*(abs(tr-2*Rt/c)Tp/2).*exp(-j*2*p

i*f0*2*Rt/c);

end

data(:,k)=record;

end

figure(1)%传统雷达图像

imagesc(abs(data))

title(传统雷达图像)

%条带SAR成像经典算法

fr=linspace(-fs/2,fs/2,Nr);%距离频率

fr=fftshift(fr);

hfr=exp(j*pi*fr.^2/kr);%滤波器，滤波器设计的三种算法之一，对于线性调频信号，直接根据公式

计算设计

hfr=transpose(hfr);

%win=hamming(Nr);%汉明窗

%win=transpose(win);

%hw=hfr.*win;%频域加窗，失配后的滤波器

fork=1:Na

data(:,k)=ifft(fft(data(:,k)).*hfr);

%data(:,k)=ifft(fft(data(:,k)).*hw);%对回波信号作fft后乘以滤波器再作逆fft

end

figure(2)

imagesc(abs(data(:,1:3:Na)))

title(距离压缩后图像)

fork=1:Nr%方位时间转化方位频率

data(k,:)=fftshift(fft(data(k,:)));

end

fig