| 查看: 2382 | 回复: 6 | |||
小木虫zb木虫 (正式写手)
|
[求助]
matlab ode15s方程中两个参数怎么调用其他函数带入求解
|
|
本人用ode15s算一个偏微分方程组，该方程组中有两个参数要通过另外一个隐函数方程组求解，不知道怎么调用，ode15s求该解偏微分方程组代码（[t,y]=ode15s(@fangcheng,tspan,y0,options);）如下： 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); q1=y(1:N);c1=y(N+1:2*N);q2=y(2*N+1:3*N);c2=y(3*N+1:4*N); %-----------定义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=[dq1dt dc1dt dq2dt dc2dt]'; 其中ce1(j)、ce2(j)、q1(j)和q2(j)的关系是q1(j)=0.45967ce1(j)^0.3839/(ce1(j)^0.0457+0.3*ce2(j)^0.0003206，q2(j)=0.17356*ce2(j)^(-1.55201)/(ce2(j)^(-2)+0.1*ce1(j)^(-1)) 如果用q1、q2表示c1、c2带入，程序是可以解出的，但是这是个隐函数方程组，不知道该怎么办，求大神相助 |
回复此楼» 猜你喜欢
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有109人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有64人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab ode求助
已经有1人回复
- 关于MATLAB用ode15s函数的问题
已经有0人回复
- matlab的ODE里面怎么写卷积积分
已经有0人回复
- Matlab ode45的运用
已经有0人回复
- matlab 中ode15s 的代码
已经有5人回复
- Matlab ode45求解非线性方程
已经有2人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- matlab编程 ODE45
已经有4人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- matlab的ode解微分方程
已经有1人回复
- 新虫求助，matlab使用ode求解模型问题
已经有0人回复
- [求助Matlab ode15s仿真步长的问题]
已经有0人回复
- 分歧理论可否用matlab的ode45求解
已经有2人回复
- 求助Matlab ode45 求解方程有时能求解，有时无法收敛，是什么情况？
已经有2人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- matlab的ode45的初值问题
已经有2人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 用matlab的ode45求解Liouville von-neumann方程的问题
已经有1人回复
2楼2013-06-04 06:22:44
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
| 偏微分方程组图片中所示，终于有人回复了，谢谢你啊！ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 360截图20130604164648097.jpg
2013-06-04 16:52:45, 28.68 K
3楼2013-06-04 16:53:49
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
4楼2013-06-04 16:57:02
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68175.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
5楼2013-06-04 19:35:26
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
|
嗯，对，就是你所说的，只是离散的方程中有个两个参数ce1(j)、ce2(j)需要通过方程q1(j)=0.45967ce1(j)^0.3839/(ce1(j)^0.0457+0.3*ce2(j)^0.0003206，q2(j)=0.17356*ce2(j)^(-1.55201)/(ce2(j)^(-2)+0.1*ce1(j)^(-1))计算得到，不知道该如何调用，http://muchong.com/bbs/viewthread.php?tid=5950673是我编的一个程序，程序没错误但是运行不出结果，你能帮下我吗，这个问题困扰我很久了，之前也收到过你的回复，你应该在这方面挺懂的 |
6楼2013-06-04 20:21:15
7楼2013-06-05 05:32:02
