匹配滤波器的实验

1、2010年 秋季学期研穷生课程考核（阅读报告、研究报告）考核科目：科学技术哲学学生所在院（系）：电气工程及自动化学院学生所在学科：仪器科学与技术学生姓名：李海洋学 号 ：10S001049学生类别：工学硕士考核结果阅卷人匹配滤波器的设计与验证实验报告实验目的：1、了解匹配滤波器的基本原理；2、掌握如何设计一个传输系统的匹配滤波器；3、深刻认识匹配滤波器的一些实际应用；实验原理：设线性滤波器输入端输入的信号与噪声的混合波形为并假定噪声为白噪声，其功率谱密度 月=%（2 ,而信号s（t）的频谱函数为3（叫 即式今a E。我们要求线性滤波器在某时刻 上。上有最大的信号瞬时功率与噪声平均功率 的比值。

2、现在就来确定在上述最大输出信噪比准则下的最佳线性滤波器的传输特性。/=总畀巳一这就是最佳线性滤波器的传输特性。式中，（田）即为sg）的复共轲。在白噪声干扰的背景下，按式 （8.7-3）设计的线性滤波器，将能在给定时刻上获得最大的输出信噪比2E/。这种滤波器就是最大信噪比意义下的最佳线性滤波器。由于它的传输 特性与信号频谱的复共轲相一致 （除相乘因子 处7咻外）,故又称其为匹配滤波器。'da NCdr匹配滤波器的传输特性 沃3），当然还可用它的冲激响应心）来表示，这时有：3日3 = P KS* 1心，一对0 dsit3g由此可见，匹配滤波器的冲激响应便是信号为了获得物理可实现的匹配滤波器

3、，A1）=0。为了满足这个条件，就要满足:式力的镜像信号式T）在时间上再平移如。营8t）= q,3要求当£口时有有=Oh t >这个条件表明，物理可实现的匹配滤波器，其输入端的信号日团必须在它输出最大信噪 比的时刻为之前消失（等于零）。这就是说，若输入信号在 1瞬间消失，则只有当今 -'1时 滤波器才是物理可实现的。一般总是希望"尽量小些，故通常选择顺便指出，当我们专门关心匹配滤波器的。（）（ ） （ ）®（）（ 。 ） 输出信号波形时，它可表示为.J' .1.1 ' . - 1 ' '由此可见，匹配滤波器的输出信号波

4、形式输入信号的自相关函数的K倍。至于常数正，实际上它是可以任意选取的，因为“与反无关。因此，在分析问题时，可令 太二1。实验过程1 .产生1000点的白噪声nt,所用命令nt=randn（1,1000）依口图一）2 .产生1000点的有用信号st, st的角频率是8000pi,相位是时间的函数0.5*k*t.*t ,幅度是1的余弦函数。具体信号形式st=a*cos（w*t+0.5*k*t*t）;（如图二）3.产生和信号xt=st+nt （如图三）白噪声图图三* .4 .根据公式h（* =hm（T） =cu （t0 -'） ,求出匹配滤波器的数学表示形式。由于表达 式中有共腕计算，而信号

5、的时域函数式无法进行共腕表示，所以要把信号变道频 域中，才能出现所要的共腕形式，所以现对信号进行 fourier变换，对变换的谱取 共腕，再进行fourier反变换，再把反变换的信号即已经取完时域共腕的信号进行 时间轴的平移，即可得到匹配滤波器htlll的时域数学表示形式。这里用到 MATLAB指令有fft （离散fourier变换），conj （共腕计算），ifft （离散fourier反变5 .用conv （卷积）命令，将xt与htlll进行卷积，看输出结果。谱传递函数输出函数谱至此实验完成。(程序以附录的形式出现)实验结论输入混合信号经过匹配滤波器后，输出信号的确在t=0处出现了最大值，

6、且其它点的函数值远小于t=0时刻的值，这便证明了匹配滤波器能给出最大输出信 噪比，最有利于信号的检测。实验讨论本实验采用的有用信号是st=a*cos(w*t+0.5*k*t.*t),这是一个相位随时间连 续变化的余弦信号，若实验采用一个相位不随时间改变的正弦或余弦信号，则会 出现下列情况。以下的实验步骤与上面所述完全相同，故只呈现最后结果。4伯万'Ifl 1 1 f Aft I0.B11,'1 1III '0.61 10.4-02 -I -0.2 -0.4 -06 -0.8 -1 I,1). U v . V.I y. vy. I-1 1 1 ”-2-1.5-1-0.50

7、0.5t(s)11.52X10-3带噪声信号35 iiiiiiii-2-1.5-1-0.500.511.52陶x江-6 x 10-可见该输出信号也在t=0时刻出现了最大值即最大输出信噪比， 所以证明了匹配 滤波器的检测能力与输入信号波形无关， 只与信号能量有关；或者说，在同样的 白噪声干扰条件下，只要信号能量相同，并实现匹配滤波，则任何信号形式都能 给出相同的检测能力。 附录(实验源程序) %随即白噪声(normal) nt=randn(1,1000);白噪声);figure,plot(nt),title(xlabel('n');ylabel('x(n)');%

8、s(t)w=8000*pi;k=10000000;a=1;t=linspace(-0.002,0.002,1000); st=a*cos(w*t);信号)；figure,plot(t,st),title(xlabel('t(s)'ylabel('s(t)'%x(t) xt=nt+st;figure,plot(t,xt),title( xlabel('t(s)'ylabel('x(t)'%y(t););););带噪信号);stt1=fft(st); stt2=conj(stt1); st2=ifft(stt2);N=size(st2,2); for n=1:N;ht(n+1000)=st2(n); end ;ht111=ht(1,1001:2000);yt=conv(ht111,xt);%?i ?yt1=0.5*10A(-5)*(-999:999)/1000;figure,plot(t1,yt),title( %H(W) HW=fft(ht111);f=10A6*(1:1000)/1000;figure,plot(f,HW),title( xlabel('f(Hz)');ylabel('H(W)');%Y(W)YW=fft(yt);f=10A6*(-999:999)/