雷达后向投影算法MATLAB实现

一、原始回波数据

二、时域脉冲压缩

三、后向投影

四、部分代码（完整代码见下载链接）

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);

下载链接https://download.csdn.net/download/qq_41186819/11966747