| 查看: 238 | 回复: 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 ] |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
wuguocheng
荣誉版主 (职业作家)
- 应助: 2 (幼儿园)
- 贵宾: 7.049
- 金币: 3674.1
- 散金: 316
- 红花: 22
- 沙发: 1
- 帖子: 3032
- 在线: 352.5小时
- 虫号: 623107
- 注册: 2008-10-11
- 性别: GG
- 专业: 力学中的基本问题和方法
- 管辖: 土木建筑
2楼2009-09-26 08:29:11
eezyzyzzzee
木虫 (著名写手)
- 应助: 45 (小学生)
- 金币: 3126.5
- 散金: 110
- 红花: 17
- 帖子: 2212
- 在线: 296小时
- 虫号: 516652
- 注册: 2008-03-03
- 专业: 可再生与替代能源利用中的
3楼2009-09-26 09:05:21