RError.com

RError.com Logo RError.com Logo

RError.com Navigation

  • 主页

Mobile menu

Close
  • 主页
  • 系统&网络
    • 热门问题
    • 最新问题
    • 标签
  • Ubuntu
    • 热门问题
    • 最新问题
    • 标签
  • 帮助
主页 / 问题 / 1287310
Accepted
Fire13nyu
Fire13nyu
Asked:2022-05-29 02:19:20 +0000 UTC2022-05-29 02:19:20 +0000 UTC 2022-05-29 02:19:20 +0000 UTC

Werlet集成方法的实现

  • 772

再会!我想在不使用速度的情况下实现 Verlet 集成方法(https://en.wikipedia.org/wiki/Verlet_integration),但我偶然发现了拥有当前坐标和以前坐标的问题。我正在介绍一个具有三个实例(一维、二维和多维案例)的坐标类:

class Coord coord where  
    infixl 6 +++
    (+++) :: coord -> coord -> coord
    infixl 7 ***
    (***) :: Double -> coord -> coord
    infixl 6 -.-
    (-.-) :: coord -> coord -> coord
    (-.-) coord1 coord2 = coord1 +++ (-1.0)***coord2

instance Coord Double where
    coord1 +++ coord2 = coord1 + coord2
    coord1 *** coord2 = coord1 * coord2

instance Coord (Double, Double) where
    (coord11,coord12) +++ (coord21,coord22) = (coord11+coord21,coord12+coord22)
    num *** (coord11,coord12) = (num*coord11,num*coord12)

instance Coord [Double] where 
    list1 +++ list2 = zipWith (+) list1 list2
    num *** list = map (* num) list

然后我用输入参数做积分函数本身:加速度(函数),初始坐标,初始速度和时间步长,并在每一步之后返回一个无限的点位置列表

verlet :: Coord coord => (coord -> coord) -> coord -> coord -> Double -> [coord]
verlet a x0 v0 dt = 
    let x1 = x0 +++ (dt***v0) +++ (0.5*(dt**2))***(a x0)
        verlet' xPrev xNow = let xNext = 2***xNow -.- xPrev +++ (dt**2))***(a xNow)
                            in (xPrev,xNow) : verlet' xNow xNext
    --- in (,) : verlet'  - вот тут проблема

如果有的话,我是一个真正的初学者,所以,很可能,这很容易解决,但它对我不起作用,所以我寻求帮助)

haskell
  • 2 2 个回答
  • 10 Views

2 个回答

  • Voted
  1. Best Answer
    Fyodor Soikin
    2022-05-30T00:27:22Z2022-05-30T00:27:22Z

    您已经拥有x0,这是给您的,并且x1您根据 Wikipedia 的说明计算得出。然后,按照相同的说明,只需将它们用作 的前两个参数verlet':

    verlet :: Coord coord => (coord -> coord) -> coord -> coord -> Double -> [coord]
    verlet a x0 v0 dt = 
        let x1 = x0 +++ (dt***v0) +++ (0.5*(dt**2))***(a x0)
            verlet' xPrev xNow = let xNext = 2***xNow -.- xPrev +++ (dt**2))***(a xNow)
                                 in (xPrev,xNow) : verlet' xNow xNext
        in verlet' x0 x1
    

    添加:

    从您的评论来看,您需要一个单独的坐标列表作为输出,而不是它们的对。如果是这样,那么只需让函数verlet'完全返回此列表:

    ...
    in xPrev : verlet' xNow xNext
    ...
    
    • 1
  2. Максим Строганов
    2022-05-29T02:31:15Z2022-05-29T02:31:15Z

    verlet :: Coord coord => (coord -> coord) -> coord -> coord -> Double -> [coord] 奇怪的系统...从好的方面来说,它是这样实现的:

    struct Body
    {
        Vec3d pos { 0.0, 0.0, 0.0 };
        Vec3d vel { 2.0, 0.0, 0.0 }; // 2m/s along x-axis
        Vec3d acc { 0.0, 0.0, 0.0 }; // no acceleration at first
        double mass = 1.0; // 1kg
        double drag = 0.1; // rho*C*Area - simplified drag for this example
    
        /**
         * Update pos and vel using "Velocity Verlet" integration
         * @param dt DeltaTime / time step [eg: 0.01]
         */
        void update(double dt)
        {
            Vec3d new_pos = pos + vel*dt + acc*(dt*dt*0.5);
            Vec3d new_acc = apply_forces(); // only needed if acceleration is not constant
            Vec3d new_vel = vel + (acc+new_acc)*(dt*0.5);
            pos = new_pos;
            vel = new_vel;
            acc = new_acc;
        }
    
        Vec3d apply_forces() const
        {
            Vec3d grav_acc = Vec3d{0.0, 0.0, -9.81 }; // 9.81m/s^2 down in the Z-axis
            Vec3d drag_force = 0.5 * drag * (vel * abs(vel)); // D = 0.5 * (rho * C * Area * vel^2)
            Vec3d drag_acc = drag_force / mass; // a = F/m
            return grav_acc - drag_acc;
        }
    };
    

    也许会有所帮助)

    • -1

