试验方法与数据处理大作业.docx

大作业一FIR、IIR滤波器设计一、题目要求已知：X（t）= sin（2··30·t）+ sin（2··100·t）+ sin（2··300·t）要求：1、绘出给定信号波形2、用低通滤波器滤波：除去100Hz，300Hz信号，留下30Hz信号，注：先用FIR（加窗口和不加窗口），再用IIR。3、用带通滤波器滤波，除去100Hz，300Hz信号，留下30Hz信号，（先用FIR（加窗口和不加窗口），再用IIR）4、在原信号上加噪声：0.2·rand（1，N）注：1、用3-3,3-4的结论做2、提示：初始点有问题可忽略不管 FIR滤波器N试6-1060-100 IIR滤波器：Ｌ：6－1001、FIR滤波器，滤去100Hz、300Hz信号，留下30Hz信号（无修正窗口）：代码：clearclcsyms kfimdt=1/1000;%满足采样定理;x1=sin(2*pi*30*k*dt);x2=sin(2*pi*100*k*dt);x3=sin(2*pi*300*k*dt);x=x1+x2+x3;F=30;n1=-100;n2=100;t=0;for i=n1:1:n2;if i==0 q=2*F*dt else q=sin(2*pi*F*i*dt)/(pi*i) end p=sin(2*pi*30*(k-i)*dt)+sin(2*pi*100*(k-i)*dt)+sin(2*pi*300*(k-i)*dt) fi=symfun(p,k) t=t+fi*q;end yk=t; hold all; g1=ezplot(yk,[0,1000]); g2=ezplot(x1,[0,1000]); g3=ezplot(x,[0,1000]); axis([0 1000 -3 3]); set(g1, color, r); set(g2, Color, b); set(g3, Color, g); legend(过滤后波形,30H原波形,给定波形); title(波形图); xlabel(横坐标); ylabel(纵坐标);波形图：其中，绿色波形图为给定原信号，红色波形图为经过滤波器滤波后所得的波形，蓝色波形为原30Hz原波形，从图形中可知，当n1取-10，n2取10时，滤波效果不好，当n1取-100，n2取100时，滤波效果较好效果图如下：2、FIR滤波器，除去100Hz，300Hz信号，留下30Hz信号（加Hanning修正窗口）：代码：clearclcsyms kfimdt=1/1000;%满足采样定理;x1=sin(2*pi*30*k*dt);x2=sin(2*pi*100*k*dt);x3=sin(2*pi*300*k*dt);x=x1+x2+x3;F=30;n1=-100;n2=100;t=0;for i=n1:1:n2; w=1/2*(1+cos((pi*i)/10)) p=sin(2*pi*30*(k-i)*dt)+sin(2*pi*100*(k-i)*dt)+sin(2*pi*300*(k-i)*dt) q=sin(2*pi*F*i*dt)/(pi*i) fi=symfun(p,k); t=t+fi*q*w;end yk=t; hold all; g1=ezplot(yk,[0,1000]); g2=ezplot(x1,[0,1000]); g3=ezplot(x,[0,1000]); axis([0 1000 -3 3]); set(g1, color, r); set(g2, Color, b); set(g3, Color, g); legend(过滤后波形,30Hz波形,给定波形); title(波形图); xlabel(横坐标); ylabel(纵坐标);波形图：3、FIR滤波器，滤去100Hz，300Hz信号，留下30Hz信号（加Hamming修正窗口）：代码：clearclcsyms kfimdt=1/1000;%满足采样定理;x1=sin(2*pi*30*k*dt);x2=sin(2*pi*100*k*dt);x3=sin(2*pi*300*k*dt);x=x1+x2+x3;F=30;n1=-100;n2=100;t=0;for i=n1:1:n2; w=0.54+0.46*cos(pi*i/10) p=sin(2*pi*30*(k-i)*dt)+sin(2*pi*100*(k-i)*dt)+si