Cubic+Spline 源代码

x=0:.12:1; y=(x.^2-3*x+5).*exp(-5*x).*sin(x); plot(x,y,x,y,'o') x1=0:.02:1; y0=(x1.^2-3*x1+5).*exp(-5*x1).*sin(x1); y1=interp…

>> x=0:0.5:10; >> y=cos(x); >> y1=cos(x1); >> y2=interp1(x,y,x1,'nearest') %没有制定外插算法，估值都返回NaN >> y2=interp1(x,y,x1,'…

x=0:.12:1; y=(x.^2-3*x+5).*exp(-5*x).*sin(x); plot(x,y,x,y,'o') x1=0:.02:1; y0=(x1.^2-3*x1+5).*exp(-5*x1).*sin(x1); y1=interp…

function v = splinetx(x,y,u) %SPLINETX Textbook spline function. % v = splinetx(x,y,u) finds the piecewise cubic interpolatory %…

x=0:.12:1; y=(x.^2-3*x+5).*exp(-5*x).*sin(x); plot(x,y,x,y,'o') x1=0:.02:1; y0=(x1.^2-3*x1+5).*exp(-5*x1).*sin(x1); y1=interp…

x=0:10;y=sin(x);xi=0:.25:10; % 将插值方法定义为单元数组 strmod={'nearest','linear','spline','cubic'} % 将X轴标识定义为单元数组 strlb=…

%filename c1_1_1.m x=[0:0.5*pi:2*pi 2*pi+0.6*pi:0.4*pi:4*pi]; y=sin(x);xi=0:0.05*pi:4*pi; y_nearest=interp1(x,y,xi,'nearest');…

%例16-3 外插运算方法和误差 x = 0:10; y = sin(x); xi=5:0.25:15; yi=sin(xi); y1=interp1(x,y,xi,'nearest') y1=interp1(x,y,xi,'nearest','e…

x=0:8; y=sin(x); x1=0:.1:8; y1=interp1(x,y,x1,'nearest') ; y2=interp1(x,y,x1,'linear') ; y3=interp1(x,y,x1,'spline') ; y4=in…

function [res] = INTE(J,JJ) % This function is used by Cubic Spline Rayleigh Ritz Method % (Algorithm 11.6). res = JJ-J+3;