filter.m

用于绘制LMF波形匹配滤波器输出波形的包络

%
%本程序在5.4节中，用于绘制匹配滤波器输出波形的包络
%
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}信号卷积匹配滤波器结果','\fontsize{9}模糊函数计算结果');
xlabel('\fontsize{9}时间\times\itB  (s)');
ylabel('\fontsize{9}相对幅度dB');
title('\fontsize{9}脉冲压缩后的回波(归一化后)');