解线性方程组课件.ppt

《解线性方程组课件.ppt》由会员分享，可在线阅读，更多相关《解线性方程组课件.ppt（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于解线性方程组现在学习的是第1页，共38页一、数学理论复习数学理论复习1、线性方程组记为记为 A x=b 其中A=(aij)mn x=(x1,xn),b=(b1,bm)现在学习的是第2页，共38页若秩(A)秩(A,b)，则无解；若秩(A)=秩(A,b)=n,存在唯一解；若秩(A)=秩(A,b)n,存在无穷多解；通解是齐次线性方程组 Ax=0 的基础解系与 Ax=b 的一个特解之和。对于线性方程组 Ax=b：Ax=0 称为齐次的线性方程组现在学习的是第3页，共38页高斯消元法对于线性方程组 Ax=b（A|b）行变换（U|v ）其中U是行简化阶梯形矩阵(1)阶梯形矩阵(2)每行首个非零元素为1,

2、并且该1所在列其 它元素都为0 现在学习的是第4页，共38页2、逆矩阵方阵A称为可逆的，如果存在方阵B，使A B=B A=E，记 B=A-1方阵A可逆的充分必要条件:A0求逆矩阵方法：A-1=A*/|A|这里A*为A的伴随矩阵（A E）行变换（E A-1）现在学习的是第5页，共38页3、特征值与特征向量对于方阵A，若存在数和非零向量x 使 A x=x，则称为A的一个特征值，x 为A 的一个对应于特征值的特征向量。特征值计算归结为:特征多项式|A-E|=0的求根。对应于特征值的特征向量是齐次线性方程组(A-E)x=0的所有非零解现在学习的是第6页，共38页二、使用二、使用MATLAB det 方

3、阵的行列式 diag 对角阵inv 方阵的逆 cond 方阵的条件数trace 方阵的迹 orth 正交规范化rank 矩阵的秩 null 求基础解系rref 矩阵的行最简形eig 特征值与特征向量jordan 约当标准形分解norm 矩阵或向量范数现在学习的是第7页，共38页1、特殊矩阵生成zeros(m,n)生成m行n列的零矩阵;ones(m,n)生成m行n列的元素全为1的阵;eye(n)生成n阶单位矩阵;当A是矩阵,diag(A)返回A的对角线元素构成的向量;当X是向量,diag(X)返回由X的元素构成的对角矩阵；现在学习的是第8页，共38页rand(m,n)生成m行n列0,1上均匀分布

4、随机数矩阵；linspace(x1,x2,n)生成x1与x2间的n维等距行向量,即将x1,x2 n-1等分。2、行列式和逆矩阵det(A)返回方阵A的行列式;inv(A)返回A的逆矩阵。现在学习的是第9页，共38页3、矩阵除法左除法 AB 求解矩阵方程AX=B右除法 B/A 求解矩阵方程XA=B(1)当A为方阵，AB与inv(A)*B基本一致：(2)当A不是方阵，除法将自动检测。若方程组无解,除法给出最小二乘意义上的近似解，即使向量AXB的长度达到最小；若方程组有无穷多解，除法将给出一个具有最多零元素的特解；若为唯一解，除法将给出解。现在学习的是第10页，共38页4、特征值和特征向量D=eig

5、(A)返回方阵A的特征值构成的列向量；V,D=eig(A)返回方阵A的特征值构成的对角阵D和每个特征值对应的特征向量按列构成的矩阵V。其中每个特征向量都是模等于1的向量，并且属于同一特征值的线性无关特征向量已正交化。现在学习的是第11页，共38页例1 解下列方程组A=1 2;3-2;B=1;4;x=AB 求得唯一解A=1 2 1;3-2 1;B=1;4;x=AB 求得一特解现在学习的是第12页，共38页 A=1 2;3-2;1-1;B=1;4;2;x=A B 求得一最小二乘近似解A=1 2;-2-4;B=1;-2;x=AB 不能直接求解A=1 2;-2-4;0 0;B=1;-2;0;x=AB仍

