第6章__MATLAB解方程与最优化问题求解.pptVIP

MATLAB程序设计教程（第二版） 刘卫国 主编 中国水利水电出版社 第6章 MATLAB解方程与最优化问题求解 6.1 线性方程组求解 6.1.1 直接解法 1．利用左除运算符的直接解法 对于线性方程组可以利用左除运算符“\”求解： x=A\b； (1) 若A为n×n方阵，自动利用高斯消元法求解，若b 是n×1的列向量，则解x为n×1的列向量,若b是n×m的矩 阵，可得到m个以A为系数矩阵的线性方程组的数值解x (n×m的矩阵),即x(:,j)=A\b(:,j),j=1,2,…,m. (2) 当A不是方阵时， Ax=b称为欠定方程组或超定方程 组，MATLAB会在最小二乘意义下解 例6-1 用直接解法求解下列线性方程组。 2．利用矩阵的分解求解线性方程组 矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。常见的矩阵分解有LU分解、QR分解、Cholesky分解，以及Schur分解、Hessenberg分解、奇异分解等。 优点：运算速度快，节省存储空间 (1) LU分解 矩阵的LU分解就是将一个矩阵表示为一个交换下三角矩阵和一个上三角矩阵的乘积形式。线性代数中已经证明，只要方阵A是非奇异的，LU分解总是可以进行的。 MATLAB提供的lu函数用于对矩