MATLAB插值源代码(19个).docVIP

MATLAB插值源代码(19个) 函数名 功能 Language 求已知数据点的拉格朗日插值多项式 Atken 求已知数据点的艾特肯插值多项式 Newton 求已知数据点的均差形式的牛顿插值多项式 Newtonforward 求已知数据点的前向牛顿差分插值多项式 Newtonback 求已知数据点的后向牛顿差分插值多项式 Gauss 求已知数据点的高斯插值多项式 Hermite 求已知数据点的埃尔米特插值多项式 SubHermite 求已知数据点的分段三次埃尔米特插值多项式及其插值点处的值 SecSample 求已知数据点的二次样条插值多项式及其插值点处的值 ThrSample1 求已知数据点的第一类三次样条插值多项式及其插值点处的值 ThrSample2 求已知数据点的第二类三次样条插值多项式及其插值点处的值 ThrSample3 求已知数据点的第三类三次样条插值多项式及其插值点处的值 BSample 求已知数据点的第一类B样条的插值 DCS 用倒差商算法求已知数据点的有理分式形式的插值分式 Neville 用Neville算法求已知数据点的有理分式形式的插值分式 FCZ 用倒差商算法求已知数据点的有理分式形式的插值分式 DL 用双线性插值求已知点的插值 DTL 用二元三点拉格朗日插值求已知点的插值 DH 用分片双三次埃尔米特插值求插值点的z坐标 function f = Atken(x,y,x0) syms t; if(length(x) == length(y)) n = length(x); else disp(x和y的维数不相等！); return; end %检错 y1(1:n) = t; %符号函数数组要赋初值 for(i=1:n-1) for(j=i+1:n) y1(j) = y(j)*(t-x(i))/(x(j)-x(i))+y(i)*(t-x(j))/(x(i)-x(j)); end y = y1; simplify(y1); end if(nargin == 3) f = subs(y1(n),t,x0); %计算插值点的函数值 else simplify(y1(n)); %化简 f = collect(y1(n)); %将插值多项式展开 f = vpa(f,6); %将插值多项式的系数化成6位精度的小数 end function f0 = BSample(a,b,n,y,y_1,y_N,x0) f0 = 0.0; h = (b-a)/n; c = zeros(n+3,1); b = zeros(n+1,1); for i=0:n-1 if(a+i*h=x0) (a+i*h+h=x0) index = i; break; end end %找到x0所在区间 A = diag(4*ones(n+1,1)); I = eye(n+1,n+1); AL = [I(2:n+1,:);zeros(1,n+1)]; AU = [zeros(1,n+1);I(1:n,:)]; A = A+AL+AU; %形成系数矩阵 for i=2:n b(i,1) = 6*y(i); end b(1) = 6*y(1)+2*h*y_1; b(n+1) = 6*y(n+1)-2*h*y_N; d = followup(A,b); %用追赶法求出系数 c(2:n+2) = d; c(1) = c(2) - 2*h*y_1; %c(-1) c(n+3) = c(3)+2*h*y_N; %c(n+1) x1 = (a+index*h-h-x0)/h; m1 = c(index+1)*(-((abs(x1)