第六章 方程组求根.pdf

Matlab and Engineering Calculation Matlab and Engineering Calculation 第6章 MATLAB解方程与函数极值 6.1 线性方程组求解 6.2 非线性方程数值求解 6.3 函数极值 Xiamen University Xiamen University Matlab and Engineering Calculation Matlab and Engineering Calculation 6.1 线性方程组求解 非齐次线性方程组 + + + a x a x ... a x b ⎧11 1 12 2 1n n 1 ⎪ AX b + + + a x a x ... a x b ⎪ 21 1 22 2 2n n 2 ⎨ 齐次线性方程组 ⎪ ⎪ a x a x ... a x b + + + ⎩m1 1 m 2 2 mn n m AX 0 系数矩阵： A 非齐次项列矢量： b 增广矩阵： B=[A,b] Xiamen University Xiamen University Matlab and Engineering Calculation Matlab and Engineering Calculation 非齐次线性方程组解的判别 1、rank(A)=rank(B)=n, 方程有唯一解； 线性无关方程数等于变量数 2 、rank(A)=rank(B)n, 方程有无穷多解； 线性无关方程数小于变量数 3 、rank(A)rank(B), 方程无解； 含有矛盾方程 4 、齐次线性方程组有非零解的条件：rank(A)n Xiamen University Xiamen University Matlab and Engineering Calculation Matlab and Engineering Calculation 6.1.1 直接解法 1. 利用逆矩阵 AX b =⇒X A−1b inv(A) ⋅b 条件限制：m=n且rank(A)=rank(B)=n或det(A)0 即变量数和方程数相等且方程线性无关 例如 A=magic(3) det(A) b=[1;2;3] X=inv(A)*b ??? mn且rank(A)=rank(B)=n, 即方程数大于变量数但线性 无关的方程数与变量数相等：唯一解