Cubic+Spline 源代码

%func1p2 %求自然边界条件的三次样条 x = [27.7,28,29,30]; y = [4.1,4.3,4.1,3.0]; cs = spline(x,[0 y 0]); xx = linspace(27.7,30,101); plot(…

x0=[0 3 5 7 9 11 12 13 14 15 ]; y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 ]; x=0:0.1:15; y1=lagr1(x0,y0,x); y2=interp1(x0,y0,x…

hours=1:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; h=1:0.1:12; t=interp1(hours,temps,h,'spline'); plot(hours,temps,'+',h,…

%******************************************************** %程序：EX212.M %功能：1997-2006年北京年平均降水量 %********************************…

clear x = 0:2:24; y = [12 9 9 10 18 24 28 27 25 20 18 15 13]; x1 = 13; % y1 = interp1(x, y, x1); % y2 = interp1(x, y, x1, '…

x0=[0 3 5 7 9 11 12 13 14 15 ]; y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 ]; x=0:0.1:15; y1=lagr1(x0,y0,x); y2=interp1(x0,y0,x…

hours=1:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; h=1:0.1:12; t=interp1(hours,temps,h,'spline'); plot(hours,temps,'+',h,…

function shili14 h0=figure('toolbar','none',... 'position',[200 200 500 500],... 'name','shili14'); axis([0 10 0 10]);…

void splie2(x1a,x2a,ya,m,n,y2a) float **y2a,**ya,x1a[],x2a[]; int m,n; { void spline(); int j; for (j=1;j

void splie2(float x1a[], float x2a[], float **ya, int m, int n, float **y2a) { void spline(float x[], float y[], int n, float y…