第6章MATLAB解方程与函数极值精读.ppt

第六章 MATLAB解方程与函数极值 6.1 线性方程组求解 6.2 非线性方程数值求解 6.3 代数方程的符号解 6.4 函数极值 6.1 线性方程组求解-直接解法 1．利用左除运算符的直接解法 Ax=b 解x=A\b 2．利用矩阵的分解求解线性方程组 矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。常见的矩阵分解有LU分解、QR分解、Cholesky分解，以及Schur分解、Hessenberg分解、奇异分解等。 (1) LU分解 LU分解就是将一个矩阵表示为一个交换下三角矩阵和一个上三角矩阵的乘积形式。线性代数中已经证明，只要方阵A是非奇异的，LU分解总是可以进行的。 [L,U]=lu(X) 产生一个上三角阵U和一个变换形式的下三角阵L(行交换)，且满足X=LU。(注:这里的矩阵X必须是方阵) [L,U,P]=lu(X) 产生一个上三角阵U和一个下三角阵L以及一个置换矩阵P，使之满足PX=LU。当然矩阵X同样必须是方阵。 例6-2 用LU分解求解例6-1中的线性方程组。 命令如下： A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-