1、第 1 页 共 4 页一、目的1通过本实验加深对拉格朗日插值和牛顿插值法构造过程的理解;2能对上述两种插值法提出正确的算法描述编程实现。二、内容与设计思想自选插值问题,编制一个程序,分别用拉格朗日插值法和牛顿插值法求解某点的函数近似值。 (从课件 jsff04.ppt 或教材习题中选题)已知 y=f(x)的数据表如下,求 t=0.63 处的函数值 z=f(t) 。i xi yi i xi yi1 0.10 0.904837 6 0.57 0.5655252 0.15 0.860708 7 0.70 0.4965853 0.25 0.778801 8 0.85 0.4274154 0.40 0.
2、670320 9 0.93 0.3945545 0.50 0.606531 10 1.00 0.367879三、使用环境操作系统:windons XP软件环境:VC 6.0四、核心代码及调试过程拉格朗日插值法:#include #include #include #define N 20double lagrange(double xN,double yN,double t,int n)double s,ft;int i,j,k,m;ft=0.0;if (nn-1)m=n-1;for (i=k;ivoid main()double x=0.10,0.15,0.25,0.40,0.50,0.57
3、,0.70,0.85,0.93,1.00;double y=0.904837,0.860708,0.77801,0.670320,0.6065310,0.565525,0.496585,0.427418,0.394554,0.367879;double c1010;double b10;double t=0.63,ft;int n,i,j,k;for (i=0;i=1;k-)bk-1=ck-1k-1+bk*(t-xk-1);ft=b0;printf(“f(t=0.63)=%lfn“,ft);运行结果:总结第 4 页 共 4 页1.通过本实验加深了对拉格朗日插值和牛顿插值法构造过程的理解。2.通过编程能对上述两种插值法提出正确的算法。六、附录