pd_5chirp_com_xiugai_2.m

基于MATLAB的雷达回波成像与脉冲压缩程序。

%%=====================================================================
clear;clc;close all;
%%=====================================================================
C=3e8; 								% 光速,单位为米/秒
Wavewidth=15;						% 主波束宽度,单位为度
Wb=Wavewidth/2;
%%=====================================================================
%%扫描向参数
P=10;								% 雷达扫描速度,单位为度/秒
S=60;								% 扫描范围，单位为度
PRF=100;							% 脉冲重复频率,单位为Hz
PRT=1/PRF; 							% 脉冲周期,单位为秒
ds=PRT; 							% 扫描向采样间距
Tscan=S/P;							% 天线扫描一次所需时间，单位为秒
M=ceil(Tscan/ds); 					% 扫描向采样数
M=2^nextpow2(M);					% 使FFT速度快
Sn=linspace(-S/2,S/2,M); 			% 扫描向角度离散化
%%=====================================================================
%%点目标实际位置参数
Ntarget=3; 							% 点目标数
Ro=1000; 							% 径向距离，单位为m
Sita1=0;
R1=Ro;
Sita2=-7;
R2=Ro;
Sita3=7;
R4=Ro/2;
Rmax=Ro;
Rmin=R4;
r=Rmax-Rmin;						% 目标最大间距
%%=====================================================================
%%发射信号参数
Tr=5e-7;                         	% 脉冲持续时间,单位为秒
Br=30e6;                        	% 信号带宽:30MHz,频偏
Kr=Br/Tr;                        	% 调频斜率
Fsr=3*Br;                        	% 采样频率,应大于2倍的最高频率
dt=1/Fsr;                         	% 采样间隔
%%=====================================================================
%%图像中距离向区域参数
Xc=1000;                        	% 图像中心
Ptarget=[Sita1,		Xc              % 目标相对位置
		 Sita2,		Xc
		 Sita3,		Xc-R4];
%%将距离转化成时间
tr=r/C;
deltaT=2*(Tr+tr+Tr);
deltaX=2*(Tr+tr+Tr)*C;				% 图像宽度
Xo=deltaX/2;                        % 图像范围:[Xc-Xo,Xc+Xo]
Tc=Xc/C;
To=deltaT/2;
N=ceil(deltaT/dt); 					% 距离向采样数
N=2^nextpow2(N);
Tn=linspace(Tc-To,Tc+To, N); 		% 距离向离散化
%%=====================================================================
%%产生回波信号数据
Rtim=ones(M,1)*Tn; 					% 距离窗
Scan=Sn'*ones(1,N); 				% 扫描窗
t=Rtim;
phase=pi*Kr*t.^2;
T=Ptarget; 							% 目标位置
PDmn=zeros(M,N); 					% 定义一个全零矩阵
for	k=1:1: Ntarget
Smn=(T(k,1)-Wb<Scan&Scan<T(k,1)+Wb);
Rmn=(T(k,2)/C-Tr/2<Rtim&Rtim<T(k,2)/C+Tr/2);
PDmn=PDmn+ exp(j*phase).*Smn.*Rmn;
end
%%=====================================================================
%%回波成像
colormap(gray); 					% colormap色图
figure(1); 							% figure创建图像窗
imagesc(real(PDmn));        		% imagesc显示高度图像
xlabel('\rightarrow\it直径方向'); 	% xlabel X轴名称
ylabel('\leftarrow\it扫描方向');
title('PD雷达点目标回波');
%%=====================================================================
%%距离压缩
Refr=exp(j*phase).*Rmn; 			% 参考函数
PDcom=ifty(fty(PDmn).*conj(fty(Refr)));
colormap(gray); 					% colormap色图
figure(2); 							% figure创建图像窗
imagesc(real(PDcom));
xlabel('\rightarrow\it直径方向'); 	% xlabel X轴名称
ylabel('\leftarrow\it扫描方向');
title('PD雷达点目标图像（经过距离压缩）');
%%=====================================================================