pipeilvbaoqi.m

匹配滤波的MATLAB实现

clear;clc;
B=6e8;
w0=1e8;
t=1e-6;
k=B/(2*t);
fs=5*w0;  %采样频率
t1=-t:1/fs:t;
t2=-2*t:1/fs:2*t;
s=10*cos(w0*t1+k*(t1.*t1)/2);   %产生输入信号
N=length(s);             %统计输入信号长度
n=randn(1,N);       %产生对应的白噪声
x=s+n;                   %得到复合信号
X=fft(x);
              %对输入信号进行傅立叶变换
s1=fft(s);
S=conj(s1);  
H=S.*exp(-i*w0*t1*1.5);  %求得匹配滤波器
h=ifft(H);               %对匹配滤波器做反傅氏变换
%Y=X.*H;
%y=ifft(Y);
y=conv(x,h);             %对复合信号进行滤波处理
y=abs(y);
figure
plot(t2,y),title('输出信号y');
figure 
plot(n),title('白噪声n');
figure
plot(t1,x),title('复合信号x');
figure 
plot(t1,s),title('输入信号s');