MATLAB信号处理

(1) 资深作者编著，图书质量更有保证：一线资深工程师执笔，系统归纳和总结了智能算法的实战经验。(2) 提供配套源码，便于读者动手实践：理论必须联系实践，本书提供源代码下载，方便读者学习使用。(3) 内含丰富实例，利于读者二次开发：书中提供了十几个信号处理算法的典型实例，读者可以据此二次开发。配套学习资源下载地址为清华大学出版社网站本书页面。系列畅销图书如下：MATLAB 编程指南——计算、编程、仿真、算法及应用MATLAB/Simulink系统仿真MATLAB GUI程序设计MATLAB智能算法MATLAB数学建模MATLAB科学计算MATLAB信号处理MATLAB图像处理MATLAB优化算法 

本书面向MATLAB的初中级读者，以MATLAB R2016a版本为平台，全面讲解了MATLAB软件在信号处理中常用的知识。本书按逻辑编排，自始至终采用实例描述，内容完整且每章相对独立，是一本全面的MATLAB信号处理应用参考书。
本书分为3个部分，共12章。*部分介绍了MATLAB的基础，涵盖的内容有MATLAB基础知识； 第二部分介绍了数字信号处理基本理论和方法及其MATLAB实现，涵盖的内容有信号与系统的分析基础、信号变换、IIR滤波器的设计、FIR滤波器设计、其他滤波器、小波在信号处理中的应用； 第三部分介绍了基于MATLAB信号处理的具体应用，涵盖的内容有基于MATLAB的语音信号处理、通信信号处理、雷达信号处理以及信号处理的图形用户界面工具与设计等内容。
本书以实用为目标，深入浅出，实例引导，讲解翔实，适合作为理工科高等院校研究生、本科生教学用书，也可作为广大科研工程技术人员的参考用书。

作者简介沈再阳  毕业于清华大学，研究生学历业，就职于中国核工业集团公司。精通MATLAB系统仿真软件，熟练应用MATLAB进行图像处理、信号处理等。在国内外期刊发表论文多篇，申请并获得授权专利多项，参与的项目获得集团三等奖。

目录

部分MATLAB基础
第1章MATLAB基础知识
1.1MATLAB概述
1.1.1MATLAB的发展历程
1.1.2MATLAB系统
1.2MATLAB工作环境
1.2.1命令行窗口
1.2.2帮助系统窗口
1.2.3图形窗口
1.2.4M文件编辑窗口
1.2.5当前文件夹
1.2.6搜索路径
1.3MATLAB程序控制结构
1.3.1顺序结构
1.3.2选择结构
1.3.3循环结构
1.3.4程序流程控制语句及其他常用命令
1.4变量、数值与表达式
1.4.1变量
1.4.2数值
1.4.3表达式
1.5数组与矩阵
1.5.1数组的创建与操作
1.5.2常见的数组运算
1.5.3矩阵的表示
1.5.4MATLAB矩阵寻访
1.5.5MATLAB矩阵的运算
1.6数据分析
1.6.1平均值、中值
1.6.2数据比较
1.6.3期望
1.6.4方差
1.6.5协方差与相关系数
1.7图形的绘制
1.7.1二维图形的绘制
1.7.2图形绘制和编辑
1.7.3三维图形的绘图
本章小结
第二部分信号处理的基本理论

