求助：命令出错了，请高手指教。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

程序命令如下，但总是出错，请高手指教，是哪里出了错误，附件为命令的txt格式，谢谢！

>> f=@(a,x) (1/((2*pi)^0.5*a(1)*x))*exp(-1/2*(log(x)-a(2))^2/a(1)^2);
x=1:0.5:19;
y=[3.155 8.621 11.777 9.170 17.110 18.813 22.961 24.887 25.331 22.961 20.739 18.517 15.702 12.888 10.666 9.184 7.999 6.814 6.370 5.481 5.185 5.037...
    4.740 4.148 4.000 3.851 3.881 3.881 3.807 3.985 3.881 3.614 3.244 3.526 2.977 3.348 3.348];
a0=[4.4013 7.4999];
a=lsqcurvefit(f,a0,x,y);
??? Error using ==> mldivide
Matrix dimensions must agree.
Error in ==> @(a,x)(1/((2*pi)^0.5*a(1)*x))*exp(-1/2*(log(x)-a(2))^2/a(1)^2)

Error in ==> lsqcurvefit at 209
            initVals.F = feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:});
Caused by:
    Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue.

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Evaluation Dimensions Valuation Objective dimension 程序 Error agree using

f=@(a,x) (1./((2*pi)^0.5*a(1)*x)).*exp(-1/2*(log(x)-a(2)).^2/a(1)^2);
x=1:0.5:19;
y=[3.155 8.621 11.777 9.170 17.110 18.813 22.961 24.887 25.331 22.961 20.739 18.517 15.702 12.888 10.666 9.184 7.999 6.814 6.370 5.481 5.185 5.037...
    4.740 4.148 4.000 3.851 3.881 3.881 3.807 3.985 3.881 3.614 3.244 3.526 2.977 3.348 3.348];
a0=[4.4013 7.4999];
a=lsqcurvefit(f,a0,x,y)
%%%%%%%%%
Optimization terminated: relative function value
changing by less than OPTIONS.TolFun.
a =
    0.0043    1.3870

epoh 发表于 2011-10-12 10:15
f=@(a,x) (1./((2*pi)^0.5*a(1)*x)).*exp(-1/2*(log(x)-a(2)).^2/a(1)^2);
x=1:0.5:19;
y=[3.155 8.621 1 ...

您好，谢谢您的解答。可是我在拟合的过程中发现，利用拟合所得系数对应的方程绘图与原始数据图差别太大了，具体代码如下，请问是为什么？
>> f=@(a,x) (1./((2*pi)^0.5*a(1)*x)).*exp(-1/2*(log(x)-a(2)).^2/a(1)^2);
x=1:0.5:19;
y=[3.155 8.621 11.777 9.170 17.110 18.813 22.961 24.887 25.331 22.961 20.739 18.517 15.702 12.888 10.666 9.184 7.999 6.814 6.370 5.481 5.185 5.037...
    4.740 4.148 4.000 3.851 3.881 3.881 3.807 3.985 3.881 3.614 3.244 3.526 2.977 3.348 3.348];
a0=[4.4013 7.4999];
a=lsqcurvefit(f,a0,x,y);
y1=(1./((2*pi)^0.5*a(1)*x)).*exp(-1/2*(log(x)-a(2)).^2/a(1)^2);
subplot(211);plot(x,y);
subplot(212);plot(x,y1);

这组数据应该是你随便给的吧

不管用lsqcurvefit.m,lsqnonlin.m,ga.m

得到的都是local solution

除非能给出正确的初始值,以及设定lb,ub

向量运算用点乘撒