24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3690  |  回复: 3

啊蓝若

新虫 (小有名气)

[求助] 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] = ...
本人菜鸟初学,求高手指点怎么修改啊?多谢啊
回复此楼

» 猜你喜欢

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

文刀星

木虫 (小有名气)

还是好好检查你的拟合方程有没有错吧!这个函数时根据你给的初始值和Xdata计算初始Y值,你给的初始值和Xdata代入你写的方程会发现,得到的是NAN,可能是你的方程出现除以0的情况,所以还是好好检查你的拟合方程,看看能不能变个型之类的避免出现被0除的情况!好像你的拟合方程中有(1-a)做除数的情况,当然连初始化都不行,更不用说拟合了!
2楼2015-12-25 20:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

文刀星

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
啊蓝若: 金币+15, ★★★很有帮助 2015-12-29 18:01:52
当然,你可以采用不同不使拟合方程出现NAN情况的初值,比如3*[1,1,1,1]。
3楼2015-12-25 20:28:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啊蓝若

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 文刀星 at 2015-12-25 20:28:21
当然,你可以采用不同不使拟合方程出现NAN情况的初值,比如3*。

谢谢,初始值确实设置有问题
4楼2015-12-29 18:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 啊蓝若 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 270求调剂 +4 海苔花生饼 2026-03-28 4/200 2026-04-03 12:15 by xingguangj
[硕博家园] 求老师收留 +6 lllq123 2026-04-03 6/300 2026-04-03 12:06 by cymywx
[考研] 298分 070300求调剂 +11 zwen03 2026-04-02 11/550 2026-04-03 10:19 by macy2011
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 085600 295分求调剂 +19 W55j 2026-03-30 23/1150 2026-04-03 09:53 by 千千运气
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 385分 生物学(071000)求调剂 +6 qf626 2026-04-01 6/300 2026-04-02 10:39 by guoweigw
[考研] 324求调剂 +5 想上学求调 2026-04-01 6/300 2026-04-02 10:16 by sanrepian
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[硕博家园] 考研调剂 +5 骆驼男人 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
[考研] 材料工程085601数二英一335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-03-31 19:07 by Wang200018
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 11408总分309,一志愿东南大学求调剂,不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 085600,专业课化工原理,321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
信息提示
请填处理意见