python-2.7 - scipy odeint を使用した時間の後方統合

Question

scipy.integrate.odeint を使用して、任意の常微分方程式を時間的に後方に統合することは可能ですか? 可能であれば、「odeint.

Answer 1

odeint引数の負の値を処理しtます。特別な治療は必要ありません。

次に例を示します。

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt


def mysys(z, t):
    """A slightly damped oscillator."""
    return [z[1] - 0.02*z[0], -z[0]]


if __name__ == "__main__":
    # Note that t starts at 0 and goes "backwards"
    t = np.linspace(0, -50, 501)

    z0 = [1, 1]
    sol = odeint(mysys, z0, t)

    plt.plot(t, sol)
    plt.xlabel('t')
    plt.show()

プロット：

Answer 2

変数 を変更してs = t_0 - t、 について微分方程式を積分できsます。odeint はこれを行いません。