《实验一 系统函数幅度相位谱的绘制.doc》由会员分享，可在线阅读，更多相关《实验一 系统函数幅度相位谱的绘制.doc（5页珍藏版）》请在金锄头文库上搜索。
1、数字信号处理实验一：系统函数幅度相位谱的绘制实验目的：通过编程绘制的幅度相位谱加深理解系统的特性。实验要求：1绘制教材习题1.59各小题系统函数的幅度谱和相位谱。（1）0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi0.30.71.11.51.92.32.7电子1002王江 幅度谱0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi-3.14
2、16-2.3562-1.5708-0.785400.78541.57082.35623.1416电子1002王江 相位谱（2）0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi0.30.71.11.51.92.32.7电子1002王江幅度谱0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi-3.1416-2.3562-1.5708-0.785
3、400.78541.57082.35623.1416电子1002王江相位谱（3）0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi1.51.92.32.7电子1002王江幅度谱0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi0电子1002王江 相位谱（4）0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi
4、/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi0.30.71.11.51.92.32.7电子1002王江 幅度谱0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi-1.5708-0.785400.78541.5708电子1002王江 相位谱（5）0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi0
5、.71.11.51.9电子1002王江 幅度谱0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi-0.785400.7854电子1002王江 相位谱2. 在ITU 关于语音压缩的标准G.723.1中，提到采样语音信号序列首先经过一个高通滤波器，该滤波的系统函数,由教材1.6.7先在纸上大致描出零极点图，看看是否是高通滤波器，并编写程序绘出幅度谱进行验证。0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/
6、83*pi/213*pi/87*pi/415*pi/82*pi0.30.71.11.5电子1001 吴鹏0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi-1.5708-0.785400.7854电子1001 吴鹏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的系统幅频
7、曲线。实验报告主要关注：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均是向量而且等长；也可以直接使用
8、plot(y)函数进行绘制，此时y向量的自然下标（1,2,3,3,4）作为横坐标。Matlab中可以用w=0:0.01:2*pi来离散计算在每一点的值，结果也是向量，则绘制可用plot(w,y)进行，这里的y表示是计算结果向量。2.可在纸上先求出公式后编程或直接书写，的语句表示为exp(j*w),求模用abs(x),abs对于负数是求绝对值，对于复数求模。相角用angle(hjw)函数。3.对于H(z),只要收敛域包含单位圆，可以直接将z换成，然后分析幅度相位谱.0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/21
9、3*pi/87*pi/415*pi/82*pi0.30.71.11.51.92.32.7?0pi/8pi/43*pi/8pi/25*pi/83*pi/47*pi/8pi9*pi/810*pi/811*pi/83*pi/213*pi/87*pi/415*pi/82*pi-3.1416-2.3562-1.5708-0.785400.78541.57082.35623.1416?电子1001吴鹏10109101034.要计算超过2*pi的区间,可写成w=0:0.01:4*pi; 4可以改为其他，注意要量力而行，w实际上在离散点上，过多的点会耗内存，导致出错！5当序列较长时，可采用如下语句按公式书写表达式：hjw=0;for di=1 to length(x) hjw=hjw+x(di)*exp(-j*w*(di-1);end然后用abs( hjw)函数求得幅值。
