西京学院数学软件实验任务书16.docVIP

西京学院数学软件实验任务书 课程名称 数学软件实验 班级 数0901 学号 0912020119 姓名 王震 实验课题 等距节点插值，Hermite插值，分段插值（线性，二次，三次） 实验目的 熟悉等距节点插值，Hermite插值，分段插值（线性，二次，三次） 实验要求 运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成 实验内容 等距节点插值，Hermite插值，分段插值（线性，二次，三次） 成绩 教师 【实验课题】 等距节点插值，Hermite插值，分段插值（线性，二次，三次） 【实验目的】 熟悉等距节点插值，Hermite插值，分段插值（线性，二次，三次） 【实验内容】 1、Hermite插值 设在节点上，，, 满足条件： 这里共有个插值条件，可唯一确定一个次数不超过的多项式，其形式为： 先求出个插值基函数及，每一个基函数都是次多项式，且满足条件： ，， 即：。 2、分段插值 由于升高插值多项式的阶数有时并不能达到提高精度的效果, 所以实际中往往采用分段插值的思想。分段插值的基本思想是将插值区间划分为若干个小区间，然后在每个小区间上做满足一定条件的低阶插值。 设已知节点上的函数值记 则满足： 1、 2、 3、在每个小区间上是线性函数。 【程序】 %%%%%%%%%%% Hermite插值 function [Hc, Hk,wcgs,Cw]= hermite (X,Y,Y1) m=length(X); n=M1;s=0; H=0;q=1;c1=1; L=ones(m,m);G=ones(1,m); for k=1:n+1 V=1; for i=1:n+1 if k~=i s=s+(1/(X(k)-X(i))); V=conv(V,poly(X(i)))/(X(k)-X(i)); end h=poly(X(k)); g=(1-2*h*s); G=g*Y(k)+h*Y1(k); end H=H+conv(G,conv(V,V)); b=poly(X(k));b2=conv(b,b); q=conv(q,b2);t=2*n+2;Hc=H;Hk=poly2sym(H);Q=poly2sym(q); end for i=1:t c1=c1*i; end syms M;wcgs=M*Q/c1; Cw=q/c1; %%%%%%%%%%%%%%%%%%分段插值 function y=lagr1(x0,y0,x) n=length(x0); m=length(x); for i=1:m z=x(i);s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end