数字信号处理实验报告第四章.doc

数字信号处理实验四 专业：电子信息工程 学号：20101560063 姓名：任子翔 Q4.2使用修改后的程序P3.1，计算并画出当时的传输函数的因果线性时不变离散时间系统的频率响应。他表示那种类型的滤波器？ 解： %程序P3.1 clf; w=0:8*pi/511:pi; num=[0.15,0,-0.15]; den=[1,-0.5,0.7]; h=freqz(num,den,w); %plot the DTFT subplot(2,1,1); plot(w/pi,real(h));grid; title(H(e^{j\omega})的实部); xlabel(\omega/\pi); ylabel(振幅); subplot(2,1,2) plot(w/pi,imag(h));grid; title(H(e^{j\omega})的虚部); xlabel(\omega/\pi); ylabel(振幅); pause subplot(2,1,1); plot(w/pi,abs(h));grid; title(|H(e^{j\omega})|幅度谱); xlabel(\omega/\pi); ylabel(振幅); subplot(2,1,2) plot(w/pi,angle(h));grid; title(相位谱[H(e^{j\omega})]); xlabel(\omega/\pi); ylabel(以弧度为单位的相位谱);  由上图可看出这是一个低通滤波器。 Q4.3对下面的传输函数重做习题Q4.2：这两题中给出的两个滤波器之间的区别是什么？你将选择哪个滤波器来滤波，为什么？ 解： %程序P3.1 clf; w=0:8*pi/511:pi; num=[0.15,0,-0.15]; den=[0.7,-0.5,1]; h=freqz(num,den,w); %plot the DTFT subplot(2,1,1); plot(w/pi,real(h));grid; title(H(e^{j\omega})的实部); xlabel(\omega/\pi); ylabel(振幅); subplot(2,1,2) plot(w/pi,imag(h));grid; title(H(e^{j\omega})的虚部); xlabel(\omega/\pi); ylabel(振幅); pause subplot(2,1,1); plot(w/pi,abs(h));grid; title(|H(e^{j\omega})|幅度谱); xlabel(\omega/\pi); ylabel(振幅); subplot(2,1,2) plot(w/pi,angle(h));grid; title(相位谱[H(e^{j\omega})]); xlabel(\omega/\pi); ylabel(以弧度为单位的相位谱);  这题中的滤波器的幅度响应较Q4.2的更为尖锐， Q4.7用程序 P4. 1计算并画出近似理想低通滤波器的冲激响应。低通有限冲激响应滤波器的长度是多少？在程序 P4. l 中，哪个语句确定滤波器的长度？哪个参数控制截止 频率？ 解： %程序P4.1 %截短的理想低通滤波器 clf; fc=0.25; n=[-6.5:1:6.5] y=2*fc*sinc(2*fc*n);k=n+6.5 stem(k,y);title(N=13);axis([0 13 -0.2 0.6]); xlabel(时间序号n);ylabel(振幅);grid; n = Columns 1 through 5 -6.5000 -5.5000 -4.5000 -3.5000 -2.5000 Columns 6 through 10 -1.5000 -0.5000 0.5000 1.5000 2.5000 Columns 11 through 14 3.5000 4.5000 5.5000 6.5000 k = Columns 1 through 8 0 1 2 3 4 5 6 7 Columns 9 through 14 8 9 10 11 12 13 低通滤波器的长度为13，n=[-6.5:1:6.5]决定了滤波器的长度。fc=0.25;控制截止频率。 Q4.8修改程序 P4.l，计算并画出（4.39 所示长度为 20 ，截止角频率为 的有限 冲激响应低通滤波器的冲激响应。 解： %程序P4.1 %截短的理想低通滤波