MatLab中的矩阵.docVIP

MatLab中的矩阵 我们知道,求解线性方程组是线性代数课程中的核心内容,而矩阵又在求解线性方程组的过程中扮演着举足轻重的角色。下面我们就利用科学计算软件MATLAB来演示如何使用矩阵,同时,也使学生对线性代数的认识更加理性。 一、矩阵的构造 在MatLab中,构造矩阵的方法有两种。一种是直接法,就是通过键盘输入的方式直接构造矩阵。另一种是利用函数产生矩阵。 例1.利用pascal函数来产生一个矩阵 A=pascal(3)A=1 1 1 1 2 3 1 3 6 例2.利用magic函数来产生一个矩阵B=magic(3)B= 8 1 6 3 5 7 4 9 2 例3.还可以利用函数产生一个4*3的随机矩阵 c=rand(4,3) c= 0.9501 0.8913 0.8214 0.2311 0.7621 0.4447 0.6068 0.4565 0.6154 0.4860 0.0185 0.7919 例4.利用直接输入法可产生列矩阵、行矩阵及常数 u=[3;1;4]u=3 1 4 v=[2 0 -1]v=2 0 -1 s=7s=7 二、矩阵的基本运算 1、四则运算 例5.矩阵的加法 X=A+BX=9 2 74 7 105 12 8 例6.矩阵的减法Y=X-AY=8 1 6 3 5 7 4 9 2 注: 若二个矩阵的大小不完全相同，则会出错！ 例如,X=A+u??? Error using == plus Matrix dimensions must agree。 例7.矩阵的乘法 X=A*B X= 15 15 15 26 38 26 41 70 39 注: 若第一个矩阵的列数和第二个矩阵行数不相同,这两个矩阵就不可以相乘。 例如,X=A*v ??? Error using == mtimes Inner matrix dimensions must agree。 在MATLAB中,矩阵的除法有两个运算符号,分别为左除“\”与右除“/”,矩阵的右除运算速度要慢一点,而左除运算可以避免奇异矩阵的影响,它们的作用主要用于求解线性方程组,我们在后面会涉及到矩阵的除法。 2、矩阵的转置、逆运算及行列式运算 与线性代数中一样,矩阵的转置只需用符号“,”来表示即可。 例8.求矩阵B的转置 X=BX=8 3 4 1 5 9 6 7 2 线性代数中求矩阵逆的运算非常复杂,而在MATLAB中,矩阵的逆运算只需要函数“inv”来实现,这大大简化了计算过程。 例9.求矩阵A的逆 X=inv(A) X= 3 -3 1 -3 5 -2 1 -2 1 在MATLAB中,求矩阵的行列式大小,可用函数“det”实现。 例10.求矩阵A的行列式 X=det(A) X= 1 注: 在求矩阵的逆和行列式时,一定要求矩阵是一个方阵,否则会出错! 例如,X=inv(u) ??? Error using == inv Matrix must be square。 再如,X=det(u) ??? Error using == det Matrix must be square。 三、矩阵的常用函数运算 1.矩阵的特征值运算 在线性代数中,计算矩阵特征值及特征向量的过程相当麻烦,但在MATLAB中,矩阵特征值运算只需要函数“eig”或“eigs”即可。 例11.求矩阵A的特征值及特征向量 [b,c]=eig(A) b= -0.5438 -0.8165 0.1938 0.7812 -0.4082 0.4722 -0.3065 0.4082 0.8599 c= 0.1270 0 0 0 1.0000 0 0 0 7.8730 上例中的b、c矩阵分别为特征向量矩阵和特征值矩阵。 2.矩阵的秩运算 矩阵的秩在求解线性方程组中应用非常广泛,而在线性代数中计算矩阵的秩也非常复杂,但在MATLAB中,矩阵的秩只需要用函数“rank”即可。 例12.求矩阵A的秩 x=rank(A) x= 3 3.矩阵的正交化运算 在MATLAB中,矩阵的正交化运算可由函数“orth”计算得到。下面的例子用来求矩阵的一组正交基,有了正交基就可以对矩阵进行正交化了。 例13.求矩阵A的正交基 x=