24小时热门版块排行榜    

查看: 1426  |  回复: 6
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

superdirac

木虫 (正式写手)

[交流] 【求助】想问个velocity-corrected verlet 算法的问题

v(t)={v(t+dt/2)+v(t-dt/2}/2 - dt*{a(t+dt)-a(t-dt)}/12 + O(dt**4)

这个式子是怎么推出来的??

谢谢!

假如将

v(t+dt/2),v(t-dt/2),a(t+dt),a(t-dt)作泰勒展开,让后合并同类项:  得出的是

v(t)={v(t+dt/2)+v(t-dt/2}/2 - dt*{a(t+dt)-a(t-dt)}/16 + O(dt**4)

两者一直对不上!

找个高手指点一下! 呵呵!谢谢!

[ Last edited by lei0736 on 2009-11-25 at 12:07 ]
回复此楼
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superdirac

木虫 (正式写手)

呵呵! 搞错了!
  这个不是 预测--校正算法
  这个叫  velocity-corrected verlet算法.
  对速度的计算可以达到4阶的精度.
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
3楼2009-03-18 22:04:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

老虎大王

木虫 (著名写手)

★
xuefei06(金币+1,VIP+0):3Q! 3-19 10:44
预报-校正算法也有很多版本的。你提的这个公式我好像还没有见过。
2楼2009-03-18 20:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★
xuefei06(金币+1,VIP+0):3Q! 3-19 10:45
呵呵。我今天又看了一下,你这个叫做改进型Beeman算法。它和Verlet算法是等效的。
4楼2009-03-19 10:17:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superdirac

木虫 (正式写手)

★ ★
xuefei06(金币+2,VIP+0):3Q! 3-19 14:59
呵呵! 这个不是 beeman

http://en.wikipedia.org/wiki/Beeman's_algorithm
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
5楼2009-03-19 11:15:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 9/450 2026-11-11 04:54 by Z7qfCsdIiD1T
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:40 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 04:32 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:12 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 20:07 by Z7qfCsdIiD1T
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
信息提示
请填处理意见