主题:急急!!!
>> x=[0 11 29]; %原始数据
y=[6.0,7.2,10.3];
S=csapi(x,y); %三次样条插值函数
fnplt(S) %画出插值图
xp=0:1:29;
yp=fnval(S,xp);
>> f=inline('a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))','a','x')
ff=optimset;ff.TolFun=1e-20;ff.TolX=1e-15;%参数设定
[a,j]=lsqcurvefit(f,1,xp,yp,[],[],ff)
f =
Inline function:
f(a,x) = a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))
??? Error using ==> optim\private\lsqncommon
User supplied expression or inline function ==> a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))
failed with the following error:
Error using ==> inlineeval
Error in inline expression ==> a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))
??? Error using ==> eval
Index exceeds matrix dimensions.
Error in ==> lsqcurvefit at 149
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
哪错了?
y=[6.0,7.2,10.3];
S=csapi(x,y); %三次样条插值函数
fnplt(S) %画出插值图
xp=0:1:29;
yp=fnval(S,xp);
>> f=inline('a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))','a','x')
ff=optimset;ff.TolFun=1e-20;ff.TolX=1e-15;%参数设定
[a,j]=lsqcurvefit(f,1,xp,yp,[],[],ff)
f =
Inline function:
f(a,x) = a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))
??? Error using ==> optim\private\lsqncommon
User supplied expression or inline function ==> a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))
failed with the following error:
Error using ==> inlineeval
Error in inline expression ==> a(1)./(1+(a(1)/6-1)*exp(-a(2)*x))
??? Error using ==> eval
Index exceeds matrix dimensions.
Error in ==> lsqcurvefit at 149
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
哪错了?