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

数字信号处理实验四 专业：电子信息工程 学号：座机电话号码063 姓名：任子翔Q4.2使用修改后的程序P3.1，计算并画出当时的传输函数的因果线性时不变离散时间系统的频率响应。他表示那种类型的滤波器？ 解： %程序P3.1 clf; w0:8*pi/511:pi; num[0.15,0,-0.15]; den[1,-0.5,0.7]; hfreqznum,den,w; %plot the DTFT subplot2,1,1; plotw/pi,realh;grid; titleHe^j\omega的实部; xlabel\omega/\pi; ylabel振幅; subplot2,1,2 plotw/pi,imagh;grid; titleHe^j\omega的虚部; xlabel\omega/\pi; ylabel振幅; pause subplot2,1,1; plotw/pi,absh;grid; title|He^j\omega|幅度谱; xlabel\omega/\pi; ylabel振幅; subplot2,1,2 plotw/pi,angleh;grid; title相位谱[He^j\omega]; xlabel\omega/\pi; ylabel以弧度为单位的相位谱; 由上图可看出这是一个低通滤波器。 Q4.3对下面的传输函数重做习题Q4.2：这两题中给出的两个滤波器之间的区别是什么？你将选择哪个滤波器来滤波，为什么？ 解： %程序P3.1 clf; w0:8*pi/511:pi; num[0.15,0,-0.15]; den[0.7,-0.5,1]; hfreqznum,den,w; %plot the DTFT subplot2,1,1; plotw/pi,realh;grid; titleHe^j\omega的实部; xlabel\omega/\pi; ylabel振幅; subplot2,1,2 plotw/pi,imagh;grid; titleHe^j\omega的虚部; xlabel\omega/\pi; ylabel振幅; pause subplot2,1,1; plotw/pi,absh;grid; title|He^j\omega|幅度谱; xlabel\omega/\pi; ylabel振幅; subplot2,1,2 plotw/pi,angleh;grid; title相位谱[He^j\omega]; xlabel\omega/\pi; ylabel以弧度为单位的相位谱;这题中的滤波器的幅度响应较Q4.2的更为尖锐，Q4.7用程序 P4. 1计算并画出近似理想低通滤波器的冲激响应。低通有限冲激响应滤波器的长度是多少？在程序 P4. l 中，哪个语句确定滤波器的长度？哪个参数控制截止 频率？ 解： %程序P4.1 %截短的理想低通滤波器 clf; fc0.25; n[-6.5:1:6.5] y2*fc*sinc2*fc*n;kn+6.5 stemk,y;titleN13;axis[0 13 -0.2 0.6]; xlabel时间序号n;ylabel振幅;grid; nColumns 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 143.5000 4.5000 5.5000 6.5000 kColumns 1 through 801234567 Columns 9 through 1489 10 11 12 13 低通滤波器的长度为13，n[-6.5:1:6.5]决定了滤波器的长度。fc0.25;控制截止频率。 Q4.8修改程序 P4.l，计算并画出（4.39 所示长度为 20 ，截止角频率为 的有限 冲激响应低通滤波器的冲激响应。 解： %程序P4.1 %截短的理想低通滤波器 clf; wc0.45; fcwc/2*pi; n[-9.5:1:9.5]; y2*fc*sinc2*fc*n;kn+9.5; stemk,y;titleN20;axis[0 20 -0.2 0.6]; xlabel时间序号n;ylabel振幅;grid;Q4.19 运行程序 P4. 3，生成每一 类线性相位有限冲激响应。 每一个有限 冲激响应滤波器的 长度是多少？验证冲激响应序列的对称性。接着验证这些滤波器 的零点位置。 使用 MATIAB 计算并绘出这些滤波器的相位响应 ，