Odeint Parameter mit Grenzwert

mit matplotlib, NumPy, pandas, SciPy, SymPy und weiteren mathematischen Programmbibliotheken.
Pumpkinpy
User
Beiträge: 10
Registriert: Freitag 26. Mai 2017, 11:45

Odeint Parameter mit Grenzwert

Beitragvon Pumpkinpy » Donnerstag 15. Juni 2017, 18:33

Hallo liebes Forum,

ich beschäftige micht zur Zeit mit DGL's und deren Lösung mit odeint. Und zwar stehe ich vor dem Problem, dass ich ein DGL System habe. Die Werte nach denen Entwickelt wird haben allerdings die Grenze 0. Ich hoffe man konnte verstehen was ich ausdrücken wollte.
Habe dazu mal etwas Code das mein Problem vielleicht besser veranschaulicht. Die abfallende Gerade die gezeigt wird darf dabei minimal den Wert 0 annehmen, habe aber keine Ahnung wie ich das Realisieren soll. Danke für eure Mithilfe.

  1. import matplotlib.pyplot as plt
  2. import scipy.integrate as sci
  3. import numpy as np
  4.  
  5. def con(c, t):
  6.        
  7.     adt = b_1
  8.     bdt = -adt
  9.  
  10.     return( (adt, bdt))
  11.  
  12. b_1 = -0.5
  13. b_2 = 0.5
  14. c0 = [1, 0]
  15.  
  16. t = np.linspace(0 , 10, 1000)
  17.  
  18. sol = sci.odeint(con, c0 , t)
  19.  
  20. plt.plot(t, sol[:,0], 'r-' )
  21. plt.plot(t, sol[:,1], 'b-' )
  22. plt.show()
Zuletzt geändert von Anonymous am Donnerstag 15. Juni 2017, 20:12, insgesamt 1-mal geändert.
Grund: Quelltext in Python-Codebox-Tags gesetzt.
Benutzeravatar
Sr4l
User
Beiträge: 1091
Registriert: Donnerstag 28. Dezember 2006, 20:02
Wohnort: Kassel
Kontaktdaten:

Re: Odeint Parameter mit Grenzwert

Beitragvon Sr4l » Freitag 16. Juni 2017, 20:47

  1. def con(c, t):
  2.     adt = b_1
  3.     bdt = -adt
  4.    
  5.     if (c[0]) < 0:
  6.         adt = 0
  7.  
  8.     return( (adt, bdt))


Aber dein Beispiel ist etwas komisch dafür braucht man keine speziellen Verfahren.
Pumpkinpy
User
Beiträge: 10
Registriert: Freitag 26. Mai 2017, 11:45

Re: Odeint Parameter mit Grenzwert

Beitragvon Pumpkinpy » Montag 19. Juni 2017, 15:24

Danke dir.
An so eine einfache Lösung hab ich gar nicht gedacht :?
Naja manchmal sieht man den Wald vor Bäumen nicht.

Wer ist online?

Mitglieder in diesem Forum: 0 Mitglieder