溢出错误:(34,“结果太大”)关于使用欧拉方法的微分方程的分辨率

我试图解决一个物理问题的微分方程:高尔夫球的坠落。我正在使用欧拉解析方法,我有这个代码:


def F(Z, t):

    res = [Z[1], (-K/m)*Z[1]*sqrt((Z[1]**2) + (Z[1]**2)), Z[3], (-K/m)*Z[3]*sqrt((Z[3]**2) + (Z[3]**2))]

    return res


def reso_z_euler(liste_t, F, CI, K):

    n = len(liste_t)

    Z = CI

    pas = (liste_t[-1] - liste_t[0])/(n-1)

    listeZ = [Z]

    for i in range(1, n):

        Z = F(Z, temps[i-1])

        listeZ.append(Z)

    return listeZ


listeZ = reso_z_euler(temps, F, Z0, K)

当我运行此代码时,终端抛出一个错误:OverflowError:(34,“结果太大”)我从未见过此错误,我不知道该怎么办。你的帮助会很好..谢谢


慕丝7291255
浏览 103回答 1
1回答

波斯汪

您正在使用未知或非本地数组作为时间,而不是tempsliste_t你没有实现欧拉方法,你必须使用某种机制来实现向量算术,例如Z = Z + pas*F(Z,liste_t[i-1])Z = [ zk+pas*fk for zk, fk in zip(Z, F(Z,liste_t[i-1])) ]您忘记在导数计算中包含重力。一般来说,欧拉方法只适合学习数值ODE积分的原理,对于任何有用的结果,请使用高阶方法。
打开App,查看更多内容
随时随地看视频慕课网APP

相关分类

Python