| 5 | 1/1 | 返回列表 |
| 查看: 1003 | 回復(fù): 4 | |||
| 【懸賞金幣】回答本帖問題,作者我無恥又無賴將贈(zèng)送您 20 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
我無恥又無賴木蟲 (著名寫手)
|
[求助]
1stOpt高版本代跑 已有1人參與
|
||
| 1stopt擬合參數(shù)比較多,希望有5.0以上版本的大神幫助運(yùn)算一下,程序見附件,非常感謝! |
木蟲 (著名寫手)
送紅花一朵 |
您好,我剛剛看了下您提出的問題,之前公式是lny=后面的一堆式子,我這里是把y單獨(dú)拿了出來,所以后面多個(gè)了指數(shù)運(yùn)算。第一個(gè)exp()運(yùn)算括號(hào)里的內(nèi)容是后面的全部式子計(jì)算出來的結(jié)果,那個(gè)很大的數(shù)字只是后面計(jì)算的一小部分?赡苁俏夜嚼锢ㄌ(hào)太多讓人不太好分辨 ![]() 發(fā)自小木蟲IOS客戶端 |
鐵桿木蟲 (職業(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對(duì)于隱函數(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),對(duì)應(yīng)x2系列最。278.15)和最大值(323.15)分別為:593368.271762467,493399.16621225,這兩個(gè)數(shù)值還能進(jìn)行指數(shù)計(jì)算嗎?試下exp(493399.16621225),早就計(jì)算溢出了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎211本科同學(xué),過A區(qū)國(guó)家線,A區(qū)非偏遠(yuǎn)一本,交叉學(xué)科課題組 +13 | lisimayy 2026-03-04 | 23/1150 |
|
|---|---|---|---|---|
|
[考研] 材料277分求調(diào)劑 +11 | 飯飯星球 2026-03-04 | 11/550 |
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +19 | February_Feb 2026-03-01 | 21/1050 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 4/200 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
[考研] 11408,學(xué)碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
[考博] 誠(chéng)招農(nóng)業(yè)博士 +3 | 心欣向榮 2026-02-28 | 3/150 |
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|