试验一 系统函数幅度相位谱的绘制[精彩]

试验一 系统函数幅度相位谱的绘制[精彩] 数字信号处理实验一:系统函数幅度相位谱的绘制 jw实验目的:通过编程绘制的幅度相位谱加深理解系统的特性。H(e) jwjw实验要求:1绘制教材习题1.59各小题系统函数的幅度谱和相位谱。H(e)H(e) (1) 电子1002王江 幅度谱 2.7 2.3 1.9 1.5 1.1 0.7 0.3 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 电子1002王江 相位谱 3.1416 2.3562 1.5708 0.7854 0 -0.7854 -1.5708 -2.3562 -3.1416 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi (2) 电子1002王江幅度谱 2.7 2.3 1.9 1.5 1.1 0.7 0.3 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/ 15*pi/84 2*pi 电子1002王江相位谱 3.1416 2.3562 1.5708 0.7854 0 -0.7854 -1.5708 -2.3562 -3.1416 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi (3) 电子1002王江幅度谱 2.7 2.3 1.9 1.5 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 电子1002王江 相位谱 0 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi (4) 电子1002王江 幅度谱 2.7 2.3 1.9 1.5 1.1 0.7 0.3 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 电子1002王江 相位谱 1.5708 0.7854 0 -0.7854 -1.5708 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi (5) 电子1002王江 幅度谱 1.9 1.5 1.1 0.7 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 电子1002王江 相位谱 0.7854 0 -0.7854 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 2. 在ITU 关于语音压缩的标准G.723.1中，提到采样语音信号序列首先经过一个高通滤波 ,11z,器，该滤波的系统函数,由教材1.6.7先在纸上大致描出零极点图，看看是Hz,，，127,11z,128 否是高通滤波器，并编写程序绘出幅度谱进行验证。 电子1001 吴鹏 1.5 1.1 0.7 0.3 0 pi/8 pi/4 3*pi/8 pi/2 5*pi/83*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 电子1001 吴鹏 0.7854 0 -0.7854 -1.5708 0 pi/8 pi/4 3*pi/8 pi/2 5*pi/83*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 3.系统函数H(z)表明了系统特性，系统的单位脉冲响应h(n)同样也可以表示一个系统的特性，按照公式(1-106)编写求系统单位脉冲响应为h(0)=1,h(1)=1,h(2)=1,h(3)=1,h(4)=1,h(n)=0,n不等于0,1,2,3,4的系统幅频曲线。 实验 报告 软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载 主要关注: 1.在编程时，如果让w的变化范围超过[0，2,]，绘制出来的图形是否和[0，2,]内的周期重复?如果重复，请说明原因。 2.报告中粗略绘制出要求2的幅度谱，假设语音信号采样频率为8000Hz,那该滤波器主要滤除哪些频率信号。 3.实验要求3绘制出的曲线是高通/低通/带通/带阻中的哪一类滤波器，给出简要图形说明。 4.对于如何编程表示出H(jw)表达式，有哪些方式，请举例说明。 提示: 1.C，Matlab或其他语言均能实现，绘图方面Matlab可能更快捷方便。绘图函数plot(x,y), 该函数在二维平面坐标系中绘出以x为自变量，y为因变量的曲线。要求x,y均是向量而且等长;也可以直接使用plot(y)函数进行绘制，此时y向量的自然下标(1,2,3,3,4…)作为横 jw坐标。Matlab中可以用w=0:0.01:2*pi来离散计算在每一点的值，结果也是向量，H(e) 则绘制可用plot(w,y)进行，这里的y表示是计算结果向量。 jw2.可在纸上先求出公式后编程或直接书写，的语句表示为exp(j*w),求模用abs(x),abse 对于负数是求绝对值，对于复数求模。相角用angle(hjw)函数。 jw3.对于H(z),只要收敛域包含单位圆，可以直接将z换成，然后分析幅度相位谱.e ??? 2.7 2.3 1.9 1.5 1.1 0.7 0.3 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi ??? 3.1416 2.3562 1.5708 0.7854 0 -0.7854 -1.5708 -2.3562 -3.1416 0 pi/8 pi/43*pi/8 pi/2 5*pi/8 3*pi/4 7*pi/8 pi 9*pi/810*pi/8 11*pi/8 3*pi/2 13*pi/8 7*pi/4 15*pi/8 2*pi 电子1001吴鹏1010910103 4.要计算超过2*pi的区间,可写成 w=0:0.01:4*pi; 4可以改为其他，注意要量力而行，w实际上在离散点上，过多的点会耗内存，导致出错～ jw5当序列较长时，可采用如下语句按公式书写表达式: H(e) hjw=0; for di=1 to length(x) hjw=hjw+x(di)*exp(-j*w*(di-1)); end 然后用abs( hjw)函数求得幅值。