数值分析上机报告c精选 .doc

数值分析上机报告c 实验3 数值积分 实验目的 通过本实验理解数值积分与微分的基本原理。掌握数值积分中常见的复合求积公式的编程实现。 掌握龙贝格算法的基本思路和迭代步骤； 培养编程与上机调试能力。 算法描述 龙贝格算法基本思路 先算出他T0，从而计算 T0 以此类推，在判断| Tn- Tn-1| 龙贝格算法计算步骤 步骤1：输入区间端点a,b，精度控制值e，循环次数M，定义函数 f，取n?1，h?b?a;R1,1?h?f?f?/2 步骤2：for k?2 to M ? {Rk,1??Rk? ? for 2 k?2 ?h?k1,1 ? to 1i?1 1?f?a?2?i?k ? ?k?h/2? ?j2 if} ?R k,j ?RR?k1 k,?j1 ??R e ?k, ?Rj1? ?k1,j1 ?/?4 j?1 ??1 ? R?k,k ,?k1 ?退出循环 步骤3：数据积分近似值Rk,k。 利用Romberg方法计算函数I? ? 10 sinxx dx 实验内容 用龙贝格算法计算： I? 10 ? sinxx 实验步骤 代码 #include #include #define e double f { double y; if { return y=; } else y=sin/x; return y; } void romberg { } int n=1,k=0; double h,T2,S2=0,C2=0,R2=0,T1,C1,S1,R1; h=/2; T2=h*+f); while )e) { R1=R2; T1=T2; S1=S2; C1=C2; double sum=0; int i; for { sum=sum+f*h); } T2=T1/2+sum*h; S2=/3; C2=/15; R2=/63; n=n*2; k++; h=h/2; } cout void main { double a,b; coutab; cout cout cout } 实验结果 实验体会 ************************************************************************************************************************************* 数值分析上机报告 姓名： 学号： 专业： 2013年10月27日 第一章 舍入误差与有效数 设SN?? j?2N 311?。 ，其精确值为1?????j?12?2NN?1? 2 1 111 ，计算SN的通用程序。 ????222 2?13?1N?1 11，计算S的通用程序。 编制按从小到大的顺序SN?1????N N2?12?122?1 编制按从大到小的顺序SN? 按两种顺序分别计算S102，S104，S106，并指出有效位数。 通过本上机题，你明白了什么？ 解： 、题程序见电子版 按从大到小顺序： 按从小到大顺序： S102= 有效位数6位 S104= 有效位数3位 S106= 有效位数3位 S102= 有效位数5位 S104= 有效位数6位 S106= 有效位数6位 通过上述实验数据可以看出此次算法使用从小到大的顺序进行得到的数据相对而言更精确，可以得到这样的启示：在计算数值时，要先分析不同算法对结果的影响，避免大数吃小数的现象，找出能得到更精确的结果的算法。 第二章 Newton迭代法 给定初值x0及容许误差?，编制Newton法解方程f?0根的通用程序。 3 给定方程f?x/3?x? 0，易知其有三个根x1?，x2? 0，x3? ? ? ? 1．由Newton方法的局部收敛