三次样条插值的Matlab实现(自然边界和第一边界条件)(共2页).docx

精选优质文档-倾情为你奉上（第一边界条件）源代码：function y=yt1(x0,y0,f_0,f_n,x) _(1)%第一类边界条件下三次样条插值；%xi 所求点；%yi 所求点函数值；%x 已知插值点；%y 已知插值点函数值；%f_0左端点一次导数值；%f_n右端点一次导数值；n = length(x0);z = length(y0);h = zeros(n-1,1);k=zeros(n-2,1);l=zeros(n-2,1);S=2*eye(n);for i=1:n-1 h(i)= x0(i+1)-x0(i);endfor i=1:n-2 k(i)= h(i+1)/(h(i+1)+h(i); l(i)= 1-k(i);end%对于第一种边界条件： k = 1;k; _(2)l = l;1; _