w05矩阵操作与线性方程组求解-正式版.ppt

矩阵的行列式 矩阵的行列式运算，如 >> A =[5 -2 0; -2 6 2; 0 2 7]; >> det(A) ans = 162 迹函数 ? 矩阵所有对角线上元素的和称为矩阵的迹，在 MATLAB 中可由 trace 函数计算得出。 >> A =[5 -2 0; -2 6 2; 0 2 7]; >> trace(A) ans = 18 矩阵的加减乘除法 ? 矩阵的加减就是对应元素的加减 ? 矩阵的点乘点除（ .* ， ./ ）就是对应元素的乘除 如： c = 8 3 7 7 3 3 3 5 8 d = 6 5 6 4 4 8 7 7 10 e= 2 -2 1 3 -1 -5 -4 -2 -2 e=c-d f = 48 15 42 28 12 24 21 35 80 f=c.*d 加减与点乘点除 乘法 c = 8 3 7 7 3 3 3 5 8 d = 6 5 6 4 4 8 7 7 10 g=c*d g = 109 101 142 75 68 96 94 91 138 矩阵的除法 ? 在 Matlab 中，求解线性方程组归结为两条简 单的命令，即：“ \ ” 和“ /” ，使得求解线性 方程组变得十分容易。 对于方程组 Ax=b ，其中 A 是一个（ n × m ）阶的矩阵， 则： 当 n=m 且非奇异时，此方程称为恰定方程 恰定方程组 Ax b ? 1 3 7 2 4 3 3 7 2 A ? ? ? ? ? ? ? ? ? ? ? ? ? ? 3 7 2 2 4 3 1 3 7 2 3 x y z x y z x y z ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 2 1 3 b ? ? ? ? ? ? ? ? ? ? ? ? >> A=[1 -3 7;2 4 -3;-3 7 2]; >> b= [2 -1 3]'; >> x=A\b 左除求解 A 矩阵的列数等于 b 矩阵的行数！ 逆函数 ? 矩阵的逆可以通过函数 inv 来求解 格式 Y=inv(X) % 求方阵 X 的逆矩阵。若 X 为 奇异阵或近似奇异阵，将给出警告信息。 ? Y=A^(-1) 也可以实现矩阵求逆运算 求逆矩阵函数： inv 求逆矩阵 >> A=[1 2 3; 2 2 1; 3 4 3]; >> Y=inv(A) Y = 1.0000 3.0000 -2.0000 -1.5000 -3.0000 2.5000 1.0000 1.0000 -1.0000 超定方程组 2 3 1 5 2 2 3 y z y z y z ? ? ? ? ? ? ? ? ? ? ? >>A=[-2 3 ; -1 5 ;1 2]; >> b=[1 2 3]'; 保留方程个数大于变量个数 2. 求出一个最小二乘解 >>A\b 欠定方程组 3 1 3 3 4 4 5 9 8 2 x y z w x y z w x y z w ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1. 解的判断 >>A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 8]; >> b=[1 4 2]'; 方程数小于变量个数 2. 求出一个特解 >>A\b 求解线性方程组的数学原理 稳定、高效的线性方程组解法 ? 克莱姆法则 ? 高斯消去法 ? 高斯列主元素消元法 ? 直接法：用有限次计算求出 “准确解” ? 迭代法：以一个猜测的初始值 开始逐步迭代求满足精度的近 似解 ? 雅可比迭代 ? 高斯 - 赛德尔迭代 克莱姆法则 由 n 个线性方程组联立构成的方程组： AX=R 克莱姆法则解方程 解为： 如三阶线性方程组： 克莱姆法则的计算量 ? 计算工作量巨大，不实用 ? 对于 n 阶方程组，采用克莱姆法则直接求解需要 （ n-1 ）× n ！次乘法和加法。如果 n=20 ，则需计算 4.62 × 10 19 ，如果计算机每秒计算 10 亿次 ，需计算 146 年！ ? 通常采用其他方法，如高斯法等 矩阵操作与线性方程组求解 Matlab 第三讲 11 1 12 2 13 3 1 21 1 22 2 23 3 2 31 1 32 2 33 3 3 a x a x a x b a x a x a x b a x a x a x b ? ? ? ? ? ? ? ? ? ? ? ? ? ? Ax