| 查看: 2460 | 回复: 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论文润色/翻译怎么收费?
已经有137人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有12人回复
-
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人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +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
5楼2013-04-25 21:50:20
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
6楼2013-04-25 22:24:56
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
7楼2013-04-25 22:27:11
8楼2013-04-25 22:32:10
小木虫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
★
csgt0: 金币+1, 谢谢 2013-10-15 13:32:09
csgt0: 金币+1, 谢谢 2013-10-15 13:32:09
|
fzero可以求解无解析式解的代数方程的符号解吗？ 我个人觉得你单独使用这个命令或此类命令恐怕都不能实现你的目的，当然，你可以再多试试几个命令。 这个代数方程组就算有解，你确定只有一组解？而且具体是哪一组（如果多组解的话），那也跟初值密切相关，要在赋初值的基础上求解的。 总的来讲，我还是觉得把那两个关于Q的微分方程转成C*的可能会好些 求导的过程可以利用matlab的符号求解功能，虽然形式复杂点，但应该不影响一般的求解 你可以再试试你的方法，如果还是不理想，可以试试转成C*的方程试试看 |
» 本帖已获得的红花（最新10朵）
小木虫zb10楼2013-04-26 10:35:37
5