 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1622 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
姜玉旺木蟲 (著名寫手)
|
[求助]
關(guān)于matlab微分、及非線性擬合問(wèn)題
|
|
|
最近一直苦苦思索化學(xué)反應(yīng)速率的問(wèn)題，想用matlab來(lái)求微分?jǐn)M合非線性曲線求出常數(shù)，可我確實(shí)不懂matlab，望各位大俠不吝賜教，幫我編寫出程序，最好能有圖形結(jié)合，小弟感激涕零： 我的方程式是（dx）/（dt）=a*c*[(K-1)*x^2-32*K*x+31*K] t,x的取值為t=0,x=0;t=1,x=0.0355;t=2,x=0.1014;t=3,x=0.1394;t=4,x=0.2147;t=5,x=0.2823;t=6,x=0.3274;t=7,x=0.3717;t=8,x=0.4095. 其中a,K未知，c=1.213. 小弟已經(jīng)計(jì)算好幾天了，其在無(wú)能為力，拜托了，幫我編出程序。。。感謝了。。。 |
回復(fù)此樓» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有17人回復(fù)
-
考研調(diào)劑
已經(jīng)有7人回復(fù)
-
343求調(diào)劑085601
已經(jīng)有3人回復(fù)
-
08工科求調(diào)劑286
已經(jīng)有3人回復(fù)
-
337求調(diào)劑
已經(jīng)有3人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有10人回復(fù)
-
291求調(diào)劑
已經(jīng)有6人回復(fù)
-
299求調(diào)劑
已經(jīng)有10人回復(fù)
-
面上5B能上會(huì)嗎？
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab數(shù)據(jù)擬合問(wèn)題求助
已經(jīng)有7人回復(fù)
- matlab中關(guān)于非線性擬合，用nlinfit.程序總出錯(cuò)。謝謝各位大俠
已經(jīng)有3人回復(fù)
- 用MATLAB或者ORIGIN怎么擬合復(fù)數(shù)數(shù)據(jù)
已經(jīng)有10人回復(fù)
- 請(qǐng)問(wèn)matlab擬合二元方程的問(wèn)題
已經(jīng)有3人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較
已經(jīng)有4人回復(fù)
- matlab解微分方程
已經(jīng)有10人回復(fù)
- 在matlab中擬合曲面的表達(dá)式
已經(jīng)有29人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- matlab解微分方程組
已經(jīng)有15人回復(fù)
- matlab非線性參數(shù)擬合問(wèn)題
已經(jīng)有7人回復(fù)
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
姜玉旺
木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2747.2
- 紅花: 3
- 帖子: 2669
- 在線: 109.5小時(shí)
- 蟲號(hào): 1892949
- 注冊(cè): 2012-07-14
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
3樓2013-01-14 14:08:23
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與 2013-01-14 00:21:28
姜玉旺: 金幣+5, ★★★很有幫助 2013-01-14 08:52:30
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與 2013-01-14 00:21:28
姜玉旺: 金幣+5, ★★★很有幫助 2013-01-14 08:52:30
|
分離微分方程變量以后，如果懶得用配方法解出微分方程的話，弄個(gè)mathematica可以直接用解析法解 Integrate[1/( a*c*((k - 1)*x^2 + 32*k*x + 31*k)), x] 16 k - x + k x ArcTanh[------------------------] Sqrt[k] Sqrt[31 + 225 k] -(---------------------------------) - 66 = t a c Sqrt[k] Sqrt[31 + 225 k] 然后顯然 f(x_n) - t_n =0 喜歡用fmincon做優(yōu)化 還是用nlinfit做非線性回歸求系數(shù)也是隨便你啦。 |
2樓2013-01-13 21:06:18
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
姜玉旺: 金幣+5, ★★★很有幫助 2013-01-15 10:48:45
姜玉旺: 金幣+5, ★★★很有幫助 2013-01-15 10:48:45
|
積分我是用MATHEMATICA做的，不是matlab。 將 16 k - x + k x ArcTanh[------------------------] Sqrt[k] Sqrt[31 + 225 k] - (---------------------------------) - 66 - t = 0 a c Sqrt[k] Sqrt[31 + 225 k] 變換成matlab的代碼： 然后在matlab可以定義newfunction.m function result = newfunction(X, beta) k = beta(1); a = beta(2); x = X(1); t = X(2); result = 前面那個(gè)公式 最后用nlinfit(X, Y, beta0)算結(jié)果 其中 X = [ 0, 0; 1, 0.0355; 等等] Y = [0 ; 0 等等]（和X等高） beta0 = [ 0, 0] 是初值，你猜一個(gè)吧�？碿ode就可以知道前面是k后面是a了。 |
4樓2013-01-15 02:54:24
姜玉旺
木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2747.2
- 紅花: 3
- 帖子: 2669
- 在線: 109.5小時(shí)
- 蟲號(hào): 1892949
- 注冊(cè): 2012-07-14
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
5樓2013-01-15 10:49:21
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +5 | 彼岸t 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +11 | 溫柔的晚安 2026-03-25 | 12/600 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +12 | 阿斯蒂芬2004 2026-03-28 | 12/600 |
|
|
|
[考研] 求調(diào)劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿國(guó)科過(guò)程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|