Cubic+Spline 源代码

function r=thp(x,y) m=length(x); h=0.1; for k=1:m s1=spline(x,y,x(k)+2*h); s2=spline(x,y,x(k)+h); s3=s…

function r=nt1(x,y) m=length(x); h=0.1; for i=1:(m-1) t=x(i)+h; s=spline(x,y,t); r(i)=(s-y(i))/h…

clear all; m=1000;n = 11; x = -1:2/(m-1):1; y = 1./(1+9*x.^2); x0 = -1:2/(n-1):1; y0 = 1./(1+9*x0.^2); y1 = lagr1(x0,y0,x); …

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.5:15; y1=lagr1(x0,y0,x); y2=interp1(x0,y0,x…

n=11;m=21; x=-5:10/(m-1):5; y=1./(1+x.^2);z=0*x;plot(x,z,x,y),gtext('y=1/(1+x^2)'),pause x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y…

function [Kp,Ti,Td,key,lb]=pid_hv_tune(key1,key2,vars) key=0; Kp=[]; Ti=[]; Td=[]; lb=[]; K=vars(1); L=vars(2); T=vars(3); sw…

Ts=0.001;n=-5:1:5;nTs=n*Ts; x=exp(-1000*abs(nTs)); Dt=0.00005;t=-0.005:Dt:0.005; xa=spline(nTs,x,t); error=max(abs(xa-exp(-100…

Ts=0.0002;n=-25:1:25;nTs=n*Ts; x=exp(-1000*abs(nTs)); Dt=0.00005;t=-0.005:Dt:0.005; xa=spline(nTs,x,t); error=max(abs(xa-exp(-…

function [Kp,Ti,Td,key,lb]=pid_hv_tune(key1,key2,vars) key=0; Kp=[]; Ti=[]; Td=[]; lb=[]; K=vars(1); L=vars(2); T=vars(3); sw…

x0=0:pi/30:pi; y0=sin(x0); I=quadspln(x0,y0,0,pi) x0=0:pi/10:pi; y0=sin(x0); I1=trapz(x0,y0) I=quadspln(x0,y0,0,pi) x0=[0…