gpt4 book ai didi

python - 求解和绘制具有相关参数的 ODE

转载 作者:太空宇宙 更新时间:2023-11-04 04:43:47 25 4
gpt4 key购买 nike

我正在编写一些代码,其中我有以下方程组 here .问题是我非常想解决多个 k 值以及为每个 k 值绘制相平面/颤动图。有人可以帮帮我吗?到目前为止,这是我对求解器的了解:

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

def model(X, t):
x = X[0]
y = X[1]
dxdt = k*x - y
dydt = x + y
return [dxdt, dydt]
#Initial state
X0 = [1,1]

#Time
t = np.linspace(0,10)

X = odeint(model, X0, t)

这就是我到目前为止的绘图:

x = X[:,0]
y = X[:,1]

plt.plot(x,y)

请注意,我并不是要简单地解决系统问题!我正在尝试通过更改多个值 (k) 并绘制生成的方程式来解决它。

最佳答案

通过更多的计算工作,您可以使用

xvalues, yvalues = np.meshgrid(np.arange(-2, 2, 0.02), np.arange(-2, 2, 0.02))
Ks = np.linspace(-1, 1, 12)

plt.subplots_adjust(hspace=0.4,wspace=0.4)
for j,kix in enumerate(Ks):
xdot, ydot = model([xvalues,yvalues],0 , kix)
plt.subplot(3,4,j+1)
plt.streamplot(xvalues, yvalues, xdot, ydot, density=0.5)
plt.title("k=%.2f"%kix)

plt.show()

得到

enter image description here

关于python - 求解和绘制具有相关参数的 ODE,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/50071763/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com