周期性矩形脉冲频谱.docVIP

第8章 MATLAB程序设计语言在信号处理中的应用 8.1 概 述 8.1.1 MATLAB程序设计语言简介 MATLAB，Matrix Laboratory的缩写，是由Mathworks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。 与大家常用的Fortran和C等高级语言相比，MATLAB的语法规则更简单，更贴近人的思维方式，被称之为“草稿纸式的语言”。截至目前，MATLAB已经发展到12.1版， 适用于所有32位的Windows操作系统， 按NTFS(NT文件系统)格式下完全安装约需 850 MB。MATLAB软件主要由主包、仿真系统和工具箱三大部分组成。 8.2 基本数值运算 8.2.1 MATLAB内部特殊变量和常数 MATLAB内部有很多变量和常数， 用以表达特殊含义。常用的有： (1) 变量ans： 指示当前未定义变量名的答案。  (2) 常数eps：表示浮点相对精度， 其值是从1.0到下一个最大浮点数之间的差值。该变量值作为一些MATLAB函数计算的相对浮点精度，按IEEE标准， eps=2-52,　近似为2.2204e-016 (3) 常数Inf： 表示无穷大。 当输入或计算中有除以0时产生Inf。  (4) 虚数单位i,j： 表示复数虚部单位， 相当于 。 （5) NaN： 表示不定型值， 是由 0/0 运算产生的。  (6) 常数pi： 表示圆周率π， 其值为3.141 592 653 589 7…。 8.2.2 变量类型 1． 变量命名规则 MATLAB中对变量的命名应遵循以下规则： （1） 变量名可以由字母、　数字和下划线混合组成， 但必须以字母开头。 （2） 字符长度不能大于31。 （3） 变量命名区分大小写 2． 局部变量和全局变量 局部变量是指那些每个函数体内自己定义的，不能从其他函数和MATLAB工作空间访问的变量。  全局变量是指用关键字“global”声明的变量。 全局变量名应尽量大写，并能反映它本身的含义。如果需要在工作空间和几个函数中都能访问一个全局变量，必须在工作空间和这几个函数中都声明该变量是全局的。 8.2.3 矩阵及其运算 MATLAB具有强大的矩阵运算和数据处理功能， 对矩阵的处理必须遵从代数规则。  1． 矩阵生成 1） 一般矩阵的生成 对于一般的矩阵MATLAB的生成方法有多种。 最简单的方法是从键盘直接输入矩阵元素。直接输入矩阵元素时应注意： 各元素之间用空格或逗号隔开，用分号或回车结束矩阵行，用中括号把矩阵所有元素括起来。 2） 特殊矩阵的生成 对于特殊的矩阵可直接调用MATLAB的函数生成。  用函数zeros生成全0矩阵：格式 B=zeros(m,n)生成m×n的全0阵。  用函数ones生成全1矩阵：格式 B=ones(m,n)生成m×n的全1阵。  用函数eye生成单位阵：格式 B=eye(m,n)生成m×n矩阵， 其中对角线元素全为1，其他元素为0。 2． 矩阵的运算 矩阵的运算有基本运算和函数运算两种类型。基本运算包括矩阵的加、减、乘、除、乘方、求转置、求逆等，其主要特点是通过MATLAB提供的基本运算符＋、－、*、/（\）、^等即可完成。函数运算主要是通过调用MATLAB系统内置的运算函数来求取矩阵的行列式（det(A)）, 求秩（rank(A)）, 求特征值和特征向量（［V， D］＝eig(A)）, 求Jordan标准形（jordan(A)）和矩阵分解等。需要用时可以参阅联机帮助和相关参考书。 例 矩阵的基本运算。 A=［1, 2, 3; 4, 5, 6］; B =［6, 5, 4; 3, 2, 1］; C =A+B %计算两个矩阵的和 D =B′ %计算矩阵B的转置 E=A*D %做矩阵乘法，必须要满足矩阵乘法的基本要求  F=det(E) %求E的行列式值 G=E^(-1) %求E的逆 8.3 基本语句 8.3.1 程序控制语句 1.　循环语句 MATLAB的循环语句包括for循环和while循环两种类型。  1) for循环 语法格式：  for 循环变量 = 起始值： 步长： 终止值 循环体 end 起始值和终止值为一整形数，步长可以为整数或小数，省略步长时，默认步长为1。执行for循环时，判定循环变量的值是否大于(步长为