 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1621 | 回復: 5 | ||
姜玉旺木蟲 (著名寫手)
|
[求助]
關于matlab微分、及非線性擬合問題
|
|
最近一直苦苦思索化學反應速率的問題，想用matlab來求微分擬合非線性曲線求出常數，可我確實不懂matlab，望各位大俠不吝賜教，幫我編寫出程序，最好能有圖形結合，小弟感激涕零： 我的方程式是（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. 小弟已經計算好幾天了，其在無能為力，拜托了，幫我編出程序。。。感謝了。。。 |
回復此樓» 猜你喜歡
材料科學與工程求調劑
已經有3人回復
-
289求調劑
已經有5人回復
-
329求調劑，一志愿西北工業(yè)大學，材料工程（085601）
已經有6人回復
-
化學工程085602 305分求調劑
已經有21人回復
-
26申博自薦
已經有5人回復
-
296求調劑
已經有5人回復
-
總分293求調劑
已經有10人回復
-
材料專碩調劑
已經有3人回復
-
調劑310
已經有12人回復
-
化學0703 調劑 306分 一志愿211
已經有6人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- matlab數據擬合問題求助
已經有7人回復
- matlab中關于非線性擬合，用nlinfit.程序總出錯。謝謝各位大俠
已經有3人回復
- 用MATLAB或者ORIGIN怎么擬合復數數據
已經有10人回復
- 請問matlab擬合二元方程的問題
已經有3人回復
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應曲線？
已經有13人回復
- matlab微分方程組參數擬合，以周為單位求解，匯總后以年為單位進行數值比較
已經有4人回復
- matlab解微分方程
已經有10人回復
- 在matlab中擬合曲面的表達式
已經有29人回復
- 高金求助matlab解微分方程組
已經有12人回復
- matlab解微分方程組
已經有15人回復
- matlab非線性參數擬合問題
已經有7人回復
- 【求助】matlab 求解微分方程中的未知參數
已經有20人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經有21人回復
- 【求助】積分微分方程matlab求解
已經有6人回復
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數 +1
xiegangmai: 金幣+2, 謝謝參與 2013-01-14 00:21:28
姜玉旺: 金幣+5, ★★★很有幫助 2013-01-14 08:52:30
感謝參與，應助指數 +1
xiegangmai: 金幣+2, 謝謝參與 2013-01-14 00:21:28
姜玉旺: 金幣+5, ★★★很有幫助 2013-01-14 08:52:30
|
分離微分方程變量以后，如果懶得用配方法解出微分方程的話，弄個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做非線性回歸求系數也是隨便你啦。 |
2樓2013-01-13 21:06:18
姜玉旺
木蟲 (著名寫手)
- 應助: 0 (幼兒園)
- 金幣: 2747.2
- 紅花: 3
- 帖子: 2669
- 在線: 109.5小時
- 蟲號: 1892949
- 注冊: 2012-07-14
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
3樓2013-01-14 14:08:23
【答案】應助回帖
★ ★ ★ ★ ★
姜玉旺: 金幣+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 = 前面那個公式 最后用nlinfit(X, Y, beta0)算結果 其中 X = [ 0, 0; 1, 0.0355; 等等] Y = [0 ; 0 等等]（和X等高） beta0 = [ 0, 0] 是初值，你猜一個吧。看code就可以知道前面是k后面是a了。 |
4樓2013-01-15 02:54:24
姜玉旺
木蟲 (著名寫手)
- 應助: 0 (幼兒園)
- 金幣: 2747.2
- 紅花: 3
- 帖子: 2669
- 在線: 109.5小時
- 蟲號: 1892949
- 注冊: 2012-07-14
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
5樓2013-01-15 10:49:21
dingd
鐵桿木蟲 (職業(yè)作家)
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
姜玉旺: 金幣+10, ★★★很有幫助 2013-11-27 18:37:37
姜玉旺: 金幣+10, ★★★很有幫助 2013-11-27 18:37:37
|
用1stOpt直接微分方程擬合： 均方差(RMSE):0.0144298065185047 殘差平方和(SSE):0.00166575452929184 相關系數(R): 0.995283478470763 相關系數之平方(R^2): 0.990589202516862 確定系數(DC): 0.986895420543766 F統(tǒng)計(F-Statistic): 457.855211610462 參數 最佳估算 -------------------- ------------- a -0.101555842272903 k -0.0151477051447373 ====== 結果輸出 ====== 文件: 數據文件-1 No t 目標 x 計算 x 1 1 0.0355 0.0562847266411599 2 2 0.1014 0.110052998758528 3 3 0.1394 0.162114481657899 4 4 0.2147 0.213185538231459 5 5 0.2823 0.263930448271775 6 6 0.3274 0.314996671293445 7 7 0.3717 0.367048236670994 8 8 0.4095 0.420800978410109 |
6樓2013-11-22 12:46:12
相關版塊跳轉
我要訂閱樓主 姜玉旺 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調劑 +9 | 芋泥小鈴鐺 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 0856材料化工調劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|
|
[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 283求調劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 291求調劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 308求調劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 求調劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 315分求調劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 295求調劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 070300化學求調劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 085601 材料工程 313分 求調劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 314求調劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 材料求調劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研]
 304材料求調劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研]
 085602化學工程求調劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 271求調劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
|
[考研] 【2026考研調劑】制藥工程 284分 求相關專業(yè)調劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 282求調劑 +3 | wcq131415 2026-03-24 | 3/150 |
|