24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1712  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕322分 +8 哈哈哈吼吼吼哈 2026-04-04 8/400 2026-04-05 14:13 by brian0809
[考研] 求调剂 +9 Hll胡 2026-04-04 9/450 2026-04-05 13:31 by wwytracy
[考研] 270求调剂 +9 小杰pp 2026-03-31 11/550 2026-04-05 11:02 by 风雨无晴
[考研] 283分求调剂 +9 试试看呗 2026-04-04 9/450 2026-04-05 10:27 by 果冻大王
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 0703求调剂383分 +8 W55j 2026-04-03 8/400 2026-04-04 20:09 by xhai2011
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 一志愿华南师范361分,化学求调剂 +7 Nicole88888 2026-04-01 7/350 2026-04-04 18:28 by macy2011
[考研] 求调剂:085600材料与化工,考材科基,总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 303求调剂 +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by ls刘帅
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 348环境工程调剂 +3 吴彦祖24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
[考研] 求调剂,一志愿南京师范大学计算机专硕,初试373,六级通过, +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 285求调剂 +11 AZMK 2026-04-01 11/550 2026-04-01 22:40 by peike
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
信息提示
请填处理意见