 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3290 | 回復(fù): 21 | |||
anxiaoxin2金蟲 (小有名氣)
|
[求助]
MATLAB應(yīng)用四階龍格庫塔法求解反應(yīng)動力學(xué)問題? 已有2人參與
|
||
|
各位英雄大俠，你們好。本人在做鐵水噴吹脫硫的優(yōu)化，遇到一個問題，文獻(xiàn)中是應(yīng)用MATLAB四階龍格庫塔法分析的，本人對MATLAB不熟，希望大家能幫幫忙！ 謝謝。。。 具體問題見附件。 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 脫硫動力學(xué)微分方程.docx
2015-03-05 14:00:32, 46.03 K
» 收錄本帖的淘帖專輯推薦
 動力學(xué)擬合 |
» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有210人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab 如何求解多元二階微分方程組
已經(jīng)有7人回復(fù)
- 【緊急求救】求解一個方程組的近似解或者估計解，有何辦法
已經(jīng)有12人回復(fù)
- 動力學(xué)方程組參數(shù)求解問題
已經(jīng)有9人回復(fù)
- 求助�。τ谡{(diào)用龍格庫塔法如何對結(jié)果判斷反饋
已經(jīng)有6人回復(fù)
- 求助用matlab編一個程序，采用龍格庫塔發(fā)求求解動力學(xué)方程的參數(shù)問題
已經(jīng)有5人回復(fù)
- 急！求大神用Matlab四階龍格庫塔解個方程！
已經(jīng)有5人回復(fù)
- matlab解方程組求助攻
已經(jīng)有4人回復(fù)
- 常微分方程解析解
已經(jīng)有5人回復(fù)
- 6元方程求解----繪制等高線
已經(jīng)有5人回復(fù)
- Matlab用龍格庫塔法求二階非線性方程
已經(jīng)有13人回復(fù)
- Matlab求解二階偏微分方程組，希望能給出相應(yīng)的求解代碼
已經(jīng)有10人回復(fù)
- 分?jǐn)?shù)階微分方程的求解（用差分法，matlab編程求解）
已經(jīng)有9人回復(fù)
- 求matlab 用四階龍格-庫塔法求解微分方程
已經(jīng)有4人回復(fù)
- 龍格庫塔法 matlab
已經(jīng)有6人回復(fù)
- 動力學(xué)參數(shù)擬合
已經(jīng)有26人回復(fù)
- 如何求解描述振動的二階微分方程
已經(jīng)有21人回復(fù)
- 一個微分方程組，求解
已經(jīng)有18人回復(fù)
- LABVIEW四階龍格庫塔法求解一階微分方程組
已經(jīng)有12人回復(fù)
- 【求助】催化反應(yīng)動力學(xué)
已經(jīng)有5人回復(fù)
- 【求助】Matlab中利用四階龍格-庫塔法求解微分方程！�。。�
已經(jīng)有9人回復(fù)
- 【求助】急請微分方程高手看下這個方程組！
已經(jīng)有14人回復(fù)
- 【求助】四階龍格-庫塔方法�。。。�
已經(jīng)有17人回復(fù)
- 【求助】使用Matlab預(yù)估動力學(xué)方程問題
已經(jīng)有13人回復(fù)
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
4樓2015-03-05 21:00:42
CelestialCYJ
木蟲 (小有名氣)
- 應(yīng)助: 33 (小學(xué)生)
- 金幣: 2482.6
- 紅花: 12
- 帖子: 144
- 在線: 43.3小時
- 蟲號: 3277506
- 注冊: 2014-06-16
- 專業(yè): 結(jié)構(gòu)工程
anxiaoxin22樓2015-03-05 16:11:26
anxiaoxin2
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2572.8
- 帖子: 178
- 在線: 86小時
- 蟲號: 1042931
- 注冊: 2010-06-17
- 性別: GG
- 專業(yè): 無機非金屬類光電信息與功
3樓2015-03-05 17:00:56
anxiaoxin2
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2572.8
- 帖子: 178
- 在線: 86小時
- 蟲號: 1042931
- 注冊: 2010-06-17
- 性別: GG
- 專業(yè): 無機非金屬類光電信息與功
5樓2015-03-06 08:37:00
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
6樓2015-03-06 08:40:39
anxiaoxin2
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2572.8
- 帖子: 178
- 在線: 86小時
- 蟲號: 1042931
- 注冊: 2010-06-17
- 性別: GG
- 專業(yè): 無機非金屬類光電信息與功
7樓2015-03-06 08:46:36
anxiaoxin2
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2572.8
- 帖子: 178
- 在線: 86小時
- 蟲號: 1042931
- 注冊: 2010-06-17
- 性別: GG
- 專業(yè): 無機非金屬類光電信息與功
|
你好，按你提供的輸入到matlab中了，顯示 dydt(1)=(6*b*Ip*U1*ip*Pm*Km*([S]*exp(-Tmix/T)-[S]CaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*([S]*exp(-Tmix/T)-[S]Mge))/Wm/Db+0.4167*k*Pm*[S]*exp(-Tmix/T)*[Mg]*exp(-Tmix/T)+A*Pz*Ks*Ls*([S]*exp(-Tmix/T)-(S)/Ls)/Wm; dydt(2)=100*B*U3*Ip/Wm-0.75*(0.4167*k*Pm*[S]*exp(-Tmix/T)*[Mg]*exp(-Tmix/T)); ??? function solve_odes | Error: Function definitions are not permitted in this context. 不知道是哪里出錯了？該怎么修改？另外輸入這些就可以作出[S]對t的圖形嗎？ |
8樓2015-03-17 17:06:06
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
9樓2015-03-17 18:05:45
anxiaoxin2
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2572.8
- 帖子: 178
- 在線: 86小時
- 蟲號: 1042931
- 注冊: 2010-06-17
- 性別: GG
- 專業(yè): 無機非金屬類光電信息與功
|
謝謝。。。 你說的問題我都已修改。。。 function solve_odes clear all;clc y0=[ ]; % [%S] [%Mg]在t=0時刻的數(shù)值，即其初值，待填入 tspan=1200; %時間t的范圍，待填入 [t,y]=ode45(@my_odes],[2 0]) figure(1),plot(t,y(:,1)); % [%s]對時間t作圖 figure(2),plot(t,y(:,2); % [%Mg]對時間t作圖 這段代碼中y0=[];此處有兩個初始值，怎么填入？ 后邊的微分方程中是不是所有的這兩個[%S] [%Mg]的初始值都要改成y0? |
10樓2015-03-18 09:10:35
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 anxiaoxin2 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +3 | vv迷 2026-03-22 | 3/150 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|