24小时热门版块排行榜    

查看: 1213  |  回复: 10
当前主题已经存档。

oytxtu

铁杆木虫 (正式写手)

[交流] 【求助】如何通过龙格库塔法求解哈密顿方程(已解决)

请问:如何通过龙格库塔法求解含时的哈密顿方程
      龙格库塔的方法我知道,但是,如何通过这个方法求解哈密顿方程,还很是模糊,麻烦各位给我一些指导!

[ Last edited by oytxtu on 2008-11-10 at 16:44 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvzhu2007

荣誉版主 (著名写手)


zt970831(金币+1,VIP+0):感谢您的交流
龙格库塔只是一种解微分方程的数值算法而已
在物理上的应用要具体问题具体分析的
2楼2008-09-04 20:57:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oytxtu

铁杆木虫 (正式写手)

zt970831(金币+0,VIP+0):可以用图片以附件上传
谢谢你的回复,
小木虫里面不能编公式,我没办法给你看方程
我只想知道,它是如何进行迭代的
这就是我的困难!
3楼2008-09-04 21:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quantumoptics2

金虫 (小有名气)



zt970831(金币+1,VIP+0):感谢您的答复
用MATLAB应该可以,直接用函数
4楼2008-09-04 23:42:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvzhu2007

荣誉版主 (著名写手)

我是自己编了个程序算的,是格点模型
 
我想一般的东西可以简单一些,有很多现成的代码可以用,你可以找些数值分析的书看,不难的
引用回帖:
Originally posted by oytxtu at 2008-9-4 21:54:
谢谢你的回复,
小木虫里面不能编公式,我没办法给你看方程
我只想知道,它是如何进行迭代的
这就是我的困难!

5楼2008-09-05 11:19:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyclone985

银虫 (小有名气)

★ ★
lvzhu2007(金币+2,VIP+0):感谢交流!
可以看看 Fortran 的IMSL 链库函数。
6楼2008-09-10 14:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stereochemistry

回复的中没有一个有实用价值的,其实楼主的问题简单的很,对于做动力学的来说,其实就是解了核的运动方程,所谓的经典轨迹。不过牵扯到直角坐标和广义坐标的转化,细节也不算太难。呵呵。我估计的。
7楼2008-10-09 21:46:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamsad

至尊木虫 (著名写手)

只說一點, 這個是微分方程組, 個數比較多, 計算要小心
8楼2008-10-10 01:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvzhu2007

荣誉版主 (著名写手)

引用回帖:
Originally posted by stereochemistry at 2008-10-9 21:46:
回复的中没有一个有实用价值的,其实楼主的问题简单的很,对于做动力学的来说,其实就是解了核的运动方程,所谓的经典轨迹。不过牵扯到直角坐标和广义坐标的转化,细节也不算太难。呵呵。我估计的。

如果是单纯地在一个连续的空间解这个方程的话 显然套现成的数值代码就可以了

如果涉及到物理模型 恐怕不是那么简单吧 至少边界条件是个问题 大多数时候绝对不是细节问题

[ Last edited by lvzhu2007 on 2008-10-10 at 07:34 ]
9楼2008-10-10 07:31:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

引用回帖:
Originally posted by oytxtu at 2008-9-4 15:21:
请问:如何通过龙格库塔法求解含时的哈密顿方程
      龙格库塔的方法我知道,但是,如何通过这个方法求解哈密顿方程,还很是模糊,麻烦各位给我一些指导!

[ Last edited by oytxtu on 2008-9-4 at 21:55 [ ...

同意!如果是初值问题,只要按照龙格库塔法所要求的形式把哈密顿方程的形式改写一下的就行了~
10楼2008-11-10 13:00:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oytxtu 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 oepdf5z972 2026-06-20 3/150 2026-06-21 21:01 by zo4669e5n4
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[硕博家园] 一篇论文同时出现在两个期刊,一模一样,这算不算学术不端,请各位老师斧正。 +11 mahl005 2026-06-14 12/600 2026-06-18 17:01 by 770904134
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
请填处理意见