版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

anxiaoxin2

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2572.8
帖子: 178
在線: 86小時(shí)
蟲號: 1042931
注冊: 2010-06-17
性別: GG
專業(yè): 無機(jī)非金屬類光電信息與功

[求助] MATLAB應(yīng)用四階龍格庫塔法求解反應(yīng)動(dòng)力學(xué)問題? 已有2人參與

各位英雄大俠，你們好。本人在做鐵水噴吹脫硫的優(yōu)化，遇到一個(gè)問題，文獻(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 : 脫硫動(dòng)力學(xué)微分方程.docx

2015-03-05 14:00:32, 46.03 K

1樓 2015-03-05 14:01:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
anxiaoxin2: 金幣+300, ★★★很有幫助, 非常有用。。。 2015-03-06 08:37:41

CODE:

function solve_odes

clear all;clc

y0=[ ];    % [%s] [%Mg]在t=0時(shí)刻的數(shù)值，即其初值，待填入

tspan=; %時(shí)間t的范圍，待填入

[t,y]=ode45(@my_odes],[2 0]) 

figure(1),plot(t,y(:,1));     % [%s]對時(shí)間t作圖

figure(2),plot(t,y(:,2);      % [%Mg]對時(shí)間t作圖

function dydt=my_odes

beta= ;    %β的數(shù)值

Ip= ;    Ip的數(shù)值

......   %以下依次把已經(jīng)參數(shù)的數(shù)值填入。

dydt(1)=        ;  % d[%s]/dt的表達(dá)式

dydt(2)=        ;  % d[%Mg]/dt的表達(dá)式

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

4樓2015-03-05 21:00:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

CelestialCYJ

木蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 2482.6
紅花: 12
帖子: 144
在線: 43.3小時(shí)
蟲號: 3277506
注冊: 2014-06-16
專業(yè): 結(jié)構(gòu)工程

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

你不給數(shù)值怎么算

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

anxiaoxin2

2樓2015-03-05 16:11:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anxiaoxin2

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2572.8
帖子: 178
在線: 86小時(shí)
蟲號: 1042931
注冊: 2010-06-17
性別: GG
專業(yè): 無機(jī)非金屬類光電信息與功

送紅花一朵

引用回帖:

2樓: Originally posted by CelestialCYJ at 2015-03-05 16:11:26
你不給數(shù)值怎么算

這些都是參數(shù)，可以設(shè)定，方便以后調(diào)整，你有編寫的程序嗎？

贊一下

回復(fù)此樓

3樓2015-03-05 17:00:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anxiaoxin2

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2572.8
帖子: 178
在線: 86小時(shí)
蟲號: 1042931
注冊: 2010-06-17
性別: GG
專業(yè): 無機(jī)非金屬類光電信息與功

謝謝版主。。。這個(gè)不需要m文件之類的嗎？

贊一下

回復(fù)此樓

5樓2015-03-06 08:37:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

引用回帖:

5樓: Originally posted by anxiaoxin2 at 2015-03-06 08:37:00
謝謝版主。。。這個(gè)不需要m文件之類的嗎？

你好，因?yàn)槟銢]有各個(gè)變量或參數(shù)的具體數(shù)據(jù)，我給出的是一個(gè)m文件的框架，你把你的參數(shù)和方程填入，就能計(jì)算的。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

6樓2015-03-06 08:40:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anxiaoxin2

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2572.8
帖子: 178
在線: 86小時(shí)
蟲號: 1042931
注冊: 2010-06-17
性別: GG
專業(yè): 無機(jī)非金屬類光電信息與功

好的，我先試試，不懂的再來請教。。。

贊一下

回復(fù)此樓

7樓2015-03-06 08:46:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anxiaoxin2

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2572.8
帖子: 178
在線: 86小時(shí)
蟲號: 1042931
注冊: 2010-06-17
性別: GG
專業(yè): 無機(jī)非金屬類光電信息與功

引用回帖:

6樓: Originally posted by 月只藍(lán) at 2015-03-06 08:40:39
你好，因?yàn)槟銢]有各個(gè)變量或參數(shù)的具體數(shù)據(jù)，我給出的是一個(gè)m文件的框架，你把你的參數(shù)和方程填入，就能計(jì)算的。...

你好，按你提供的輸入到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.
不知道是哪里出錯(cuò)了？該怎么修改？另外輸入這些就可以作出[S]對t的圖形嗎？

