- 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析实验1解线性方程组的消元法
PAGE
佛山科学技术学院
实 验 报 告
课程名称 数值分析
实验项目 解线性方程组的消元法
专业班级 姓名 学号
指导教师 成 绩 日 期
实验目的
掌握程序的录入和matlab的使用和操作;
了解影响线性方程组解的精度的因素——方法与问题的性态。
学会Matlab提供的“\”的求解线性方程组。
实验要求
1、按照题目要求完成实验内容;
2、写出相应的Matlab 程序;
3、给出实验结果(可以用表格展示实验结果);
4、分析和讨论实验结果并提出可能的优化实验。
5、写出实验报告。
实验步骤
1、用调试好的程序解决如下问题:
用Gauss列主元素消去法求的解,其中
a),
解:首先 A=[1 -1 2 -1;2 -2 3 -3;1 1 1 0;1 -1 4 3];
b=[-8 -20 -2 4];
再由Gauss列主元素消去法输入 [x det index]=Gauss(A,b)
得x =
-7.0000
3.0000
2.0000
2.0000
最后得x =(-7.0000 3.0000 2.0000 2.0000)’
b) ,
解: 首先 A=[2.0000 -2.0000 3.0000 -3.0000;0.5000 2.0000 -0.5000 1.5000;0.5000 0 2.5000 4.5000;0.5000 0 0.2000 -0.4000];
b=[-14 6 4 4];再由Gauss列主元素消去法输入 [x det index]=Gauss(A,b)
得: X =
25.1613
-16.6129
-22.1290
10.3871
最后得到:x= (25.1613 -16.6129 -22.1290 10.3871)’
2、调用matlab中的“\”解上述算例
解:a)在上述1操作后由x=a\b得x=( -7.0000 3.0000 2.0000 2.0000)’
b) 同理由x=a\b得x=(25.1613 -16.6129 -22.1290 10.3871)’
3、分别用上述程序和matlab中的“\”求线性方程组
当时的数值解,求解精度为
解:当n=10时, %用函数hilb(n)产生矩阵H
H=hilb(10);
b=1:10;
b=b;
x =H\b
sprintf(%.8f\n,x)
得到:-1000 98014 -232857323304119 -121084912 359469840-632329125651135750 -36232494484065880
当n=20时,H=hilb(20);
b=1:n;
b=b;
x =H\b
sprintf(%.8f\t,x)
得到: -549674495265711069 311289655 -1862367427 41687986268820832828 -76086257686 189177306613-221654283869145321340241-215291605101406559902000-155076143081-468652323929537668585562 27477567610-351461886109210995747353-40393197950当n=30时,H=hilb(30);
b=1:n;
b=b;
x =H\b
sprintf(%.8f\t,x)
得到: -5050756876.28
文档评论(0)