- 1、本文档共13页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
本科实验报告
课程名称: 计算机数值方法B
实验项目: 实验2 线性方程组的直接解法
实验地点: 致远楼B501
专业班级:软件1229班学号:2012005829
学生姓名: 田亚鹏
指导教师: 王华
2014年 05月 19日
一、实验目的和要求:
培养理论联系实际的思路,训练综合运用已经学过的理论和实际知识去分析和解决实际问题的能力。
帮助学生全面消化已学的相关课程内容,深刻理解计算数值方法课程的内涵,培养使用电子计算机进行科学计算和解决问题的能力。
进行基本技能训练和巩固。使学生得到选择算法、编写程序、分析数值结果、写数值试验报告、课堂讨论等环节的综合训练。
二、实验内容和原理:
合理利用Gauss消元法、LU分解法或追赶法求解下列方程组:
?1 2 3??x
? ?14?
? ?? ? ? ?①?0 1 2??x1?
? ?? ? ? ?
2
3??2 4 1????x?? ??13??
3
?0.3?10?15 59.14 3 1??x? ?59.17?
? 5.291 ?6.130 ?1 2??x1? ?46.78?
②? ??2??? ?
? 11.2
9 5 2??x
? ? 1 ?
? 1 2 1 1??x3? ? 2 ?
4??
4
?4 2 1
5??x
??? ? ? ?
? ??2?
?? ? ? ???8 7 2 10??x1?
?? ? ? ?
?
③ 2 ?
?4 8 3 6??x? ??7?
4? ??x3? ? ?
4
?12
6 11 20??
? ??3?
?2 1 ??x
?1 2 1 ??x1
? ??7?
? ??5?
? ?? 2? ? ?
④? ? ? ? ?? ? ?????(n=5,10,100,…)
? 1 2 1??x
? ??5?
? ??n?1? ? ?
n?? 1 2????x ?? ???5??
n
二、主要仪器设备:
笔记本电脑;软件vc++6.0;
三、操作方法与实验步骤:
1,Guass消元法:
#includeStdio.h#defineN3
voidGauss_eliminant(doublea[N][N+1])
{
doubleb[N];doubles;
inti,j,k;for(k=0;kN;k++)
{
for(i=k;iN-1;i++)
{
for(j=N;j=0;j--)
{
a[i+1][j]=a[i+1][j]-a[k][j]*a[i+1][k]/a[k][k];
}
}
}
printf(消元后得到上三角矩阵:);printf(\n);
for(i=0;iN;i++)
{
for(j=0;j=N;j++)
{
printf(%lf ,a[i][j]);
}
printf(\n);
}
printf(回代求得原方程组的解为:);printf(\n);
for(i=N-1;i=0;i--)
{
s=0;
for(j=N-1;ji;j--)
{
s=a[i][j]*b[j]+s;
}
b[i]=(a[i][N]-s)/a[i][i];
}
for(i=0;iN;i++)
{
printf(x(%d)=%lf\n,i+1,b[i]);
}
}
main()
{
inti,j;
doublea[N][N+1];
printf(输入系数及常数项矩阵:);printf(\n);
for(i=0;iN;i++)
{
for(j=0;j=N;j++)
{
scanf(%lf,a[i][j]);
}
}
Gauss_eliminant(a);getch();
}
2.Gauss_列主元消元法;
#includeStdio.h#defineN2
voidGauss_eliminant(doublea[N][N+1])
{
doubleb[N];doubles;
inti,j,k,m,n;for(k=0;kN;k++)
{
for(m=N-1;m0;m--)
{
for(n=N-m;nm-1;n++)
{
if(a[n][k]a[n+1][k])
{
s=a[n][k];a[n][k]=a[n+1][k];a[n+1][k]=s;
}
}
}
for(i=k;iN-1;i++)
{
for(j=N;j=0;j--)
{
a[i+1][j]=a[i+1][j]-a[k][j]*a[i+1][k]/a[k][k];
}
}
}
printf(消元后得到上三角矩阵:);printf(\n);
for(i=0;iN;i++)
{
for(j=0;j=N;
文档评论(0)