矩阵在MATLAB中的运算与应用.docxVIP

矩阵在MATLAB中的运算与应用摘要:介绍了Matlab在矩阵运算方而的功能。演示了用Matlab构造矩阵，获取矩阵的相关信息，进行矩阵运算的方法，对矩阵运算进行了分析，对矩阵作图进行了研究。关键词:矩阵； Matlab1 引言Matlab的含义是矩阵实验室( Matrix Laboratory) ,是由美国Mathwork公司于1984年推出的一套高性能的数值计算和可视化软件［1］。现在，它己发展为国际上最优秀的科技应用软件。如果能将它用到相关学科课的学习上无疑是非常有意义的。Matlab赋予学习者一个可实验的环境，一个强大的数值计算和分析及可视化(图形)工具。矩阵论是高等院校理、工科研究生的一门重要基础课程。有人认为“科学计算，归根结底就是矩阵的计算”［2］。因此，对于将来从事科学技术工作的研究生来说，矩阵理论和方法是必不可少的数学工具。矩阵的理论和方法在数学和其他学科中都具有重要的意义，但许多学生无法克服矩阵庞大的计算量带来的恐惧，从而丧失了学习的兴趣和动力。本文展示了如何方便地用Matlab构造矩阵，获取矩阵的相关信息以及完成矩阵的运算，展示了矩阵的结构和运算，以此来说明在机器计算环境中，庞大复杂的计算不再是令人头疼的事情。矩阵及其运算矩阵是进行数据处理和运算的基本元素。在MATLAB中：a、通常意义上的数量（标量）可看成是“1*1”的矩阵；b、n维矢量可看成是“n*1”的矩阵；c、多项式可由它的系数矩阵完全确定。2.1 矩阵的创建在MATLAB中创建矩阵有以下规则：a、矩阵元素必须在“[ ]”内；b、矩阵的同行元素之间用空格（或“,”）隔开；c、矩阵的行与行之间用“;”（或回车符）隔开；d、矩阵的元素可以是数值、变量、表达式或函数；e、矩阵的尺寸不必预先定义。2.1.1直接输入法最简单的建立矩阵的方法是从键盘直接输入矩阵的元素，输入的方法按照上面的规则。建立向量的时候可以利用冒号表达式，冒号表达式可以产生一个行向量，一般格式是： e1:e2:e3，其中e1为初始值，e2为步长，e3为终止值。还可以用linspace函数产生行向量，其调用格式为：linspace(a,b,n) ，其中a和b是生成向量的第一个和最后一个元素，n是元素总数。可以看出来linspace(a,b,n)与a:(b-a)/(n-1):b等价。2.1.2利用MATLAB函数创建矩阵基本矩阵函数如下：(1) ones()函数：产生全为1的矩阵，ones(n)：产生n*n维的全1矩阵，ones(m,n)：产生m*n维的全1矩阵；(2) zeros()函数：产生全为0的矩阵；(3) rand()函数：产生在（0，1）区间均匀分布的随机阵；(4) eye()函数：产生单位阵；(5) randn()函数：产生均值为0，方差为1的标准正态分布随机矩阵。2.1.3利用文件建立矩阵当矩阵尺寸较大或为经常使用的数据矩阵，则可以将此矩阵保存为文件，在需要时直接将文件利用load命令调入工作环境中使用即可。同时可以利用命令reshape对调入的矩阵进行重排。reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m*n的二维矩阵。2.2 矩阵的拆分2.2.1矩阵元素可以通过下标（行列索引）引用矩阵的元素，如 Matrix(m,n)。也可以采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依次类推。序号(Index)与下标(Subscript)是一一对应的，以m*n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)*m+i。其相互转换关系也可利用sub2ind和ind2sub函数求得。2.2.2矩阵拆分1、利用冒号表达式获得子矩阵：(1) A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列的元素。(2) A(i:i+m,:)表示取A矩阵第i~i+m行的全部元素；A(:,k:k+m)表示取A矩阵第k~k+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第i~i+m行内，并在第k~k+m列中的所有元素。此外，还可利用一般向量和end运算符来表示矩阵下标，从而获得子矩阵。end表示某一维的末尾元素下标。2、利用空矩阵删除矩阵的元素：在MATLAB中，定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。注意，X=[]与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为0。2.2.3特殊矩阵(1) 魔方矩阵：魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵，其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方