| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 1651 | 回復(fù): 16 | |||||
hanmei98232木蟲(chóng) (小有名氣)
|
[求助]
求助:以下非線性擬合方程可以用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)該怎么解? 求助求助。。!非常感謝 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è)樣編寫(xiě)對(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)附件。 |

木蟲(chóng) (小有名氣)
主管區(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 |
木蟲(chóng) (小有名氣)
木蟲(chóng) (小有名氣)
木蟲(chóng) (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
pk3725069 2026-03-19 | 17/850 |
|
|---|---|---|---|---|
|
[考博] 申博26年 +4 | 八6八68 2026-03-19 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5 | niko- 2026-03-22 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +3 | SYA! 2026-03-23 | 3/150 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚(yú) 2026-03-18 | 7/350 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|