贊一下

回復(fù)此樓

8樓2015-03-17 17:06:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

引用回帖:

8樓: Originally posted by anxiaoxin2 at 2015-03-17 17:06:06
你好，按你提供的輸入到matlab中了，顯示
dydt(1)=(6*b*Ip*U1*ip*Pm*Km*(*exp(-Tmix/T)-CaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*(*exp(-Tmix/T)-Mge))/Wm/Db+0.4167*k*Pm**exp(-Tmix/T)**exp(-Tmix/T)+A*Pz*Ks*Ls*(*ex ...

在MATLAB中，[S]這種代碼寫法不對的，用單個(gè)字母代替，如果字母有重復(fù)，可以用字母+數(shù)字的屬性，比如S1 S2 S_1 S_2

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

9樓2015-03-17 18:05:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anxiaoxin2

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2572.8
帖子: 178
在線: 86小時(shí)
蟲號: 1042931
注冊: 2010-06-17
性別: GG
專業(yè): 無機(jī)非金屬類光電信息與功

引用回帖:

9樓: Originally posted by 月只藍(lán) at 2015-03-17 18:05:45
在MATLAB中，這種代碼寫法不對的，用單個(gè)字母代替，如果字母有重復(fù)，可以用字母+數(shù)字的屬性，比如S1 S2 S_1 S_2...

謝謝。。。
你說的問題我都已修改。。。
function solve_odes
clear all;clc
y0=[ ]; % [%S] [%Mg]在t=0時(shí)刻的數(shù)值，即其初值，待填入
tspan=1200; %時(shí)間t的范圍，待填入
[t,y]=ode45(@my_odes],[2 0])
figure(1),plot(t,y(:,1)); % [%s]對時(shí)間t作圖
figure(2),plot(t,y(:,2); % [%Mg]對時(shí)間t作圖
這段代碼中y0=[];此處有兩個(gè)初始值，怎么填入？
后邊的微分方程中是不是所有的這兩個(gè)[%S] [%Mg]的初始值都要改成y0?

贊一下

回復(fù)此樓

10樓2015-03-18 09:10:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 anxiaoxin2 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-04-02	7/350	2026-04-02 16:20 by chran16
[考研] 330分求調(diào)劑 +12	qzenlc 2026-03-29	12/600	2026-04-02 15:38 by 不吃魚的貓
[考研] 材料調(diào)劑 +6	懶羊羊輕置玉臀 2026-04-02	6/300	2026-04-02 14:47 by 5896
[考研] 材料求調(diào)劑一志愿哈工大324 +10	閆旭東 2026-03-28	12/600	2026-04-02 14:37 by olim
[考研] 能源動(dòng)力調(diào)劑 +3	不破不立0 2026-04-02	3/150	2026-04-02 12:46 by ffffjjjj
[考研] 材料調(diào)劑 +8	一樣YWY 2026-04-02	8/400	2026-04-02 12:30 by 173785543
[考研] 329求調(diào)劑 +8	miaodesi 2026-04-02	10/500	2026-04-02 10:24 by 不吃魚的貓
[考研] 337求調(diào)劑 +11	《樹》 2026-03-29	11/550	2026-04-02 10:20 by 不吃魚的貓
[考研] 求生物學(xué)調(diào)劑 +8	15172915737 2026-04-01	8/400	2026-04-02 06:49 by ilovexiaobin
[考研] 310分求調(diào)劑 +4	成功上岸wang 2026-04-01	4/200	2026-04-01 20:35 by liu823948201
[考研] 英一數(shù)二生物信息學(xué)287分，本科生物科學(xué)，求調(diào)劑 +3	碧水xyz 2026-03-29	4/200	2026-04-01 15:35 by plum
[考研] 08工科275求調(diào)劑，可跨考。 +5	AaAa7420 2026-03-31	5/250	2026-04-01 15:21 by 159357hjz
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑六級521 +4	愛吃大鴨梨 2026-03-31	4/200	2026-03-31 11:51 by 搏擊518
[考研] 08工科求調(diào)劑286 +5	tgs_001 2026-03-28	5/250	2026-03-31 08:18 by 一只好果子?
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan

24小時(shí)熱門版塊排行榜

[求助] MATLAB應(yīng)用四階龍格庫塔法求解反應(yīng)動(dòng)力學(xué)問題? 已有2人參與

» 本帖附件資源列表

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）