相关问题

  • 在 Haskell 中实现自动机算法

  • Haskell 中的工作表

  • 如何旋转 Haskell 文件的内容

Sidebar

Stats

  • 问题 10021
  • Answers 30001
  • 最佳答案 8000
  • 用户 6900
  • 常问
  • 回答
  • Marko Smith

    表格填充不起作用

    • 2 个回答
  • Marko Smith

    提示 50/50,有两个,其中一个是正确的

    • 1 个回答
  • Marko Smith

    在 PyQt5 中停止进程

    • 1 个回答
  • Marko Smith

    我的脚本不工作

    • 1 个回答
  • Marko Smith

    在文本文件中写入和读取列表

    • 2 个回答
  • Marko Smith

    如何像屏幕截图中那样并排排列这些块?

    • 1 个回答
  • Marko Smith

    确定文本文件中每一行的字符数

    • 2 个回答
  • Marko Smith

    将接口对象传递给 JAVA 构造函数

    • 1 个回答
  • Marko Smith

    正确更新数据库中的数据

    • 1 个回答
  • Marko Smith

    Python解析不是css

    • 1 个回答
  • Martin Hope
    Alexandr_TT 2020年新年大赛! 2020-12-20 18:20:21 +0000 UTC
  • Martin Hope
    Alexandr_TT 圣诞树动画 2020-12-23 00:38:08 +0000 UTC
  • Martin Hope
    Air 究竟是什么标识了网站访问者? 2020-11-03 15:49:20 +0000 UTC
  • Martin Hope
    Qwertiy 号码显示 9223372036854775807 2020-07-11 18:16:49 +0000 UTC
  • Martin Hope
    user216109 如何为黑客设下陷阱,或充分击退攻击? 2020-05-10 02:22:52 +0000 UTC
  • Martin Hope
    Qwertiy 并变成3个无穷大 2020-11-06 07:15:57 +0000 UTC
  • Martin Hope
    koks_rs 什么是样板代码? 2020-10-27 15:43:19 +0000 UTC
  • Martin Hope
    Sirop4ik 向 git 提交发布的正确方法是什么? 2020-10-05 00:02:00 +0000 UTC
  • Martin Hope
    faoxis 为什么在这么多示例中函数都称为 foo? 2020-08-15 04:42:49 +0000 UTC
  • Martin Hope
    Pavel Mayorov 如何从事件或回调函数中返回值?或者至少等他们完成。 2020-08-11 16:49:28 +0000 UTC

热门标签

javascript python java php c# c++ html android jquery mysql

Explore

  • 主页
  • 问题
    • 热门问题
    • 最新问题
  • 标签
  • 帮助

Footer

RError.com

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

帮助

© 2023 RError.com All Rights Reserve   沪ICP备12040472号-5