第2章信号与系统的分析基础
2.1离散时间信号的概念
2.2采样定理
2.3离散时间序列
2.3.1单位采样序列
2.3.2单位阶跃序列
2.3.3正弦序列
2.3.4实指数序列
2.3.5复指数序列
2.3.6周期序列
2.4信号的基本运算
2.4.1序列相加与相乘
2.4.2序列累加与序列值乘积
2.4.3序列翻转与序列移位
2.4.4常用连续时间信号的尺度变换
2.4.5常用连续时间信号的奇偶分解
2.4.6信号的积分和微分
2.4.7卷积运算
2.5信号波形的产生
2.5.1线性调频函数与方波函数
2.5.2随机函数与三角波函数
2.5.3rectpuls函数与diric函数
2.5.4sinc函数与tripuls函数
2.5.5gauspuls函数与pulstran函数
2.6连续时间系统的时域分析
2.6.1连续时间系统的零状态与零输入响应的求解分析
2.6.2连续时间系统数值求解
2.6.3连续时间系统冲激响应和阶跃响应分析
2.6.4连续时间系统卷积求解
2.7离散时间信号在MATLAB中的运算
2.7.1离散时间系统
2.7.2离散时间系统响应
2.7.3离散时间系统的冲激响应和阶跃响应
2.7.4离散时间信号的卷积和运算
本章小结
第3章信号的变换
3.1Z变换概述
3.1.1Z变换的定义
3.1.2Z变换的收敛域
3.2Z变换的性质
3.2.1线性性质
3.2.2时域的移位
3.2.3时域扩展性
3.2.4时域卷积性质
3.2.5微分性
3.2.6积分性
3.2.7时域求和
3.2.8初值定理
3.2.9终值定理
3.3Z反变换
3.4离散系统中的Z域描述
3.4.1离散系统函数频域分析
3.4.2离散系统函数零点分析
3.4.3离散系统差分函数求解
3.5傅里叶级数和傅里叶变换
3.6周期序列的离散傅里叶级数
3.7离散的傅里叶变换
3.8离散傅里叶变换的性质
3.8.1线性
3.8.2循环移位
3.8.3循环卷积定理
3.8.4共轭对称性
3.9频率域采样
3.9.1频率响应的混叠失真
3.9.2频谱泄漏
3.9.3栅栏效应
3.9.4频率分辨率
3.10快速傅里叶变换
3.10.1直接计算DFT的问题及改进途径
3.10.2基2时分的FFT算法
3.10.3基2频分的FFT算法
3.10.4快速傅里叶变换的MATLAB实现
3.11离散余弦变换
3.11.1一维离散余弦变换
3.11.2二维离散余弦变换
3.11.3离散余弦函数
3.12Chirp Z变换
3.13Gabor函数
3.13.1Gabor函数定义
3.13.2Gabor函数的一般求法与解析理论
3.13.3Gabor展开
本章小结
第4章IIR滤波器的设计
4.1IIR滤波器结构
4.1.1直接型
4.1.2级联型
4.1.3并联型
4.2模拟滤波器的基础知识与原型设计
4.2.1巴特沃斯滤波器设计
4.2.2切比雪夫Ⅰ型滤波器设计
4.2.3切比雪夫Ⅱ型滤波器设计
4.2.4椭圆滤波器设计
4.3频带变换
4.3.1低通到低通的频带变换
4.3.2低通到高通的频带变换
4.3.3低通到带通的频带变换
4.3.4低通到带阻的频带变换
4.4冲激响应不变法与双线性变换法
4.5滤波器小阶数选择
4.6滤波器设计
4.6.1滤波器设计步骤
4.6.2经典滤波器设计
本章小结
第5章FIR滤波器设计
5.1FIR滤波器的结构
5.1.1直接型结构
5.1.2级联型结构
5.1.3频率采样型结构
5.1.4快速卷积型结构
5.2线性相位FIR滤波器的特性
5.2.1相位条件
5.2.2线性相位FIR滤波器频率响应的特点
5.2.3线性相位FIR滤波器的零点特性
5.3常用的窗函数法FIR滤波器设计
5.3.1窗函数的基本原理
5.3.2矩形窗
5.3.3汉宁窗
5.3.4海明窗
5.3.5布莱克曼窗
5.3.6巴特窗
5.3.7凯塞窗
5.3.8窗函数设计法
5.4频率采样的FIR滤波器的设计
5.4.1设计的思路与约束条件
5.4.2误差设计
5.5FIR数字滤波器的设计
5.5.1均方误差小化准则
5.5.2误差小化准则
5.5.3切比雪夫一致逼近
本章小结
第6章其他滤波器
6.1维纳滤波器
6.2卡尔曼滤波器
6.3自适应滤波器
6.3.1自适应滤波器简介
6.3.2自适应滤波器在MATLAB中的应用
6.4Lattice滤波器
6.4.1全零点Lattice滤波器
6.4.2全极点Lattice滤波器
6.4.3零极点的Lattice结构
6.5线性预测滤波器
6.5.1AR模型
6.5.2MA模型
6.5.3ARMA模型
本章小结
第7章随机信号处理
7.1随机信号处理基础
7.1.1随机信号的简介与时域统计描述
7.1.2平稳随机序列及其数字特征
7.1.3平稳随机序列的功率谱
7.1.4基于随机信号处理的MATLAB函数
7.2随机信号的功率谱分析
7.2.1非参量类方法
7.2.2参数法
7.2.3子空间法
本章小结
第8章小波在信号处理中的应用
8．1小波分析概述
8．1．1傅里叶变换与小波变换的比较
8．1．2多分辨分析
8．2小波变换
8．2．1一维连续小波变换
8．2．2高维连续小波变换
8．2．3离散小波变换
8．3小波包分析
8．3．1小波包的定义
8．3．2小波包的性质
8．3．3几种常用的小波
8．4小波工具箱介绍
8．4．1启动小波工具箱
8．4．2一维连续小波分析工具
8．5信号的重构
8．5．1idwt函数
8．5．2wavedec函数
8．5．3upcoef函数
8．5．4upwlev函数
8．5．5wrcoef函数
8．5．6wprec函数
8．5．7wprcoef函数
8．6提升小波变换用于信号处理
8．7信号去噪
8．7．1信号阈值去噪
8．7．2常用的去噪函数
8．8小波变换在信号处理中的应用
8．8．1分离信号的不同成分
8．8．2识别信号的频率区间与发展趋势
8．8．3基于小波变换的图像信号的局部压缩
8．8．4小波在数字图像信号水印压缩方面的应用
本章小结

