Tag: scipy

Erhalten eines Fehlers "konnte das Eingabefeld aus der Form (252,4) nicht in Form (4) in der Optimierung übertragen

Ich eine relativ neue python scipy bibliothek Ich habe versucht, die scipy.optimize verwenden, um den maximalen Wert der sharpe () – Funktion in dem folgenden Code zu finden def sharpe(dr, wts): portfolio=np.ones(dr.shape[0]) dr[0:]=wts*dr[0:] sharpe_ratio=-np.mean(np.sum(dr, axis=1))/np.std(np.sum(dr, axis=1)) return sharpe_ratio def wts_con(wts): return wts[0]+wts[1]+wts[2]+wts[3]-1 def sharpe_optimize(dr, sharpe_func): wts_guess=np.array([0.25,0.25,0.25,0.25]) con=[{"type":"eq", "fun":wts_con}] bnds=((0,1),(0,1),(0,1),(0,1)) result=spo.minimize(sharpe_func, wts_guess, args=(dr,), method="SLSQP", constraints=con, bounds=bnds, […]

Numpy / Scipy Broadcast Berechnung Skalar Produkt für ein bestimmtes Element

Ich habe eine spärliche Matrix wie A Und ein dataframe (df) mit Zeilen, die genommen werden sollten, um Skalarprodukt zu berechnen. Row1 Row2 Value 2 147 scalar product of vectors at Row1 and Raw2 in matrix A Kann ich es im Rundfunk ohne Looping machen? In meinem Fall A wie 1m * 100k Größe und […]

Interpolation basierend auf einem Array-Wert

Ich habe zwei Arrays mit Werten: x = np.array([100, 123, 123, 118, 123]) y = np.array([12, 1, 14, 13]) Ich möchte zB die Funktion auswerten: def func(a, b): return a*0.8 * (b/2) Also, ich möchte die fehlenden Werte füllen. Ich benutze: import numpy as np from scipy import interpolate def func(a, b): return a*0.8 * […]

Extrahieren Sie benötigte Bytes aus einer Datei in Python

Ich habe hier eine Binärdatei: ftp://n5eil01u.ecs.nsidc.org/SAN/GLAS/GLA06.034/2003.02.21/GLA06_634_1102_001_0079_3_01_0001.DAT Ich muss folgende Daten aus dieser Datei extrahieren: Byte Offset: 176 Data type: 4-byte (long) integer Total bytes: 160 Ich habe wie folgt versucht: import numpy as np fname = 'GLA06_634_1102_001_0079_3_01_0001.DAT' with open(fname,'rb') as fi: fi.seek (176,0) data= np.fromfile(fi,dtype='long',count=160) print data Kein Erfolg, was ist los mit meiner Idee?

Wie man Coroutine IPython <-> einen Rückruf ()

In einem IPython-Terminal möchte ich eine Funktion tief in main() zurück zu IPython, wo ich drucken kann, setze … wie üblich, dann lauf weiter main() : IPython run main.py … def callback( *args ): … try: back_to_ipython() # <– how to do this ? In[]: print, set *args … … except KeyboardInterrupt: # or IPython […]

Python: Berechnung der Kosinus Ähnlichkeit von zwei Dicts schneller

Ich habe zwei dicts: d1 = {1234: 4, 125: 7, …} d2 = {1234: 8, 1288: 5, …} Die Längen der Dicts variieren von 10 bis 40000. Um die Cosinus Ähnlichkeit zu berechnen, benutze ich diese Funktion: from scipy.linalg import norm def simple_cosine_sim(a, b): if len(b) < len(a): a, b = b, a res = […]

Einschränkungen bei Kurvenanpassungsparametern

Ich versuche, in meinem quadratischen Kurvenanpassungsprozess Grenzen und Beschränkungen aufzuerlegen. Ziel ist es, Koeffizienten a,b und c . Imposing Einschränkung auf b: delta-2*a*x ist mein Zweifel. Wie kann ich eine Variable hinzufügen, x in meine Einschränkungen. Bearbeitbarer Code: from lmfit import Model, Parameters #create x and y data to be used for curve fitting xip=[ […]

Deblur ein Bild mit scikit-Bild

Ich versuche, skimage.restoration.wiener zu benutzen, aber ich komme immer mit einem Bild mit einem Haufen von 1 (oder -1), was mache ich falsch? Das Originalbild kommt von Uni von Waterloo . import numpy as np from scipy.misc import imread from skimage import color, data, restoration from scipy.signal import convolve2d as conv2 def main(): image = […]

Angabe von Einschränkungen für fmin_cobyla in scipy

Ich benutze Python 2.5. Ich gehe an die Cobyla-Optimierung vorbei: import numpy from numpy import asarray Initial = numpy.asarray [2, 4, 5, 3] # Initial values to start with #bounding limits (lower,upper) – for visualizing #bounds = [(1, 5000), (1, 6000), (2, 100000), (1, 50000)] # actual passed bounds b1 = lambda x: 5000 – […]

Scipy.integrate.trapz und diskontinuierliche Funktionen

Die Funktion scipy.integrate.trapz verwendet die Newton-Cotes-Formel von Ordnung 1, wie sie in der scipy Dokumentation steht. Bei der Ableitung dieser Formel wird jedoch in der Regel davon ausgegangen Der Integrand ist eine stetige Funktion und Die Punkte, in denen der Wert des Integrandes bekannt ist, sind verschieden. Allerdings habe ich versucht, das Integral der Funktion […]

Python ist die beste Programmiersprache der Welt.