第2章-MATLAB的数值计算.ppt

第2章 MATLAB的数值计算 所谓数值计算，就是指计算的表达式、变量中不得包含未经定义的自由变量。 本章主要内容： 2．1 变量及其赋值 2．2 矩阵和数组的算术运算 2．3 线性方程组的求解 2．4 矩阵函数 2．5 多项式运算 2. 6 练习 2.1 变量及其赋值 在命令窗口可直接输入变量并赋值，无须声明其类型和维数，MATLAB将自动处理。 MATLAB 的变量名称由英文大小写52个字母、数字和下划线等符号组成，并且第一个字符必须是英文字母，变量和常量最长允许31个字符。须注意的是，在默认状态下，MATLAB对字母大小写敏感 。 MATLAB数据格式与显示形式 MATLAB既可用传统的十进制数表达数值，也可以用科学计数表达数值，用e代表10的指数形式，用i和j来代表虚数。 MATLAB内部数据格式只有一种，是IEEE浮点标准的双精度二进制（64位），相应于十进制的16位有效数，范围为10-308 ~ 10+308。 为了人机交互的友好性，数据输出显示格式有8种 。可用菜单选项或format命令选择。 MATLAB定义的特别变量及其意义 2．2 矩阵和数组的算术运算 2.2.1 矩阵的定义 MATLAB中的变量或常量都代表矩阵； 矩阵变量中的元素用“( )”中的数字来注明其下标。维数不同，数字个数也不同，数字之间用“ ，”隔开 使用下标，用户可以单独给元素赋值，如x(2)=1.7321，a(2,3)=6等 实例：访问矩阵中的元素 x=(0:pi/4:2*pi); y=sin(x); A=[x,y] %先构造矩阵A 上面构造的A由列向量x和在x基础上生成的列向量y合并组成。由定义，x是一长度为9的列向量，因此这样构造的A是一个9 ×2的矩阵。 b=A(2,1) 此命令访问了A的第二行第一列元素。 创建特殊矩阵的函数 zeros() 矩阵元素都是零 ones() 矩阵元素都是1 注意： () 中的数字用于定义矩阵维数， ()内不是数字的情况参看Help。 rand() 矩阵元素是0到1之间均匀分布的随机数 randn() 矩阵元素是正态分布的随机数 2.2.2 矩阵运算 MATLAB提供矩阵算术运算有： 加 (+)、减 (-)、乘 (*)、除 (/)、幂次方 (^)、转置() 注意，对于矩阵来说，乘法有左乘右乘之分，A *B时要求A的列数等于B的行数，即满足矩阵相乘的条件。 MATLAB中，为求解线性方程组，引入了矩阵除法。 MATLAB中的矩阵除法有左除(\)和右除(/)两种。 以下是矩阵运算的合法命令形式： A+B 、A-B、A*B、A/B、A\B、A^p（此处p只能是标量，A必须为方阵）、A 矩阵的左除和右除 Case 1：方程为 A *x = B 若A非奇异，用左除解得 x = A \ B 如A为方阵，则A \ 形式上相当于A的逆阵。但须注意，MATLAB中，求解线性方程组并不采用求逆的方法。 Case 2：方程为 x *A = B 若A非奇异，用右除解得 x = B / A 如A为方阵， / A形式上相当于A的逆矩阵。 2.2.3 数组运算 和矩阵运算符比较，除了加减符号外，其余的数组运算符号均须加上一个 “ . ”符号，即： 加 (+)、减 (-)、乘 (.*)、左除 (./)、右除(.\)、幂次方 (.^) 数组运算是矩阵中元素对元素的操作，因而进行数组运算时，加、减、乘、除均要求A、B有相同的维数。对加减运算而言，A和B中有一个是标量也可以。 数组相除也有左除和右除之分，其中A./B表示A的元素为被除数，而A.\B 则表示A的元素是除数。 数组运算A .^p中，如p是标量，则此运算是指对A中每个元素按幂次p进行乘方运算。但若数组运算A.^B中，A为标量，B为矩阵，则A.^B的含义为A B(i)。 范数 范数：信号处理问题中，信号常表示为线性空间中的一个向量，而向量的长度是一个非常重要的量。利用上述的数组运算A.^p，向量X的长度的平方很容易表示为 sum(X.^2) 范数（norm）概念可容易地用来度量向量长度。 MATLAB 中设置了专用的norm函数 norm(X，p)= sum(X.^p)^(1/p