第三部分信号处理的综合实例
第9章基于语音信号处理
9．1语音产生的过程
9．2语音信号产生的数学模型
9．2．1激励模型
9．2．2声道模型
9．2．3辐射模型
9．2．4语音信号的数字化和预处理
9．3语音信号分析和滤波处理
9．3．1语音信号的采集
9．3．2语音信号的读入与打开
9．3．3语音信号分析
9．3．4含噪语音信号的合成
9．3．5滤波器的设计
9．4小波变换在语音信号处理中的应用
9．4．1小波在语音信号增强中的应用
9．4．2小波变换在语音信号压缩上的应用
本章小结
第10章基于通信信号处理
10．1幅度调制
10．1．1DSBAM调制
10．1．2普通AM调制
10．1．3SSBAM调制
10．1．4残留边带幅度调制
10．2角度调制
10．3数字调制
10．3．1FSK调制
10．3．2PSK调制
10．3．3QAM调制
10．4自适应均衡
10．4．1递归小二乘算法(RLS)
10．4．2盲均衡算法
本章小结
第11章基于雷达信号处理
11．1雷达的基本原理
11．2雷达的用途
11．2．1双/多基达
11．2．2相控阵雷达
11．2．3宽带／超宽带雷达
11．2．4合成孔径雷达
11．2．5毫米波雷达
11．2．6激光雷达
11．3线性调频脉冲压缩雷达仿真
11．3．1匹配滤波器
11．3．2线性调频信号(LFM)
11．3．3相位编码信号
11．3．4噪声和杂波的产生
11．3．5杂波建模与MATLAB实现
11．4动目标的显示与检测
本章小结
第12章信号处理的图形用户界面工具与设计
12．1SPTool工具
12．1．1主窗口
12．1．2信号浏览器
12．1．3滤波浏览器
12．1．4频谱浏览器
12．1．5滤波器设计器
12．2图形用户界面(GUI)简介
12．2．1GUI的设计原则及步骤
12．2．2GUI模板与设计窗口
12．3控制框对象及属性
12．3．1按钮
12．3．2滑块
12．3．3单选按钮
12．3．4复选框
12．3．5静态文本
12．3．6可编辑文本框
12．3．7弹出式菜单
12．3．8列表框
12．3．9切换按钮
12．3．10面板
12．3．11按钮组
12．3．12轴
12．4MATLAB专用对话框
12．5GUI的设计工具
12．5．1布局编辑器
12．5．2对象浏览器
12．5．3用属性查看器设置控制框属性
12．5．4对齐对象
12．5．5Tab键顺序编辑器
12．5．6菜单编辑器
12．5．7编辑器
12．6回调函数
12．6．1Callback程序基本操作
12．6．2CreateFcn 
12．7脉搏信号处理的GUI设计
本章小结
参考文献

