雷达回波信号的脉冲压缩matlab仿真

 

%LFM脉冲的匹配滤波
clear all;
T=10e-6;                                   %发射脉宽10us
B=30e6;                                   %调频带宽30MHz
K=B/T;                                    %频率调制斜率
Fs=10*B;Ts=1/Fs;                           %计算机仿真的采样频率和采样周期
N=T/Ts;
t=linspace(-T/2,T/2,N);
St=exp(j*pi*K*t.^2);                         %产生线性调频信号
Ht=exp(-j*pi*K*t.^2);                        %匹配滤波器单位冲激响应
Sot=conv(St,Ht);                            %匹配滤波输出
figure(1)
L=2*N-1;
t1=linspace(-T,T,L);
Z=abs(Sot);Z=Z/max(Z);                     %仿真信号对数归一化
Z=20*log10(Z+1e-6);
Z1=abs(sinc(B.*t1));                        %产生理伦输出信号并归一化
Z1=20*log10(Z1+1e-6);
t1=t1*B;                                  %时间轴与1/B归一化
plot(t1,Z,t1,Z1,'r.');
axis([-15,15,-50,inf]);grid on;
legend('\fontsize{9}仿真结果','理论结果');
xlabel('\fontsize{9}时间\times\itB  (s)');
ylabel('\fontsize{9}相对幅度dB');
title('\fontsize{9}脉冲压缩后的回波(归一化后)');
figure(2)                                   %仿真结果局部放大
N0=3*Fs/B;
t2=-N0*Ts:Ts:N0*Ts;
t2=B*t2;
plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');
axis([-inf,inf,-50,inf]);grid on;
set(gca,'Ytick',[-13.4,-4,0],'Xtick',[-3,-2,-1,-0.5,0,0.5,1,2,3]);
legend('\fontsize{9}仿真结果','理论结果');
xlabel('\fontsize{9}时间\times\itB  (s)');
ylabel('\fontsize{9}相对幅度dB');
title('\fontsize{9}脉冲压缩后的回波(归一化后局部图)');
D154