欢迎来到天天文库
浏览记录
ID:12146328
大小:561.00 KB
页数:17页
时间:2018-07-15
《数值分析实验报告1》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数值分析实验报告姓名学号班级软件0801指导教师实验名称数值分析实验报告开设学期2010年5月评定成绩评定人签字评定日期东北大学软件学院2010年Ø课题一迭代格式的比较一、实验目的1、通过实验进一步了解方程求根的算法;2、认识选择计算格式的重要性;3、掌握迭代算法和精度控制;4、明确迭代收敛性与初值选取的关系。一、程序设计流程图二、源程序代码#include#includeusingnamespacestd;doublex1=1.8;doublex2=-0.3;doublee=0.001;inti=0;voidfun1(doublea
2、){i++;doubleb;b=(3*a+1)/(a*a);if(b-a>100
3、
4、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<100
5、
6、b-a<-100){cout<<"发散"<7、代:"<-e)){cout<<"X(k+1)="<1008、9、b-a<-100){cout<<"发散"<10、-e)){cout<<"X(k+1)="<10011、12、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<13、10014、15、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<16、l;i=0;}elsefun5(b);}b=a-((a*a*a-3*a-1)/(a*a-1))/3;if(b-a>10017、18、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<19、n4;2.有的迭代格式局部收敛,如fun5,在1.8附近收敛于x1,但在-0.3附近发散;3.好的迭带格式整体收敛,如fun6;Ø课题二线性方程组的直接算法一、流程图Gauss顺序消去法Gauss列主元消去法二、源代码和运行结果Gauss顺序消去法:#include#include#includeusingnamespacestd;classCGAUSSSOLVEEQU{private:vector>m_equset;vectorm_answe
7、代:"<-e)){cout<<"X(k+1)="<100
8、
9、b-a<-100){cout<<"发散"<10、-e)){cout<<"X(k+1)="<10011、12、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<13、10014、15、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<16、l;i=0;}elsefun5(b);}b=a-((a*a*a-3*a-1)/(a*a-1))/3;if(b-a>10017、18、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<19、n4;2.有的迭代格式局部收敛,如fun5,在1.8附近收敛于x1,但在-0.3附近发散;3.好的迭带格式整体收敛,如fun6;Ø课题二线性方程组的直接算法一、流程图Gauss顺序消去法Gauss列主元消去法二、源代码和运行结果Gauss顺序消去法:#include#include#includeusingnamespacestd;classCGAUSSSOLVEEQU{private:vector>m_equset;vectorm_answe
10、-e)){cout<<"X(k+1)="<100
11、
12、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<13、10014、15、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<16、l;i=0;}elsefun5(b);}b=a-((a*a*a-3*a-1)/(a*a-1))/3;if(b-a>10017、18、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<19、n4;2.有的迭代格式局部收敛,如fun5,在1.8附近收敛于x1,但在-0.3附近发散;3.好的迭带格式整体收敛,如fun6;Ø课题二线性方程组的直接算法一、流程图Gauss顺序消去法Gauss列主元消去法二、源代码和运行结果Gauss顺序消去法:#include#include#includeusingnamespacestd;classCGAUSSSOLVEEQU{private:vector>m_equset;vectorm_answe
13、100
14、
15、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<16、l;i=0;}elsefun5(b);}b=a-((a*a*a-3*a-1)/(a*a-1))/3;if(b-a>10017、18、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<19、n4;2.有的迭代格式局部收敛,如fun5,在1.8附近收敛于x1,但在-0.3附近发散;3.好的迭带格式整体收敛,如fun6;Ø课题二线性方程组的直接算法一、流程图Gauss顺序消去法Gauss列主元消去法二、源代码和运行结果Gauss顺序消去法:#include#include#includeusingnamespacestd;classCGAUSSSOLVEEQU{private:vector>m_equset;vectorm_answe
16、l;i=0;}elsefun5(b);}b=a-((a*a*a-3*a-1)/(a*a-1))/3;if(b-a>100
17、
18、b-a<-100){cout<<"发散"<-e)){cout<<"X(k+1)="<19、n4;2.有的迭代格式局部收敛,如fun5,在1.8附近收敛于x1,但在-0.3附近发散;3.好的迭带格式整体收敛,如fun6;Ø课题二线性方程组的直接算法一、流程图Gauss顺序消去法Gauss列主元消去法二、源代码和运行结果Gauss顺序消去法:#include#include#includeusingnamespacestd;classCGAUSSSOLVEEQU{private:vector>m_equset;vectorm_answe
19、n4;2.有的迭代格式局部收敛,如fun5,在1.8附近收敛于x1,但在-0.3附近发散;3.好的迭带格式整体收敛,如fun6;Ø课题二线性方程组的直接算法一、流程图Gauss顺序消去法Gauss列主元消去法二、源代码和运行结果Gauss顺序消去法:#include#include#includeusingnamespacestd;classCGAUSSSOLVEEQU{private:vector>m_equset;vectorm_answe
此文档下载收益归作者所有