| 查看: 316 | 回复: 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 ] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
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
