一、原始回波数据
二、时域脉冲压缩
三、后向投影
四、部分代码(完整代码见下载链接)
clc;
clear all;
close all;
%% 参数定义
%=========================================
c=3e8;
j=sqrt(-1);
pi=3.1416;
fc=5.3e9;
lamda=c/fc;
D=4; % 方位向天线尺寸
Va=150; % 等效雷达速度
Kr=20e12; % 距离调频率
Tr=2.5e-6; % 发射脉冲时宽
sq_ang=3.5/180*pi; % 波束斜视角 (增大波束斜视角,方向维合成孔径长度增加)
Br=Kr*Tr; % 带宽
Frfactor=1.2; % 距离向过采样率
Fr=Br*Frfactor; % 距离向采样频率
Ba=0.886*2*Va*cos(sq_ang)/D; % 多普勒带宽
Fafactor=1.2; % 方位向过采样率
Fa=Ba*Fafactor; % 方位向采样频率
R_near=2e4; % 场景近端(零多普勒位置斜距)
R_far=R_near+1000; % 场景远端(零多普勒位置斜距)
La_near=0.886*R_near*lamda/(cos(sq_ang).^2)/D; % 场景近端合成孔径长度
La_far=0.886*R_far*lamda/(cos(sq_ang).^2)/D; % 场景远端合成孔径长度
Tc_near=-R_near*tan(sq_ang)/Va; % 场景近端波束中心穿过时间(以零多普勒位置为参考时间)
Tc_far=-R_far*tan(sq_ang)/Va; % 场景远端波束中心穿过时间(以零多普勒位置为参考时间)
fdc=2*Va*sin(sq_ang)/lamda; % 多普勒中心频率
Rmin=sqrt(R_near^2+(Tc_near*Va+La_near/2)^2); % 雷达和该条带的最短斜距
Rmax=sqrt(R_far^2+(Tc_far*Va-La_far/2)^2); % 雷达和该条带的最远斜距
disp('参数');
disp('雷达的实波束宽度');disp(0.886*lamda/D);
disp('雷达的斜视角');disp(sq_ang/pi*180);
disp('雷达和场景间的最短斜距:');disp(Rmin);
disp('雷达和场景间的最远斜距:');disp(Rmax);
disp('距离分辨率:');disp(0.886*(c/2/Br));
disp('方位分辨率:');disp(0.886*Va/Ba);
disp('多普勒中心:');disp(fdc);
%% 回波模型
%====================================================
Nr=(2*Rmax/c+Tr-2*Rmin/c)*Fr;
Nr=2^nextpow2(Nr); % 距离向采样点数
tr=linspace(-Tr/2+2*Rmin/c,Tr/2+2*Rmax/c,Nr); % 距离向采样时间(以每个发射脉冲时刻为参考时间)
Fr=(Nr-1)/(Tr/2+2*Rmax/c-(-Tr/2+2*Rmin/c)); % 距离向采样频率
Na=((Tc_near+La_near/2/Va)-(Tc_far-La_far/2/Va))*Fa;
Na=2^nextpow2(Na); % 方位向采样点数
ta=linspace(Tc_far-La_far/2/Va,Tc_near+La_near/2/Va,Na); % 方位向采样时间(以零多普勒位置为参考时间)
Fa=(Na-1)/(Tc_near+La_near/2/Va-(Tc_far-La_far/2/Va)); % 方位向采样频率
disp('由方向维采样频率限制的雷达与场景任意一点的最远斜距:');disp(1/Fa*c/2);