帮我用python写一个单个FHN模型的代码

我很乐意为你解答，请参考以下代码：import numpy as np# 参数
#w：振荡频率
#a：激励阈值
#b：反馈阈值
#dt：时间步长
#x：状态变量
#v：激励
#u：反馈w = 1.0
a = 0.2
b = 0.3
dt = 0.001x = np.array([[0.0],[0.0]])
v = np.zeros(x.shape)
u = np.zeros(x.shape)for t in range(10000):
v = v + dt*( -w*w*x - 2*a*v + u )
x = x + dt*v
u = u + dt*(b*x)
if t%1000 == 0:
print('t=',t,'x=',x[0],'v=',v[0])