三次样条matlab程序,含详细注释.docVIP

clear; f = input(请输入函数表达式:f(x) = , s);%注s，表明允许用户输入一个字符串 a = input(请输入区间左端值a：); b = input(请输入区间右端值b：); n = input(请输入区间等分值n：); for i=1:n+1 x(i) = a + (b-a)/n*(i-1); y(i) = eval(subs(f,x(i),x)); end n=n+1; lamda(1)=1; %构造向量 miu(n)=1; %构造向量 h=diff(x); %若X为向量，Y = diff（X）= [X(2)-X(1),X(3)-X(2),...,X(n)-X(n-1)] df=diff(y)/diff(x); d(1)=6*(df(1)-1/2)/h(1); d(n)=6*(0.5*81^-0.5-df(n-1))/h(n-1); for j=2:n-1 lamda(j)=h(j)/(h(j-1)+h(j)); miu(j)=h(j-1)/(h(j-1)+h(j)); d(j)=6*(df(j)-df(j-1))/(h(j-1)+h(j)); end m