6、可求一近似特解增加方程 0 x+0y=0现在学习的是第13页，共38页例2 线性方程组的通解解 在无穷多解情况下可用三种方法求通解，用rref化为行最简形以后求解；用除法求出一个特解，再用null求得一个齐次组的基础解系；用符号工具箱中的solve求解。现在学习的是第14页，共38页a=1-1 1-1;-1 1 1-1;2-2-1 1;b=1;1;-1;r=rank(a),rank(a,b);x0=ab,xx=null(a);%x0为一特解,xx为对应齐次组的基础解系运行后得:r=(2,2)说明系数矩阵秩和增广矩阵秩相等,自由未知量为4-2=2个0010 x0=-0.7071 0-0.7071

7、 0-0.0000 0.7071-0.0000 0.7071xx=方法一方法一:方程组的解=特解+对应齐次组的通解其中c1和c2为任意实数结果为结果为:现在学习的是第15页，共38页t=1 -1 0 0 00 0 1 -1 10 0 0 0 0a=1-1 1-1;-1 1 1-1;2-2-1 1;b=1;1;-1;r=rank(a),rank(a,b);t=rref(a,b);%此时得出一个行简化阶梯形矩阵 解法二:运行后得:从而知原方程组等价于从而知原方程组等价于现在学习的是第16页，共38页结果为结果为:其中c1和c2为任意实数现在学习的是第17页，共38页例3 判定下列线性方程组是否有解

8、?若有解,求出其解现在学习的是第18页，共38页(1)a=2-2 3;-1 1-2;1-1 1;(2)b=5;3;4;(3)r1=rank(a);r2=rank(a,b)r1 r2无解唯一解(2)a=2-2 3;-1 1-2;2-3 1;b=5;3;0;r1=rank(a);r2=rank(a,b)r1=r2=3x=ab 或x=inv(a)*b 现在学习的是第19页，共38页(3)a=2-2 3;-1 1-2;1-1 1;b=5;3;8;r1=rank(a);r2=rank(a,b)r1=r2=23x0=ab x=null(a1)%运行后得基础解x=(0.7071,0.7071,0)无穷解a1

9、=2-2 3;-1 1-2;1-1 1;0 0 0;b1=5;3;8;0;x1=a1b1;%经运行后可得出一个特解x1=(0,-19,-11)结果为结果为:其中c为任意实数现在学习的是第20页，共38页三、国民经济投入产出分析三、国民经济投入产出分析 设有n个经济部门，xi为部门i的总产出，cij为部门j单位产品对部门i产品的消耗，di为外部对部门i的需求，fj为部门j新创造的价值。那么各经济部门总产出应满足下列关系式：消耗平衡方程组j=1,2,n现在学习的是第21页，共38页令 C=（cij），X=(x1,xn),D=(d1,dn)，F=(f1,fn)则 X=CX+D令 A=EC，E为单位矩

10、阵，则 AX=DC称为直接消耗矩阵，A称为列昂杰夫（Leontief）矩阵。分配平衡方程组i=1,2,n现在学习的是第22页，共38页Y=1,1,1 BY表示各部门的总投入，称为投入向量。新创造价值向量 F=X Y B=CB表示各部门间的投入产出关系，称为投入产出矩阵。现在学习的是第23页，共38页四、实验例题 例4 某地有三个产业，一个煤矿，一个发电厂和一条铁路，开采一元钱的煤，煤矿要支付0.25元的电费及0.25元的运输费;生产一元钱的电力，发电厂要支付0.65元的煤费，0.05元的电费及0.05元的运输费;创收一元钱的运输费,铁路要支付0.55元的煤费和0.10元的电费，在某一周内煤矿接

11、到外地金额50000元定货，发电厂接到外地金额25000元定货，外界对地方铁路没有需求。现在学习的是第24页，共38页解：这是一个投入产出分析问题。设x1为本周内煤矿总产值，x2为电厂总产值,x3为铁路总产值,则问三个企业间一周内总产值多少才能满足自身及外界需求？三个企业间相互支付多少金额？三个企业各创造多少新价值?现在学习的是第25页，共38页直接消耗矩阵C=外界需求向量 D=产出向量X=则原方程为则原方程为(E-C)X=D 投入产出矩阵为 B=C*diag(X)总投入向量 Y=ones(1,3)*B 新创造价值向量 F=X-Y现在学习的是第26页，共38页Matlab程序:C=0 0.65

