| 查看: 2738 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
小木虫zb木虫 (正式写手)
|
[求助]
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
» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有173人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有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人回复
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
|
我修改了下程序，但是还是有很多错误， 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))'); |
9楼2013-04-26 09:09:40
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-04-25 16:25:30
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
3楼2013-04-25 16:42:42
4楼2013-04-25 21:48:19
100