基于Matlab求解线性方程组、非线性方程组.docVIP

Matlab求解线性方程组、非线性方程组 姓名：罗宝晶 学号：1012208015 专业：材料学院高分子系 第一部分 数值计算 Matlab求解线性方程组AX=B或XA=B在MATLAB中，求解线性方程组时，主要采用除法运算符“/”和“\”。如：X=A\B表示求矩阵方程AX＝B的解；X＝B/A表示矩阵方程XA=B的解。对方程组X＝A\B，要求A和B用相同的行数，X和B有相同的列数，它的行数等于矩阵A的列数，方程X＝B/A同理。如果矩阵A不是方阵，其维数是m×n，则有：m＝n 恰定方程，求解精确解；mn 超定方程，寻求最小二乘解；mn 不定方程，寻求基本解，其中至多有m个非零元素。针对不同的情况，MATLAB将采用不同的算法来求解。恰定方程组恰定方程组由n个未知数的n个方程构成，方程有唯一的一组解，其一般形式可用矩阵，向量写成如下形式：Ax=b 其中A是方阵，b是一个列向量；在线性代数中，最常用的方程组解法有：（1）利用ramer公式来求解法；（2）利用矩阵求逆解法，即x=A-1b；（3）利用aussian消去法；（4）利用u法求解。一般来说，对维数不高，条件数不大的矩阵，上面四种解法所得的结果差别不大。前三种解法的真正意义是在其理论上，而不是实际的数值计算。MATLAB中，出于对算法稳定性的考虑，行列式及逆的计算大都在u分解的基础上进行。在MATLAB中，求解这类方程组的命令