12、 0.55;0.25 0.05 0.1;0.25 0.05 0;D=50000;25000;0;A=eye(3)-C;X=AD;%总产出矩阵向量B=C*diag(X);%投入产出矩阵Y=ones(1,3)*B;%总投入向量F=X-Y%新创造价值向量现在学习的是第27页，共38页 消耗部门外界需求煤矿电厂铁路生产部门煤矿0365061558250000电厂255222808283325000铁路25522280800新创造的价值 51044140419915总产出1020885616328330投入产出分析表现在学习的是第28页，共38页例4（隐性病遗传）染色体遗传中，后代是从父母体的基因对中各

13、继承一个基因，形成自己的基因型。如果所考虑的遗传特征是由两个基因A和a控制，那么就有三种基因型，上表给出父母基因型的所有可能组合使其后代形成每种基因对的概率。现在学习的是第29页，共38页设金鱼某种遗传病染色体的正常基因为A，不正常基因为a,那么AA，Aa，aa分别表示正常金鱼，隐性患者，显性患者。设初始分布为90%正常金鱼，10%的隐性患者，无显性患者。考虑下列两种配种方案对后代该遗传病基因型分布的影响方案一：同类基因结合，均可繁殖；方案二：显性患者不允许繁殖，隐性患者必须与正常金鱼结合繁殖现在学习的是第30页，共38页解 设初始分布X(1)=(0.9 0.1 0),第n代分布为X(n)=A

14、=B=则 X(n)=An-1X(1)X(n)=Bn-1X(1)分别是 两种情况下第n代的基因型分布AAAaaa现在学习的是第31页，共38页Matlab程序:方案一:A=1 1/4 0;0 1/2 0;0 1/4 1;x=0.9 0.1 0;for i=2:20 x=A*x;endx20=x方案二:clear;B=1 1/2 0;0 1/2 0;0 0 0;y=0.9 0.1 0;for i=2:20 y=B*y;endy20=y现在学习的是第32页，共38页运行程序后得结果x20=(0.9500,0.0000,0.0500)y20=(1.0000,0.0000,0.0000)可见按方案：很多

15、代以后将出现5%的稳定显性患者按方案：很多代以后显性患者将趋于消失方案体现了杂交的优势现在学习的是第33页，共38页补充内容 解的误差分析解的误差分析现在学习的是第34页，共38页u 解的误差分析解的误差分析对于实际问题导出的方程组对于实际问题导出的方程组 Ax=b,系数矩阵系数矩阵A与向量与向量b往往往往带有误差（扰动），下面讨论带有误差（扰动），下面讨论A或或b的微小变化对解的微小变化对解x的影的影响。响。解线性方程解线性方程Ax=b 即求解线性方程组即求解线性方程组例例:可得出解为可得出解为若方程右端变为若方程右端变为,则方程的解变为则方程的解变为 可见可见x对对b的的扰动敏感扰动敏感现

16、在学习的是第35页，共38页从图可以看出，原方从图可以看出，原方程组对应的两条直线程组对应的两条直线（红与黑）交于（红与黑）交于（2，0）点，但由于两直线）点，但由于两直线几近平行，所以当第二几近平行，所以当第二个方程有微小变化（从个方程有微小变化（从2到到2.01）时，交点）时，交点变（变（1，1），变化很），变化很大。大。对对Ax=b=b，如果解，如果解x 对对b 或或A 的扰的扰动敏感，就称动敏感，就称方程组是病态的方程组是病态的，也，也称称系数矩阵系数矩阵A 是病态的是病态的。现在学习的是第36页，共38页为了定量地估计为了定量地估计x对对b或或A的扰动敏感的程度的扰动敏感的程度,需要度量需要度量向量或矩阵向量或矩阵“大小大小”的数量指标。的数量指标。向量范数向量范数或或矩阵范数矩阵范数正是这样的指标，它们分别用正是这样的指标，它们分别用 来表示。来表示。向量范数向量范数:设设 ,范数记作范数记作 常见的向量范数：常见的向量范数：现在学习的是第37页，共38页感感谢谢大大家家观观看看现在学习的是第38页，共38页