插值法的程序实现.doc

插值法的程序实现 专业班级：信息与计算科学061 姓名：袁金龙 学号 一 实验目的Matlab编程； 2.学习插值方法及程序设计算法 二 实验题目 三 实验原理与理论基础… ②由Lk(x)=((x-x0)*…*(x-x(k-1))*(x-x(k+1))…(x-xn))/ ((xk-x0)*…*(xk-x(k-1))*(xk-x(k+1))…(xk-xn))得出Li(x); ③由Y=y1* L1(x)+…+yn*Ln(x)得出Y关于x的表达式。 ④带值计算即可。 2. 牛顿插值算法设计 ①利用已知条件得到xi，yi，i=0，1，2，… ②利用差商公式f[x0,…xk]=(f[x0,…,x(k-2),xk]-f[x0,…,x(k-1)])/(xk-x(k-1))各阶差商。 ③利用牛顿插值公式f(x)=f(x0)-f[x0,x1]*(x-x0)+…f[x0,x1,…xn]*(x-x0)*…(x-x(n-1)). ④带值计算即可。 3. 自然样条函数算法设计 ①利用已知条件得到xi，yi，i=0，1，2，… ②利用公式求出h(i) u(i) k(i);di; ③利用h(i) u(i) k(i);di解出向量M; ④ 将相关变量带入自然样条表达式中即可。 ⑤带值计算即可。 四 实验内容1920年到1970年的人口表： 年份 1920 1930 1940 1950 1960 1970 人口 (千人) 105711 123203 131669 150697 179323 203212 1.用表中数据构造一个5次拉格朗日插值多项式，并用此估计1910，1965和2002年的人口。在1910年的实际人口约请判断插值计算得到的1965年和2002年的人口数据准确性是多少？ 2 用牛顿插值估计： （1）1965年的人口数； （2）2002年的人口数。? 用自然样条函数估计在1910，1965和2002年的人口数。 请比较以上三种方法所求值的效果。那一种方法最优？ （二）实验代码： 1.①用matlab编写的拉格朗日插值M文件如下： function [m,y]=cz5(x) y1=105711;y2=123203;y3=131669;y4=150697;y5=179323;y6=203212; x1=1920; x2=1930; x3=1940; x4=1950; x5=1960; x6=1970; a1=((x-x2)*(x-x3)*(x-x4)*(x-x5)*(x-x6))/((x1-x2)*(x1-x3)*(x1-x4)*(x1-x5)*(x1-x6)); a2=((x-x1)*(x-x3)*(x-x4)*(x-x5)*(x-x6))/((x2-x1)*(x2-x3)*(x2-x4)*(x2-x5)*(x2-x6)); a3=((x-x1)*(x-x2)*(x-x4)*(x-x5)*(x-x6))/((x3-x1)*(x3-x2)*(x3-x4)*(x3-x5)*(x3-x6)); a4=((x-x1)*(x-x2)*(x-x3)*(x-x5)*(x-x6))/((x4-x1)*(x4-x2)*(x4-x3)*(x4-x5)*(x4-x6)); a5=((x-x2)*(x-x3)*(x-x4)*(x-x1)*(x-x6))/((x5-x1)*(x5-x3)*(x5-x2)*(x5-x4)*(x5-x6)); a6=((x-x2)*(x-x3)*(x-x4)*(x-x5)*(x-x1))/((x6-x1)*(x6-x2)*(x6-x3)*(x6-x4)*(x6-x5)); 实际值大约为： y=y1*a1+y2*a2+y3*a3+y4*a4+y5*a5+y6*a6; if x==1910 m=91772; 估计值大约为： else m=(91772-31872)/31872*y+y; 估计值大约为： End ②用c++编写程序如下： #includeiostream.h #includestdlib.h #includestdio.h #includemath.h double H(double x,double a,double p[6]) //求得插值基函数 { double h=1.0,f=1.0; for(int i=0;i6;i