条带sar成像实验报告

clear all 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); %二维数组 for k=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)

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

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

if(abs(cita)

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

Rt=sqrt((xr(k)-xm3)^2+(yr-ym3)^2); cita=asin((xm3-xr(k))/Rt); if(abs(cita)

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

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; %频域加窗，失配后的滤波器 for k=1:Na

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

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

figure(2)

imagesc(abs(data(:,1:3:Na))) title('距离压缩后图像')

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

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

figure(3) imagesc(abs(data)) title('实现方位校正后图像')

fa=linspace(-PRF/2,PRF/2,Na);%方位频率 det_r=lamd^2*R0/(8*V^2)*fa.^2; t0=2*det_r/c;%关于方位频率的变量

for k=1:Na

r_shift=exp(j*2*pi*fr*t0(k)); %计算每个频率k的距离徙动

r_shift=r_shift.';%转置，向量乘法

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

figure(4)

imagesc(abs(data))

title('校正过距离徙动后的图像')

R=c*tr/2; %目标距离 for k=1:Nr

ka=-2*V^2/lamd/R(k); %对于每个目标计算线性调频参数 hfa=exp(j*pi*fa.^2/ka); %方位压缩滤波器设计 data(k,:)=ifft(data(k,:).*hfa);%同距离压缩 end

figure(5) imagesc(abs(data)) title('方位压缩后图像')

figure(6)

imagesc(abs(data)) plot(abs(data(128,:)))

title('未加窗目标（0,0），（50,0）SAR成像')

figure(7)

imagesc(abs(data)) plot(abs(data(88,:)))

title('未加窗目标（0,50）SAR成像')

% figure(8)

% imagesc(abs(data)) % plot(abs(data(128,:)))

% title('加窗后目标（0,0），（50,0）SAR成像')

% figure(9)

% imagesc(abs(data)) % plot(abs(data(88,:)))

% title('加窗后目标（0,50）SAR成像')

本文来源：https://www.bwwdw.com/article/4vbd.html