资源描述:
《数值计算实验3》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、本科实验报告课程名称:计算机数值方法B实验项目:线性方程组的迭代解法实验地点:逸夫楼302专业班级:学号:学生姓名:指导教师:2013年4月20日一、实验目的和要求学习使用雅可比迭代法或高斯-赛德尔迭代法,进行基木技能训练和巩固。二、实验内容和原理使用雅可比迭代法或高斯■赛德尔迭代法对下列方程组进行求解。10%!-x2-2x3=7.2*—兀]+10x9—2花=8.3x2+5心=4.21、雅克比迭代法:设线性方程组:Ax=b的系数矩阵A可逆且主对角元索Qii,a22,・・・,均不为零,令并将A分解成从而线性方程组可写成则有迭代公式D=diag(an,a22,
2、…,ann)A=(A-D)+1)Dx二(I)-A)x+b其中,BlI-LA,f]二Lb。■1i-ln工严)=J-(一工S'『-为JM+bJaay=iy=i+i(心1,2,…加2、高斯-塞德尔迭代兀严=丄(-曾-«11(R+1)一1(nv(k+l)Xn^~anX1一ann缩写为:£2)==于s-乞SW"1〉一土%乂产]auy=ij=i^i2n三、主要仪器设备惠普ProBook6460bVC++运行环境雅可比迭代法#inelude"stdio.h"#include"math,h"#inelude"stdafx.h"#include"iostream"usin
3、gnamespacestd;intmain(){floata[3][3]={{10,-1,-2),{-1,10,—2},{—1,—1,5}},b[3]二{7.2,8.3,4.2};floatx[3]={0,0,0},sum;inti,j,k,n=3;printf(z,ttX[l]ttX[2]ttX[3]");for(k=0;k<8;k++){for(i=0;i<3;i++){sum=0;for(j=0;j4、];}printf(*第%d次迭代:t",k+1);for(i=0;i5、m2;inti,j,k,n=3;for(k=0;k<10;k++){for(i=0;i6、uDocumentsVisualStudio2008ProjectsprojectlDebugproject1^:^:^:^:^:^:^:^:制次次次次次次次次轨1234567第第第第第第第第请1.0430801・1671881.2820541.0931301・1957241・2977711.0991261・1994671・2997191・0998901.1999331.2999651.0999861.1999921.2999961.0999981.1999991.2999991・1000001・2000001・300000雅可比迭代法3BC:
7、Windowssystem32cmd.exe720000^[2]=0・902000,X[3]=1.164400,dl]=1.043080,x[2]=l.167188,x[31=1.282054,dl]=1.093130,x[2]=1.195724,x[31=1.297771,cLl]=1.099126,xL2]=l.199467,xL3]=1.299719,c[1]=1.099890>x[21=1.199933>x[31=1.299965>dl]=l.099986,x[21=l.199992,x[31=1.299996,dl]=l.099998,x[
8、21=l.199999,x[31=1.299999,([1]=1•