 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 575 | 回復(fù): 2 | ||||
zhang66163銀蟲 (正式寫手)
|
[求助]
matlab程序編寫 已有2人參與
|
|
按照教程編寫了一個(gè)程序：如下 function Batchreactor clear all;clc T=224.6+273 R=8.314 k0=[5.78052E+10, 3.92317E+12, 1.64254E+4, 6.264E+8] Ea=[124670, 150386, 77954, 111528] C0=[1, 0, 0, 0] tspan=[0, 1E+4] [t,C]=ode45(@Massequations,tspan,C0,[],k0,Ea,T,R); plot(t,C(:,1),'r-',t,C(:,2),'O-',t,C(:,3),'b-',t,C(:,4),'g-'); function dCdt=Massequations(t,C,k0,Ea,T,R); k=k0*exp(-Ea/(R*T)); k(5)=2.1667E-4; rA=-(k(1)+k(2))*C(1); rB=k(1)*C(1)-k(3)*C(2); rC=k(2)*C(1)-k(4)*C(3); rD=k(3)*C(2)-k(5)*C(4); rE=k(4)*C(3)+k(5)*C(4); 可是運(yùn)行后，出現(xiàn)了錯(cuò)誤；哪位大俠可以幫忙解釋一下。謝謝 Error using * Inner matrix dimensions must agree. Error in concentration>Massequations (line 16) k=k0*exp(-Ea/(R*T)); k(5)=2.1667E-4; Error in odearguments (line 88) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode45 (line 114) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in concentration (line 10) [t,C]=ode45(@Massequations,tspan,C0,[],k0,Ea,T,R); dCdt=[rA;rB;rC;rD;rE] |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有218人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 計(jì)算強(qiáng)帖: 2
- 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時(shí)
- 蟲號(hào): 449046
- 注冊(cè): 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計(jì)算模擬
2樓2015-08-11 09:08:23
月只藍(lán)
主管區(qū)長(zhǎng) (職業(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í)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
zhang66163: 金幣+40, ★★★很有幫助 2015-08-11 11:02:11
感謝參與，應(yīng)助指數(shù) +1
zhang66163: 金幣+40, ★★★很有幫助 2015-08-11 11:02:11
|
少一個(gè)初值條件，由C0=[1, 0, 0, 0]改為：C0=[1, 0, 0, 0,0]，即默認(rèn)組分E初始濃度為0；如果初值條件不是這樣的，請(qǐng)自行需要改，須填滿5個(gè)，因?yàn)橛?個(gè)方程。 附圖1.png |
3樓2015-08-11 09:27:15
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhang66163 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +7 | 自由煎餅果子 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 招08考數(shù)學(xué) +4 | laoshidan 2026-03-20 | 8/400 |
|
|
|
[考研] 306求調(diào)劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +4 | Sxy112 2026-03-21 | 6/300 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研]
 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|