数值分析报告-迭代格式的比较.docVIP

《数值分析》上机实验报告 课题：迭代格式的比较 姓 名： 周生刚 学 号： 2006202033 专 业： 材料物理与化学 授课教师： 胡 杰 昆 明 理 工 大 学 06工科硕士 《数值分析》上机实验报告 专业： 材料物理与化学 姓名： 周生刚 学号： 2006202033 任课教师： 胡 杰 作业完成实验室： 实验成绩： 理论描述(10) 数值公式(10) 程序流程图和程序结构(20) 数据和结果(20) 讨论(20) 源程序(20) 总分(100) 实验题目：迭代格式的比较 一、问题提出 设方程f(x)=x- 3x –1=0 有三个实根 x=1.8793 , x=-0.34727 ,x=-1.53209现采用下面六种不同计算格式，求 f(x)=0的根 x 或x 。 x = x = x = x = x = x = x - 要 求 1、编制一个程序进行运算，最后打印出每种迭代格式的敛散情况； 2、用事后误差估计来控制迭代次数，并且打印出迭代的次数； 3、初始值的选取对迭代收敛有何影响； 4、分析迭代收敛和发散的原因。 二、目的和意义 1、通过实验进一步了解方程求根的算法； 2、认识选择计算格式的重要性； 3、掌握迭代算法和精度控制； 4、明确迭代收敛性与初值选取的关系。 本实验通过C语言编程的源程序编写及TC环境运行，完成了对数值分析题目迭代格式的比较的算法选取、程序设计、上机调试和结果分析，并进一步了解方程求根的算法； 三、计算公式 1、 2、 3、 4、 5、 6、 四、结构程序设计 1、程序源代码： #includeiostream #includemath.h int k=1,k1=1,k2=1,k3=1,k4=1,k5=1,k6=1;; float x1,x2,x3,x4,x5,x6; float one(float x0) { x1=(3*x0+1)/(x0*x0); return(x1); } float two(float x0) { x2=(pow(x0,3)-1)/3; return(x2); } float three(float x0) { x3=pow(3*x0+1,0.33333); return(x3); } float four(float x0) { x4=1/(x0*x0-3); return(x4); } float five(float x0) { x5=sqrt(3+1/x0); return(x5); } float six(float x0) { x6=x0-(pow(x0,3)-3*x0-1)/(3*(x0*x0-1)); return(x6); } main() { float x,x0; printf(please input x0=); scanf(%f,x); x0=x; x1=one(x0); printf(第一个公式迭代结果：\n); while(fabs(x0-x1)1e-5) { printf(x1=%6.5f\n,x1); x0=x1; x1=one(x0); k1++; } printf(x1=%6.5f\n,x1); printf(k1=%i\n,k1); x0=x; x2=two(x0); printf(第二个公式迭代结果：\n); while(fabs(x0-x2)1e-5) { printf(x2=%6.5f\n,x2); x0=x2; x2=two(x0); k2++; } printf(x2=%6.5f\n,x2); printf(k2=%i\n,k2); x0=x; x3=three(x0); printf(第三个公式迭代结果：\n); while(fabs(x0-x3)1e-5) { printf(x3=%6.5f\n,x3); x0=x3; x3=three(x0); k3++; } printf(x3=%6.5f\n,x3); printf(k3=%i\n,k3); x0=x; x4=four(x0); printf(第四个公式迭代结果：\