24小时热门版块排行榜    

查看: 1682  |  回复: 5

ghw_nit

铁杆木虫 (正式写手)

[求助] 一个偏微分方程组的求解问题 已有1人参与

在一篇文献中看到的某系统的方程,及相应的一些条件如图所示,文章中说用pdepe就可以求解,但是我就是解不出来呀,请大家看看有没有什么问题。
一个偏微分方程组的求解问题
360截图-3516750.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

这是我编写的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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

这个方程应该是没有什么问题的,这是一个温度与浓度的传播方程,好多篇文献中都以这个方程为例进行讨论,尽管写成了不同的形式,但都是这个方程,不管什么形式我都求解不出来呀,怎么回事,大家看看。如果大家有兴趣我也可以相关文章贴出来。
3楼2014-01-23 20:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

自己顶一下
4楼2014-01-25 07:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfg25aw

新虫 (初入文坛)

【答案】应助回帖

最好用迎风格式去解。
5楼2014-02-06 10:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

有人给我解答了,模型本身存在问题,因此求解不出,只是我看到的这个模型在很多的文献中都被作为例子进行讨论,而且不同的论文我都找了一下,没有一篇我是能够仿出来,真是奇怪,不知道这些作者是怎么弄的,我是无论如何也搞不出来。
这些论文中的参数没有差别,只是在最后的常数项前面的正负号有差别,我都试过了,还是不行呀。哎放弃了。
6楼2014-02-08 14:58:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ghw_nit 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +6 墨墨漠 2026-03-20 6/300 2026-03-21 09:23 by 天青等烟雨
[考研] 268求调剂 +7 简单点0 2026-03-17 7/350 2026-03-21 08:32 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见