求助:关于lsqcurvefit函数的问题

shuchi2001

function F=myfun(x,xdata)
F=[xdata(:,4)*0.5+(0.005*0.12*x^0.63*550/xdata(:,1))]/[xdata(:,2)+x*xdata(:,3)+0.352*exp(0.0013*(440-xdata(:,1)))+3.665*exp(0.01*(440-xdata(:,1)))+xdata(:,4)*0.5+(0.005*0.12*x^0.63*550/xdata(:,1))];

通过以下的调用
>> ydata = [0.037471 0.034791 0.032611 0.031011 0.029912 0.029635 0.030214 0.030515 0.030613 0.029495 0.026599];
>>xdata = [600 0.245 0.052466152 0.0014;605 0.268 0.076288668 0.0013;610 0.290 0.108099298 0.0013;615 0.300 0.125276835 0.0012;620 0.31 0.107250119 0.0012;625 0.315 0.67450492 0.0011;630 0.32 0.36982112 0.0011;635 0.325 0.020722678 0.0010;640 033 0.13441035 0.0010;645 0.34 0.009481882 0.001;650 0.35 0.007643467 0.0010];
>>x0=40;
>> [x,resnorm] = lsqcurvefit(@myfun,x0,xdata,ydata);
出现问题
Error using ==> c:/matlab6p5p1/toolbox/optim/private/lsqncommon
User supplied function failed with the following error:

Error using ==> +
Matrix dimensions must agree.

Error in ==> C:\MATLAB6p5p1\toolbox\optim\lsqcurvefit.m
On line 129  ==> [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
请高手指教

eight

原帖由 shuchi2001 于 2007-3-12 15:23 发表
function F=myfun(x,xdata)
F=/;

通过以下的调用
>> ydata = ;
>>xdata = ;
>>x0=40;
>>  = lsqcurvefit(@myfun,x0,xdata,ydata);
出现问题
Error using ==> c:/matlab6 ...

请查看*和.*、^和.^之间的区别

shuchi2001

感谢提醒,经过修改,无法得出结果:
Optimization terminated successfully:
First-order optimality less than OPTIONS.TolFun, and no negative/zero curvature detected
请问440-xdata(:,1)这种写法有没有问题?我想用440分别减去数组的第一列的值

eight

原帖由 shuchi2001 于 2007-3-12 15:53 发表
感谢提醒,经过修改,无法得出结果:
Optimization terminated successfully:
First-order optimality less than OPTIONS.TolFun, and no negative/zero curvature detected
请问440-xdata(:,1)这种写法有没有 ...

可以，没有问题

shuchi2001

那为什么会无法得出结果呢?是因为函数过于复杂吗?

eight

原帖由 shuchi2001 于 2007-3-12 16:01 发表
那为什么会无法得出结果呢?是因为函数过于复杂吗?

拟合这东西我不懂，你搜索一下版面吧。

ps：按照你的结果，不是成功吗？Optimization terminated successfully

shuchi2001

啊,看错了,的确有结果,谢谢了:@)