应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 原子、分子&光学 » iteration method for density equations_help_me
162/2返回列表上一页12
查看: 924  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的回帖
相关版块跳转 我要订阅楼主 十点钟的咖啡 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