24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

小木虫zb

木虫 (正式写手)

应助: 11 (小学生)
金币: 1560.9
红花: 5
帖子: 345
在线: 142.1小时
虫号: 1227944
注册: 2011-03-10
性别: GG
专业: 分离过程

[求助] matlab解偏微分方程求助

本人用matlab中的线上法（ODE15s）解一偏微分方程模型，单组份已经成功解出，模型如单组份模型.jpg所示，相应matlab程序如附件textcuvre.m。
现在想用matlab解双组份模型，模型如图双组份模型.jpg，不知道怎么在单组份matlab程序中修改，主要是在双组份模型的平衡吸附模型的方程为隐函数方程，不知道怎么加入程序中，特向大神求助
单组份模型.jpg

双组份模型.jpg

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : textcurve.m

2013-04-25 15:32:52, 1.15 K

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有84人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】matlab求解偏微分方程已经有4人回复
求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？已经有5人回复
matlab求解偏微分方程组已经有3人回复
matlab有限差分法求解复杂的偏微分方程组已经有19人回复
Matlab求解偏微分方程组已经有22人回复
偏微分方程数值解有基本理论和一些例子（化工类）MATLAB源程序已经有296人回复
高金求助matlab解微分方程组已经有12人回复
【求助】用pdepe求解偏微分方程已经有10人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】如何求解下列偏微分方程已经有13人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】用MATLAB求解偏微分方程的数值解悬赏20 已经有8人回复
【求助】有限元法求解二阶偏微分方程组已经有6人回复
【求助】求助用matlab解偏微分方程已经有15人回复
【求助】用matlab编过有限差分法的请进已经有9人回复

1楼 2013-04-25 15:37:51

已阅关注TA 给TA发消息送TA红花 TA的回帖

小木虫zb

木虫 (正式写手)

应助: 11 (小学生)
金币: 1560.9
红花: 5
帖子: 345
在线: 142.1小时
虫号: 1227944
注册: 2011-03-10
性别: GG
专业: 分离过程

引用回帖:

8楼: Originally posted by SEEDEN at 2013-04-25 22:32:10
呵呵，但我还是觉得那两个较复杂的代数式并不太好处理
如果你的方法可行（在原有程序上的修改），可以把你的方法大概说下

...

我修改了下程序，但是还是有很多错误，
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);%颗粒表观密度
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))');