原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
追赶法解线性方程组
一．实验目的
1.熟练运用已学计算方法求解方程组
2.加深对计算方法技巧，选择正确的计算方法来求解各种方程组
3.培养使用电子计算机进行科学计算和解决问题的能力
二．实验内容
合理利用追赶法法求解下列方程组：
三．实验环境
Visual C++ C语言
四．实验公式
（1）LU分解：
l11= a11
u12= a12/ l11
lii-1= aii-1
lii= aii- lii-1ui-1i
uii+1= aii+1/ lii i=2,3,…,n
（2）求解：
“追” y1=b1/l11
yi=(bi-lii-1yi-1)/lii i=2,3,…,n
“赶” xn= yn
xi= yi-uii+1xi+1 i=n-1,…,2,1
五．源程序
#includestdio.h
#includemath.h
#includeconio.h
void main()
{
//追赶法解三对角方程组
float u[10],y[10],a[10],b[10],c[10],f[10],l[10],x[10],exi=1e-5f;
int n,i,exit;
printf(输入方程组维数:\n);
scanf(%d,n);
printf(输入三对角线性方程组:按ai,bi,ci,fi的顺序分别输入。\n);
printf(输入ai:\n);
for(i=2;i=n;i++)
scanf(%f,a[i]);
printf(输入bi:\n);
for(i=1;i=n;i++)
scanf(%f,b[i]);
printf(输入ci:\n);
for(i=1;i=n-1;i++)
scanf(%f,c[i]);
printf(输入fi:\n);
for(i=1;i=n;i++)
scanf(%f,f[i]);
u[0]=0;y[0]=0;
for(i=1;i=n;i++)
{
l[i]=b[i]-a[i]*u[i-1];
if(fabs(l[i])=exi) exit=1;
else
{
y[i]=(f[i]-a[i]*y[i-1])/l[i];
if(i!=n)
u[i]=c[i]/l[i];
}
}
x[n]=y[n];
for(i=n-1;i=1;i--)
x[i]=y[i]-u[i]*x[i+1];
if(exit==1)
printf(li=0,停止计算\n);
else
{printf(方程组得解为:\n);
for(i=1;i=n;i++)
printf(%-14f,x[i]);
}
_getch();
}
六．运行结果
文档评论（0）