24小时热门版块排行榜    

查看: 3241  |  回复: 9

hitgysh

银虫 (正式写手)

[交流] 【求助】用matlab编过有限差分法的请进 已有7人参与

有限差分法在matlab中求解偏微分方程是求解的图示结果中有很多尖或者是刺,向上向下都有,非常诡异 完全不是我们所需,结果如图示,目前已经排除了精度问题。请高手告诉一下,谢谢!!!!!!!


回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitgysh

银虫 (正式写手)

是不是有限差分不稳定造成的啊
引用回帖:
Originally posted by hitgysh at 2010-04-27 16:49:02:
有限差分法在matlab中求解偏微分方程是求解的图示结果中有很多尖或者是刺,向上向下都有,非常诡异 完全不是我们所需,结果如图示,目前已经排除了精度问题。请高手告诉一下,谢谢!!!!!!!


[img]http://pic.e ...

[ Last edited by hitgysh on 2010-4-28 at 10:54 ]
2楼2010-04-27 16:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hunter2002

金虫 (正式写手)


nono2009(金币+1):鼓励应助!欢迎常来仿真模拟版交流。 2010-04-30 19:11
参数没有设好,不收敛,可以调整一下,空间,时间步长的比率
3楼2010-04-30 08:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitgysh

银虫 (正式写手)

谢谢你的回帖 我用的是CN那种差分格式 不是无条件收敛吗 只有显式才需要设置网比,是吗?
引用回帖:
Originally posted by hunter2002 at 2010-04-30 08:09:40:
参数没有设好,不收敛,可以调整一下,空间,时间步长的比率

4楼2010-05-01 00:00:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hefbingo

捐助贵宾 (知名作家)


nono2009(金币+1):鼓励应助!欢迎常来仿真模拟版交流。 2010-05-01 11:26
hitgysh(金币+3):谢谢您的指导,那就是说CN的网比应该是R<0.25,是吗? 2010-05-03 20:02
C-N不是无条件收敛的,n时间参数状态和n+1时间参数状态各取了一半,收敛步长只比显性多允许一倍,隐性才是无条件收敛
5楼2010-05-01 01:37:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hefbingo

捐助贵宾 (知名作家)


adu886886(金币+1):谢谢提供意见 2010-05-05 09:28:26
hitgysh(金币+2):谢谢 2010-05-08 10:44:30
隐形,C-N,显性都是指的在时间上积分的不同处理方式。显性收敛的时间步长条件为△t<Γ(△x)^2/2,C-N为△t<Γ(△x)^2,你可以看看《An Introduction to Computational Fluid Dynamics the Finite Volume Method》这本书,作者是Malalasekera,里面讲的挺清楚,这本书网上电子版很好找
6楼2010-05-04 15:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
7楼2011-07-12 11:23:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomeizifly

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by hefbingo at 2010-05-04 15:02:28
隐形,C-N,显性都是指的在时间上积分的不同处理方式。显性收敛的时间步长条件为△t<Γ(△x)^2/2,C-N为△t<Γ(△x)^2,你可以看看《An Introduction to Computational Fluid Dynamics the Finite Volume Method》这 ...

请问,可不可以帮我看看我的问题呢?
http://muchong.com/bbs/viewthread.php?tid=7740033&pid=2#pid2
非常感谢!
8楼2014-08-07 02:51:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

-四-叶-草-

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我最近也在看这方面的知识,求推荐学习资料
9楼2015-07-01 10:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

201110140244

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
请问有MS的DMol板块的手册吗?以及各种输出结果怎么读的资料,谢谢啦!
10楼2016-03-08 20:30:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hitgysh 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见