- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Matlab求解方程组.doc
在求解线性方程组时,会遇到以下几种情形:定解方程组、不定方程组、超定方程组、奇异方程组。作为示例,首先以定解线性方程组为例:
在分析如上方程组时,需要知道,方程中有3个未知数,而方程也有3个,所以可以求出(x,y,z)值,转化为矩阵即为: AX = B,其中A为系数矩阵,B为右边值向量。而X即为未知数构成的向量,转化后即为:
A = [2,3,1;
? ? ? ? ? ? ?4,2,3;
? ? ? ? ? ? ?7,1,-1];
如上为系数矩阵;
B = [4;17;1];
如上为右边值矩阵;
利用矩阵除法:
X = A\B?
求得结果如下图红色箭头所示:
x = 1.0000 ; y = -1.0000 ; z = 5.0000;
不定方程组求解
在不定方程组求解时,遇到的方程组常如下所示:
数学上分析可以知道,未知数多于方程式数目,所以解有无数个。但是,可以利用matlab求解一个特定的解(特定解),如下所示输入:
A = [4,5,1;
? ? ? ? ? ? ?1,2,4];
B = [3;15];
X = A\B
可以求得一个特解,如下所示:
超定方程组求解
数学分析上可以知道,当方程数目多于未知数数目时,可以知道该方程组无法求出准确解。如下方程所示:
在求解时,如下操作即可:
A = [4,5;
? ? ? ? ? ? ?1,2;
? ? ? ? ? ? ?3,1];
B = [3;15;12];
X = A\B
求解输出如下图所示,需要说明时,求得结果是以一最小二乘近似解。
在求解奇异方程组,可以发现多个方程之间有重复,如下示例:
在分析时,如上述求解方式:
A = [2,3;
? ? ? ? ? ? -4,-6];
B = [1;-2];
X = A\B
那么,可以看到matlab不能求出解,如下所示:
此时,可以做同解异构,如下所示:
A = [2,3;
? ? ? ? ? ? -4,-6;
? ? ? ? ? ? ?0,0];
B = [1;-2;0];
X = A\B
输出的一个特解如下所示:
HYPERLINK /album/636f38bb3caa88d6b846109d.html?picindex=9 \t _self 步骤阅读
文档评论(0)