数字信号处理 中国通信学会普通高等教育“十二五”规划教材立项项目 作者 王艳芬 张晓光 刘卫东 王刚 第7章 数字信号处理MATLAB实现及综合应用1.pptVIP

第7章 数字信号处理MATLAB实现及综合应用 7.1 数字信号处理MATLAB实现 7.1.1 MATLAB简介 7.1.2 离散时间信号与系统分析应用举例 7.1.3 离散傅里叶变换应用举例 7.1.4 IIR滤波器设计常用函数及举例 7.1.5 FIR滤波器设计常用函数及举例 7.1.1 MATLAB简介 7.1.2 离散时间信号与系统分析应用举例 7.1.4 IIR滤波器设计常用函数及举例 7.2 在双音拨号系统中的应用 7.2.1 电话的双音拨号系统 7.2.2 Goertzel算法 7.2.3 检测DTMF信号的DFT参数选择 7.2.3 检测DTMF信号的DFT参数选择 7.3 在wav信号分析方面的应用 7.3.1 wav文件的一次性傅里叶 变换 7.3.2 wav文件的分段傅里叶分析 7.4　在自适应滤波器方面的应用 7.4.1　自适应滤波原理 7.4.2 系统辨识与MATLAB仿真 7.4.3 噪声抵消与MATLAB仿真 7.4.2 系统辨识与MATLAB仿真 sound（w，fs，bits）——和前面的waveread一样的参数表示，它将序列的数据通过声卡转化为声音。 下面先介绍wav文件的一次性傅里叶变换。 7.3.1 wav文件的一次性傅里叶变换 1．声波主要频率的分析 图7-23 信号的时域波形 图7-24 信号频域的幅值 2．声音信号的重构 前面介绍了对wav文件直接进行一次性傅里叶变换的分析，但存在两个问题：第一，对分析的wav文件直接一次性做FFT，点数太大，而在实际运用中DSP的FFT的点数是有限的，一般只能达到千点；第二，语音是分音节的，应该对它分段分析，所以下面引入分段处理数据，并连接每个段的频域结果。 7.3.2 wav文件的分段傅里叶分析 图7-25 分段后频域幅值变化图 图7-26 合成波形及线性差值平滑波形 图7-27 采用傅里叶反变换得到的波形 （1）巴特沃思滤波器：[b,a] = butter(n,wn,ftype,s) 其中，n为滤波器阶数；wn为滤波器截止频率；s为模拟滤波器，缺省时为数字滤波器。 （2）切比雪夫滤波器：[b,a] = cheby1(n,Rp,wn,ftype,s)或[b,a] = cheby2(n,Rs,wn,ftype,s) （3）椭圆滤波器：[b,a] = ellip(n,Rp,Rs,wn,ftype,s) [bz,az]=impinvar/ bilinear (b,a,fs) 5．脉冲响应和双线性变换函数 其中，b、a分别为模拟滤波器的分子和分母多项式系数向量；fs为采样频率，单位为Hz，缺省时fs默认为1Hz；bz，az分别为数字滤波器分子和分母多项式系数向量。 图7-10 不同采样频率对所设计数字滤波器频率响应的影响 图7-11 ChebyshevⅠ型数字高通滤波器幅频特性 图7-12 7阶ChebyshevⅡ型数字低通滤波器幅度和相位图 7.1.5 FIR滤波器设计常用函数及举例 1．窗函数法设计FIR滤波器 （2）计算截止频率为wc的理想低通滤波器的单位脉冲响应hd(n)，语句如下 hd=ideal_lp(wc,N) （3）h(n)为h=hd.* win。 （4）已知h(n)，为了得到频率响应的幅度、相位及群时延，需用语句如下 [db,mag,pha,grd,w]=freqz_m(h,1) 图7-13 设计的实际低通滤波器结果 （1）h(n)为h=real(ifft(H,N))。 （2）MATLAB程序中，可以用freqz计算频率响应，但不能从中确定振幅响应。 这是因为在MATLAB函数中，不存在与abs函数功能相当而能找到振幅的函数。 但在MATLAB中，提供了4种扩展函数，可以计算FIR 4种类型滤波器的振幅响应。 ● Hr_type1 ● Hr_type2 ● Hr_type3 ● Hr_type4 2．频率采样法设计FIR滤波器 前面介绍的窗函数法设计FIR滤波器，完全是根据窗函数法的设计思路和设计流图一步一步进行的，这有助于掌握方法的原理。 在实际应用中，还可以采用另外一种方法——直接利用MATLAB中的内部函数来进行FIR滤波器的设计。 MATLAB信号处理工具箱提供了两个用窗函数法设计FIR滤波器的函数，见表7-2。 值得一提的是，FIR滤波器的窗函数设计法不仅可以设计标准频率响应（低通、高通、带通及带阻）滤波器，还可设计任意频率响