二MATLAB数值运算.docVIP

第二章 MATLAB数值运算 MATLAB中的数值运算是MATLAB最基本的也是最重要的部分。掌握矩阵、向量、数组和多项式的构造和运算方法对于解决众多实际问题将提供更多的便利。 矩阵 MATLAB是基于矩阵运算的处理工具。因此矩阵也是MATLAB的核心。系统默认的操作对象并不是一个个数字，而是矩阵。在线性代数（高等代数）中矩阵是基本运算单元，矩阵可以进行多种运算，如加、减、内积、逆矩阵、转置、线性方程式、特征值、特征向量、矩阵的分解与分块等。MATLAB支持线性代数所定义的全部矩阵运算，用户可通过MATLAB方便地处理线性代数的运算，即省时又简单。 作为MATLAB本身，它最初的目的也是为了解决矩阵问题而开发的，甚至连其英文名称都是matrix（矩阵）laboratory（实验室）两个英文单词缩写的组合。可见矩阵在MATLAB中的重要地位。 2.1.1 矩阵的构造 通常矩阵是指含有m行n列数字的矩形结构。其定义如下： 由个数排成行列的数表称为矩阵，记为 其中可以是实数或复数，由此可以将矩阵分为实数矩阵或复数矩阵。在MATLAB中输入矩阵的具体步骤如下： 1）用中括号[]把所有矩阵元素括起来； 2）同一行中的不同元素用空格或逗号隔开； 3）用分号指定一行结束；（也可以分成几行输入，用回车符代替分号） 4）数据的元素可以是表达式，系统将自动计算结果。 例1 请输入矩阵。 根据上述矩阵的输入方法，直接在命令窗口中输入即可，具体方法如下： a=[1,2,3;4,5,6] a = 1 2 3 4 5 6 b=[2 -3 4 -2;7 6 -5 8;0 4 9 -1] b = 2 -3 4 -2 7 6 -5 8 0 4 9 -1 c=[1-i,2;3,7+3i] c = 1.0000 - 1.0000i 2.0000 3.0000 7.0000 + 3.0000i d=[sqrt(3),pi 5^(1/4),exp(1)] d = 1.7321 3.1416 1.4953 2.7183 矩阵不仅可以通过直接在命令窗口中输入的方法产生，还可以通过以下方法来产生： 通过m文件来产生 点击MATLAB主窗口中file(New(M-File将出现下述界面 在此窗口中输入矩阵即可。输入方法和MATLAB命令窗口的输入方法是一样的。输入完毕后保存（CTRL+S），并给m文件起名。起完名后在MATLAB命令窗口中输入m文件名就可以显示矩阵信息了。 利用外部数据文件导入到指定的矩阵 作为例子在这里简单介绍如何导入Excel中的数据，其它数据的导入将在后面讲到。打开Excel选择菜单栏中的工具(加载宏(在该窗口中点击浏览按钮(找到MATLAB安装目录底下的toolbox并点击exlink目录下的excellink，此时excel界面中多出现startmatlab putmatrix文字，如下图所示 此后选定即将导入的Excel数据，并点击putmatrix，届时将出现Variable name in MATLAB，给变量取好名字并选择确定键。此时回到MATLAB命令窗口中输入变量名即可访问该变量的详细信息了。需要注意的是如果选定的单元格是空的，则系统将此单元格自动赋值为0。 利用已经构造好的矩阵来构成更大的矩阵 例如: x1=[1,2,3,4] x1 = 1 2 3 4 x2=[5,6,7,8] x2 = 5 6 7 8 x3=[0,x1,5] x3 = 0 1 2 3 4 5 x4=[x1,x2] x4 = 1 2 3 4 5 6 7 8 2.1.2 矩阵元素的访问与子矩阵提取 和线性代数中对矩阵元素的表示方法一样，在MATLAB中也可以采用下标来定位矩阵元素。而且，MATLAB还进一步扩充了矩阵下标的表达功能。通过矩阵的下标，不仅能提取某一个矩阵元素的值，也可以提取某个子块。以下列出了访问矩阵A的相关内容。 A(m,n) %提取第m行n列的元素 A(:,n) %提取第n列元素 A(m,:) %提取第m行元素 A(m1:m2,n1:n2) %提取m1行到m2行和第n1列到n2列的所有元素。 A(m:end,n) %提取第n列的第m行到最末行的子块 A(m,n:end) %提