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

第7章 数字信号处理MATLAB实现及综合应用 ; 随着信息学科和计算机学科的飞速发展，数字信号处理的重要性日益显著，数字信号的处理已成为电子信息学科中的工程技术人员必不可少的知识。 由于数字信号处理的概念比较抽象，再加上其数值计算又比较繁琐，因此，国外很早就开始把MATLAB用于数字信号处理的教学过程，并取得了很好的经验，现在MATLAB已经成为解决数字信号处理问题的公认的标准软件。; 为了使读者更好地掌握数字信号处理的基本概念和基本理论，本章首先在7.1节给出前面各章的MATLAB典型例题。 通过这些例题，帮助读者弄清一些重要的概念，加深对课程基本概念和基本原理的理解，学会利用计算机软件解决在理论学习中不易解决的问题。 然后在7.2节~7.4节分别给出了数字信号处理技术在通信、信号分析、自适应滤波等方面的实际应用，并给出了MATLAB仿真程序和??果分析，这些内容可以扩宽读者数字信号处理技术应用的知识面，促进理论与实际相结合。; 本章和下一章上机实验都是课程学习的重要实践环节，它不仅能帮助同学们理解和掌握课程的重要概念和基本内容，而且也有益于读者深入学习和应用信号处理知识来解决实际问题。;7.1 数字信号处理MATLAB实现 ;7.1.1 MATLAB简介; MATLAB采用开放式的环境，你可以读到它的原码、了解它的算法，并能改变当前的函数或增添自己编写的函数。 目前，在国内外大学和研究机构中，MATLAB是一种非常流行的计算机语言，许多重要的学术刊物上发表的论文均是用MATLAB来分析计算和绘制各种图形。 它还是一种有力的教学工具，在大学的线性代数、自动控制理论、数字信号处理、动态系统仿真等课程的教学中，已成为标准的教学工具。;1．变量 2．基本操作符 3．输入、输出语句 4．库函数 5．绘图操作; 实际的图形显示有多种情况，下面分别加以说明。 （1）需要用多个图形窗显示多个图形时，应在图形显示语句之前加上图形窗指令figure； （2）在一个图形窗中显示多个图形时，应在图形显示语句之前调用图形窗分块函数subplot( )； （3）在一个图形坐标系中显示多条曲线时，要在plot( )语句前先执行指令hold on用以保持当前曲线，显示完成之后，再执行指令hold off； （4）需要选择坐标轴尺寸时，在plot( )语句之后调用axis( )函数。; 信号处理工具箱覆盖了经典信号处理的大多数内容，是一个非常优秀的算法研究、辅助设计工具。 它提供了大量的工具函数、滤波器分析和设计工具，常用的函数有以下几种。 （1）设计IIR滤波器相关函数 buttord、butter、cheblord、cheb2ord、ellipord、cheby1、cheby2等。;（2）设计FIR滤波器相关函数 fir1、fir2、firls、remezord、remez等； 窗函数：boxcar、hanning、hamming、blackman、kaiser。 （3）傅里叶变换相关函数 FFT（一维快速傅里叶变换）、IFFT（一维快速傅里叶逆变换）、FFT2（二维快速傅里叶变换）、abs（绝对值或模）、angle（相角）等。;（4）频率变换相关函数 freqz、freqs、lp2lp、lp2hp、lp2bp、lp2bs。 （5）频率响应相关函数 freqz、freqs。;7.1.2 离散时间信号与系统分析应用举例;图7-1 零极点图;图7-2 频率响应图;;图7-3 仿真波形;7.1.3 离散傅里叶变换应用举例;图7-4 4点DTFT和DFT图;图7-5 DFT图（N=8，N=16）;;图7-6 x(n)的前10点数据对应的x(n)、X(ejw)、X(k) ;图7-7 x(n)补零至100点对应的x(n)、X(ejw)、X(k) ;;图7-9 128点x(n)对应的x(n)、X(ejw)、X(k);7.1.4 IIR滤波器设计常用函数及举例; （1）低通到低通的频率变换 [b1,a1]=lp2lp(b,a,w0)，其中，w0为低通滤波器的截止频率（rad/s）。 （2）低通到高通的频率变换 [b1,a1]=lp2hp(b,a,w0)，其中，w0为高通滤波器的截止频率（rad/s）。; （3）低通到带通的频率变换 [b1,a1]=lp2bp(b,a,w0,Bw) ，其中，w0为带通滤波器的中心频率，Bw为带通滤波器的带宽。 当滤波器通带的下截止频率为w1，上截止频率为w2时，w0=sqrt(w1*w2)，Bw=w2-w1。