MATLAB编程与工程应用——第7章_MATLAB解方程与函数极值.pptVIP

MATLAB解方程与函数极值 MATLAB解方程与函数极值 第7章 MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7.4 函数极值 MATLAB解方程与函数极值 7.1 线性方程组求解 一、直接解法 对于线性方程组Ax=b，可以使用下面方法直接求解： 1．利用左除运算符 x=A\b 2. 利用矩阵的逆 x=inv(A)*b 或 x=A^(-1)*b 例7.1 用直接法求解线性方程组 exp7_1.m MATLAB解方程与函数极值 一、直接解法 3. 利用矩阵分解求解 矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。常见的矩阵分解有LU分解、QR分解、Cholesky分解、 Schur分解、奇异分解等。 对系数矩阵分解后再求解，可大大提高求解速度 (1) LU分解 矩阵的LU分解就是将一个矩阵表示为一个交换下三角矩阵和一个上三角矩阵的乘积形式。自变量A要求是方阵。 [L,U]=lu(A) x=U\(L\b)=inv(L)*inv(U)*b A=L*U，其中L是一个下三角阵，U是一个上三角阵 [L,U,P]=lu(A) x=U\(L\P*b)=inv(L)*inv(U)*P*b P*A=L*U，其中P是置换矩阵 例7.2 用LU分解求解例7.1 exp7_2.m MATLAB解方程与函数