第四讲 方程求解.pptVIP

线性方程组的求解 1. 求线性方程组的唯一解或特解（第一类问题） 这类问题的求法分为两类：一类主要用于解低阶稠密矩 阵 —— 直接法；另一类是解大型稀疏矩阵 —— 迭代法。 1．利用矩阵除法求线性方程组的特解（或一个解） 方程：AX=b 解法：X=A\b 求方程组下列的解。 A=[5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5]; B=[1 0 0 0 1]; R_A=rank(A) %求秩 X=A\B %求解 用函数rref求解： C=[A,B] %由系数矩阵和常数列构成增广矩阵C R=rref(C) %将C化成行最简行 R = 1.0000 0 0 0 0 2.2662 0 1.0000 0 0 0 -1.7218 0 0 1.0000 0 0 1.0571 0 0 0 1.0000 0 -0.5940 0 0 0 0 1.0000 0.3188 则R的最后一列元素就是所求之解。 求方程组 的一个特解。 解： A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; B=[1 4 0]; X=A\B %由于系数矩阵不满秩，该解法可能存在误差。 X =[ 0 0 -0.5333 0.6000]’（一个特解近似值）。 若用rref求解，则比较精确： A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; B=[1 4 0]; C=[A,B]; %构成增广矩阵 R=rref(C) R = 1.0000 0 -1.5000 0.7500 1.2500 0 1.0000 -1.5000 -1.7500 -0.2500 0 0 0 0 0 由此得解向量X=[1.2500 – 0.2500 0 0]’（一个特解）。 2. 求线性齐次方程组的通解 在Matlab中，函数null用来求解零空间，即满足A·X=0的解空间，实际上是求出解空间的一组基（基础解系）。 格式 z = null % z的列向量为方程组的正交规范基，满足 。 % z的列向量是方程AX=0的有理基 求解方程组的通解： 解： A=[1 2 2 1;2 1 -2 -2;1 -1 -4 -3]; format rat %指定有理式格式输出 B=null(A,r) %求解空间的有理基 运行后显示结果如下： B = 2 5/3 -2 -4/3 1 0 0 1 syms k1 k2 X=k1*B(:,1)+k2*B(:,2) %写出方程组的通解 运行后结果如下： X = [ 2*k1+5/3*k2] [ -2*k1-4/3*k2] [ k1] [ k2] 3. 求非齐次线性方程组的通解 非齐次线性方程组需要先判断方程组是否有解， 若有解，再去求通解。 因此，步骤为： 第一步：判断AX=b是否有解，若有解则进行第二步 第二步：求AX=b的一个特解 第三步：求AX=0的通解 第四步：AX=b的通解= AX=0的通解+AX=b的一个