Lösen einer Gleichung mit scipy ist fsolve

Ich versuche, lösen Sie die Gleichung f(x) = x-sin(x) -n*t -m0

In dieser Gleichung n und m0 sind Attribute, die definiert werden in meiner Klasse. Weiter t ist eine Konstante ganze Zahl in der Gleichung, aber es hat jedes mal ändert.

Hab ich gelöst, die Gleichung so bekomme ich eine "neue Gleichung". Ich habe importiert scipy.optimieren

def f(x, self):
    return (x - math.sin(x) -self.M0 - self.n*t)

def test(self,t):
    return fsolve(self.f, 1, args=(t))

Korrekturen und Vorschläge zu machen, damit es funktioniert?

  • Können Sie uns ein vollständiges Beispiel, auch das, was über dieses "funktioniert nicht"?
  • Ich glaube nicht, dass ich je gesehen habe f(x, self) vor.. das erste argument erhalten die Instanz, unabhängig davon, was Sie es nennen.
  • es ist mit fsolve aus scipy.optimieren
  • das Hauptproblem ist, dass es al ist, getan zu werden, in eine Klasse. das ist, warum ich brauchen, um zu definieren die Gleichung in f(x, self)
  • Der erste Fehler, die ich brauche, um loszuwerden, ist TypeError: func() takes exactly 2 arguments (3) also hat es wahrscheinlich etwas damit zu tun, args=(t), weil vor dem habe ich Hinzugefügt, dass es gar nicht geben, dass Fehler
InformationsquelleAutor Saybia | 2012-05-18
  1. 4

    Können, sehe ich mindestens zwei Probleme: das haben Sie verwechselt die Reihenfolge der Argumente zu f, und du bist nicht zu geben f Zugang zu t. So etwas wie dies funktionieren sollte:

    import math
from scipy.optimize import fsolve

class Fred(object):
    M0 = 5.0
    n = 5

    def f(self, x, t):
        return (x - math.sin(x) -self.M0 - self.n*t)

    def test(self, t):
        return fsolve(self.f, 1, args=(t))

    [Hinweis: ich war faul und machte M0 und n Klasse Mitglieder]

    gibt:

    >>> fred = Fred()
>>> fred.test(10)
array([ 54.25204733])
>>> import numpy
>>> [fred.f(x, 10) for x in numpy.linspace(54, 55, 10)]
[-0.44121095114838482, -0.24158955381855662, -0.049951288133726734,
 0.13271070588400136, 0.30551399241764443, 0.46769772292130796, 
 0.61863201965219616, 0.75782574394219182, 0.88493255340251409, 
 0.99975517335862207]
    • Es hat funktioniert !!! Danke Mann, ich wirklich zu schätzen, dass!
    InformationsquelleAutor DSM
  2. 1

    Müssen Sie definieren f() etwa so:

      def f(self, x, t):
    return (x - math.sin(x) - self.M0 - self.n * t)

    In anderen Worten:

    1. self kommt zuerst (immer);
    2. dann kommt der aktuelle Wert der x;
    3. dann kommen die Argumente, die Sie liefern, zu fsolve().
    InformationsquelleAutor NPE
  3. 0

    Bist du mit einem root-finding-Algorithmus in irgendeiner Form. Es gibt mehrere, die die gemeinsame Nutzung, so wäre es hilfreich zu wissen, was man.

    Müssen Sie wissen, drei Dinge:

    1. Der Algorithmus, die Sie verwenden
    2. Die Gleichung, die Sie sind Suche nach den Wurzeln für
    3. Die erste Vermutung und der Bereich, in dem Sie suchen

    Müssen Sie wissen, dass einige Kombinationen können keine Wurzeln.

    Visualisierung der Funktionen, die von Interesse sein können, hilfreich. Sie haben zwei: eine lineare Funktion und eine sinusiod. Wenn Sie waren, um eine Darstellung der beiden, die Sätze von Konstanten geben würde, Sie Schnittmengen? Der Schnittpunkt ist die Wurzel, die Sie suchen.

    • das Hauptproblem ist, dass es al ist, getan zu werden, in eine Klasse. das ist, warum ich brauchen, um zu definieren die Gleichung in f(x selbst).
    InformationsquelleAutor duffymo
Schreibe einen Kommentar