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

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * 课堂练习 写出用MATLAB求解以下方程组解的命令 把以上方程系数矩阵的第3至第6个元素，赋值给变量c 求以上方程系数矩阵中大于5小于10的元素位置 程序： function demo A=[4,3;5,8] b=[5;18] X=A\b 结果：-0.8235 2.7647 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 矩阵的行列式 矩阵的行列式运算，如 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且非奇异时，此方程称为恰定方程 恰定方程组 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 超定方程组 A=[-2 3 ; -1 5 ;1 2]; b=[1 2 3]’; 保留方程个数大于变量个数 2.求出一个最小二乘解 A\b 欠定方程组 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×1019，如果计算机每秒计算10亿次，需计算146年！ 通常采用其他方法，如高斯法等 高斯消元法 基本思想 将系数矩阵化为上三角矩阵 回代求解方程 基本过程 消去过程 回代过程 用高斯消去法解方程组 对第一列进行消元，(2)-2/0.5*(1)和(3)-6/0.5*(1)，使a21和a31为0。得 用高斯消去法解方程组 对第二列进行消元， (3)-16.43/4.9*(2)，使a32为0。得 高斯消去