| 5 | 1/1 | 返回列表 |
| 查看: 998 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題,作者我無恥又無賴將贈(zèng)送您 20 個(gè)金幣 | ||
我無恥又無賴木蟲 (著名寫手)
|
[求助]
1stOpt高版本代跑 已有1人參與
|
|
| 1stopt擬合參數(shù)比較多,希望有5.0以上版本的大神幫助運(yùn)算一下,程序見附件,非常感謝! |
鐵桿木蟲 (職業(yè)作家)
木蟲 (著名寫手)
送紅花一朵 |
非常抱歉,當(dāng)時(shí)發(fā)完貼以為附件上傳了,沒有仔細(xì)看 。我這是個(gè)2個(gè)因變量一個(gè)自變量的非線性擬合,函數(shù)是個(gè)隱函數(shù)。關(guān)系式是y=exp(-(48900/8.314*x2)*(x2/436.45-1)+((a4*(1-y)*x1*exp(-a7*a4)+a5*(1-y)*(1-x1)*exp(-a8*a5))*((1-y)*x1*exp(-a7*a4)+(1-y)*(1-x1)*exp(-a8*a5)))/(y+(1-y)*x1*exp(-a7*a4)+(1-y)*(1-x1)*exp(-a8*a5))^2+((1-y)*x1*exp(-a7*a1)*(a1*(1-y)*x1+a1*(1-y)*(1-x1)*exp(-a9*a6)-a6*(1-y)*(1-x1)*exp(-a9*a6)))/(y*exp(-a7*a1)+(1-y)*x1+(1-y)*(1-x1)*exp(-a9*a6))^2+((1-y)*(1-x1)*exp(-a8*a2)*(a2*(1-y)*(1-x1)+a2*(1-y)*x1*exp(-a9*a3)-a3*(1-y)*x1*exp(-a9*a3)))/(y*exp(-a8*a2)+(1-y)*x1*exp(-a9*a3)+(1-y)*(1-x1))^2) x1,x2是自變量,y是因變量。a1-a9是參數(shù),其中a7,a8,a9這三個(gè)參數(shù)取值是0-1,在論壇上看matlab對于隱函數(shù)的擬合好像不是太友好,然后我自己也用低版本的1stopt試了試,并不能出結(jié)果。不清楚高版本的1stopt能不能運(yùn)行出來,求大佬幫忙指點(diǎn)一下 |
鐵桿木蟲 (職業(yè)作家)
| 感覺樓主的公式有問題:exp(-(48900/8.314*x2)*(x2/436.45-1)+....,其中“-(48900/8.314*x2)*(x2/436.45-1)”項(xiàng),對應(yīng)x2系列最。278.15)和最大值(323.15)分別為:593368.271762467,493399.16621225,這兩個(gè)數(shù)值還能進(jìn)行指數(shù)計(jì)算嗎?試下exp(493399.16621225),早就計(jì)算溢出了。 |
木蟲 (著名寫手)
送紅花一朵 |
您好,我剛剛看了下您提出的問題,之前公式是lny=后面的一堆式子,我這里是把y單獨(dú)拿了出來,所以后面多個(gè)了指數(shù)運(yùn)算。第一個(gè)exp()運(yùn)算括號里的內(nèi)容是后面的全部式子計(jì)算出來的結(jié)果,那個(gè)很大的數(shù)字只是后面計(jì)算的一小部分?赡苁俏夜嚼锢ㄌ柼嘧屓瞬惶梅直 ![]() 發(fā)自小木蟲IOS客戶端 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 哇塞王帥 2026-03-03 | 5/250 |
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
[考研] 0856材料工程,初試313調(diào)劑 +7 | 賣個(gè)關(guān)子吧 2026-03-03 | 7/350 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 環(huán)境工程專碩307求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
[考研] 化學(xué)0703求調(diào)劑 學(xué)碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考博] 26超級電容器申博 +3 | dhdjdjend 2026-02-25 | 3/150 |
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
[考研] 一志愿鄭大材料學(xué)碩298分,求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
[考研] 材料學(xué)調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|