24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3692  |  回复: 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的回帖

文刀星

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖
查看全部 4 个回答

文刀星

木虫 (小有名气)

还是好好检查你的拟合方程有没有错吧!这个函数时根据你给的初始值和Xdata计算初始Y值,你给的初始值和Xdata代入你写的方程会发现,得到的是NAN,可能是你的方程出现除以0的情况,所以还是好好检查你的拟合方程,看看能不能变个型之类的避免出现被0除的情况!好像你的拟合方程中有(1-a)做除数的情况,当然连初始化都不行,更不用说拟合了!
2楼2015-12-25 20:10:49
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 一志愿哈工大 材料与化工 +38 洛神哥哥 2026-03-31 40/2000 2026-04-03 18:00 by Jimmyandyou
[考研] 一志愿0817化学工程与技术,求调剂 +17 我不是只因 2026-04-02 19/950 2026-04-03 17:08 by lijunpoly
[考研] 调剂 +6 JLLLLLLLLLL 2026-04-03 6/300 2026-04-03 15:47 by 晴空210210
[考研] 368求调剂 +4 今华习 2026-04-03 6/300 2026-04-03 13:46 by 今华习
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 312求调剂 +6 小小墨123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
[考研] 302求调剂一志愿华中师范大学 +7 小江小江江江 2026-04-02 7/350 2026-04-02 21:48 by macy2011
[考研] 环境科学与工程334分求调剂 +7 王一一依依 2026-03-30 9/450 2026-04-02 21:15 by 1104338198
[考研] +4 雾与海 2026-04-02 5/250 2026-04-02 19:16 by 土木硕士招生
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[考研] 【求调剂】新能源材料本科,一志愿211,初试321 +6 求调剂学校, 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 各位老师好,我的一志愿为北京科技大学085601材料专硕 +13 Koxui 2026-03-28 13/650 2026-04-02 09:35 by 哒哒哒呱呱呱
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 350求调剂 +7 阿佳~ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 352分-085602-一志愿985 +6 海纳百川Ly 2026-03-29 6/300 2026-03-31 21:06 by yuq
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
信息提示
请填处理意见