序言
前言

数字信号处理是从20世纪60年代以来，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科，它的重要性日益在各个领域的应用中表现出来。简言之，数字信号处理是把信号用数字或符号表示的序列，通过计算机或信号处理设备，用数字的数值计算方法处理，以达到提取有用信息、便于应用的目的。MATLAB 是一个功能强大的数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。目前，MATLAB已成为信号处理、图像处理、通信原理、自动控制等专业的重要基础课程的实验平台，而对于学生而言有效的学习途径是结合某一专业课程的学习掌握该软件的使用与编程。1. 本书特点由浅入深，循序渐进： 本书以初中级读者为对象，内容安排上考虑到MATLAB进行仿真和运算分析时的基础知识和实践操作，从基础开始，由浅入深地帮助读者掌握MATLAB的分析方法。步骤详尽，内容新颖： 本书结合作者多年MATLAB使用经验与实际应用案例，将MATLAB软件的使用方法与技巧详细地讲解给读者，使读者在阅读时能够快速掌握书中所讲内容。实例典型，轻松易学： 学习实际工程应用案例的具体操作是掌握MATLAB好的方式。本书通过综合应用案例，透彻详尽地讲解了MATLAB在各方面的应用。2. 本书内容本书结合多年MATLAB使用经验与实际工程应用案例，将MATLAB软件的使用方法与技巧详细地讲解给读者。本书在讲解过程中步骤详尽、内容新颖，讲解过程辅以相应的图片，使读者在阅读时一目了然，从而快速掌握书中所讲内容。本书对数字信号处理的基本理论、算法及MATLAB实现进行系统的论述。全书分为3部分，共12个章节，具体内容如下： 部分： MATLAB基础。介绍MATLAB的基础知识、发展史及基本运算等内容，让读者对MATLAB有一个概要性的认识。具体的章节安排如下： 第1章MATLAB基础知识第二部分： 信号处理的基本理论。介绍了数字信号处理基本理论和方法及其MATLAB实现，向读者展示了MATLAB在处理数字信号方面的方法及技巧。具体的章节安排如下： 第2章信号与系统的分析基础第3章信号变换第4章IIR滤波器的设计第5章FIR滤波器设计第6章其他滤波器第7章随机信号处理第8章小波在信号处理中的应用 第三部分： 信号处理的综合实例。介绍了MATLAB基于图像信号处理、语音信号处理、通信信号处理、雷达信号处理等在实际中应用，让读者进一步领略到MATLAB的强大功能和广泛的应用范围。具体的章节安排如下： 第9章基于语音信号处理第10章基于通信信号处理第11章基于雷达信号处理第12章信号处理的图形用户界面工具与设计3. 读者对象本书适合于MATLAB初学者和期望提高应用MATLAB进行信号处理能力的读者，具体说明如下： ★ 初学MATLAB的技术人员★ 广大从事信号处理的科研工作人员★ 大中专院校的教师和在校生★ 相关培训机构的教师和学员★ 参加工作实习的“菜鸟”★ MATLAB爱好者4. 读者服务为了方便解决本书疑难问题，读者朋友在学习过程中遇到与本书有关的技术问题，可以发邮件到邮箱[email protected]，或者访问博客http://blog.sina.com.cn/caxart，编者会尽快给予解答，我们将竭诚为您服务。5. 本书作者本书主要由沈再阳编著。此外，付文利、王广、张岩、温正、林晓阳、任艳芳、唐家鹏、孙国强、高飞等也参与了本书部分内容的编写工作，在此表示感谢。虽然作者在本书的编写过程中力求叙述准确、完善，但由于水平有限，书中欠妥之处在所难免，希望读者和同仁能够及时指出，共同促进本书质量的提高。后再次希望本书能为读者的学习和工作提供帮助！编著者2017年5月

