避免绘制Odes intgenet解ODEint

问题描述

我试图绘制一个相平面，我希望它看起来不错，但是由于初始条件，方程组的某些解决方案有所不同。有什么办法可以使try / except链顺序排列，以便在解决方案分歧时不会绘制它。这是我的代码：

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

def aux_func(x):  
  y = x[0]-x[1]
  if (np.abs(y) <= 1):
    f = y**3 + 0.5*y 
  else: 
    f = 2*y - np.sign(y)
  return f 

def function(x,t): 
  x1_dot = x[1]  
  x2_dot = -x[1] - aux_func(x)  
  return [x1_dot,x2_dot]

ts = np.linspace(0,20,300)  
ic_1 = np.linspace(-1,1,10)
ic_2 = np.linspace(-1,10) 

for r1 in ic_1:
  for r2 in ic_2:
    x0 = (r1,r2)
    try:
      xs = odeint(function,x0,ts)
      plt.plot(xs[:,0],xs[:,1],"r-",linewidth=.8) 
    except: 
      pass

# Nombre de los ejes,limites,plt.xlabel("$x_1$",fontsize=12)
plt.ylabel("$x_2$",fontsize=12)
# plt.tick_params(labelsize=10)
# plt.xticks(np.linspace(0,11))
# plt.yticks(np.linspace(0,11))
plt.xlim(-1,1)
plt.ylim(-1,1)
# Grafica el campo vectorial
X1,X2 = np.mgrid[-1:1:20j,-1:1:20j]
u=X2
d= X1-X2

t = np.zeros(np.shape(d))
for i in range(len(d)):
    for j in range(len(d[0])):
        if np.abs(d[i][j]) > 1:
            t[i][j]= 2*d[i][j]-0.5*np.sign(d[i][j])
        else:
            t[i][j] =d[i][j]**3 + 0.5*d[i][j]
v=-X2-t

pl.quiver(X1,X2,u,v,color = 'b',width = .002)
plt.grid()
plt.title('Plano de Fase Punto 1')
#plt.savefig('FasePunto4.png')
plt.show()

该代码正在绘制以下内容：

enter image description here

感谢帮助。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

numpy ode python