24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3691  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 啊蓝若 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +9 yelck 2026-04-03 10/500 2026-04-03 13:20 by yelck
[考研] 282求调剂 +5 呼吸都是减肥 2026-03-31 5/250 2026-04-03 12:03 by 1753564080
[考研] 生物学硕341求调剂 +4 你笑起来像云朵 2026-04-03 4/200 2026-04-03 10:32 by macy2011
[考研] 化工调剂303分,过四级 +25 栖梧待风 2026-04-02 25/1250 2026-04-03 10:22 by 哦哦嗯哈
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 282求调剂 +3 aaa车辆 2026-04-02 3/150 2026-04-02 21:55 by zllcz
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 282求调剂 +18 ycy1201 2026-04-01 20/1000 2026-04-02 16:39 by liumengping
[考研] 301求调剂 +4 细胞相关蛋白 2026-04-02 8/400 2026-04-02 16:30 by 271179835
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 311求调剂 +16 zchqwer 2026-04-01 18/900 2026-04-02 10:15 by 不吃魚的貓
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 食品学硕362求调剂 +3 xuanxianxian 2026-04-01 3/150 2026-04-01 21:05 by 啊李999
[考研] 085601材料工程找调剂 +20 oatmealR 2026-03-29 21/1050 2026-04-01 21:00 by lijunpoly
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +7 @taotao 2026-03-30 7/350 2026-04-01 14:30 by chenqifeng666
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 抱歉 +3 田洪有 2026-03-30 3/150 2026-03-30 19:11 by 迷糊CCPs
[考研] 085600,专业课化工原理,321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
信息提示
请填处理意见