原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
班级：071211 学号 姓名：袁玉凡
PAGE1 / NUMPAGES7
MATLAB高斯列主元消去
高斯列主元消去法的算法过程
对于线性方程组AX=b，A为n*n矩阵：
（1）Step1——在增广矩阵（A，b）第一列中找到绝对值最大的元素，将其所在 行与第一行交换，再对（A，b）做初等行变换使原方程组的第一行第一列的系数 为1，且第一列其它系数都为0.
Step2——在增广矩阵（A，b）中的第二列中（从第二行开始）找到绝对值 最大的元素，将其所在行与第二行交换，再对（A，b）做初等行变换使原方程组 的第二行第二列的系数为1，且第二列中第二行以下的系数都为0.
Step3——在增广矩阵（A，b）中的第三列中（从第三行开始）找到绝对值 最大的元素，将其所在行与第三行交换，再对（A，b）做初等行变换使原方程组 的第三行第三列的系数为1，且第三列中第二行以下的系数都为0。
重复此过程…….
（2）如果增广矩阵的秩为n，则可以将增广矩阵的前n列化为主对角线都为1的矩阵； 如果增广矩阵的秩为mn，则可以将增广矩阵的前m行m列矩阵化为主对角线 都是1的矩阵，m行后的系数均为0；如果增广矩阵的前n列的秩不等于增广矩 阵的秩，此时方程组无解。
（3）接下来，通过初等变换把对角线上每个1所对应的列上的其他元素变为0，就可 以得到增广矩阵的最简阶梯型。这时可以容易得到最简阶梯型对于方程组的解， 即得到方程组的解。
算法的流程图
开始
生成6*6和6*1的矩阵A和b，要求A、b均为随机整数，且值介于[-100, 100]之间，且数据的正负大致各半得到方程组的增广矩阵C
C的秩是否与A的秩相等
是 否
C的秩是否等于n
否 是
高斯列主元消元法求解方程组
方程组有无穷解
方程组无解
结束
运行结果
(1)
随机生成增广矩阵(6*7)：
C =
-94 -89 -41 88 83 -57 -14
41 70 -27 -30 -82 93 43
-67 -24 -31 -62 -49 -14 83
18 -21 26 83 -15 56 78
21 -35 32 -43 15 5 -74
54 11 98 10 79 -34 -77
方程组的解为：
x =
-811/542
554/1015
1329/1117
521/628
-1763/1053
-673/4428
检查结果：Right result!
(2)
随机生成增广矩阵(6*7)：
C =
78 -42 24 86 68 -43 8
30 -54 60 -81 11 -97 -44
-92
文档评论（0）