专家评论★“科学与工程计算技术丛书”系列主题反映了MATLAB和Simulink汇集的领域——大规模编程、机器学习、科学计算、机器人等。我们高兴地看到“科学与工程计算技术丛书”支持MathWorks一直以来追求的目标：助您加速工程技术和科学研究。——Jim Tung (MathWorks Fellow)★ 《MATLAB信号处理》面向MATLAB的入门级读者由浅入深地全面论述了MATLAB软件在信号处理中常用的知识。按逻辑编排，自始至终采用实例描述；内容完整且每章相对独立，是一本全面的MATLAB信号处理应用参考书。——张  岩（北京航空航天大学）  ★ 《MATLAB信号处理》的主要内容来源于作者多年工作经验的总结，无论是从案例的选择、讲解的详细程度、书中涉及的专业知识等各方面都充分考虑读者的喜好，力求使得读者尽快掌握MATLAB信号处理应用。——付文利（烽火通信集团）   ★ 《MATLAB信号处理》一书是在MathWorks公司技术人员的大力支持下，由中科院、清华、北理工等相关科研院校机构的一线专家联合编著，内容专业性强、实践性强，书中很多项目实例拿来可用，二次开发，值得读者认真学习。——刘成龙（中国科学院）★ 《MATLAB信号处理》一书通过包含信号处理学科多个重要知识点的简单明了的实例介绍，使读者在较短的时间内掌握利用MATLAB工具解决信号处理相关问题的基本方法。——林晓阳（中国科学院）★ 《MATLAB信号处理》一书构思新颖，实践性强，内容丰富，深入浅出，所有程序及应用实例均可通过MATLAB上机调试。无论是对数字信号处理的学习，还是利用MATLAB开发实践性课题，或是培养创新精神和解决实际问题，本书都具有十分重要的意义。        ——王  广（北京化工大学）

第3章信号的变换
信号与系统的分析方法中，除了时域分析方法外，还有变换域分析的方法。连续时间信号与系统的变换域分析方法主要是傅里叶变换和拉普拉斯变换。离散时间信号的Z变换是分析线性时不变离散时间系统问题的重要工具，在数字信号处理、计算机控制系统等领域有着广泛的应用。学习目标： （1） 了解、熟悉Z变换的概念与性质； （2） 理解Z反变换的相关内容； （3） 掌握离散系统中的Z域描述方法； （4） 了解、熟悉傅里叶级数与变换； （5） 理解离散傅里叶变换及其性质； （6） 实现频率域采样和快速傅里叶变换； （7） 熟悉实现离散余弦变换、Chirp Z变换和Gabor函数。3.1Z变换概述连续系统一般使用微分方程、拉普拉斯变换的传递函数和频率特性等概念进行研究。一个连续信号f(t)的拉普拉斯变换F(s)是复变量s的有理分式函数，而微分方程通过拉普拉斯变换后也可以转换为s的代数方程，从而可以大大简化微分方程的求解，从传递函数可以很容易地得到系统的频率特征。因此，拉普拉斯变换作为基本工具将连续系统研究中的各种方法联系在一起。计算机控制系统中的采样信号也可以进行拉普拉斯变换，从中找到简化运算的方法，引入了Z变换。3.1.1Z变换的定义序列x(n)的Z变换(简称ZT)定义为

X(z)=∑ ∞n=－∞x(n)z－n

上式称为双边Z变换。如果x(n)的非零值区间为(－∞,0］或者［0, ∞),则上式可变为

X(z)=∑0n=－∞x(n)z－n
X(z)=∑ ∞n=0x(n)z－n

此时，称为序列x(n)的单边Z变换。序列的ZT存在的条件为

|X(z)|=∑ ∞n=－∞x(n)z－n≤∑ ∞n=－∞x(n)z－n=∑ ∞n=－∞|x(n)||z－n|< ∞

