数值计算实验报告

《数值计算实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、2012级6班###（学号）计算机数值方法实验报告成绩册姓名：安元龙学号：2012060501成绩：实验编号实验名称实验学时实验成绩实验一Newton插值2实验二Romberg算法2实验三求拟合算法2实验四弦截断法2实验五Courant求解方程组2实验六超松弛迭代求解方程组2实验七幕法和反幕法2实验八R-K方法求解初值问题2t计算方法与算法实验报告学期：2014至2015第_1_学期2014年月26日课程名称:数值计算方法与算法_专业:信息与计算科学12级LW.实验编号：1实验项口Neton插值多项式指导教师__孙崎怀姓名：安元龙学号：2012060501实验成绩：一、实验目的及要求实验目

2、的：掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。实验要求：1.给出Newton插值算法2.用C语言实现算法二、实验内容三、实验步骤(该部分不够填写.请填写附页)1.算法分析：卜*面用伪码描述Newton插值多项式的算法：Stepl输入插值节点数n，插值点序列｛x(i)，f(i)｝，i=l，2,，n，耍计算的插值点x.Step2形成差商表fori=0tonforj=ntoif(j)=((f(j)-f(j-l)/(x(j)-x(j-l-i));Stcp3置初始值tcmp=l,ncwton=f(0)Step4fori=ltonte

3、mp:(x-x(i-1))*temp/*由temp(k)=(x-x(k-1))*temp(k-1)形成(x~x(0)(x~x(i~l)*/Newton=newton+temp*f(i);Step5输出f(x)的近似数值newton(x)=newton.#include#defineMAX_N20typedefstructtagPOINT{doublex;doubley;}POINT;intmain(){intn;inti，j;POINTpointslMAX_N+l];doublediff[MAX_N+1】;doublex，tmp，newton=0;printf(nln

4、putnvalue:");scanf(’’％dn，&n);if(n>MAX_N){printfC’TheinputnislargerthenMAX_N,pleaseredefinetheMAX_N.Mreturn1;}if(n<=0){printf("Pleaseinputanumberbetween1and%d."，MAX_N);return1;}printf("Nowinputthe(x_i，y_i)，i=O，..”％d:"，n);for(i=0;i<=n;i++)scanf("%lf%lf’，&points[ij.x，&points⑴.y);printf("Nowi叩ut

5、theXvalue:");scanf("%lf*，&x);for(i=0;i<=n;i++)diff[i]=points[i].y;for(i=0;ii;j-){diff[j]=(diff[j]-diff[j-1])/(points[j].x-points[j-1-i].x);}}tmp=1;newton=diff[()l;for(i=0;i

6、TCTC.EXEInputnvalue:2Nowinputthe(x_i，y_i〉，i=0”••:110.190809120.207912130224951NowinputtheXvalue:11.5newton<11.500000>=-28118.647513Inputnvalue:2Nowinputthe=0.199369数值计算方法与算法实验报告学期：2014至2015第_J_学期2014年月26日课程名称：

7、数值计算方法与算法专业:倍息与计算科孕11级5班实验编号：2实验项目Romberg算'法指导教师孙峪怀姓名：安元龙学号：2012060501实验成绩：一、实验目的及要求实验目的：熟练掌握数值积分的基本思想和动手推导积分公式的能力，并写出自已设计的变步长算法，用此算法解决实际问题。主要锻炼分析问题及对数值方法产生的误差在实际执行算法中的变化，体会事前误差和事后误差的不同之处。实验要求.•写出Romberg算法的描述，并用此