在信号与系统中的应用一总结概括.pptVIP

例7.19 用DFT求连续信号频谱 加大截取长度Tp可提高频率分辨率;选择合适的窗函数可降低谱间干扰;而频谱混叠失真要通过提高采样频率fs和(或)预滤波(在采样之前滤除折叠频率以外的频率成分)来改善。 编写程序q719.m验证截断效应及加窗的改善作用,先选取以下参数: ? 采样频率fs = 400Hz,T = 1/fs ? 采样信号序列 ? 对x(n)作4096点DFT作为的近似频谱Xa(jf)。 例7.19 用DFT求连续信号频谱 如图7.19所示。图中X1(jf),X4(jf)和X8(jf)分别表示Tp = 0.04s,0.04*4s和0.04*8s时的谱分析结果。由图可见,由于截断使原频谱中的单频谱线展宽(又称之为泄漏),Tp越大泄漏越小,频率分辨率越高。Tp = 0.04s时,25Hz与50Hz两根谱线已分辨不清了。所以实际谱分析的截取时间Tp是由频率分辨率决定的。 另外,在本应为零的频段上出现了一些参差不齐的小谱包(称为谱间干扰)。谱间干扰的大小取决于加窗的类型。用矩形窗比用Hamming窗的频率分辨率高(泄漏小),但谱间干扰刚好相反。 例7.20 IIR滤波器直接型的转换 程序调用了信号处理工具箱函数tf2sos和扩展函数dir2par, ? [sos, g] = tf2sos(B, A) 实现从直接型到级联型(二阶分割形式)的转换。g为式中的增益,sos为L×6阶矩阵,表示式中的系数。 例7.20 IIR滤波器直接型的转换 ? [Cp, Bp, Ap] = dir2par(B, A) 实现从直接型到并联型的转换。B为直接型H(z)的分子多项式系数向量,A为直接型H(z)的分母多项式系数向量;Cp,Bp,Ap的含义与扩展函数dir2par中的C,B,A相同。 dir2par中又调用了复共轭对比较函数cplxcomp。由于dir2par和cplxcomp是文献[7]中开发的,不属于MATLAB工具箱函数,所以将其M文件清单附在程序q720.m之后。 例7.20 IIR滤波器直接型的转换 根据计算结果, 级联型H(z)表达式: 级联型结构图。 例7.20 IIR滤波器直接型的转换 并联型结构H(z)表达式 并联型结构图 例7.21 直接型结构到格型梯形结构 ? tf2latc函数实现直接型到格型转换 ? [K, C] = tf2latc (B, A) 求出零-极点IIR系统格型梯形结构的格型参数向量K和梯形参数向量C(用A(1)归一化)。注意,当系统函数在单位圆上有极点时发生错误。 ? K = tf2latc (1, A) 求出全极点IIR系统的格型结构参数向量K。如果使用格式[K, C] = tf2latc (1, A),返回的系数C为标量。 ? K = tf2latc (B) 求出FIR系统的格型梯形结构参数(反射系数)向量K(用H(z)的常数项B(1)归一化)。 例7.21 直接型结构到格型梯形结构 直接型系统函数 转换为格型梯形结构 例7.22 FIR滤波器直接型到其他型 系统函数为 调用信号处理工具箱函数tf2sos和tf2latc,给变元A赋值1,B=[2,13/12,5/4,2/3]即可. 级联型结构系数 sos = 1.0000 0.5360 0 1.0000 0 0 1.0000 0.0057 0.6219 1.0000 0 0 g = 2 格型结构系数(反射系数): K = 0.2500 0.5000 0.3333 例7.22 FIR滤波器直接型到其他型 得出级联型为 格型结构为 例7.23 FIR格型到直接型转换 给定K = [2,1/4,1/2,1/3] ,用函数 latc2tf即可 由B=latc2tf(K)得到B,写出直接型结构 例7.24 系统函数的计算机推导 数字滤波器的网络结构图实际上也是一种信号流图。因此【例6.20】中介绍的方法和公式同样可以用来求离散域的数字滤波器的系统函数。不同的地方仅仅在于节点方程中出现了作为系数的符号变量z?1,它将出现在系数矩阵中。MATLAB是不能处理上标变量的,因此在程序中设q=z ?1,在计算完成后再人工地把结果中的q恢复为z?1。 例7.24的结构图与方程 例7.24的方程的矩阵形式 由此可以求出系统函数 例7.24解出的系统函数 程序运行的结果 如果加入一个激励x(n)?A?(n),则得出 7.5 FIR数字滤波器设计 滤波器的特性指标 用绝对值δ1,δ2表示; 用分贝Rp,Rs表示 (1)窗函数法设计FIR滤波器 先根据?c和N求