版塊導(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í)熱門版塊排行榜

返回列表

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
金幣: 68127.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
金幣: 68127.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
金幣: 68127.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ā)表

[考研] 招08考數(shù)學(xué) +3	laoshidan 2026-03-20	7/350	2026-03-22 12:54 by 在風(fēng)落中
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	16/800	2026-03-22 11:20 by 脫穎而出
[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 資源與環(huán)境調(diào)劑申請(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 求助 +5	夢里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 298求調(diào)劑 +4	上岸6666@ 2026-03-20	4/200	2026-03-21 17:14 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 268求調(diào)劑 +9	簡單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 310求調(diào)劑 +3	baibai1314 2026-03-16	3/150	2026-03-21 03:56 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[論文投稿] 申請回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl

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

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

» 本帖附件資源列表

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

» 猜你喜歡

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

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

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

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

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