Matlab做三次样条插值，怎么设定边界条件 第一边界条件first boundary condition第二边界条件second boundary condition周期边界条件periodic boundary condition自然边界条件natural boundary condition
matlab 怎么进行三次样条插值 spline函数可以实现三次样条 插值： x=0：10；y=sin(x)；xx=0：.25：10；yy=spline(x，y，xx)；plot(x，y，'o'，xx，yy) （另外fnplt csapi这两个函数也是三次样条插值函数） 。
三次样条插值的Matlab实现(自然边界和第一边界条件) 最低0.27元开通文库会员，查看完整内容>；原发布者：snowwolfjay（第一边界条件）源代码：functiony=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)；fori=1：n-1h(i)=x0(i+1)-x0(i)；endfori=1：n-2k(i)=h(i+1)/(h(i+1)+h(i))；l(i)=1-k(i)；end%对于第一种边界条件：k=[1；k]；(2)l=[l；1]；(3)%构建系数矩阵S：fori=1：n-1S(i，i+1)=k(i)；S(i+1，i)=l(i)；end%建立均差表：F=zeros(n-1，2)；fori=1：n-1F(i，1)=(y0(i+1)-y0(i))/(x0(i+1)-x0(i))；endD=zeros(n-2，1)；fori=1：n-2F(i，2)=(F(i+1，1)-F(i，1))/(x0(i+2)-x0(i))；D(i，1)=6*F(i，2)；end%构建函数D：d0=6*(F(1，2)-f_0)/h(1)；(4)dn=6*(f_n-F(n-1，2))/h(n-1)；(5)D=[d0；D；dn]；(6)m=S\\D；寻找x所在位置，并求出对应插值：fori=1：length(x)forj=1：n-1if(x(i)=x0(j))y(i)=(m(j)*(x0(j+1)-x(i))^3)/(6*h(j))+.(m(j+1)*(x(i)-x0(j))^3)/(6*h(j))+.(y0(j)-(m(j)
