| 查看: 310 | 回复: 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 ] |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
帮忙看看fileCode
已经有5人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
eezyzyzzzee
木虫 (著名写手)
- 应助: 45 (小学生)
- 金币: 3127.3
- 散金: 110
- 红花: 17
- 帖子: 2212
- 在线: 296小时
- 虫号: 516652
- 注册: 2008-03-03
- 专业: 可再生与替代能源利用中的
3楼2009-09-26 09:05:21
wuguocheng
荣誉版主 (职业作家)
- 应助: 2 (幼儿园)
- 贵宾: 7.049
- 金币: 3676.6
- 散金: 316
- 红花: 22
- 沙发: 1
- 帖子: 3032
- 在线: 353.1小时
- 虫号: 623107
- 注册: 2008-10-11
- 性别: GG
- 专业: 力学中的基本问题和方法
- 管辖: 土木建筑
2楼2009-09-26 08:29:11
100