匹配滤波器原理

数字通信课程设计匹配滤波器摘要在通信系统中，滤波器是重要的部件之一，滤波器特征的选择直接影响数字信号的恢复。在数字信号接收中，滤波器的作用有两个方面，使滤波器输出有用信号成分尽可能强；抑制信号带外噪声，使滤波器输出噪声成分尽可能小，减少噪声对信号判决的影响。对最佳线性滤波器的设计有一种准则是使滤波器输出信噪比在特定时刻到达最大，由此导出的最佳线性滤波器称为匹配滤波器。在数字通信中，匹配滤波器具有广泛的应用。因此匹配滤波器是指滤波器的性能与信号的特征取得某种一致，使滤波器输出端的信号瞬时功率与噪声平均功率的比值最大。本文设计并仿真了一种数字基带通信系统接收端的匹配滤波器。一、课程设计的目的通过本次对匹配滤波器的设计，让我们对匹配滤波器的原理有更深一步的理解，掌握具体的匹配滤波器的设计方法与算法。二、课程设计的原理设接收滤波器的传输函数为)(fH,冲击响应为)(th,滤波器输入码元)(ts的持续时间为sT,信号和噪声之和)(tr为)()()(tntstrsTt0式中，)(ts为信号码元，)(tn为白噪声。并设信号码元)(ts的频谱密度函数为)(fS,噪声)(tn的双边功率谱密度为2/0nPn,0n为噪声单边功率谱密度。假定滤波器是线性的，根据叠加定理，当滤波器输入信号和噪声两部分时，滤波器的输出也包含相应的输出信号和输出噪声两部分，即)()()(00tntsty由于：)()()()()()(2*fPfHfPfHfHfPRRY)(fPR为输出功率谱密度，)(fPR为输入功率谱密度，2/)(0nfPR这时的输出噪声功率0N等于dffHndfnfHN20020)(22)(在抽样时刻0t上，输出信号瞬时功率与噪声平均功率之比为dffHndfefSfHNtsrftj20220000)(2)()()(0为了求出0r的最大值，利用施瓦兹不等式求0r的最大值02022202202)(2)()()(2)()(0nEdffHndffSdffHdffHndfefSfHrftj且当时02*)()(ftjefkSfH，等式成立，即得到知道的信噪比为02nE在白噪声干扰的背景下，按上式的设计的线性滤波器，将能在给定时刻0t上获得最大输出信噪比02nE。是输出信噪比最大时刻。这种滤波器就是最大信噪比意义下的最佳线性滤波器，由于它的传输特性与信号频谱的复共轭一致，称为匹配滤波器。匹配滤波器的特性还可以用冲激响应函数)(th:022*2()()()jftjftjfthtHfedfkSfeedf**()()()00ksttdkstt可见，匹配滤波器的冲激响应)(th就是信号)(ts的镜像)(0tts，但在时间轴上平移了0t。0t是输出信噪比最大时刻。三、课程设计的内容设计仿真了数字基带通信系统在接收端接收到模拟信号后通过匹配滤波器的效果，其中匹配滤波器的输入信号我们假定为正弦信号，经过匹配滤波器的输出信号作为抽样判决器的输入信号。我们通过比较匹配滤波器的输入输出信号，验证得出其性能特性。如下图所示的仿真结果，我们假定匹配滤波器的输入信号为正弦信号与白噪声的叠加，其频率,采样频率,延时时间,。由实验的结果，我们可以看出匹配滤波器的单位冲激响应幅度谱与输入信号的幅度谱一致，符合上面我们推导的公式；由滤波器的输出信号我们可以得出，在抽样时刻处，信号有最大值，在输出信号的幅度谱中，干扰噪声白噪声明显减小，信号的有用成分被保留下来，体现出了匹配滤波器抑制噪声，保证在抽样时刻有最大值的特性，从而进一步保证了输出信噪比在抽样时刻最大，为下一步抽样判决提供了有利的条件。由此可以看出匹配滤波器对数字基带通信系统的通信可靠性的提高有着很重要的作用。四、附录实验程序：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];00.511.52-10010t/sS(t)输入信号+噪声-5005000.51f(Hz)|S(W)|输入信号幅度谱00.511.52-10010t/sh(t)冲激响应-5005000.51f(Hz)|H1(W)|单位冲激响应幅度谱01234-100001000t/sSo(t)匹配滤波器输出信号-50050050100f(Hz)|So1(W)|输出信号幅度谱freq1=linspace(-fs/2,fs/2,length(t));freq2=linspace(-fs/2,fs/2,length(t1));Noise=randn(1,length(t));%白噪声signal=3*cos(2*pi*fc*t+pi*k*t.^2)+Noise;%信号h1=3*cos(2*pi*fc*(to-t)+pi*k*(to-t).^2)+Noise;%匹配滤波器So1=conv2(signal,h1);%信号卷积subplot(321)plot(t,signal);axis([0,2,-10,10]);xlabel('t/s');ylabel('S(t)');title('输入信号+噪声');subplot(322)plot(freq1,fftshift(abs(fft(signal))/length(t))),xlabel('f(Hz)'),ylabel('|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('f(Hz)'),ylabel('|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('f(Hz)'),ylabel('|So1(W)|'),title('输出信号幅度谱');