哈工大实验方法第一次大作业.pdf

《哈工大实验方法第一次大作业.pdf》由会员分享，可在线阅读，更多相关《哈工大实验方法第一次大作业.pdf（23页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、实验方法及数字信号处理分析第一次大作业哈尔滨工业大学哈尔滨工业大学试验方法及数字信号处理分析试验方法及数字信号处理分析第一次大作业第一次大作业数字滤波器设计指导老师：学生姓名：学生学号：15S1实验方法及数字信号处理分析第一次大作业第一部分第一部分 作业题目作业题目一、设计题目一、设计题目1、杂波信号：x(t)sin(210t)sin(280t)sin(2200t)2、要求：（1）绘出杂波信号波形。（2）分别用 FIR IIR 滤波器设计低通和带通滤波器，保留 10Hz，80Hz频率。绘出滤波后波形，并与理想波形比较。（3）在原信号加上白噪声信号，再比较分析。第二部分第二部分 具体设计内容具体

2、设计内容第一节第一节 卷积滤波器的设计卷积滤波器的设计一、低通滤波一、低通滤波1、低通滤波器参数计算（1）FIR 滤波频率响应：H(f)iN1N2fie j2fit1;0 f F（2）低通期望频率响应：H(f)（3）通过、计算滤波因子当i 0时，f 2Ft当i 0时，f 取f f sin(2Fi t)i0;f 0,f F可得近似理想低通滤波器：ykN2fixkiiN1（4）由于题目x（t）的最高频率fmax=200。基于采样定理，2实验方法及数字信号处理分析第一次大作业f2fmax=400。本例取 f=5fmax=1000。故 t=0.001s。2、设计程序程序参数：t1t；F=低通截止频率；

3、t=0:10-3:0.5;t1=10-3;F=20;x=sin(2*pi*10*t)+sin(2*pi*80*t)+sin(2*pi*200*t);x1=sin(2*pi*10*t);f(1)=2*F*t1;i=2;while i2fmax=400。本例取 f=5fmax=1000。故 t=0.001s。2、设计程序t=0:10-3:0.5;t1=10-3;F1=50;F2=110;F0=(F1+F2)/2;Ft=(F2-F1)/2;x=sin(2*pi*10*t)+sin(2*pi*80*t)+sin(2*pi*200*t);x1=sin(2*pi*80*t);f(1)=4*Ft*t1;i=

4、2;while i2fmax=400。本例取 f=5fmax=1000。故 t=0.001s。2、设计程序t=0:10-3:0.5;t1=10-3;F=25;x=sin(2*pi*10*t)+sin(2*pi*80*t)+sin(2*pi*200*t);x1=sin(2*pi*10*t);w=tan(pi*F*t1);f(1)=w2/(1+sqrt(2)*w+w2);f(2)=2*w2/(1+sqrt(2)*w+w2);f(3)=f(1);g(1)=-2*(1-w2)/(1+sqrt(2)*w+w2);1 3实验方法及数字信号处理分析第一次大作业1 4g(2)=(1-sqrt(2)*w+w2)

