 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4046 | 回復: 20 | |||||
kingfisher8288銅蟲 (小有名氣)
|
[交流]
【求助】matlab 求解微分方程中的未知參數(shù) 已有4人參與
|
||||
|
求助，本人剛學matlab，不太熟悉�，F(xiàn)在碰到一個問題。 一個反應體系會同時進行n個化學反應，彼此反應物和產(chǎn)物互相關聯(lián)，我列了一系列的反應速率方程，當所有的反應速率常數(shù)已知時，用ode15s解微分方程可以得到任何反應物濃度隨時間的變化。 現(xiàn)在其中有幾個速率常數(shù)是未知的，但是我知道一種反應物的濃度隨時間變化的曲線（反應時間及其對應的反應物濃度），用什么方法可以求出這幾個未知的速率常數(shù)？ ps：已知的反應曲線數(shù)量是大于未知數(shù)的，我覺得理論上應該可以求出那幾個速率常數(shù)。 各位大俠幫幫忙，感激不盡！ |
回復此樓» 收錄本帖的淘帖專輯推薦
 matlab | DFT模擬 |
» 猜你喜歡
289求調劑
已經(jīng)有4人回復
-
292求調劑
已經(jīng)有4人回復
-
080500求調劑
已經(jīng)有3人回復
-
291求調劑
已經(jīng)有3人回復
-
材料292調劑
已經(jīng)有8人回復
-
一志愿211 初試270分 求調劑
已經(jīng)有6人回復
-
求調劑
已經(jīng)有3人回復
-
274求調劑
已經(jīng)有4人回復
-
300求調劑，材料科學英一數(shù)二
已經(jīng)有5人回復
-
招08考數(shù)學
已經(jīng)有15人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- matlab解微分方程
已經(jīng)有10人回復
- Matlab求解偏微分方程組
已經(jīng)有22人回復
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復
- 高金求助matlab解微分方程組
已經(jīng)有12人回復
- 幫幫忙看看我的難題
已經(jīng)有15人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復
★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+1): 謝謝應助！可否請蟲友介紹一下forcal�。� 2011-04-11 21:05:29
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+1): 謝謝應助！可否請蟲友介紹一下forcal�。� 2011-04-11 21:05:29
|
本帖內容被屏蔽 |
3樓2011-04-11 16:47:03
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應助！ 2011-04-15 21:57:50
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應助！ 2011-04-15 21:57:50
|
本帖內容被屏蔽 |
9樓2011-04-15 08:24:52
★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzuwangshilei(金幣+1): 尋優(yōu)很不方便 2011-04-15 09:49:53
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzuwangshilei(金幣+1): 尋優(yōu)很不方便 2011-04-15 09:49:53
|
本帖內容被屏蔽 |
10樓2011-04-15 08:37:47
2樓2011-04-10 16:21:49
kingfisher8288
銅蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 197.2
- 散金: 47
- 帖子: 56
- 在線: 39.5小時
- 蟲號: 708266
- 注冊: 2009-02-25
- 性別: MM
- 專業(yè): 污染控制化學
4樓2011-04-11 22:31:34
★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
|
本帖內容被屏蔽 |
5樓2011-04-14 03:11:13
kingfisher8288
銅蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 197.2
- 散金: 47
- 帖子: 56
- 在線: 39.5小時
- 蟲號: 708266
- 注冊: 2009-02-25
- 性別: MM
- 專業(yè): 污染控制化學
★
xiegangmai(金幣+1): 鼓勵自己動手，積極討論交流！ 2011-04-14 23:24:32
xiegangmai(金幣+1): 鼓勵自己動手，積極討論交流！ 2011-04-14 23:24:32
|
function ydot = f(t, y) k1 = 3.5*10^4; k2 = 1.0*10^3; k3 = k=a*ka+b*kb+c*kc; r1 = k1*y(1)*y(2); r2 = k2*y(1)*y(3); r3 = k3*y(1)*y(4); z1 = -r1-r2-r3; z2 = -r1; z3 = -r2+r1; z4 = -r3+r2; ydot = [z1; z2; z3; z4;]; y0 = [30.0*10^-6;10.2*10^-6;4.1*10^-6;6.7] 當[a b c]=[100 5 20]時 t=[0, 10, 20, 30, 40] y(1)=[30.0*10^-6, 17.5*10^-6, 15*10^-6, 14*10^-6, 13.5*10^-6 ] 當[a, b, c]=[120, 15, 25]時 t=[0 10 20 30 40] y(1)=[30.0*10^-6, 12.5*10^-6, 10*10^-6, 8*10^-6, 7*10^-6 ] 求ka，kb，kc 我自己的模型東西太多了，看著費勁。上面的例子和數(shù)據(jù)是我自己編的，但基本能描述我所要的東西，請幫忙解答，謝謝。 [ Last edited by kingfisher8288 on 2011-4-14 at 04:26 ] |
6樓2011-04-14 04:22:10
★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+1): 謝謝建議！對照著數(shù)學模型，方便理解。 2011-04-14 23:25:24
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+1): 謝謝建議！對照著數(shù)學模型，方便理解。 2011-04-14 23:25:24
|
本帖內容被屏蔽 |
7樓2011-04-14 19:46:12
kingfisher8288
銅蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 197.2
- 散金: 47
- 帖子: 56
- 在線: 39.5小時
- 蟲號: 708266
- 注冊: 2009-02-25
- 性別: MM
- 專業(yè): 污染控制化學
|
謝謝forcal關注，我不知道什么是數(shù)學語言，我就大概解釋一下吧，k1，k2，k3等是反應速率常數(shù)；r1，r2，r3 等是基元反應速率；y（1），y（2），y（3）等代表反應物的濃度；基元反應速率等于反應速率常數(shù)乘以相應的反應物濃度；z1，z2，z3 分別代表三種反應物的反應速率，它是由各個基元反應加和而成的，因為對于每一種物質來說有的基元反應是消耗該物質，有的是生成該物質。對于【y（1），y（2），y（3），y(4)】 有一組初始值 [30.0*10^-6;10.2*10^-6;4.1*10^-6;6.7]。k3是未知的，但是它和a，b，c 三個物性常數(shù)有關，k3 = k=a*ka+b*kb+c*kc; 當給定一組物性常數(shù)【a，b，c 】=[100， 5， 20]時，獲得一組y（1）（一種反應物的濃度）對時間的數(shù)據(jù)t=[0, 10, 20, 30, 40] y(1)=[30.0*10^-6, 17.5*10^-6, 15*10^-6, 14*10^-6, 13.5*10^-6 ]； 當改變物性常數(shù)的值為[a, b, c]=[120, 15, 25]時， t=[0 10 20 30 40] y(1)=[30.0*10^-6, 12.5*10^-6, 10*10^-6, 8*10^-6, 7*10^-6 ]；當然還可以通過改變物性常數(shù)獲得n組數(shù)據(jù)。我需要確定的是反應速率常數(shù)對于物性常數(shù)的相關系數(shù) ka，kb，kc。ka，kb，kc 有可能為正值，也可能為負值，也就是說a，b，c可能正作用于k3，也可能負作用于k3. 不知道我的描述是不是能有助于你的理解，謝謝！ |
8樓2011-04-14 23:26:48
相關版塊跳轉
我要訂閱樓主 kingfisher8288 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277分求調劑，跨調材料 +3 | 考研調劑lxh 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 生物學學碩求調劑 +6 | 小羊睡著了? 2026-03-23 | 6/300 |
|
|
|
[考研] 收08調劑生 +3 | komorebi69 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調劑 +3 | gdlf9999 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿中國石油大學（華東） 本科齊魯工業(yè)大學 +4 | 石能偉 2026-03-17 | 4/200 |
|
|
[考研] 280分求調劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 石河子大學（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 298求調劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿北京化工大學070300 學碩336求調劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 354求調劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考！ +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
|
[考研]
 一志愿東華大學控制學碩320求調劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 0805 316求調劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 268求調劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 085700資源與環(huán)境308求調劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 一志愿華南師大 070300（化學）304分求調劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 中南大學化學學碩337求調劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|