用matplotlib绘制三常数系统的相图

import numpy as np
import matplotlib.pyplot as plt

def system(X, t):
    x, y, z = X
    dx_dt = a * (y - x)
    dy_dt = x * (b - z) - y
    dz_dt = x * y - c * z
    return [dx_dt, dy_dt, dz_dt]

a = 1.0
b = 2.0
c = 0.5
t = np.linspace(0, 10, 1000)  # 时间变量的范围和步长
x0 = 0.1
y0 = 0.2
z0 = 0.3

from scipy.integrate import odeint

X0 = [x0, y0, z0]
X = odeint(system, X0, t)

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot(X[:, 0], X[:, 1], X[:, 2])
ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Z')
plt.show()