| 查看: 2453 | 回复: 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
» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有299人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】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
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 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
