第十个.txt

来自「计算机端口扫描端口程序应用和开发的作用,供参考.」· 文本 代码 · 共 47 行

FIR:
xn=[4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];
N=21;
wn=0.2*pi/pi;
hn=fir1(N-1,wn);
yn=conv(xn,hn);
subplot(1,1,1);
plot(yn);


IIR:
xn=[4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];
N=21;
a1=[1,0.4164,-0.6006];
a2=[1,0.3093,-0.1888];
a3=[1,0.2692,-0.0349];
b=[0.2318,0.4636,0.2318];
yn1=filter(b,a1,xn);
yn2=filter(b,a2,yn1);
yn3=filter(b,a3,yn2);
subplot(1,1,1);
plot(yn3); 


比较：
xn=[4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];
N=21;
wn=0.2*pi/pi;
hn=fir1(N-1,wn);
yn=conv(xn,hn);
subplot(3,1,1);
plot(xn);
subplot(3,1,2);
plot(yn);
a1=[1,0.4164,-0.6006];
a2=[1,0.3093,-0.1888];
a3=[1,0.2692,-0.0349];
b=[0.2318,0.4636,0.2318];
yn1=filter(b,a1,xn);
yn2=filter(b,a2,yn1);
yn3=filter(b,a3,yn2);
subplot(3,1,3);
plot(yn3);