溫馨提示×

SciPy中怎么進(jìn)行動態(tài)系統(tǒng)分析

小億
82
2024-05-24 11:51:13
欄目: 編程語言

在SciPy中進(jìn)行動態(tài)系統(tǒng)分析可以使用scipy.integrate模塊中的odeint函數(shù)來求解微分方程組。以下是一個(gè)簡單的示例：

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

# 定義微分方程組
def system(y, t):
    dydt = np.zeros_like(y)
    dydt[0] = -0.1*y[0] + 0.2*y[1]
    dydt[1] = 0.1*y[0] - 0.2*y[1]
    return dydt

# 初始條件
y0 = [1.0, 0.0]

# 時(shí)間點(diǎn)
t = np.linspace(0, 10, 100)

# 求解微分方程組
sol = odeint(system, y0, t)

# 繪制結(jié)果
plt.plot(t, sol[:, 0], label='y1')
plt.plot(t, sol[:, 1], label='y2')
plt.legend()
plt.xlabel('Time')
plt.ylabel('Values')
plt.show()

在這個(gè)示例中，我們定義了一個(gè)簡單的微分方程組，然后使用odeint函數(shù)求解該微分方程組，并繪制了結(jié)果。您可以根據(jù)自己的動態(tài)系統(tǒng)模型來修改微分方程組的定義和初始條件，以實(shí)現(xiàn)自己的動態(tài)系統(tǒng)分析。

0 踩