原创力文档- 1、本文档共14页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
6.2线性方程组求解
线性方程组包括齐次线性方程组和非齐次线性方程组.非齐次线性方程组的通解等于对应的齐次方程的通解加上非齐次方程的一个特解.在MATLAB中，可以用null(A)得到齐次线性方程组的基础解系；可以用inv、rank、null、左除(\)等命令求解非齐次线性方程组.
例6.6求解齐次线性方程组.解方法1先求出系数矩阵A的行最简形矩阵，再求解.clearA=[1221;21-2-2;1-1-4-3];formatratB=rref(A)B=10-2-5/30124/30000
即得与原方程组同解的方程组由此即得写出向量形式，得到通解
方法2先求出齐次线性方程组的基础解系，再求解.clearA=[1221;21-2-2;1-1-4-3];null(A,r)ans=25/3-2-4/31001即得方程组的基础解系得到方程组的通解
例6.7?求解方程组.解首先计算系数矩阵和增广矩阵的秩，判断方程组解的结构，clear;a=[1-11-1;-111-1;2-2-11];b=[1;1;-1];r1=rank(a)%系数矩阵的秩2r2=rank([a,b])?%增广矩阵的秩2计算表明，系数矩阵和增广矩阵的秩都为2，小于变量的个数4，说明原方程组有无穷组解.有几种方法求原方程组的通解.
方法1用rref命令化为行最简形式求解.clear;a=[1-11-1;-111-1;2-2-11];b=[1;1;-1];rref([a,b])ans=??1???-1????0????0????0????0????0????1???-1????1????0????0????0????0????0
由上述行最简形式得到方程组，即可知原方程组的通解为其中为任意常数.
方法2由于非齐次方程的通解等于齐次方程的通解加非齐次方程的一个特解，可以用null命令求对应的齐次方程组的一个基础解系.clear;a=[1-11-1;-111-1;2-2-11];b=[1;1;-1];x0=a\b?%非齐次方程的一个特解x1=null(a,’r’)?%齐次方程的通解结果为x0=????0????0????1????0
x1=10100101原方程组的通解为其中为任意常数.
例6.8求解线性方程组.解该方程组的系数矩阵A是方阵，可以先计算|A|，A=[1-1-1;2-1-3;32-5]；D=det(A)D=3.0000可得系数行列式|A|=3，由克拉默法则可知方程组由唯一解；也可利用逆矩阵求解.
方法1克拉默法则法formatratb=[2;1;1];A1=[b,A(:,[23])];%b代替A中第1列A2=[A(:,1),b,A(:,3)];%b代替A中第2列A3=[A(:,[12]),b];%b代替A中第3列x1=det(A1)/Dx1=17/3x2=det(A2)/Dx2=1/3x3=det(A3)/Dx3=10/3
方法2逆矩阵法x=inv(A)*bx=17/3
文档评论（0）