5、/(1+sqrt(2)*w+w2);y(2)=0;for k=3:500y(k)=f(1)*x(k)+f(2)*x(k-1)+f(3)*x(k-2)-g(1)*y(k-1)-g(2)*y(k-2);endy(length(t)=0;m(2)=0;for k=3:500m(k)=f(1)*y(k)+f(2)*y(k-1)+f(3)*y(k-2)-g(1)*m(k-1)-g(2)*m(k-2);endm(length(t)=0;n(2)=0;for k=3:500n(k)=f(1)*m(k)+f(2)*m(k-1)+f(3)*m(k-2)-g(1)*n(k-1)-g(2)*n(k-2);endn(

6、length(t)=0;plot(t,x,r,t,x1,b,t,y,y,t,m,k,t,n,g);legend(原图,理想图,2 阶滤波图,4 阶滤波图,6 阶滤波图);title(IIF 低通 F=25);xlabel(x);ylabel(y);实验方法及数字信号处理分析第一次大作业3、结果分析F=15放大图1 5实验方法及数字信号处理分析第一次大作业F=20放大图1 6实验方法及数字信号处理分析第一次大作业分析：上图展示了用 IIR 滤波在 2 阶、4 阶、六阶不同情况下的对比，同时也比较了 F 取不同值时的滤波情况。综合分析可发现，在其他参数不变的情况下，F=15 时，其幅值比理想波形小

7、，而 F=20 时，幅值基本可以吻合。而在 F 不变的情况下，滤波的阶数越高，波形的光滑程度越好。但同时波形有迟滞现象，阶数越大，迟滞现象越严重。这是 IIR 低通滤波器与 FIR 滤波器最大的区别。二、带通滤波二、带通滤波1、带通滤波器参数计算（1）IIR 滤波频率响应：H(f)Nfie j2fi t1g1e j2fii1i0Lt（2）带通期望频率响应：H(f)1;F1 f F20;the,others（3）基于巴特沃斯函数设计的滤波器巴特沃斯函数：Bn(s)n(ss1)(ss2)(ss3)1 7实验方法及数字信号处理分析第一次大作业（5）带通 IIR 滤波器系数z22z 1映射满足要求s

8、T(z)z21H(Z)B2(T(z)1()22()1ss由、可得对应系数如下：2f0212f1 022f2 212f3 02f4 f0122g1(42 2)2124222g2122g32 24212122g4212可得近似理想低通滤波器：yk f0 xk f1xk1 f2xk2 f3xk3 f4xk4 g1yk1 g2yk2 g3yk3 g4yk4（6）由于题目x（t）的最高频率fmax=200。基于采样定理，f2fmax=400。本例取 f=5fmax=1000。故 t=0.001s。1 8实验方法及数字信号处理分析第一次大作业1、程序设计clear;t=0:10-3:0.5;t1=10-3

9、;F1=60;F2=100;x=sin(2*pi*10*t)+sin(2*pi*80*t)+sin(2*pi*200*t);x1=sin(2*pi*80*t);w=tan(pi*(F2-F1)*t1);b=cos(pi*(F2+F1)*t1)/cos(pi*(F2-F1)*t1);f(1)=w2/(1+sqrt(2)*w+w2);f(2)=0;f(3)=-2*w2/(1+sqrt(2)*w+w2);f(4)=0;f(5)=f(1);g(1)=-(4*b+2*sqrt(2)*w*b)/(1+sqrt(2)*w+w2);g(2)=(4*b2+2-2*w2)/(1+sqrt(2)*w+w2);g(3

10、)=(2*sqrt(2)*w*b-4*b)/(1+sqrt(2)*w+w2);g(4)=(1-sqrt(2)*w+w2)/(1+sqrt(2)*w+w2);y(4)=0;for k=5:500y(k)=f(1)*x(k)+f(2)*x(k-1)+f(3)*x(k-2)+f(4)*x(k-3)+f(5)*x(k-4)-g(1)*y(k-1)-g(2)*y(k-2)-g(3)*y(k-3)-g(4)*y(k-4);endy(length(t)=0;m(4)=0;for k=5:500m(k)=f(1)*y(k)+f(2)*y(k-1)+f(3)*y(k-2)+f(4)*y(k-3)+f(5)*y(

11、k-4)-g(1)*m(k-1)-g(2)*m(k-2)-g(3)*m(k-3)-g(4)*m(k-4);endm(length(t)=0;plot(t,x,r,t,x1,b,t,y,y,t,m,k);%plot(t,x,r,t,x1,b,t,y,g);legend(原图,理想图,4 阶滤波图,8 阶滤波图);title(IIF 带通 F=25);xlabel(x);ylabel(y);1 9实验方法及数字信号处理分析第一次大作业3、实验分析放大图2 0实验方法及数字信号处理分析第一次大作业分析：由于之前低通分析，FIR 带通滤波器的 F 对波形影响不大，所以本处没有罗列不同 F 的波形情况。

12、同时计算出来的带通是 4 阶，所以进一步做了 8 阶的实验。分析可知，带通同样有迟滞现象，且阶数越高，迟滞越严重，但是波形也越光滑。总体光滑程度比 FIR 效果好。三、加入白噪声后，图形结果三、加入白噪声后，图形结果1、原理：用 0.5*rand（1）函数，添加一维的随机数，放在 x（t）后面，产生噪声。2、程序：t=0:10-3:0.5;t1=10-3;F=25;x=sin(2*pi*10*t)+sin(2*pi*80*t)+sin(2*pi*200*t)+0.5*rand(1);x1=sin(2*pi*10*t);w=tan(pi*F*t1);f(1)=w2/(1+sqrt(2)*w+w2

13、);f(2)=2*w2/(1+sqrt(2)*w+w2);f(3)=f(1);g(1)=-2*(1-w2)/(1+sqrt(2)*w+w2);g(2)=(1-sqrt(2)*w+w2)/(1+sqrt(2)*w+w2);y(2)=0;for k=3:500y(k)=f(1)*x(k)+f(2)*x(k-1)+f(3)*x(k-2)-g(1)*y(k-1)-g(2)*y(k-2);endy(length(t)=0;m(2)=0;for k=3:500m(k)=f(1)*y(k)+f(2)*y(k-1)+f(3)*y(k-2)-g(1)*m(k-1)-g(2)*m(k-2)2 1实验方法及数字信号

14、处理分析第一次大作业;endm(length(t)=0;n(2)=0;for k=3:500n(k)=f(1)*m(k)+f(2)*m(k-1)+f(3)*m(k-2)-g(1)*n(k-1)-g(2)*n(k-2);endn(length(t)=0;plot(t,x,r,t,x1,b,t,y,y,t,m,k,t,n,g);legend(原图,理想图,2 阶滤波图,4 阶滤波图,6 阶滤波图);title(含有白噪声的 IIF 低通);xlabel(x);ylabel(y);3、结果分析放大图2 2实验方法及数字信号处理分析第一次大作业可见加入白噪声后，图形无明显变化，该滤波器可以滤出白噪声。

15、第三节第三节 试验总结试验总结本次实验总体比较成功，滤波效果较好。总结的规律如下。1.比较发现，低通滤波器对 F 比带通敏感。2.感官上看，IIF 的滤波效果比同阶的 RIF 好得多，但是带来的问题是，IIF 滤波有波形迟滞现象，且其阶数越高，迟滞现象越严重。实际运用时要综合考虑。3.RIF 没有迟滞现象，IIF 有迟滞现象，其程度与阶数有关。4.截止频率 F 与通波频率不能太接近，过于接近，滤波幅值会减小。对于白噪声，滤波器可以很容易去除。实验使我加深了对滤波器设计的理解，也大概了解了滤波器的设计过程，有助于我以后对滤波器的学习和认识。实验进行过程也使我对 MATLAB 编程又熟悉了一般，编程能力得到进一步提高。总的来说，这次大作业对我的学习帮助很大，期望第二次的作业能让我有更多的收获。2 3