Articles of entrometido

Producto puntual de un vector en SciPy / NumPy (obteniendo ValueError: los objetos no están alineados)

Acabo de empezar a aprender SciPy y estoy luchando con las funciones más básicas. Considere el siguiente vector estándar: In [6]: W=array([[1],[2]]) In [7]: print W [[1] [2]] Si lo comprendo correctamente, esta debería ser la representación en SciPy de un vector matemático estándar de 2×1, como este: (1) (2) El producto puntual de este […]

¿Cómo leer / recorrer / cortar matrices dispersas Scipy (LIL, CSR, COO, DOK) más rápido?

Para manipular matrices Scipy, normalmente, se utilizan los métodos incorporados. Pero a veces es necesario leer los datos de la matriz para asignarlos a tipos de datos no dispersos. Por el bien de la demostración, creé una matriz dispersa aleatoria de LIL y la convertí a una matriz Numpy (¡los tipos de datos de python […]

ValueError: objeto demasiado profundo para la matriz deseada

“”” ___ “”” from scipy.optimize import root import numpy as np LENGTH = 3 def process(x): return x[0, 0] + x[0, 1] * 5 def draw(process, length): “”” “”” X = np.matrix(np.random.normal(0, 10, (length, 2))) y = np.matrix([process(x) for x in X]) y += np.random.normal(3, 1, len(y)) return yT, XT def maximum_likelyhood(y, X): def objective(b): […]

Etiqueta erosionante erosión

¿Cómo puedo mantener un anillo de píxeles alrededor de las regiones etiquetadas en una matriz numpy? En un caso simple, restaría la erosión. Ese enfoque no funciona cuando las tags se tocan. ¿Cómo puedo obtener B de A ? A = array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, […]

¿Cómo encontrar todos los ceros de una función usando numpy (y scipy)?

Supongamos que tengo una función f(x) definida entre b . Esta función puede tener muchos ceros, pero también muchas asíntotas. Necesito recuperar todos los ceros de esta función. ¿Cuál es la mejor manera de hacerlo? En realidad, mi estrategia es la siguiente: Evalúo mi función en un número dado de puntos Detecto si hay un […]

Muestreo de flotadores aleatorios en un rango en números

¿Cómo puedo muestrear flotantes aleatorios en un intervalo [a, b] en números? No solo enteros, sino cualquier número real. Por ejemplo, random_float(5, 10) devolvería números aleatorios entre [5, 10] . Gracias.

¿Cómo se implementa la función numpy.cov ()?

Tengo mi propia implementación de la función de covarianza basada en la ecuación: ”’ Calculate the covariance coefficient between two variables. ”’ import numpy as np X = np.array([171, 184, 210, 198, 166, 167]) Y = np.array([78, 77, 98, 110, 80, 69]) # Expected value function. def E(X, P): expectedValue = 0 for i in […]

Cómo incluir la restricción a la solución de la función NNLS de Scipy para que sume 1

Tengo el siguiente código para resolver el mínimo cuadrado no negativo. Utilizando scipy.nnls. import numpy as np from scipy.optimize import nnls A = np.array([[60, 90, 120], [30, 120, 90]]) b = np.array([67.5, 60]) x, rnorm = nnls(A,b) print x #[ 0. 0.17857143 0.42857143] # Now need to have this array sum to 1. Lo que […]

¿Por qué NumPy y SciPy tienen muchas de las mismas funciones? ¿Cuál debería preferir?

Posible duplicado: Relación entre scipy y numpy Por ejemplo, NumPy tiene las funciones de ventana bartlett , blackman , hamming , hanning , kaiser , mientras que SciPy tiene estas y varias más , pero parecen producir resultados idénticos. NumPy tiene numpy.fft.fft2(a, s=None, axes=(-2, -1)) . SciPy tiene scipy.fftpack.fft2(x, shape=None, axes=(-2, -1), overwrite_x=0) . ¿Por […]

Cómo obtener los valores de una matriz NumPy usando múltiples índices

Tengo una matriz NumPy que se ve así: arr = [100.10, 200.42, 4.14, 89.00, 34.55, 1.12] ¿Cómo puedo obtener múltiples valores de esta matriz por índice? Por ejemplo, ¿cómo puedo obtener los valores en las posiciones de índice 1, 4 y 5? Estaba intentando algo como esto, lo cual es incorrecto: arr[1, 4, 5]