 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1624 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
姜玉旺木蟲(chóng) (著名寫(xiě)手)
|
[求助]
關(guān)于matlab微分、及非線性擬合問(wèn)題
|
||
|
最近一直苦苦思索化學(xué)反應(yīng)速率的問(wèn)題，想用matlab來(lái)求微分?jǐn)M合非線性曲線求出常數(shù)，可我確實(shí)不懂matlab，望各位大俠不吝賜教，幫我編寫(xiě)出程序，最好能有圖形結(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ù)此樓» 猜你喜歡
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有11人回復(fù)
-
332求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料專碩 085600求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有14人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有14人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有11人回復(fù)
-
環(huán)境科學(xué)與工程334分求調(diào)劑
已經(jīng)有8人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(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)助回帖
★ ★ ★ ★ ★
姜玉旺: 金幣+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è)吧。看code就可以知道前面是k后面是a了。 |
4樓2013-01-15 02:54:24
【答案】應(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
姜玉旺
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2747.2
- 紅花: 3
- 帖子: 2669
- 在線: 109.5小時(shí)
- 蟲(chóng)號(hào): 1892949
- 注冊(cè): 2012-07-14
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
3樓2013-01-14 14:08:23
姜玉旺
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2747.2
- 紅花: 3
- 帖子: 2669
- 在線: 109.5小時(shí)
- 蟲(chóng)號(hào): 1892949
- 注冊(cè): 2012-07-14
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
5樓2013-01-15 10:49:21
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程專碩求調(diào)劑 +6 | hyl3153942 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085601材料工程找調(diào)劑 +16 | oatmealR 2026-03-29 | 17/850 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +10 | 麻辣魷魚(yú) 2026-03-27 | 11/550 |
|
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
|
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +8 | 舉個(gè)栗子oi 2026-03-24 | 9/450 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我！��！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書(shū) 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|