第4章节MATLAB解方程与函数极值.pptVIP

第4章 MATLAB解方程与函数极值 4.1 线性方程组求解 4.2 非线性方程数值求解 4.3 常微分方程初值问题的数值解法 4.4 无约束优化问题 4.5 约束优化问题 4.1 线性方程组求解 4.1.1 直接解法 1．利用左除运算符的直接解法 对于线性方程组Ax=b，可以利用左除运算符“\”求解： x=A\b 对于线性方程组xA=b，可以利用右除运算符“/”求解： x=A/b 例4-1：用直接解法求解下列线性方程组。 命令如下： A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; b=[13,-9,6,0]; x=A\b 结果： x = -66.5556 25.6667 -18.7778 26.5556 2．利用矩阵的分解求解线性方程组 矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。常见的矩阵分解有LU分解、QR分解、Cholesky分解，以及Schur分解、Hessenberg分解、奇异分解等。 (1) LU分解 矩阵的LU分解就是将一个矩阵表示为一个交换下三角矩阵和一个上三角矩阵的乘积形式。线性代数中已经证明，只要方阵A是非奇异的，LU分解总