高等数学上教学课件作者张圣勤黄勇林姜玉娟MATLAB实验ex77课件.pptVIP

机械工业出版社 Jixie gongye chubanshe 高等数学实验 MATLAB7.0实验七 矩阵的基本运算 求解线性方程组 方程的数值解 线性代数 上海电机学院 张圣勤 一、矩阵的基本运算 简单矩阵函数 函 数 名 功 能 与 含 义 Det(A) 求方阵A的行列式值。即求|A| A’ 将矩阵A转置。即把A中相应行写成相应的列。 Rank(A) 求矩阵A的秩。即求A中线性无关的行数或列数。 Inv(A) 求方阵A的逆矩阵。如果A是奇异矩阵或近似奇异矩阵，则给出错误信息。 Pinv(A) 求矩阵A的伪逆。若A是m×n矩阵，则伪逆的大小为n×m。若A是非奇异矩阵，则pinv(A)=inv(A)。 Trace(A) 求矩阵A的迹。即求A中对角线元素之和。 Eig(A) 求矩阵A的特征值。 例1 设有矩阵A、B 求它们的乘积AB、BA、A/B及它们各自的秩和A的迹。 具体操作如下： A=[1 3 5;4 5 6;5 6 7]; %输入矩阵A B=[2 4 6;9 6 3;4 3 7]; %输入矩阵B A*B %计算矩阵乘积AB ans = 49 37 50 77 64 81 92 77 97 B*A %计算矩阵乘积BA ans = 48 62 76 48 75 102 51 69 87 rank(A) %求矩阵A的秩 ans = 2 rank(B) %求矩阵B的秩 ans = 3 trace(A) %求矩阵A的迹 ans = 13 A/B %求A除以B的商 ans = 0.8750 -0.0833 -0.0000 0.8750 0.2500 0 1.0000 0.3333 -0.0000 二、求解线性方程组 线性方程组的求解用矩阵除法计算，只要输入系数矩阵A和常数矩阵B， 方程组就成为矩阵方程AX=B，则该矩阵方程的解为X=A\B。 例2、求线性方程组 的解。 由线性方程组可知： 系统矩阵为 ，常数矩阵为 。操作如下： A=[2 4 6;9 6 3;4 3 7]; %输入系数矩阵A B=[7;9;6]; %输入常数矩阵B X=A\B %求未知矩阵X X = 0.0250 1.3250 0.2750 如果系数矩阵A是非奇异方阵，也可以先用inv命令求出A的逆阵NA， 再用X=NA*B求出方程组的解。 NA=inv(A) %求矩阵A的逆阵 NA = -0.2750 0.0833 0.2000 0.4250 0.0833 -0.4000 -0.0250 -0.0833 0.2000 X=NA*B %求未知矩阵X X = 0.0250 1.3250 0.2750 三、方程（组）的数值解 方程与方程组分为三大类，一类是线性方程，一类是非线性方程，还有 一类就是微分方程。不管是哪一类方程都是给出未知数（或函数）之间（或 函数与自变量之间）的关系，根据这些关系求出适合方程的数或函数，这在 数学中叫解方程。 在MATLAB中求解方程的函数有roots、fsolve和fzero等。 1 、方程的求解 单个方程按其系数的性质分为线性方程和非线性方程。这些方程的求解 一般有二种方法，一种是直接用MATLAB中的解方程的命令函数求解，一种 是根据求解方程的方法编写M文件，然后用自己编写的M文件求解。这里只 介绍直接求解法。 解方程函数调用