| 查看: 1984 | 回复: 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
» 猜你喜欢
25级博士招生
已经有9人回复
-
求助大佬: 多晶薄膜xrd显示织构, 但chi扫和摇摆曲线无峰
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有55人回复
-
求ICSD file no. 35-6351标准卡片
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】matlab求解偏微分方程
已经有4人回复
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- matlab求解偏微分方程组
已经有3人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- 偏微分方程数值解有基本理论和一些例子（化工类）MATLAB源程序
已经有294人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】Matlab解分数阶微分方程问题
已经有18人回复
- 【求助】用pdepe求解偏微分方程
已经有10人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】如何求解下列偏微分方程
已经有13人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
- 【求助】求助用matlab解偏微分方程
已经有15人回复
- 【求助】偏微分方程的matlab计算
已经有3人回复
- 【求助】用matlab编过有限差分法的请进
已经有9人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68035.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3290.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-04-25 16:25:30
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.1
- 红花: 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.1
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
6楼2013-04-25 22:24:56
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.1
- 红花: 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.1
- 红花: 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
50