| 查看: 308 | 回复: 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 ] |
回复此楼» 猜你喜欢
据悉今年马上要出结果了
已经有4人回复
-
2026国自然放榜时间
已经有6人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有7人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有15人回复
-
面上项目filecode邪修
已经有7人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
wuguocheng
荣誉版主 (职业作家)
- 应助: 2 (幼儿园)
- 贵宾: 7.049
- 金币: 3676.6
- 散金: 316
- 红花: 22
- 沙发: 1
- 帖子: 3032
- 在线: 353小时
- 虫号: 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
20