 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2422 | 回復(fù): 12 | |||
小木蟲(chóng)zb木蟲(chóng) (正式寫(xiě)手)
|
[求助]
matlab解偏微分方程求助
|
||
|
本人用matlab中的線(xiàn)上法（ODE15s）解一偏微分方程模型，單組份已經(jīng)成功解出，模型如單組份模型.jpg所示，相應(yīng)matlab程序如附件textcuvre.m。 現(xiàn)在想用matlab解雙組份模型，模型如圖雙組份模型.jpg，不知道怎么在單組份matlab程序中修改，主要是在雙組份模型的平衡吸附模型的方程為隱函數(shù)方程，不知道怎么加入程序中，特向大神求助 單組份模型.jpg 雙組份模型.jpg |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : textcurve.m
2013-04-25 15:32:52, 1.15 K
» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有128人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專(zhuān)博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】matlab求解偏微分方程
已經(jīng)有4人回復(fù)
- 求助：如何用MATLAB的偏微分方程工具箱畫(huà)出不同區(qū)域內(nèi)的溫度分布圖？
已經(jīng)有5人回復(fù)
- matlab求解偏微分方程組
已經(jīng)有3人回復(fù)
- matlab有限差分法求解復(fù)雜的偏微分方程組
已經(jīng)有19人回復(fù)
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- 偏微分方程數(shù)值解有基本理論和一些例子（化工類(lèi)）MATLAB源程序
已經(jīng)有296人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- 【求助】用pdepe求解偏微分方程
已經(jīng)有10人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】如何求解下列偏微分方程
已經(jīng)有13人回復(fù)
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復(fù)
- 【求助】用MATLAB求解偏微分方程的數(shù)值解 懸賞20
已經(jīng)有8人回復(fù)
- 【求助】有限元法求解二階偏微分方程組
已經(jīng)有6人回復(fù)
- 【求助】求助用matlab解偏微分方程
已經(jīng)有15人回復(fù)
- 【求助】用matlab編過(guò)有限差分法的請(qǐng)進(jìn)
已經(jīng)有9人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專(zhuān)家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線(xiàn): 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專(zhuān)業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
2樓2013-04-25 16:25:30
小木蟲(chóng)zb
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1560.9
- 紅花: 5
- 帖子: 345
- 在線(xiàn): 142.1小時(shí)
- 蟲(chóng)號(hào): 1227944
- 注冊(cè): 2011-03-10
- 性別: GG
- 專(zhuān)業(yè): 分離過(guò)程
3樓2013-04-25 16:42:42
4樓2013-04-25 21:48:19
5樓2013-04-25 21:50:20
小木蟲(chóng)zb
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1560.9
- 紅花: 5
- 帖子: 345
- 在線(xiàn): 142.1小時(shí)
- 蟲(chóng)號(hào): 1227944
- 注冊(cè): 2011-03-10
- 性別: GG
- 專(zhuān)業(yè): 分離過(guò)程
6樓2013-04-25 22:24:56
小木蟲(chóng)zb
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1560.9
- 紅花: 5
- 帖子: 345
- 在線(xiàn): 142.1小時(shí)
- 蟲(chóng)號(hào): 1227944
- 注冊(cè): 2011-03-10
- 性別: GG
- 專(zhuān)業(yè): 分離過(guò)程
7樓2013-04-25 22:27:11
8樓2013-04-25 22:32:10
小木蟲(chóng)zb
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1560.9
- 紅花: 5
- 帖子: 345
- 在線(xiàn): 142.1小時(shí)
- 蟲(chóng)號(hào): 1227944
- 注冊(cè): 2011-03-10
- 性別: GG
- 專(zhuān)業(yè): 分離過(guò)程
|
我修改了下程序，但是還是有很多錯(cuò)誤， function dydt=fangcheng(t,y) kf=200; a=0.363;%空隙率 u=0.361032/a;%床層顆粒間隙中流體流速 c0=7.03; L=8.31e-2; N=20; dz=L/N; p=436.8/(1-0.363);%顆粒表觀(guān)密度 q=y(1:N);c=y(N+1:2*N); [ce1,ce2]=equilibrium(q1,q2); %-----------定義dq/dt----------------------------- for j=1:N dq1dt(j)=kf/p*(c1(j)-ce1(j)); dq2dt(j)=kf/p*(c2(j)-ce2(j)); end %---------定義dC/dt----------------------------------- dc1dt(1)=-u*(c1(1)-c0)/dz-(1-a)/a*p*dq1dt(1); dc2dt(1)=-u*(c2(1)-c0)/dz-(1-a)/a*p*dq2dt(1); for j=2:N dc1dt(j)=-u*(c1(j)-c1(j-1))/dz-(1-a)/a*p*dq1dt(j); dc2dt(j)=-u*(c2(j)-c2(j-1))/dz-(1-a)/a*p*dq2dt(j); end dydt=[dqdt dcdt]'; function [ce1,ce2]=equilibrium(q1,q2) [ce1,ce2]=fzero('q1-0.45967ce1^0.3839/(ce1^0.0457+0.3*ce2^0.0003206',... 'q2-0.17356*ce2^(-1.55201)/(ce2^(-2)+0.1*ce1^(-1))'); |
9樓2013-04-26 09:09:40
★
csgt0: 金幣+1, 謝謝 2013-10-15 13:32:09
csgt0: 金幣+1, 謝謝 2013-10-15 13:32:09
|
fzero可以求解無(wú)解析式解的代數(shù)方程的符號(hào)解嗎？ 我個(gè)人覺(jué)得你單獨(dú)使用這個(gè)命令或此類(lèi)命令恐怕都不能實(shí)現(xiàn)你的目的，當(dāng)然，你可以再多試試幾個(gè)命令。 這個(gè)代數(shù)方程組就算有解，你確定只有一組解？而且具體是哪一組（如果多組解的話(huà)），那也跟初值密切相關(guān)，要在賦初值的基礎(chǔ)上求解的。 總的來(lái)講，我還是覺(jué)得把那兩個(gè)關(guān)于Q的微分方程轉(zhuǎn)成C*的可能會(huì)好些 求導(dǎo)的過(guò)程可以利用matlab的符號(hào)求解功能，雖然形式復(fù)雜點(diǎn)，但應(yīng)該不影響一般的求解 你可以再試試你的方法，如果還是不理想，可以試試轉(zhuǎn)成C*的方程試試看 |
» 本帖已獲得的紅花（最新10朵）
小木蟲(chóng)zb10樓2013-04-26 10:35:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 小木蟲(chóng)zb 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛(ài)吃番茄的旭 2026-03-24 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +4 | 星野? 2026-03-26 | 4/200 |
|
|
|
[考研] 材料求調(diào)劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話(huà) 2026-03-24 | 7/350 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|