24小时热门版块排行榜

返回列表

啊蓝若

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 121.3
散金: 50
帖子: 77
在线: 74.6小时
虫号: 1274322
注册: 2011-04-22
专业: 高分子合成化学

[求助] lsqcurvefit最小二乘提示错误，望高手指点已有1人参与

x =[770,850,860,960,990,1050,100,1060,1090,1090];
y =[0,0.3,0.7,1.0,1.3,1.7,2.0,2.3,0.7,3.0];
f=@(a, x) a(1)./a(2).*(a(3)./(1700+850.*a(1).*(1-a(4))) .*log((1700.*a(4)+a(1).*(1-a(4)).*(x-770))./(a(4).*a(3).*850-a(3).*(x-770)))+1./(a(1).*(1-a(4))) .*log(1700+(a(1).*(1-a(4)).*(x-770))./a(4))-a(3)./(1700+850.*a(1).*(1-a(4)) ).* log(1700/(850.*a(3) ))-1./(a(1)*(1-a(4))) .*log(1700))
[a,res]=lsqcurvefit(f,[1,1,1,1], x,y);
y1=f(a, x);
plot(x,y, x,y1,’x’)

f =

@(a,x)a(1)./a(2).*(a(3)./(1700+850.*a(1).*(1-a(4))).*log((1700.*a(4)+a(1).*(1-a(4)).*(x-770))./(a(4).*a(3).*850-a(3).*(x-770)))+1./(a(1).*(1-a(4))).*log(1700+(a(1).*(1-a(4)).*(x-770))./a(4))-a(3)./(1700+850.*a(1).*(1-a(4))).*log(1700/(850.*a(3)))-1./(a(1)*(1-a(4))).*log(1700))

错误使用 snls (line 47)
Objective function is returning undefined values at initial point. lsqcurvefit cannot continue.

出错 lsqncommon (line 149)
[xC,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msgData]=...

出错 lsqcurvefit (line 252)
[xCurrent,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
本人菜鸟初学，求高手指点怎么修改啊？多谢啊

回复此楼

» 猜你喜欢

1楼 2015-12-23 11:08:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

文刀星

木虫 (小有名气)

应助: 13 (小学生)
金币: 3651.2
散金: 500
红花: 5
帖子: 60
在线: 100.9小时
虫号: 3693767
注册: 2015-02-22
性别: GG
专业: 导航、制导与传感技术

还是好好检查你的拟合方程有没有错吧！这个函数时根据你给的初始值和Xdata计算初始Y值，你给的初始值和Xdata代入你写的方程会发现，得到的是ＮＡＮ，可能是你的方程出现除以０的情况，所以还是好好检查你的拟合方程，看看能不能变个型之类的避免出现被０除的情况！好像你的拟合方程中有（１－ａ）做除数的情况，当然连初始化都不行，更不用说拟合了！