- 1、本文档共21页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
矩阵与数值分析上机实验题及程序
1.给定n阶方程组,其中
,
则方程组有解。对和,分别用Gauss消去法和列主元消去法解方程组,并比较计算结果。
Gauss消去法:
Matlab编程function x=GS(n)
A=[];b=[];
for i=1:n-1
A(i,i)=6;
A(i,i+1)=1;
A(i+1,i)=8;
b(i)=15;
end
A(n,n)=6;b(1)=7;b(n)=14;b=b;
for k=1:n-1
for i=k+1:n
m(i,k)=A(i,k)/A(k,k);
A(i,k:n)=A(i,k:n)-m(i,k)*A(k,k:n);
b(i)=b(i)-m(i,k)*b(k);
end
end
b(n)=b(n)/A(n,n);
for i=n-1:-1:1
b(i)=(b(i)-sum(A(i,i+1:n).*b(i+1:n)))/A(i,i);
end
clear x;
x=b;
disp( AX=b的解x是)
end
计算结果:
在matlab命令框里输出GS(10)得:
GS(10)
AX=b的解x是
ans =
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
在matlab命令框里输出GS(84)得:
GS(84)
AX=b的解x是
ans =
1.0e+008 *
0.0000
…
…
…
0.0000
-0.0000
0.0000
-0.0000
0.0000
-0.0000
0.0000
-0.0000
0.0000
-0.0000
0.0000
-0.0000
0.0000
-0.0001
0.0002
-0.0003
0.0007
-0.0013
0.0026
-0.0052
0.0105
-0.0209
0.0419
-0.0836
0.1665
-0.3303
0.6501
-1.2582
2.3487
-4.0263
5.3684
列主元消去法:
Matlab编程function x=GLZX(n)
A=[];b=[];
eps=10^-2;
for i=1:n-1
A(i,i)=6;
A(i,i+1)=1;
A(i+1,i)=8;
b(i)=15;
end
A(n,n)=6;b(1)=7;b(n)=14;b=b;
for k=1:n-1
[mainElement,index]=max(abs(A(k:n,k)));
index=index+k-1;%index
if abs(mainElement)eps
disp(列元素太小!!);
break;
elseif indexk
temp=A(k,:);temp1=b(k);
A(k,:)=A(index,:);b(k)=b(index);
A(index,:)=temp;b(index)=temp1;
end
for i=k+1:n
m(i,k)=A(i,k)/A(k,k);%A(k,k) ;
A(i,k:n)=A(i,k:n)-m(i,k)*A(k,k:n);
b(i)=b(i)-m(i,k)*b(k);
end
end
b(n)=b(n)/A(n,n);
for i=n-1:-1:1
b(i)=(b(i)-sum(A(i,i+1:n).*b(i+1:n)))/A(i,i);
end
clear x;
x=b;
disp(AX=b的解x是)
end
计算结果:
在matlab命令框里输出GLZX(10)得:
GLZX(10)
AX=b的解x是
ans =
1
1
1
1
1
1
1
1
1
1
在matlab命令框里输出GLZX(84)得:
GLZX(84)
AX=b的解x是
ans =
1.0000
1.0000
1.0000
1.0000
1.0000
1.00
您可能关注的文档
- 普儿科值班(远东版).doc
- 普通物理-2008.pdf
- 易研教育:2014金融专硕考研30所重点推荐院校.pdf
- 更换11盘区3煤无极绳绞车钢丝绳的安全技术措施.doc
- 智能过程控制仿真系统.pdf
- 曾雪麟盼圣火点亮“足球之乡”复兴路.doc
- 最具创意婚宴布置的几个完美细节.doc
- 最好的公司福利逐样看.docx
- 最实用最强的禁止客户端上网.doc
- 普通法兰管道式流量计项目商业计划书(2013年融资成功案例范文)-协助企业融资-投资专家免费咨询.pdf
- 2025年休闲食品健康化转型下的市场拓展策略与渠道创新报告[001].docx
- 农村电商服务站品牌建设:2025年运营模式创新与品牌传播报告.docx
- 2025年互联网金融平台合规整改与金融科技监管的合规监督与可持续发展策略.docx
- 2025年短视频平台内容监管与平台社会责任关系研究报告.docx
- 农产品溯源体系在农产品质量安全监管中的信用体系构建与实施策略研究.docx
- 2025年工业互联网平台数据备份与恢复策略创新方案探讨.docx
- 2025年线上法律咨询服务平台法律咨询行业竞争策略研究.docx
- 农业产业集群农产品流通体系构建与发展机制报告.docx
- 智能制造领航者:2025年3D打印技术在规模生产中的应用效果分析报告[001].docx
- 农村金融服务创新与农村金融消费者权益保护法规实施效果分析报告.docx
文档评论(0)