| 查看: 2060 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
ghw_nit铁杆木虫 (正式写手)
|
[求助]
一个偏微分方程组的求解问题 已有1人参与
|
|
在一篇文献中看到的某系统的方程，及相应的一些条件如图所示，文章中说用pdepe就可以求解，但是我就是解不出来呀，请大家看看有没有什么问题。 360截图-3516750.jpg |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 偏微分方程求解！
已经有4人回复
- 求推荐求解偏微分方程组的书
已经有3人回复
- 一道方程组的求解
已经有16人回复
- 哪位大神给指导下如果求解一个偏微分方程组
已经有11人回复
- 使用matlab怎么解微分方程组啊？
已经有5人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 一个耦合的偏微分方程组怎么数值解？
已经有5人回复
- 求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 抛物型偏微分方程的数值解求解问题，谢谢！
已经有14人回复
- 二阶非线性微分方程求解
已经有14人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 一个微分方程组，求解
已经有18人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 求解一个含字母系数的微分方程
已经有5人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】这个微分方程组simulink求解到底错在哪里？
已经有3人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】偏微分方程的基本解
已经有5人回复
- 【分享】PDECOL-经典数值求解非线性偏微分方程组源程序
已经有142人回复
ghw_nit
铁杆木虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 7901.9
- 散金: 147
- 红花: 6
- 帖子: 653
- 在线: 175.8小时
- 虫号: 774851
- 注册: 2009-05-20
- 性别: GG
- 专业: 控制理论与方法
6楼2014-02-08 14:58:06
ghw_nit
铁杆木虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 7901.9
- 散金: 147
- 红花: 6
- 帖子: 653
- 在线: 175.8小时
- 虫号: 774851
- 注册: 2009-05-20
- 性别: GG
- 专业: 控制理论与方法
|
这是我编写的matlab程序，大家都可以调一下，看看怎么样，总是说步长不够小，我已经调到很小了，比文献中说的要小好多，都是不能求解的，不知道是怎么回事。 function [c,f,s]=pdefun(x,t,u,dudx) Pe_C=1; Pe_T=1; B_C=1; B_T=2; gamma=20; beta_T=1; u_c=0; y=u(2)/(1+u(2)/gamma); R_11=B_T*B_C*(1-u(1))*exp(y); R_12=-beta_T*(u_c-u(2)); R_13=-B_C*(1-u(1))*exp(y); c=[1 1]'; f=[1/Pe_C 1/Pe_T]'.*dudx; s=[R_13-dudx(1) R_11+R_12-dudx(2)]'; function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t) pl=[-ul(1),-ul(2)]'; ql=[1 1]'; pr=[0 0]'; qr=[1 1]'; function u0=icfun(x) u0=[0 0]'; m=0; t=linspace(0,1,500); x=linspace(0,1,100); sol=pdepe(m,@pdefun,@icfun,@bcfun,x,t); |
2楼2014-01-23 19:57:37
ghw_nit
铁杆木虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 7901.9
- 散金: 147
- 红花: 6
- 帖子: 653
- 在线: 175.8小时
- 虫号: 774851
- 注册: 2009-05-20
- 性别: GG
- 专业: 控制理论与方法
3楼2014-01-23 20:01:05
5楼2014-02-06 10:57:55
20