在笛卡尔坐标系中绘制受简单圆形轨道-双曲吸引子影响的流场- Python

import numpy as np
import matplotlib.pyplot as plt

def flow_field(x, y):
    # 定义圆形轨道的半径和圆心坐标
    radius = 1.0
    center_x = 0.0
    center_y = 0.0
    
    # 计算流场的速度分量
    vx = -y / np.sqrt(x**2 + y**2) + (x - center_x) / np.sqrt((x - center_x)**2 + (y - center_y)**2)
    vy = x / np.sqrt(x**2 + y**2) + (y - center_y) / np.sqrt((x - center_x)**2 + (y - center_y)**2)
    
    return vx, vy

def plot_flow_field():
    # 定义绘图区域的范围
    x = np.linspace(-3, 3, 100)
    y = np.linspace(-3, 3, 100)
    X, Y = np.meshgrid(x, y)
    
    # 计算流场的速度分量
    U, V = flow_field(X, Y)
    
    # 绘制流场
    plt.figure(figsize=(8, 8))
    plt.quiver(X, Y, U, V)
    plt.xlabel('x')
    plt.ylabel('y')
    plt.title('Flow Field')
    plt.show()

plot_flow_field()