希望哪位MATLAB高手,帮我看看这个程序哪里错了?急。

来源:百度知道 编辑:UC知道 时间:2024/05/21 19:22:29
function S=csfit(x,y,dx0,dxn)
n=length(x)-1;
h=diff(x);
d=diff(y)./h;
a=h(2:n-1);
b=2*(h(1:n-1)+h(2:n));
c=h(2:n);
u=6*diff(d);
b(1)=b(1)-h(1)/2;
u(1)=u(1)-3*(d(1)-dx0);
b(n-1)=b(n-1)-h(n)/2;
u(n-1)=u(n-1)-3*(dxn-d(n));
for k=2:n-1
temp=a(k-1)/b(k-1);
b(k)=b(k)-temp*c(k-1);
u(k)=u(k)-temp*u(k-1);
end
m(n)=u(n-1)/b(n-1);
for k=n-2:-1:1
m(k+1)=(u(k)-c(k)*m(k+2))/b(k);
end
m(1)=3*(d(1)-dx0)/h(1)-m(2)/2;
m(n+1)=3*(dxn-d(n))/h(n)-m(n)/2;
for k=0:n-1
s(k+1,1)=(m(k+2)-m(k+1))/(6*h(k+1));
s(k+1,2)=m(k+1)/2;
s(k+1,3)=d(k+1)-h(k+1)*(2*m(k+1)+m(k+2))/6;
s(k+1,4)=y(k+1);
end

输入数据后:
>> x=[1,12,4,5];y=[1,4,6,4];dx0=1;dxn=0;
>> S=csfit(x,y,dx0,dxn)
Warning: One or more output arguments not assigned during call to 'csfit'.
望哪位高手指出是哪里错了,我是新手。这个语言我才学。这个程序是算三次样条插值的,当然也希望高手给出个更好的算三次样条插值的程序。

大小写的问题
有两种改法
第一:S=csfit(x,y,dx0,dxn) ->s=csfit(x,y,dx0,dxn)
将大写的S改成小写的s
第二:在程序的最后加上这句 S=s;

这样程序就没有问题了

去买个新的