| 5 | 1/1 | 返回列表 |
| 查看: 1728 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
tobetheone銀蟲 (小有名氣)
|
[求助]
怎樣同時擬合三條曲線,得到一套參數(shù) 已有1人參與
|
||
|
各位大俠,小弟由于研究需要,要進行參數(shù)估計,遇到了一個難題:怎樣把涉及到2條曲線的同一組參數(shù)估計出來,請各位大俠鼎力相助。 據(jù)我所知,1stopt可以很方便實現(xiàn)這個功能,但是目前手上沒有正版軟件,實現(xiàn)不了計算功能,特將代碼貼此如下: Variables x,y1,y2,y3; Parameters p1,p2,p3,p4,p5,p6; SharedModel; Function y1=p1*exp(-1/2.6*x)+p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)+p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x); y2=-1/2*(-5658500*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*p4-108250*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-95745393*p6^2*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5-95745393*p6^2*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5+108250*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)+102621*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-100997250*p6^2*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4-102621*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-108309*p5*p4*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-9368500*p6*p5*p4*p1*exp(-1/2.6*x)-2668146*p5^2*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)+103595250*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4^2+103651713*p5*p4^2*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-102621000*p5^2*p6*p1*exp(-1/2.6*x)+103595250*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4^2-100997250*p6^2*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4-2816034*p5^2*p4*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-2816034*p5^2*p4*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-2668146*p5^2*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-108309000*p5^2*p4*p1*exp(-1/2.6*x)+93744500*p5*p6^2*p1*exp(-1/2.6*x)+103651713*p5*p4^2*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-5658500*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*p4+98940500*p4*p6^2*p1*exp(-1/2.6*x)+108309*p5*p4*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))/(-433*p6+474*p5)/(102621*p6*p5+108309*p5*p4+108250*p4*p6); y3=1/2*(-14483906*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*p4-108250*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-6875607*p6^2*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5-6875607*p6^2*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5+108250*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)+102621*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-7252750*p6^2*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4-102621*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)-108309*p5*p4*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2)+10250000*p6*p5*p4*p1*exp(-1/2.6*x)-99952854*p5^2*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)+103595250*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4^2+103651713*p5*p4^2*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)+112338000*p5^2*p6*p1*exp(-1/2.6*x)+103595250*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4^2-7252750*p6^2*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p4-105492966*p5^2*p4*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-105492966*p5^2*p4*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-99952854*p5^2*p6*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)+118500000*p5^2*p4*p1*exp(-1/2.6*x)-102621000*p5*p6^2*p1*exp(-1/2.6*x)+103651713*p5*p4^2*p2*exp(-1/1000*(933*p6+957*p4+1000*1/2.6+974*p5-(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)-14483906*p6*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*p5*p4-108250000*p4*p6^2*p1*exp(-1/2.6*x)+108309*p5*p4*p3*exp((-933/1000*p6-957/1000*p4-1/2.6-487/500*p5-1/1000*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))*x)*(870489*p6^2-994238*p4*p6-817484*p6*p5+915849*p4^2-916236*p5*p4+948676*p5^2)^(1/2))/(-433*p6+474*p5)/(102621*p6*p5+108309*p5*p4+108250*p4*p6); Data; 0 6.1616 10.767 3.704 0.05 5.8161 10.003 3.5127 0.1 5.5093 9.2265 3.158 0.15 5.0878 8.3416 2.8289 0.2 4.2967 7.0175 2.2858 0.25 3.7063 6.425 2.0596 0.3 3.3479 5.7608 1.7427 0.35 2.921 5.083 1.4903 0.4 2.6975 4.7799 1.5039 0.5 2.2055 4.1941 1.3195 0.6 1.916 3.6889 1.1094 0.7 1.5855 3.1527 0.96846 0.8 1.5025 3.0746 0.95926 0.9 1.3666 2.7722 0.88038 1.1 1.3352 2.5721 0.83152 1.3 1.1575 2.3093 0.66482 1.5 1.2741 2.5158 0.75388 1.7 1.2544 2.3014 0.59305 1.9 1.0744 2.0033 0.53928 2.2 0.95854 1.7051 0.4712 2.5 0.82586 1.4852 0.41313 2.8 0.7686 1.2639 0.34796 3.1 0.60757 1.1147 0.26026 3.9 0.54272 0.8547 0.21022 4.9 0.30678 0.5852 0.1671 方程比較長,我試過單獨一條擬合是可以在1stopt擬合的,煩請各位大神幫我代算或者給些建議唄? 2.matlab也可以實現(xiàn)相關功能,但是比較復雜,能不能給些建議謝謝 ![]() ![]() ![]() ![]() ![]() ![]() |

鐵桿木蟲 (著名寫手)
專家顧問 (正式寫手)
![]() |
專家經驗: +65 |

專家顧問 (正式寫手)
![]() |
專家經驗: +65 |

銀蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 考研調劑 +3 | 來好運來來來 2026-03-21 | 3/150 |
|
|
[考研] 08工科 320總分 求調劑 +7 | 梨花珞晚風 2026-03-17 | 7/350 |
|
|
[考研] 297求調劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 0703化學297求調劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 279分求調劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 22408 344分 求調劑 一志愿 華電計算機技術 +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 330求調劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 一志愿蘇州大學材料求調劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 290求調劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 一志愿武理材料工程348求調劑 +3 |  ̄^ ̄゜汗 2026-03-19 | 4/200 |
|
|
[考研] 261求B區(qū)調劑,科研經歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 08工學調劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
[考研] 085410人工智能專碩317求調劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 333求調劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 0856求調劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|