实验四线性时不变离散时间系统的频域分析.docxVIP

Q4.1修改程序，取三个不同的M值，当时0=w=2*pi时计算并画出式(2.13)所示滑动平 均滤波器的幅度和相位谱。证明由幅度和相位谱表现出的对称类型。它表示了哪能种类型 的滤波器？你现 修改后的程序： w=0:pi/511:2*pi; M=input(M=); num=ones(l,M)/M; h=freqz (nu m,l,w); subplot(2,l,l); plot(w/pi,abs(h));grid; title(-H(eA{j\omega})幅度谱*); xlabel(、omega/\pi);ylabel(振幅)； subplot(乙 1,2); plot(w/piza ngle(h));grid; title(相位谱 H(eA{j\omega})); xlabel(\omega/\pij;ylabel(以弧度为单位的相位 当 M=5： 当 M“0： H(#)幅度谓相位谓H(0)080 02ags^a丸MKZ518 2H(皿度语16160s￡￡M^feJmD8 1 12o/i相位谱 H(#)幅度谓 相位谓H(0) 08 0 02 ags^a丸MKZ5 18 2 H(皿度语 16 16 0s￡￡M^feJm D8 1 12 o/i 相位谱H((鬥 08 当 M=20： H(护度谱0 02 04 06 08 1 1.2 H(护度谱 0 02 04 06 08 1 1.2 1 4 1.6 1 8 co/x 结论：由图可以看出，它表示了一个低通滤波器。 Q4.2使用修改后的程序，计算并画出当0=w=pi时给出的传输函数的因果线性时不变离散 时间第统的频率响应，它表示哪种类型的滤波器? 程序： w=0:pi/511:pi; num=[0.15 0 -0.15];den=[l-0.5 0.刀； h=freqz(num,de n,w);03 04 05 06 07 0.8 09o/i相位谱H(产)subplot(2,l,l)； plot(w/pi,abs(h));grid; titleCH(eA{j\omega})幅度谱); xlabel(、omega/\pi);ylabel(‘振幅); 程序： w=0:pi/511:pi; num=[0.15 0 -0.15]; den=[l-0.5 0.刀； h=freqz(num,de n,w); 03 04 05 06 07 0.8 09 o/i 相位谱H(产) subplot(2,l,l)； plot(w/pi,abs(h));grid; titleCH(eA{j\omega})幅度谱); xlabel(、omega/\pi);ylabel(‘振幅); subplot(2,l,2); plot(w/pi,a ngle(h));grid; title(湘位谱 H(eA{j\omega}))； xlabel(\omega/\p r); 01 02 03 04 05 06 07 08 09 a/x 1 ylabef以弧度为单位的相位，); 得右上图： 结论：由图可以得出它表示带通滤波器。 Q4.3对下面给出的传输函数求因果线性时不变离散时间系统的频率响应，它表法哪种类型 的滤波器？ 在上题程序中修改den=[0.7-0.5 1] 得图： 结论： 比较上题得出的图可知，幅度谱是 相同的，相位谱线，上题得出的是 连续的，而本题的相位谱线存在跳 变，我将选择上题那个滤波器来滤 波。 相位馆H(尸) Q4.7运行程序计算并画出近似理想 低通滤波器的冲激响应。 程序： elf; fc=0.25; n 二卜 6.5:165]; y=2*fc*sinc(2*fc*n);k=n+6.5; 二 ■ ■ Jo 2 0 2 4 ? stem(k,y);title(N=13,);axis([0 13 -0.2 0.6]); xlabel(时间序号 n);ylabel(J5幅];grid 低通有限冲激响应滤波器的长度是14, 程序中，n=[-6.5:l:6.5];ig句确定滤波器的长度, fc参数控制截止频率。 得图： Q4.8修改程序，计算并画出式(4.39)所示长度为20,截止角频率为wc=0.45的有限冲激 响应低通滤波器的冲激响应。 N?20M?15时间序号nM=2; ??滑动平均低通滤波器的增益响应程序 N?20 M?15 时间序号n M=2; ??滑动平均低通滤波器的增益响应程序 修改后程序： elf; fc=0.45/(2*pi); n=[-9.5:l:9.5]; y=2*fc*sinc(2*fc*n);k=n+9.5; stem(k,y);title(,N=20,);axis([0 20 -0.2 0.6]); xlabel(时间序号 rV);ylabelC振幅];grid 得右图： Q4.9