追赶法解线性方程组.doc

  1. 1、本文档共3页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 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)

xingyuxiaxiang + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档