版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

找呀找錢

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 190.3
帖子: 51
在線: 30.4小時(shí)
蟲號(hào): 1612167
注冊(cè): 2012-02-11
專業(yè): 流體力學(xué)

引用回帖:

9樓: Originally posted by optics_ at 2014-08-27 10:21:47
恩恩，1stopt我試了，可以解出來正確的參數(shù)，擬合度也挺高，但是把1stopt解出來的值用matlab再解一遍，matlab就解不出來了，不知道為什么，求指導(dǎo)！...

計(jì)算量大，matlab就死了

贊一下

回復(fù)此樓

11樓2014-08-27 10:23:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

optics_

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10.5
帖子: 22
在線: 8.4小時(shí)
蟲號(hào): 3371702
注冊(cè): 2014-08-18
性別: MM
專業(yè): 激光

引用回帖:

11樓: Originally posted by 找呀找錢 at 2014-08-27 10:23:41
計(jì)算量大，matlab就死了...

謝謝，那matlab做非線性擬合或者解非線性方程有沒有好一點(diǎn)的方法呢？

贊一下

回復(fù)此樓

12樓2014-08-27 11:12:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

NicoleLas

木蟲 (正式寫手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時(shí)
蟲號(hào): 2872321
注冊(cè): 2013-12-15
性別: GG
專業(yè): 計(jì)算機(jī)軟件

引用回帖:

12樓: Originally posted by optics_ at 2014-08-27 11:12:26
謝謝，那matlab做非線性擬合或者解非線性方程有沒有好一點(diǎn)的方法呢？...

lsqcurvefit+MultiStart+Parallel Computing Toolbox，我也在思考這類問題，不過還未真的去做，初值選取是個(gè)麻煩，或許可以以1stOpt求解的結(jié)果作為初值再代入MATLAB迭代，或者就用1stOpt吧，其在這方面還是有優(yōu)勢(shì)的。

贊一下

回復(fù)此樓

Be water, my friend.

13樓2014-08-27 12:54:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

optics_

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10.5
帖子: 22
在線: 8.4小時(shí)
蟲號(hào): 3371702
注冊(cè): 2014-08-18
性別: MM
專業(yè): 激光

引用回帖:

13樓: Originally posted by NicoleLas at 2014-08-27 12:54:31
lsqcurvefit+MultiStart+Parallel Computing Toolbox，我也在思考這類問題，不過還未真的去做，初值選取是個(gè)麻煩，或許可以以1stOpt求解的結(jié)果作為初值再代入MATLAB迭代，或者就用1stOpt吧，其在這方面還是有優(yōu)勢(shì)的 ...

我的想法跟你一樣，用1stopt求解的結(jié)果作為初值再代入MATLAB迭代，換了好幾組不同的數(shù)據(jù)都發(fā)現(xiàn)MATLAB還會(huì)在1stopt求解的結(jié)果的基礎(chǔ)上再迭代3,4次左右，最后得到的最終解也就變了，不懂為什么會(huì)這樣，1stopt的收斂極限是10^-10，而MATLAB的收斂極限是10^-6，不應(yīng)該是1stopt得到的結(jié)果更精確么，為什么MATLAB還要將這個(gè)解迭代？求指導(dǎo)！還有繼續(xù)感謝，Isqcurvefit和nlinfit擬合我都試試剛剛，依然困惑于初值，所以得到的解還是跟正確的解偏差很大，MultiStart和Parallel Computing Toolbox我沒了解過之前，再去看看，謝謝！

贊一下(1人)

回復(fù)此樓

14樓2014-08-27 14:53:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hytao2012

鐵桿木蟲 (正式寫手)

木頭蟲子

應(yīng)助: 53 (初中生)
金幣: 6341.8
散金: 115
紅花: 11
帖子: 479
在線: 206.6小時(shí)
蟲號(hào): 2050091
注冊(cè): 2012-10-08
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

8樓: Originally posted by optics_ at 2014-08-27 10:19:22
恩恩，謝謝，我再去試試！我還想問，如果我每個(gè)參數(shù)值都設(shè)置了擬合上下限，還用設(shè)置初始值么？就是說一般情況下，同樣的擬合上下限不同的初始值，擬合出來的結(jié)果一樣么？...

這個(gè)可以試試，反正也挺簡(jiǎn)單的。不過真的懷疑沒有你要的解。1stOpt解出來的解畫出圖形和數(shù)據(jù)吻合好嗎？

贊一下

回復(fù)此樓

15樓2014-08-27 15:06:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

NicoleLas

木蟲 (正式寫手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時(shí)
蟲號(hào): 2872321
注冊(cè): 2013-12-15
性別: GG
專業(yè): 計(jì)算機(jī)軟件

引用回帖:

14樓: Originally posted by optics_ at 2014-08-27 14:53:35
我的想法跟你一樣，用1stopt求解的結(jié)果作為初值再代入MATLAB迭代，換了好幾組不同的數(shù)據(jù)都發(fā)現(xiàn)MATLAB還會(huì)在1stopt求解的結(jié)果的基礎(chǔ)上再迭代3,4次左右，最后得到的最終解也就變了，不懂為什么會(huì)這樣，1stopt的收斂極 ...

或許MATLAB的算法可能在求解某些問題上效果不太好吧，專業(yè)優(yōu)化軟件一般會(huì)有優(yōu)勢(shì)一些，至于為什么還迭代，這個(gè)我也不很好說，我一般看結(jié)果，R-Square，F(xiàn)-Statistic，SSE什么的，然后殘差圖，以及參數(shù)的Coeff Interval。也許你的擬合模型還可以調(diào)整一下或者因?yàn)槭终`哪個(gè)地方弄錯(cuò)了。
我也不保證MultiStart+Parallel Computing Toolbox能提高結(jié)果精度，其具體操作時(shí)還有一堆策略..............
不過，祝好運(yùn)啦！

贊一下

回復(fù)此樓

Be water, my friend.

16樓2014-08-27 15:28:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

optics_

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10.5
帖子: 22
在線: 8.4小時(shí)
蟲號(hào): 3371702
注冊(cè): 2014-08-18
性別: MM
專業(yè): 激光

引用回帖:

15樓: Originally posted by hytao2012 at 2014-08-27 15:06:41
這個(gè)可以試試，反正也挺簡(jiǎn)單的。不過真的懷疑沒有你要的解。1stOpt解出來的解畫出圖形和數(shù)據(jù)吻合好嗎？...

還是不行，試了很多組初值，R還是負(fù)的，擬合的一點(diǎn)也不好，真心不懂了，還用origin的自定義擬合函數(shù)試了試，擬合出來的R也是負(fù)的，也調(diào)了初值還是不行，然后樓下有人建議用用lsqcurvefit和nlinfit，我試了試解出來的值都和我想要的解差的很多，無語了啊啊啊。1stopt解出來的解還是挺吻合的，圖我就不貼了，下面是報(bào)告。但是現(xiàn)在我都開始懷疑1stopt靠不靠譜了，為什么只有它能解出來。

迭代數(shù): 112
計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:18:30:379
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn)
均方差(RMSE): 0.000261221902678085
殘差平方和(SSE): 0.000142683321179445
相關(guān)系數(shù)(R): 0.999999999997606
相關(guān)系數(shù)之平方(R^2): 0.999999999995212
決定系數(shù)(DC): -2.04562986764419
卡方系數(shù)(Chi-Square): 0.148455468079492
F統(tǒng)計(jì)(F-Statistic): 436303836271944

參數(shù) 最佳估算
---------- -------------
a       0.800337165807115
b       0.000475434654389602
c       0.989996861754167
d       1.00265650022597
e       0.0534333541169663
f       2.34756185300648E-11
g       0.0254246658821726

贊一下

回復(fù)此樓

17樓2014-08-27 18:08:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hytao2012

鐵桿木蟲 (正式寫手)

木頭蟲子

應(yīng)助: 53 (初中生)
金幣: 6341.8
散金: 115
紅花: 11
帖子: 479
在線: 206.6小時(shí)
蟲號(hào): 2050091
注冊(cè): 2012-10-08
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

17樓: Originally posted by optics_ at 2014-08-27 18:08:27
還是不行，試了很多組初值，R還是負(fù)的，擬合的一點(diǎn)也不好，真心不懂了，還用origin的自定義擬合函數(shù)試了試，擬合出來的R也是負(fù)的，也調(diào)了初值還是不行，然后樓下有人建議用用lsqcurvefit和nlinfit，我試了試解出來 ...

不知道為什么，用你給的這組數(shù)據(jù)做出的圖像近似一條直線。而你給的那組數(shù)是一條彎曲的啊

贊一下

回復(fù)此樓

18樓2014-08-27 19:00:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

這是非線性回歸或者方程擬合問題。
要說是求解方程組問題的話，由于你的x~y的數(shù)據(jù)有幾百組，而未知數(shù)僅有幾個(gè)，是典型的方程數(shù)大于未知數(shù)個(gè)數(shù)的問題，極可能出現(xiàn)無解的情況。
所以理解是非線性回歸才是合理的。
推薦用1stopt軟件來解決該問題，該問題參數(shù)個(gè)數(shù)多，初值不易確定。就擬合而言，通常1stopt軟件比MATLAB強(qiáng)大。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

19樓2014-08-27 20:38:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

optics_

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10.5
帖子: 22
在線: 8.4小時(shí)
蟲號(hào): 3371702
注冊(cè): 2014-08-18
性別: MM
專業(yè): 激光

引用回帖:

16樓: Originally posted by NicoleLas at 2014-08-27 15:28:01
或許MATLAB的算法可能在求解某些問題上效果不太好吧，專業(yè)優(yōu)化軟件一般會(huì)有優(yōu)勢(shì)一些，至于為什么還迭代，這個(gè)我也不很好說，我一般看結(jié)果，R-Square，F(xiàn)-Statistic，SSE什么的，然后殘差圖，以及參數(shù)的Coeff Interv ...

嗯��！O(∩_∩)O謝謝。

贊一下

回復(fù)此樓

20樓2014-08-28 10:12:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 optics_ 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085700資源與環(huán)境308求調(diào)劑 +7	墨墨漠 2026-03-18	7/350	2026-03-20 05:55 by laoshidan
[考研] 279分求調(diào)劑一志愿211 +3	chaojifeixia 2026-03-19	3/150	2026-03-20 05:47 by laoshidan
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +3	sbdksD 2026-03-19	3/150	2026-03-19 23:21 by fmesaito
[考研] 288求調(diào)劑 +15	于海海海海 2026-03-19	15/750	2026-03-19 22:41 by 學(xué)員8dgXkO
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	12/600	2026-03-19 16:17 by 脫穎而出
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	10/500	2026-03-19 16:05 by 余麟余
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 268求調(diào)劑 +6	簡(jiǎn)單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[基金申請(qǐng)] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike