匹配滤波器原理

1、数字通信课程设计匹配滤波器摘要在通信系统中，滤波器是重要的部件之一，滤波器特征的选择直接影响数字 信号的恢复。在数字信号接收中，滤波器的作用有两个方面，使滤波器输出有用 信号成分尽可能强；抑制信号带外噪声，使滤波器输出噪声成分尽可能小，减少噪声对信号判决的影响。对最佳线性滤波器的设计有一种准则是使滤波器输出信 噪比在特定时刻到达最大，由此导出的最佳线性滤波器称为匹配滤波器。在数字通信中，匹配滤波器具有广泛的应用。因此匹配滤波器是指滤波器的性能与信号 的特征取得某种一致，使滤波器输出端的信号瞬时功率与噪声平均功率的比值最 大。本文设计并仿真了一种数字基带通信系统接收端的匹配滤波器。一、课程设计的

2、目的通过本次对匹配滤波器的设计，让我们对匹配滤波器的原理有更深一步的理 解，掌握具体的匹配滤波器的设计方法与算法。二、课程设计的原理设接收滤波器的传输函数为H(f),冲击响应为h(t),滤波器输入码元s(t)的持续时间为Ts,信号和噪声之和r(t)为式中，s(t)为信号码元，n(t)为白噪声。并设信号码元s(t)的频谱密度函数为S(f),噪声n(t)的双边功率谱密度为Pn no/2 , no为噪声单边功率谱密度。假定滤波器是线性的，根据叠加定理，当滤波器输入信号和噪声两部分 时，滤波器的输出也包含相应的输出信号和输出噪声两部分，即r*2由于：PY(f) H (f)H(f)PR(f) H(f)

3、f(f)PJf)为输出功率谱密度，PR(f)为输入功率谱密度，PR(f) n。/2这时的输出噪声功率No等于在抽样时刻t0上，输出信号瞬时功率与噪声平均功率之比为为了求出ro的最大值，利用施瓦兹不等式求ro的最大值且当时 H(f) kS*(f)ej2 fto，等式成立，即得到知道的信噪比为2En。在白噪声干扰的背景下，按上式的设计的线性滤波器，将能在给定时刻to上获得最大输出信噪比2Eno。是输出信噪比最大时刻。这种滤波器就是最大信 噪比意义下的最佳线性滤波器，由于它的传输特性与信号频谱的复共轭一致，称为匹配滤波器。匹配滤波器的特性还可以用冲激响应函数 h(t):可见，匹配滤波器的冲激响应h(

4、t)就是信号s(t)的镜像s(to t)，但在时间轴上平移了 to o to是输出信噪比最大时刻。三、课程设计的内容设计仿真了数字基带通信系统在接收端接收到模拟信号后通过匹配滤波 器的效果，其中匹配滤波器的输入信号我们假定为正弦信号，经过匹配滤波 器的输出信号作为抽样判决器的输入信号。我们通过比较匹配滤波器的输入 输出信号，验证得出其性能特性。如下图所示的仿真结果，我们假定匹配滤波器的输入信号为正弦信号与白 噪声的叠加，其频率，采样频率,延时时间，。由实验的结果，我们可以看出匹配滤波器的单位冲激响应幅度谱与输入信号 的幅度谱一致，符合上面我们推导的公式；由滤波器的输出信号我们可以得 出，在抽样

5、时刻处，信号有最大值，在输出信号的幅度谱中，干扰噪声白噪 声明显减小，信号的有用成分被保留下来，体现出了匹配滤波器抑制噪声， 保证在抽样时刻有最大值的特性，从而进一步保证了输出信噪比在抽样时刻 最大，为下一步抽样判决提供了有利的条件。由此可以看出匹配滤波器对数 字基带通信系统的通信可靠性的提高有着很重要的作用。四、附录实验程序：clearall ;closeall ;fc=20;%余弦信号的频率fs=5*fc;%采样频率ts=1/fs;k=1;to=2;%采样时刻t=0:ts:2-ts;t1=0:ts:4-2*ts; freq1=linspace(-fs/2,fs/2,length(t); f

6、req2=linspace(-fs/2,fs/2,length(t1);Noise=randn(1,length(t);%白噪声sig nal=3*cos(2*pi*fc*t+pi*k*t.A2)+Noise; h1=3*cos(2*pi*fc*(to-t)+pi*k*(to-t).A2)+Noise;So1=conv2(signal,h1);%信号卷积subplot(321)plot(t,signal); axis(0,2,-10,10); xlabel( 't/s');ylabel( 'S(t)');title( ' 输入信号 +噪声 ' )

7、; subplot(322) plot(freq1,fftshift(abs(fft(signal)/length(t),xlabel( el( '|S(W)|' ),title( ' 输入信号幅度谱 ' ); subplot(323)plot(t,h1); axis(0,2,-10,10); xlabel( 't/s');ylabel( 'h(t)');title( ' 冲激响应 ' );subplot(324) plot(freq1,fftshift(abs(fft(h1)/length(t),xlabel( |H1(W)|' ),title( ' 单位冲激响应幅度谱 ' );subplot(325) plot(t1,So1); xlabel( 't/s');ylabel( 'So(t)');title( ' 匹配滤波器输出信号 ' );subplot(326) plot(freq2,fftshift(abs(fft(So1)/length(t1),xlabel( '|So