资源描述:
《最优滤波器的设计《生物医学信号处理》实习报告4》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、《生物医学信号处理》实习报告学生姓名:学号:实验室名称:项目名称:最优滤波器的设计项目内容:编写函数解W-H方程,寻找最优的滤波器,并检验该程序的准确性和掌握用法。1)运行MATLAB,已知信号的自相关函数和噪声的能量,来编写函数解W-H方程,寻找最优滤波器。2)编写程序,仿真出真实信号,噪声,观察信号,然后把观察信号通过滤波器得到的输出信号与原来的真实信号进行比较,观察是否达到了去噪目的。原理(写出具体的计算公式)有限脉冲响应法求解维纳-霍夫方程:设h(n)是一个因果序列且可以用有限长(N点长)的序列去逼近它,最后得到:,j=0,1,2,…,N-1。简化形式为:。若信号s(n)与噪声w(
2、n)互不相关,即则有最后化为:,j=0,1,2,…,N-1。编写的源程序:function[H,Emin]=WH(Rs,Rw,N)L1=(length(Rs)+1)/2;Rss=zeros(1,L1);fork=1:L1Rss(k)=Rs(k+L1-1);%原信号的自相关endL2=(length(Rw)+1)/2;Rww=zeros(1,L2);fork=1:L2Rww(k)=Rw(k+L2-1);%噪声的自相关endRx=zeros(1,N);fork=1:NRx(k)=Rss(k)+Rww(k);endRxx=zeros(N,N);fori=1:Nforj=1:Nif(i<=j)Rx
3、x(i,j)=Rx(j-i+1);elseRxx(i,j)=Rx(i-j+1);endendendH=pinv(Rxx)*Rss';Emin=Rss(1)-sum(H*Rss);编写的源程序:clc;t=-50:0.1:50;w=0.1*pi;s=cos(w.*t);w=randn(1,length(t));x=s+w;Rs=xcorr(s);Rw=xcorr(w);[h,E]=WH(Rs,Rw,length(t));se=filter(h,1,x);subplot(4,1,1);plot(t,s);title('原序列');subplot(4,1,2);plot(t,x);title('
4、混入噪声的信号');subplot(4,1,3);plot(t,se);title('经维纳滤波后的信号');e=s-se;subplot(4,1,4);plot(t,e);title('误差');结论(画出要求的图形)总结有限脉冲响应法求解维纳-霍夫方程,设h(n)是一个因果序列且可以用有限长(N点长)的序列去逼近它,最后得到:,j=0,1,2,…,N-1。简化形式为:。若信号s(n)与噪声w(n)互不相关,即则有最后化为:,j=0,1,2,…,N-1。思考题:1.维纳滤波器在什么情况下能明显去噪?答:对于平稳随机过程,且已知信号和噪声的自相关函数,可通过求最优的维纳霍夫滤波器来滤波,有
5、较显著的效果。2.有时仿真去噪的效果不好,考虑原因。答:可能是数据量太少和滤波阶数过小,信号不是平稳随机过程,或信号和噪声的自相关函数未知。实习报告分数:指导教师: