| 查看: 1650 | 回復(fù): 16 | ||||
hanmei98232木蟲 (小有名氣)
|
[求助]
求助:以下非線性擬合方程可以用1stopt求解么?用什么方法
|
|
以下非線性擬合方程可以用1stopt求解么?用什么方法 0.07272703=-(m*0.17809/(n*0.17809^0.5+1))*8.314*298.15*(ln(0.17809)-0.784*(1.7321*0.17809^0.5/(1+2.0785*0.17809^0.5)+1.6667*ln(1+2.0785*0.17809^0.5))+2.6667*0.17809*(c+d*(1+(-1-3.464*0.17809^0.5+6*0.17809)*e0.17809p(-3.464*0.17809^0.5))/(12*0.17809))+2.8289*0.17809^2*e)+0.07214; 0.07319177=-(m*0.36061/(n*0.36061^0.5+1))*8.314*298.15*(ln(0.36061)-0.784*(1.7321*0.36061^0.5/(1+2.0785*0.36061^0.5)+1.6667*ln(1+2.0785*0.36061^0.5))+2.6667*0.36061*(c+d*(1+(-1-3.464*0.36061^0.5+6*0.36061)*e0.36061p(-3.464*0.36061^0.5))/(12*0.36061))+2.8289*0.36061^2*e)+0.07214; 0.07391741=-(m*0.53672/(n*0.53672^0.5+1))*8.314*298.15*(ln(0.53672)-0.784*(1.7321*0.53672^0.5/(1+2.0785*0.53672^0.5)+1.6667*ln(1+2.0785*0.53672^0.5))+2.6667*0.53672*(c+d*(1+(-1-3.464*0.53672^0.5+6*0.53672)*e0.53672p(-3.464*0.53672^0.5))/(12*0.53672))+2.8289*0.53672^2*e)+0.07214; 0.07451667=-(m*0.72049/(n*0.72049^0.5+1))*8.314*298.15*(ln(0.72049)-0.784*(1.7321*0.72049^0.5/(1+2.0785*0.72049^0.5)+1.6667*ln(1+2.0785*0.72049^0.5))+2.6667*0.72049*(c+d*(1+(-1-3.464*0.72049^0.5+6*0.72049)*e0.72049p(-3.464*0.72049^0.5))/(12*0.72049))+2.8289*0.72049^2*e)+0.07214; 0.07521378=-(m*0.92046/(n*0.92046^0.5+1))*8.314*298.15*(ln(0.92046)-0.784*(1.7321*0.92046^0.5/(1+2.0785*0.92046^0.5)+1.6667*ln(1+2.0785*0.92046^0.5))+2.6667*0.92046*(c+d*(1+(-1-3.464*0.92046^0.5+6*0.92046)*e0.92046p(-3.464*0.92046^0.5))/(12*0.92046))+2.8289*0.92046^2*e)+0.07214; 0.07661342=(m*1.3003/(n*1.3003^0.5+1))*8.314*298.15*(ln(1.3003)-0.784*(1.7321*1.3003^0.5/(1+2.0785*1.3003^0.5)+1.6667*ln(1+2.0785*1.3003^0.5))+2.6667*1.3003*(c+d*(1+(-1-3.464*1.3003^0.5+6*1.3003)*e1.3003p(-3.464*1.3003^0.5))/(12*1.3003))+2.8289*1.3003^2*e)+0.07214; 0.07871152=(m*1.67471/(n*1.67471^0.5+1))*8.314*298.15*(ln(1.67471)-0.784*(1.7321*1.67471^0.5/(1+2.0785*1.67471^0.5)+1.6667*ln(1+2.0785*1.67471^0.5))+2.6667*1.67471*(c+d*(1+(-1-3.464*1.67471^0.5+6*1.67471)*e1.67471p(-3.464*1.67471^0.5))/(12*1.67471))+2.8289*1.67471^2*e)+0.07214; 0.07934748=(m*1.86291/(n*1.86291^0.5+1))*8.314*298.15*(ln(1.86291)-0.784*(1.7321*1.86291^0.5/(1+2.0785*1.86291^0.5)+1.6667*ln(1+2.0785*1.86291^0.5))+2.6667*1.86291*(c+d*(1+(-1-3.464*1.86291^0.5+6*1.86291)*e1.86291p(-3.464*1.86291^0.5))/(12*1.86291))+2.8289*1.86291^2*e)+0.07214; 0.08119555=(m*2.42891/(n*2.42891^0.5+1))*8.314*298.15*(ln(2.42891)-0.784*(1.7321*2.42891^0.5/(1+2.0785*2.42891^0.5)+1.6667*ln(1+2.0785*2.42891^0.5))+2.6667*2.42891*(c+d*(1+(-1-3.464*2.42891^0.5+6*2.42891)*e2.42891p(-3.464*2.42891^0.5))/(12*2.42891))+2.8289*2.42891^2*e)+0.07214; 0.08388612=(m*2.94245/(n*2.94245^0.5+1))*8.314*298.15*(ln(2.94245)-0.784*(1.7321*2.94245^0.5/(1+2.0785*2.94245^0.5)+1.6667*ln(1+2.0785*2.94245^0.5))+2.6667*2.94245*(c+d*(1+(-1-3.464*2.94245^0.5+6*2.94245)*e2.94245p(-3.464*2.94245^0.5))/(12*2.94245))+2.8289*2.94245^2*e)+0.07214; 0.08544611=(m*3.49978/(n*3.49978^0.5+1))*8.314*298.15*(ln(3.49978)-0.784*(1.7321*3.49978^0.5/(1+2.0785*3.49978^0.5)+1.6667*ln(1+2.0785*3.49978^0.5))+2.6667*3.49978*(c+d*(1+(-1-3.464*3.49978^0.5+6*3.49978)*e3.49978p(-3.464*3.49978^0.5))/(12*3.49978))+2.8289*3.49978^2*e)+0.07214; 0.08798993=(m*4.05405/(n*4.05405^0.5+1))*8.314*298.15*(ln(4.05405)-0.784*(1.7321*4.05405^0.5/(1+2.0785*4.05405^0.5)+1.6667*ln(1+2.0785*4.05405^0.5))+2.6667*4.05405*(c+d*(1+(-1-3.464*4.05405^0.5+6*4.05405)*e4.05405p(-3.464*4.05405^0.5))/(12*4.05405))+2.8289*4.05405^2*e)+0.07214; 即x<1 function y=-(m*X/(n*X^0.5+1))*8.314*298.15*(ln(x)-0.784*(1.7321*x^0.5/(1+2.0785*x^0.5)+1.6667*ln(1+2.0785*x^0.5))+2.6667*x*(c+d*(1+(-1-3.464*x^0.5+6*x)*exp(-3.464*x^0.5))/(12*x))+2.8289*x^2*e)+0.07214; x>1 y=(m*X/(n*X^0.5+1))*8.314*298.15*(ln(x)-0.784*(1.7321*x^0.5/(1+2.0785*x^0.5)+1.6667*ln(1+2.0785*x^0.5))+2.6667*x*(c+d*(1+(-1-3.464*x^0.5+6*x)*exp(-3.464*x^0.5))/(12*x))+2.8289*x^2*e)+0.07214; data; 0.17809 0.07272703 0.36061 0.07319177 0.53672 0.07391741 0.72049 0.07451667 0.92046 0.07521378 1.3003 0.07661342 1.67471 0.07871152 1.86291 0.07934748 2.42891 0.08119555 2.94245 0.08388612 3.49978 0.08544611 4.05405 0.08798993 如果1stopt解不了用matlab應(yīng)該怎么解? 求助求助。。。》浅8兄x Parameters c,d,e; Variable x1,x2,y1,y2 ; function y1=(m*X1/(n*X1^0.5+1))*8.314*298.15*(ln(x1)-0.784*(1.7321*x1^0.5/(1+2.0785*x1^0.5)+1.6667*ln(1+2.0785*x1^0.5))+2.6667*x1*(c+d*(1+(-1-3.464*x1^0.5+6*x1)*exp(-3.464*x1^0.5))/(12*x1))+2.8289*x1^2*e)+0.07214; y2=-(m*X2/(n*X2^0.5+1))*8.314*298.15*(ln(x2)-0.784*(1.7321*x2^0.5/(1+2.0785*x2^0.5)+1.6667*ln(1+2.0785*x2^0.5))+2.6667*x2*(c+d*(1+(-1-3.464*x2^0.5+6*x2)*exp(-3.464*x2^0.5))/(12*x2))+2.8289*x2^2*e)+0.07214; Data; //x1 y1 x2 y2 0.17809 0.07272703 1.3003 0.07661342 0.36061 0.07319177 1.67471 0.07871152 0.53672 0.07391741 1.86291 0.07934748 0.72049 0.07451667 2.42891 0.08119555 0.92046 0.07521378 2.94245 0.08388612 3.49978 0.08544611 4.05405 0.08798993 求教 這個(gè)樣編寫對(duì)么?為什么出來(lái)的圖和y對(duì)不上呢? [ Last edited by hanmei98232 on 2014-3-15 at 13:46 ] |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
| 1stopt軟件擬合、解方程的標(biāo)準(zhǔn)代碼參見(jiàn)附件。 |

木蟲 (小有名氣)
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
|
0.07272703=-(m*0.17809/(n*0.17809^0.5+1))*8.314*298.15*(ln(0.17809)-0.784*(1.7321*0.17809^0.5/(1+2.0785*0.17809^0.5)+1.6667*ln(1+2.0785*0.17809^0.5))+2.6667*0.17809*(c+d*(1+(-1-3.464*0.17809^0.5+6*0.17809)*e0.17809p(-3.464*0.17809^0.5))/(12*0.17809))+2.8289*0.17809^2*e)+0.07214; 0.07319177=-(m*0.36061/(n*0.36061^0.5+1))*8.314*298.15*(ln(0.36061)-0.784*(1.7321*0.36061^0.5/(1+2.0785*0.36061^0.5)+1.6667*ln(1+2.0785*0.36061^0.5))+2.6667*0.36061*(c+d*(1+(-1-3.464*0.36061^0.5+6*0.36061)*e0.36061p(-3.464*0.36061^0.5))/(12*0.36061))+2.8289*0.36061^2*e)+0.07214; 0.07391741=-(m*0.53672/(n*0.53672^0.5+1))*8.314*298.15*(ln(0.53672)-0.784*(1.7321*0.53672^0.5/(1+2.0785*0.53672^0.5)+1.6667*ln(1+2.0785*0.53672^0.5))+2.6667*0.53672*(c+d*(1+(-1-3.464*0.53672^0.5+6*0.53672)*e0.53672p(-3.464*0.53672^0.5))/(12*0.53672))+2.8289*0.53672^2*e)+0.07214; 0.07451667=-(m*0.72049/(n*0.72049^0.5+1))*8.314*298.15*(ln(0.72049)-0.784*(1.7321*0.72049^0.5/(1+2.0785*0.72049^0.5)+1.6667*ln(1+2.0785*0.72049^0.5))+2.6667*0.72049*(c+d*(1+(-1-3.464*0.72049^0.5+6*0.72049)*e0.72049p(-3.464*0.72049^0.5))/(12*0.72049))+2.8289*0.72049^2*e)+0.07214; 0.07521378=-(m*0.92046/(n*0.92046^0.5+1))*8.314*298.15*(ln(0.92046)-0.784*(1.7321*0.92046^0.5/(1+2.0785*0.92046^0.5)+1.6667*ln(1+2.0785*0.92046^0.5))+2.6667*0.92046*(c+d*(1+(-1-3.464*0.92046^0.5+6*0.92046)*e0.92046p(-3.464*0.92046^0.5))/(12*0.92046))+2.8289*0.92046^2*e)+0.07214; 0.07661342=(m*1.3003/(n*1.3003^0.5+1))*8.314*298.15*(ln(1.3003)-0.784*(1.7321*1.3003^0.5/(1+2.0785*1.3003^0.5)+1.6667*ln(1+2.0785*1.3003^0.5))+2.6667*1.3003*(c+d*(1+(-1-3.464*1.3003^0.5+6*1.3003)*e1.3003p(-3.464*1.3003^0.5))/(12*1.3003))+2.8289*1.3003^2*e)+0.07214; 0.07871152=(m*1.67471/(n*1.67471^0.5+1))*8.314*298.15*(ln(1.67471)-0.784*(1.7321*1.67471^0.5/(1+2.0785*1.67471^0.5)+1.6667*ln(1+2.0785*1.67471^0.5))+2.6667*1.67471*(c+d*(1+(-1-3.464*1.67471^0.5+6*1.67471)*e1.67471p(-3.464*1.67471^0.5))/(12*1.67471))+2.8289*1.67471^2*e)+0.07214; 0.07934748=(m*1.86291/(n*1.86291^0.5+1))*8.314*298.15*(ln(1.86291)-0.784*(1.7321*1.86291^0.5/(1+2.0785*1.86291^0.5)+1.6667*ln(1+2.0785*1.86291^0.5))+2.6667*1.86291*(c+d*(1+(-1-3.464*1.86291^0.5+6*1.86291)*e1.86291p(-3.464*1.86291^0.5))/(12*1.86291))+2.8289*1.86291^2*e)+0.07214; 0.08119555=(m*2.42891/(n*2.42891^0.5+1))*8.314*298.15*(ln(2.42891)-0.784*(1.7321*2.42891^0.5/(1+2.0785*2.42891^0.5)+1.6667*ln(1+2.0785*2.42891^0.5))+2.6667*2.42891*(c+d*(1+(-1-3.464*2.42891^0.5+6*2.42891)*e2.42891p(-3.464*2.42891^0.5))/(12*2.42891))+2.8289*2.42891^2*e)+0.07214; 0.08388612=(m*2.94245/(n*2.94245^0.5+1))*8.314*298.15*(ln(2.94245)-0.784*(1.7321*2.94245^0.5/(1+2.0785*2.94245^0.5)+1.6667*ln(1+2.0785*2.94245^0.5))+2.6667*2.94245*(c+d*(1+(-1-3.464*2.94245^0.5+6*2.94245)*e2.94245p(-3.464*2.94245^0.5))/(12*2.94245))+2.8289*2.94245^2*e)+0.07214; 0.08544611=(m*3.49978/(n*3.49978^0.5+1))*8.314*298.15*(ln(3.49978)-0.784*(1.7321*3.49978^0.5/(1+2.0785*3.49978^0.5)+1.6667*ln(1+2.0785*3.49978^0.5))+2.6667*3.49978*(c+d*(1+(-1-3.464*3.49978^0.5+6*3.49978)*e3.49978p(-3.464*3.49978^0.5))/(12*3.49978))+2.8289*3.49978^2*e)+0.07214; 0.08798993=(m*4.05405/(n*4.05405^0.5+1))*8.314*298.15*(ln(4.05405)-0.784*(1.7321*4.05405^0.5/(1+2.0785*4.05405^0.5)+1.6667*ln(1+2.0785*4.05405^0.5))+2.6667*4.05405*(c+d*(1+(-1-3.464*4.05405^0.5+6*4.05405)*e4.05405p(-3.464*4.05405^0.5))/(12*4.05405))+2.8289*4.05405^2*e)+0.07214; 以上放到function中,qiuyicicanshu 將data 分成兩段,針對(duì)不同的擬合函數(shù),分別做一次擬合。 將一個(gè)程序分成三個(gè)小程序處理,低版本1stop就可以勝任了。 |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
木蟲 (小有名氣)
木蟲 (小有名氣)
木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料292調(diào)劑 +8 | 橘頌思美人 2026-03-23 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +4 | 顧九笙要謙虛 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 7/350 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|