满足上式的z的取值范围称为Z变换的收敛域(Region of Convergence,ROC),它通常为z平面上的一个环状域，即

Rx－
3.1.2Z变换的收敛域序列Z变换的收敛域与序列的形态有关。反之，同一个Z变换的表达式，不同的收敛域，确定了不同序列形态。下面根据序列形态不同，分别讨论其收敛域。对于任意给定的序列x(n)，能使X(z)=∑∞n=－∞x(n)z－n收敛的所有z值集合为收敛域。即满足

∑∞n=－∞|x(n)z－n|
不同的x(n)的Z变换，由于收敛域不同，可能对应于相同的Z变换，故在确定Z变换时，必须指明收敛域。1. 有限长序列有限序列的描述函数是

x(n)=x(n)n1≤n≤n2
0其他

其Z变换为

X(z)=∑n2n=n1x(n)z－n

因此Z变换式是有限项之和，故只要级数的每一项有界，则级数就收敛。收敛域为

0
2. 右边序列右边序列的描述函数是

x(n)=x(n)n≥n1
0其他

其Z变换为

X(z)=∑∞n=n1x(n)z－n

因此Z变换样式是无限项之和，当n1≥0时，由根值判别法有

limn→∞n|x(n)z－n|<1

所以此时收敛域为

|z|>limn→∞n|x(n)|=R1

当n1<0时，此时级数全收敛，所以右边序列的收敛域为R1
x(n)=x(n)n≤n2
0其他

其Z变换为

X(z)=∑n2n=－∞x(n)z－n=∑∞n=－n2x(－n)zn

当n2<0时，由根值判别法有

limn→∞n|x(－n)zn|<1

由此求得的收敛域为

|z|
当n2>0时，此时相当于增加了一个n2>0的有限长序列，还应除去原点，左边序列的收敛域为

0
4. 双边序列双边序列的描述函数为

x(n)=x(n)［u(－n－1) u(n)］

其Z变换为

X(z)=∑∞n=－∞x(n)z－n=∑－1n=－∞x(n)z－n ∑∞n=0x(n)z－n

因为∑∞n=0x(n)z－n的收敛域为|z|>R1，∑－1n=－∞x(n)z－n的收敛域为|z|
R1
3.2Z变换的性质3.2.1线性性质
假设

Z［x1(k)］=X1(z)（|z|>Rx1)
Z［x2(k)］=X2(z)(|z|>Rx2)

则有

Z［ax1(k) bx2(k)］=aX1(z) bX2(z)

其中，a、b为任意常数。3.2.2时域的移位假设Z［f(t)］=F(z)，那么有

Z［f(t nT)］=znF(z)－∑n－1k=0f(kT)z－k

假设Z［f(t)］=F(z)，那么有

Z［f(t－nT)］=z－nF(z)

3.2.3时域扩展性若函数f(t)有Z变换F(z)，则

Z［eatf(t)］=F(ze±aT)

根据Z变换定义有

Z［eatf(t)］=∑∞k=0f(kT)eakTz－k

令z1=ze±aT，则上式可写成

Z［eatf(t)］=∑∞k=0f(kT)z－k1=F(z1)

代入z1=ze±aT，得

Z［eatf(t)］=F(ze±aT)

3.2.4时域卷积性质已知

x(k)X(z)(α1h(k)H(z)(α2
则有

x(k)*h(k)X(z)H(z)

3.2.5微分性如果有

x(k)X(z)α
那么有

kx(k)－zdX(z)dzα
3.2.6积分性已知

x(k)X(z)α
则有

x(k)k mzm∫∞zX(η)ηm 1dηα
3.2.7时域求和如果有

x(k)X(z)α
那么有

f(k)=∑ki=－∞x(i)zz－1X(z)max(α,1)
3.2.8初值定理如果函数f(t)的Z变换为F(z)，并存在极限limz→∞F(z)，则

limk→0f(kT)=limz→∞F(z)

3.2.9终值定理假定f(t)的Z变换为F(z)，并假定函数(1－z－1)F(z)在z平面的单位圆上或圆外没有极点，则

limk→∞f(kT)=limz→1(1－z－1)F(z)