24小时热门版块排行榜    

查看: 913  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

土木小和尚

木虫 (正式写手)

[求助] 求高手帮忙解偏微分方程

偏微分方程如下:

偏微分方程.jpg
好几个月了,一直弄不出来!matlab跟maple都试过了!
下面是用maple写的,但是出错了,麻烦帮忙看下:

Maple.jpg



[ Last edited by 土木小和尚 on 2013-4-8 at 16:58 ]
回复此楼
佛说:三千繁华,弹指刹那,百年过后,不过一捧黄沙。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

定解问题:
Cb*P^2u/Pz^2=u^5*Pu/Pt                                       (1)
    t=0, 0≤z≤H  u=u0                                         (2)
0<t<∞, z=0    u=0                                          (3)
0<t<∞, z=H   Pu/Pz=0                                       (4)
t=∞, 0≤z≤H    u=0                                          (5)

求解过程:
使用分离变量发求解,令u(z,t)=X(z)*T(t),代入原方程得到:
Cb*T*d^2X/dz^2=X^5*T^5*X*dT/dt
Cb*X^(-6)*d^2X/dz^2=T^4*dT/dt                                  (6)
由于上式左边是z的函数,右边是t的函数,两者衡等意味着双方都等于某一常数A。
故          Cb*X^(-6)*d^2X/dz^2= A                         (7)
  和          T^4*dT/dt = A                                  (8)
由(8)   T=(5*A*t+B)^(1/5)                                    (9)
由(7),令Y=dX/dz
    则 d^2X/dz^2=dY/dz=dy/dX*dX/dz=Y*dy/dX                   (10)
(10)代入(7)得到:
Cb *X^(-6)*Y*dY/dX= A
故 Y^2=2*A/(7*Cb)*X^7 + D
即dX/dz=±sqrt[2*A/(7*Cb)*X^7 + D]
Integral{dX/{ sqrt[2*A/(7*Cb)*X^7 + D]}} = ±z + F

以上各个式子中的A、B、D、F均为积分常数。
由于t趋向无穷 u=0,故A=0 故 T=B^(1/5)=G
故X=D*z+F
以上各个式子中的A、B、D、F、G均为积分常数。
故u(z,t)=G*(D*z+F)
由0<t<∞, z=0    u=0 和0<t<∞, z=H  Pu/Pz=0,得:
u(z,t)≡0

这个定解问题的描述有问题吧?怎么只有零解?
15楼2013-12-01 19:37:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

SEEDEN

银虫 (正式写手)

随便算了一下,不知道咱算的对不对。。。。

一阶向后差分,将[0,H]离散100段,时长是你提的864000s
当然,没试其他的离散格式。。。。
而且很明显可以看出在初始阶段,速度增加很快。。。

计算的结果.jpg

2楼2013-04-08 20:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

土木小和尚

木虫 (正式写手)

引用回帖:
2楼: Originally posted by SEEDEN at 2013-04-08 20:22:44
随便算了一下,不知道咱算的对不对。。。。

一阶向后差分,将离散100段,时长是你提的864000s
当然,没试其他的离散格式。。。。
而且很明显可以看出在初始阶段,速度增加很快。。。

计算的结果.jpg
...

首先表示感谢!但是图像的变化趋势不对!两边都应该是下降
佛说:三千繁华,弹指刹那,百年过后,不过一捧黄沙。
4楼2013-04-09 10:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

土木小和尚

木虫 (正式写手)

引用回帖:
3楼: Originally posted by hfj1988 at 2013-04-08 20:43:09
非线性方程啊

是啊
佛说:三千繁华,弹指刹那,百年过后,不过一捧黄沙。
5楼2013-04-09 10:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见