24小时热门版块排行榜    

查看: 2499  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

有机物

至尊木虫 (著名写手)

[交流] 【求助】高维PDE的求解 已有6人参与

请教:如下的高维偏微分方程如何解?如果用有限差分法,电脑算起来太慢。现在想寻求近似解析解,以前没解过高维PDE,有没有朋友解过类似的方程,请提宝贵经验,谢谢。



[ Last edited by javeey on 2010-7-1 at 23:14 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

★ ★ ★ ★
有机物(金币+1):谢谢参与
javeey(金币+3):谢谢参与交流 2010-07-02 07:27:38
实际上你的问题FDM比较快了,可能一般情况下没有简单的近似方法。

如果你的精度不是非常高的话,可以分步差分

也可以采用交替迭代

总之不会太难。也可以按照楼上的说法,不过做起来更麻烦一些。


不过我觉得采用3阶TVD时间推进格式,应该很容易做,而且快的不行。就不贴代码了阿。


另外,楼主时计算生物的阿?我可以合作做很多研究阿,怕你觉得我水平不够~




[ Last edited by onesupeng on 2010-7-2 at 02:11 ]
长期招收博士生,参见http://fsl-unsw.com
3楼2010-07-02 01:29:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

saladin983

铁杆木虫 (正式写手)

★ ★ ★
有机物(金币+1):谢谢参与
javeey(金币+2):谢谢参与交流。我能看到图片 2010-07-02 07:27:14
这两天有时候看不到小木虫图片服务器的图,真是怪事!

偏微分方程的数值方法是很大的领域,看不到图,就只能泛泛而谈了。

如果算法能求解,那么最简单的解决办法还是提高计算速度。

从硬件上来说就是用更快的机器算。

从算法本身来看,preconditioning是非常有效的提速办法。比如子问题里用了Krylov子空间迭代法的情况下。

从算法实现的角度来看,优化代码有可能会有惊喜。比如用matlab的时候降低关键for循环的使用对提速的效果很明显。

先mark一下,等我什么时候能看到图了再回头看看。
2楼2010-07-02 00:18:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

补充一个,做研究之道是先能实现,在讲求优化和效率,否则觉得别人不怎么样却做不出来。
长期招收博士生,参见http://fsl-unsw.com
4楼2010-07-02 01:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

有机物

至尊木虫 (著名写手)

引用回帖:
Originally posted by onesupeng at 2010-07-02 01:29:45:
实际上你的问题FDM比较快了,可能一般情况下没有简单的近似方法。

如果你的精度不是非常高的话,可以分步差分

也可以采用交替迭代

总之不会太难。也可以按照楼上的说法,不过做起来更麻烦一些。


不 ...

非常感谢。惭愧!PDE解得不多。我曾经用Crank-Nicolson型有限差分法解过3维的PDE,已经很慢了。也许是我的电脑性能比较差,当时好像是1G的内存,2.8G的CPU。

我要解的这个高维PDE,精度可以不用很高,但希望算法能有一定的稳定性,最好是绝对稳定的。我查阅了一些解PDE的资料,都没有介绍高维PDE。您在图中给出的计算方法哪本书有介绍,或者哪篇文献有介绍?谢谢。

生物是我打算做的另一个方向,很惭愧,刚开始入门。主要想做神经性疾病,以后可以多多交流,
5楼2010-07-03 12:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:30 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:37 by hZiFeudZyoGR
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:35 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-30 19:26 by hZiFeudZyoGR
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:17 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:07 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 17:34 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
请填处理意见