(9)--5.2数学软件数学软件.pptVIP

第五章数值运算1

本章内容5.1多项式5.2线性方程组求解5.3数值微积分5.4插值和拟合5.5微分方程求解2

5.2求解线性方程组5.2.1齐次线性方程组的解法对于齐次线性方程组AX=0而言，可以通过求系数矩阵A的秩来判断解的情况：1、如果系数矩阵的秩=n（方程组中未知数的个数），则方程组只有零解。2、如果系数矩阵的秩n，则方程组有无穷多解。可以利用MATLAB函数null(A)，求它的一个基本解。3

5.2.1齐次线性方程组的解法例5-10用matlab求解方程组A=[1111-3-11;1000110;-200-10-1-2];r=rank(A);%求矩阵A的秩x=null(A)

x的列向量为Ax=0的一个基本解。4

5.2.2非齐次线性方程组的解法对于非齐次线性方程组AX=b而言，则要根据系数矩阵A的秩和增广矩阵B=[Ab]的秩和未知数个数n的关系，才能判断方程组AX=b的解的情况。（1）如果系数矩阵的秩=增广矩阵的秩=n，则方程组有唯一解。（2）如果系数矩阵的秩=增广矩阵的秩n，则方程组有无穷多解。（3）如果系数矩阵的秩增广矩阵的秩，则方程组无解。5

5.2.2非齐次线性方程组的解法求非齐次线性方程组（A*X=b）的通解时，需要先判断方程组是否有解，若有解，再去求通解。求非齐次线性方程组（A*X=b）的通解的步骤为：第一步：判断AX=b是否有解，若有解则进行第二步；第二步：求AX=b的一个特解；第三步：求AX=0的通解；第四步：AX=b的通解为：AX=0的通解加上AX=b的一个特解。6

5.2.2非齐次线性方程组的解法用matlab求解时，求Ax=b对应的齐次方程组Ax=0的通解，可以利用函数null；求Ax=b的特解，根据方程组中方程的个数m和未知数的个数n，可以把方程组Ax=b分为：恰定方程组（m=n），超定方程组（mn），欠定方程组（mn）。（1）m=n，恰定方程组，可以尝试计算精确解；（2）mn,超定方程组，可以尝试计算最小二乘解；（3）mn，欠定方程组，可以尝试计算含有至少m个解的基解。7

5.2.2非齐次线性方程组的解法1、恰定方程组的求特解方程Ax=b（A为非奇异）x=A-1b两种方法:x=inv(A)?b—采用求逆运算解方程x=A\b—采用左除运算解方程若A为奇异矩阵，则A\b给出出错信息8

5.2.2非齐次线性方程组的解法恰定方程组的求特解9例:x1+2x2=82x1+3x2=13A=[12;23]b=[8;13]x=inv(A)*bx=A\b

5.2.2非齐次线性方程组的解法2、超定方程组的求特解——一般求最小二乘解x=A\b——matlab用最小二乘法找一个准确地基本解。10方程个数大于未知量个数的方程组，超定方程一般是不存在解的矛盾方程。例如，如果给定的三点不在一条直线上，我们将无法得到这样一条直线，使得这条直线同时经过给定这三个点。也就是说给定的条件（限制）过于严格，导致解不存在。比较常用的方法是最小二乘法。形象的说，就是在无法完全满足给定的这些条件的情况下，求一个最接近的解。

5.2.2非齐次线性方程组的解法超定方程组的求特解11例:x1+2x2=12x1+3x2=23x1+4x2=3a=[12;23;34]b=[1;2;3]x=a\bx=1.00000.0000

5.2.2非齐次线性方程组的解法3、欠定方程组的求特解当方程数少于未知量个数时（mn）,有无穷多个解存在。matlab可求出两个解：用除法求的解x是具有最多零元素的解基于伪逆pinv求得的是具有最小长度或范数的解。12

5.2.2非齐次线性方程组的解法欠定方程组的求特解13x1+2x2+3x3=12x1+3x2+4x3=2x=a\bx=pinv(a)?bx=x=1.000.8300.330-0.17

5.2.2非齐次线性方程组的解法?14