24小时热门版块排行榜    

查看: 1676  |  回复: 11

有机物

至尊木虫 (著名写手)

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

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



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

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

金虫 (职业作家)

★ ★ ★ ★
有机物(金币+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的回帖

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的回帖

onesupeng

金虫 (职业作家)

你要绝对稳定,你就按照绝对稳定的格式去做就行了,不过求解代数方程稍微麻烦一点而已了,这个每个人做都有自己的平衡点

四维和二维、三维没有本质区别,你会做三维的,我没有明白你四维的难度何在?

我一贯坚持:先能实现,再讲求优化和效率

[ Last edited by onesupeng on 2010-7-3 at 13:12 ]
长期招收博士生,参见http://fsl-unsw.com
6楼2010-07-03 13:09:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiahl

金虫 (小有名气)


有机物(金币+1):谢谢参与
可以试试有限元,其实从理论角度来讲,你的这个方程就是一个热传导方程。系数如果退化的话可能会稍微有点困难。解PDE的话边界条件和初值对最终的解影响很大。所以你要多朝这方面考虑问题。
业精于勤 荒于嬉 行成于思 毁于随
7楼2010-07-03 15:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

引用回帖:
Originally posted by jiahl at 2010-07-03 15:19:11:
可以试试有限元,其实从理论角度来讲,你的这个方程就是一个热传导方程。系数如果退化的话可能会稍微有点困难。解PDE的话边界条件和初值对最终的解影响很大。所以你要多朝这方面考虑问题。

人家连FDM都嫌慢,再去做FEM那更慢。FDM对边界处理没有那么弱
长期招收博士生,参见http://fsl-unsw.com
8楼2010-07-04 01:16:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)

有机物(金币+3):谢谢回复 2010-12-21 20:28:33
显式格式的稳定性可以通过调节时间步和空间步得到。看问题似乎离散之后不会太特殊。如果计算代价在求解线性子问题上,我还是建议看看矩阵结构,看看能否优化条件加速。这方面的研究比较活跃,可以用的知识比较多,而且实现起来比较简单。
9楼2010-07-04 05:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

有机物(金币+14):多谢! 2010-12-21 20:28:50
其实我觉得用显示推进完全满足他的应用需求,但不知道为什么要要求一个绝对稳定的要求

前面我列出的那个格式就是三阶,在求解对流扩散方程上屡试不爽,就连有间断的情形,也有很多人用这个格式,非常快捷编程也很容易。(这个根本无须迭代)

非要用隐式格式也可以,比如交替迭代等。搞什么一个矩阵Ax=B的求解多麻烦,而且我不相信这个方程会比我们遇到的压力泊松方程难解。
长期招收博士生,参见http://fsl-unsw.com
10楼2010-07-04 05:34:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 有机物 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 材料投稿速度快的二三区期刊有哪些! +3 加油努力7 2024-06-24 4/200 2024-06-27 18:09 by keyaner23
[论文投稿] 关于论文第一通讯作者的问题? +7 winsaint 2024-06-26 7/350 2024-06-27 17:23 by 莱茵润色
[基金申请] 75批面上出结果了么?系统卡的进不去了 +4 wziguang33 2024-06-27 4/200 2024-06-27 15:53 by xter9665
[基金申请] 国基在研影响申请结果吗 +14 WOWO159357 2024-06-26 22/1100 2024-06-27 15:36 by lucky_my2010
[基金申请] 请问三个字的中文名字,在文中到底怎么表示? 5+3 zhongyuan764 2024-06-22 7/350 2024-06-27 11:28 by f1yfish
[基金申请] 博后面上今天下午会公布吗?有无消息? +26 hajkdfdf 2024-06-24 45/2250 2024-06-27 09:19 by kyukitu
[基金申请] 博后网上办公系统维护 +14 取名字烦人 2024-06-25 22/1100 2024-06-26 22:15 by tuanbao2023
[基金申请] 自然基金变化?? +15 今晚推荐22 2024-06-21 25/1250 2024-06-26 19:40 by 朋友别哭010
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[硕博家园] 博士该不该读? +8 L1009225316 2024-06-25 8/400 2024-06-26 16:26 by huixiong0627
[教师之家] 神奇的中医 +8 水冰月月野兔 2024-06-24 10/500 2024-06-26 12:10 by lyfbangong
[基金申请] 今天能不能出来名单 +8 地理学1995 2024-06-25 10/500 2024-06-26 09:46 by msjy
[有机交流] 过柱子,产品,杂质在是 pe:ea=100:1 也一起出来? +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[硕博家园] 申博雅思求助 +6 刘嘉研究究 2024-06-24 7/350 2024-06-26 07:08 by 小李醒yy
[考博] 申博好难 +6 自强不息a?a 2024-06-24 6/300 2024-06-25 23:02 by 考研吒儿
[教师之家] 复旦夏同学退学理由说明,哪儿可以下载? +7 苏东坡二世 2024-06-21 8/400 2024-06-25 12:42 by 药嘿1233
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[教师之家] 复旦夏同学提出高校成年人学生退学不应该让家长审核,大家认同吗? +10 苏东坡二世 2024-06-22 17/850 2024-06-24 16:52 by wanghuawei
[基金申请] 国自然青年基金,1A4B能上会吗?青年和面上的上会标准是一样的吗? +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[基金申请] 教育部基金 +5 m1393 2024-06-21 5/250 2024-06-21 21:13 by odes
信息提示
请填处理意见