24小时热门版块排行榜    

查看: 924  |  回复: 15

mshwangg

至尊木虫 (正式写手)

★ ★ ★ ★
十点钟的咖啡(leongoall代发): 金币+4, 多谢参与讨论! 2014-06-28 19:29:45
引用回帖:
10楼: Originally posted by walk1997 at 2014-06-28 05:27:17
不确定 因为我也没实践过
你可以动手试试
(1,2维的解起来应该效率还好)
我猜测:
   后者也选用合适的基展开的话  使积分能解析解出来
   这样  就只需要做矩阵对角化之类的工作了

另外 我感觉 任何的 ...

理论上同意你得看法
要数值求解,还得积分精度够高才行
walk有好一点的积分算法推荐一下?
我用复合辛普森积分,感觉不怎么样,速度慢不说,精度也不太满意
11楼2014-06-28 11:25:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
十点钟的咖啡(leongoall代发): 金币+10, 多谢应助! 2014-06-28 19:29:21
引用回帖:
11楼: Originally posted by mshwangg at 2014-06-28 11:25:58
理论上同意你得看法
要数值求解,还得积分精度够高才行
walk有好一点的积分算法推荐一下?
我用复合辛普森积分,感觉不怎么样,速度慢不说,精度也不太满意...

我算法差不多一点都不会呀  只会用mathematica做数值积分 : (
我感觉新版本的mathematica数值积分速度已经快追上c的包了
以前是调用cernlib中的积分库(fortran下)
现在mathematica下也可以调用c写的cuba包  和自带的感觉差不多
(这个包应该也提供c和fortran的调用  cuba提供了4种算法)

很高维的积分估计一般还只能用MC (如果被积函数积分结果没真实奇点的话)
(原生的不好的话,自己产生特定的分布试试)
(粒子中有些积分带真实奇点的多维积分维的柯西主值积分)

我见过比较强大的计算方法算是:
   利用多维"复平面"上的柯西定理 结合MC做积分 (自己产生特定权重的分布)
---- 我重复不出来: (
12楼2014-06-28 12:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)

引用回帖:
12楼: Originally posted by walk1997 at 2014-06-28 12:44:00
我算法差不多一点都不会呀  只会用mathematica做数值积分 : (
我感觉新版本的mathematica数值积分速度已经快追上c的包了
以前是调用cernlib中的积分库(fortran下)
现在mathematica下也可以调用c写的cuba包  和自 ...

你太猛,跟不上你的节奏
自己写的算法真心比不上那些函数库和软件
找时间研究下你说的
13楼2014-06-28 15:49:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mahong123

禁虫 (正式写手)

感谢参与,应助指数 +1
本帖内容被屏蔽

14楼2014-06-28 16:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)

引用回帖:
13楼: Originally posted by mshwangg at 2014-06-28 15:49:42
你太猛,跟不上你的节奏
自己写的算法真心比不上那些函数库和软件
找时间研究下你说的...

缪赞了
我只会动嘴皮子(还可能都不着边)
真动手不行 呵呵
15楼2014-06-28 19:43:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

十点钟的咖啡

金虫 (正式写手)

引用回帖:
5楼: Originally posted by walk1997 at 2014-06-27 20:58:31
把Eq.2中t+\Delta t 换成t', t'标记成t''
这样就得到\rho(t')的积分表达式
把这1表达式代入Eq.2中的右边  就得到3项
第一项移左边  第三项交换下积分和[]次序 就得到eq.3

还真的是真么回事,原来这是微扰。下面的同志们,也谢谢你们的回答,不过我不是要做数值分析。只是形式求解而已。
16楼2014-06-29 13:41:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 十点钟的咖啡 的主题更新
信息提示
请填处理意见