| 查看: 313 | 回复: 2 | |||
| 当前主题已经存档。 | |||
monitor2885至尊木虫 (知名作家)
队长
|
[交流]
【求助】解偏微分方程的方法
|
||
|
偏微分方程用差分方法和特征线法，分别用matlab的什么函数？谢谢。贴图是我的例子，我根据pdepe编写的，按理说是波动方程，就像正弦曲线一样，往下传播才对，但是根本不传递这种波动性。不知道哪里写错了。（左侧边界是一个随着时间变化的函数，我采用了分段函数输入的。） function plotproblem x=0:1000:80000; t=0:0.007:1; m=0; options=odeset('AbsTol',1e-16); sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t,options); Q= sol(:,:,1); surf(x,t,Q) function [c,f,s]=pdefun(x,t,Q,dQdx) %目标pde函数 c=1.884*Q^(-0.4); f=-Q; s=0; function Q0 = pdeic(x) %初始条件 Q0=2.5; function [pl,ql,pr,qr] = pdebc(xl,Ql,xr,Qr,t) %边界条件 if (0<=t&t<0.252) p1=-162.19;p2=40.61;p3=2.4005; y=p1*t^2+p2*t+p3; else y=2.5; end left=y; pl=left; ql=1; pr=2.5; qr=1; [ Last edited by monitor2885 on 2009-9-26 at 14:40 ] |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
wuguocheng
荣誉版主 (职业作家)
- 应助: 2 (幼儿园)
- 贵宾: 7.049
- 金币: 3793.9
- 散金: 316
- 红花: 22
- 沙发: 1
- 帖子: 3060
- 在线: 360.9小时
- 虫号: 623107
- 注册: 2008-10-11
- 性别: GG
- 专业: 力学中的基本问题和方法
- 管辖: 土木建筑
2楼2009-09-26 08:29:11
eezyzyzzzee
木虫 (著名写手)
- 应助: 45 (小学生)
- 金币: 3127.3
- 散金: 110
- 红花: 17
- 帖子: 2212
- 在线: 296小时
- 虫号: 516652
- 注册: 2008-03-03
- 专业: 可再生与替代能源利用中的
3楼2009-09-26 09:05:21
30