韶关学院级数字信号处理.doc

第3、4页miss 三、1．已知线性因果网络用下面差分方程描述 y(n)=0.9y(n－1)+x(n)+0.9x(n－1) （1）求网络的系统函数H(z)及单位脉冲响应h(n)； （2）写出网络频率响应函数H(ejω)的表达式， 并定性画出其幅频特性曲线； （3）设输入x(n)=ejω0n， 求输出y(n)。 解： （1）y(n)=0.9y(n－1)+x(n)+0.9x(n－1) Y(z)=0.9Y(z)z－1+X(z)+0.9X(z)z－1 令 n≥1时，c内有极点0.9， n=0时，c内有极点0.9，0， 最后得到 h(n)=2·0.9nu(n－1)+δ(n) （2） 极点为z1=0.9，零点为z2=－0.9。极零点图如题24解图(a)所示。按照极零点图定性画出的幅度特性如题24解图(b)所示。 （3） 题24解图 2． 用微处理机对实数序列作谱分析， 要求谱分辨率F≤50 Hz， 信号最高频率为 1 kHz， 试确定以下各参数： (1) 最小记录时间Tp min； (2) 最大取样间隔Tmax； (3) 最少采样点数Nmin； (4) 在频带宽度不变的情况下， 使频率分辨率提高1倍（即F缩小一半）的N值。 解：（1）已知F=50 Hz，因而 （2） （3） （4） 频带宽度不变就意味着采样间隔T不变， 应该使记录时间扩大1倍， 即为0.04 s， 实现频率分辨率提高1倍（F变为原来的1/2）。 3．分别画出8点基2DIT-FFT和DIF-FFT运算流图，并计算其复数乘次数，。 4．设计低通数字滤波器，要求通带内频率低于0.2π rad时，容许幅度误差在1 dB之内；频率在0.3π到π之间的阻带衰减大于10 dB。试采用巴特沃斯型模拟滤波器进行设计，用脉冲响应不变法进行转换，采样间隔T=1 ms。 解: 本题要求用巴特沃斯型模拟滤波器设计，所以，由巴特沃斯滤波器的单调下降特性，数字滤波器指标描述如下： 　　　　ωp=0.2 π rad, αp=1 dB 　　　　ωs=0.3 π rad, αs=10 dB 采用脉冲响应不变法转换，所以，相应的模拟低通巴特沃斯滤波器指标为 (1) 求滤波器阶数N及归一化系统函数G(p)： 取N=5。查教材6.1节的表6.2.1(第157页)，可知模拟滤波器系统函数的归一化低通原型为 将G(p)部分分式展开： 其中，系数为 A0=－0.1382+j0.4253，A1=－0.8091－j1.1135，A2=1.8947，A3=－0.8091+j1.1135， A4=－0.1382－j0.4253 (2) 去归一化求得相应的模拟滤波器系统函数Ha(s)。 我们希望阻带指标刚好， 让通带指标留有富裕量， 所以按教材（6.2.20）式求3 dB截止频率Ωc为 其中, Bk=ΩcAk, sk=Ωcpk。 (3) 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z): 我们知道，脉冲响应不变法的主要缺点是存在的频率混叠失真，使设计的滤波器阻带指标变差。另外，由该题的设计过程可见。当N较大时，部分分式展开求解系数Ak或Bk相当困难，所以实际工作中用得很少，主要采用双线性变换法设计，见第10题。 10． 要求同题9，试采用双线性变换法设计数字低通滤波器 解：已知条件如下： 数字滤波器指标： ωp=0.2π rad, αp=1 dB ωs=0.3π rad, αs=10 dB 采用双线性变换法，所以要进行预畸变校正，确定相应的模拟滤波器指标（为了计算方便，取T=1 s）： (1) 求相应模拟滤波器阶数N： 其中， ksp与题9相同（因为αp、 αs相同）， 即 (2)查教材表6.2.1，得 (3) 去归一化，求出Ha(s)： (4) 用双线性变换法将Ha(s)转换成H(z)： 请读者按T=1 ms进行设计，比较设计结果。 四、1.设计一个工作于采样频率80 kHz的巴特沃斯低通数字滤波器，要求通带边界频率为4 kHz，通带最大衰减为0.5 dB，阻带边界频率为20 kHz，阻带最小衰减为45 dB。调用MATLAB工具箱函数buttord和butter设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。 题14*解图 解： 本题以模拟频率给定滤波器指标，所以，程序中先要计算出对应的数字边界频率，然后调用MATLAB工具箱函数buttord和butter来设计数字滤波器。设计程序为ex614.m。 　　% 程序ex614.m 　　Fs=80000； T=1/Fs； 　　wp=2*pi*4000/Fs； ws=2*pi*20000/Fs； rp=0.5； rs=45； 　　［N， wc］