| 查看: 329 | 回复: 2 | |||
| 当前主题已经存档。 | |||
monitor2885至尊木虫 (知名作家)
队长
|
[交流]
【讨论】pdepe函数的边界条件
|
||
|
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) %边界条件 abc=[1,5,9,8,0,4,6,8,7,0] pl=Ql-abc; ql=0; pr=Qr-2.5; qr=0; x=0:10:1000; t=0:0.1:1; m=0; sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t); Q= sol(:,:,1); surf(x,t,Q) 我的问题是，左侧边界条件要每一个时间步长取一个数值，但是运行结果程序只是取了一个值，后面9个值都被忽略了。怎么才能让左侧边界的10个数值都能取到呢？ [ Last edited by monitor2885 on 2009-9-20 at 16:24 ] |
回复此楼» 猜你喜欢
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
-
国自然结果
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有9人回复
-
好奇怪的filecode
已经有5人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有21人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有10人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
liang_2004
木虫 (正式写手)
始作俑者
- 应助: 16 (小学生)
- 金币: 1803.3
- 散金: 56
- 红花: 5
- 帖子: 377
- 在线: 124.4小时
- 虫号: 848455
- 注册: 2009-09-15
- 专业: 无机非金属类光电信息与功
2楼2009-09-19 19:38:42
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★
kuhailangyu(金币+3,VIP+0):非常感谢参与讨论解答 9-21 18:52
nono2009(金币+3,VIP+0):谢谢专家！ 9-21 18:53
nono2009(金币+0,VIP+0):和K版撞分了 :) 9-21 18:54
kuhailangyu(金币+3,VIP+0):非常感谢参与讨论解答 9-21 18:52
nono2009(金币+3,VIP+0):谢谢专家！ 9-21 18:53
nono2009(金币+0,VIP+0):和K版撞分了 :) 9-21 18:54
|
function [pl,ql,pr,qr] = pdebc(xl,Ql,xr,Qr,t) %边界条件 abc=[1,5,9,8,0,4,6,8,7,0]; switch t case ... pl=... case... pl=... .... end ql=0; pr=Qr-2.5; qr=0; |
3楼2009-09-21 18:38:04
100