1stOpt高版本代跑
1stopt擬合參數(shù)比較多，希望有5.0以上版本的大神幫助運算一下，程序見附件，非常感謝！
返回小木蟲查看更多
今日熱帖
板塊導(dǎo)航
- 網(wǎng)絡(luò)生活
- 資源共享
- 化學(xué)化工
- 專業(yè)學(xué)科
- 科研生活
- 生物醫(yī)藥
- 材料
- 計算模擬
- 學(xué)術(shù)交流
- 出國留學(xué)
- 注冊執(zhí)考
24小時熱帖
換一批
- 下載小木蟲APP
- 與700萬科研達(dá)人隨時交流
-
二維碼
-
IOS
-
安卓
京公網(wǎng)安備 11010802022153號
自己發(fā)完貼都不帶再看一眼的�。扛郊谀膬�？
非常抱歉，當(dāng)時發(fā)完貼以為附件上傳了，沒有仔細(xì)看。我這是個2個因變量一個自變量的非線性擬合，函數(shù)是個隱函數(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這三個參數(shù)取值是0-1，在論壇上看matlab對于隱函數(shù)的擬合好像不是太友好，然后我自己也用低版本的1stopt試了試，并不能出結(jié)果。不清楚高版本的1stopt能不能運行出來，求大佬幫忙指點一下
，
感覺樓主的公式有問題：exp(-(48900/8.314*x2)*(x2/436.45-1)+....，其中“-(48900/8.314*x2)*(x2/436.45-1)”項，對應(yīng)x2系列最�。�278.15）和最大值（323.15）分別為：593368.271762467，493399.16621225，這兩個數(shù)值還能進行指數(shù)計算嗎？試下exp(493399.16621225)，早就計算溢出了。