Articles of scipy

Normalizar scipy.ndimage.filters.correlate

¿Alguien tiene una idea de cómo normalizar el scipy.ndimage.filters.correlate función para obtener: XCM = 1/N(xc(a-mu_a,b-mu_b)/(sig_a*sig_b)) ¿Qué es N para la correlación? Por lo general, es el # de puntos de datos / píxeles para las imágenes. ¿Qué valor elegiré para scipy.ndimage.filters.correlate ? Mis imágenes difieren en tamaño. Supongo que la función de correlación scipy rellena […]

python: fsolve con unknown dentro del límite superior de una integral

¿Es posible aplicar el método fsolve para una integral con lo desconocido en el límite superior y en el integrando? Estoy usando el método cuádruple para la integración en Python. ¡¡¡Gracias por adelantado!!!

El diseño del filtro de paso alto de Python scipy.signal.remez produce una función de transferencia extraña

Estoy tratando de diseñar filtros de paso alto equivalentes utilizando la función scipy.signal.remez de python. Sin embargo, las funciones de transferencia resultantes me parecen muy extrañas con un máximo de ~ 15 db en la banda de paso y solo una atenuación de la banda de parada de 6 dB. El diseño de paso bajo […]

Directamente “trazar” segmentos de línea para numpy array

Uno de mis primeros proyectos realizados en python es la simulación de percolación de palo de Monte Carlo. El código creció continuamente. La primera parte fue la visualización de la percolación del palo. En un área de ancho * largo, una densidad definida (palos / área) de palos rectos con una cierta longitud se traza […]

Distancia de Hausdorff entre grillas 3D

Tengo varias cuadrículas (matrices numpy [Nk, Ny, Nx]) y me gustaría usar la distancia de Hausdorff como una métrica de similitud de estas cuadrículas. Hay varios módulos en scipy (scipy.spatial.distance.cdist, scipy.spatial.distance.pdist) que permiten calcular la distancia euclidiana entre matrices 2D. Ahora, para comparar las cuadrículas, tengo que elegir una sección transversal (por ejemplo, cuadrícula1 [0 […]

Cómo redondear adecuadamente las matrices de flotadores numpy

Tengo algún tipo de problema de redondeo al redondear flotantes. x = np.array([[1.234793487329877,2.37432987432],[1.348732847,8437.328737874]]) np.round(x,2) array([[ 1.23000000e+00, 2.37000000e+00], [ 1.35000000e+00, 8.43733000e+03]]) ¿Hay una manera de mostrar estos números sin las extensiones cero?

Error al intentar instalar sklearn desde Pycharm | arrayobject.h no puede ser absoluto

Declaración de error completa: ValueError: la ruta ‘/home/andy/anaconda3/lib/python3.5/sitepackages/numpy/core/include/numpy/arrayobject.h’ no puede ser absoluta He instalado la versión Scipy y numpy (mkl) descargando los archivos de la rueda comstackdos desde este enlace y luego instalándolos desde el cmd usando pip . Soy más de un Python Noob, mi enfoque es principalmente el aprendizaje automático. Por favor ayúdame.

Resolviendo un sistema de ecuaciones trascendentales con python.

Suponiendo que tengo las siguientes cuatro ecuaciones: cos (x) / x = a cos (y) / y = b a + b = 1 c sinc (x) = d sinc (y) para variables desconocidas x, y, a b . Tenga en cuenta que cos(x)/x=a tiene múltiples soluciones. Lo mismo ocurre con la variable y . […]

scipy ode update set_f_params dentro de la función establecida como set_solout

Al integrar una oda con scipy, ode acepta una función con más argumentos que t y y. Por ejemplo: def fun(t, y, param1, param2): y el valor de estos argumentos se puede establecer utilizando el método set_f_params . Sin embargo, al usar también el método set_solout e intentar actualizar los parámetros con set_f_params dentro de […]

¿Resolver un sistema de odas (con cambio de constante) usando scipy.integrate.odeint?

Actualmente tengo un sistema de odas con una constante dependiente del tiempo. P.ej def fun(u, t, a, b, c): x = u[0] y = u[1] z = u[2] dx_dt = a * x + y * z dy_dt = b * (yz) dz_dt = -x*y+c*yz return [dx_dt, dy_dt, dz_dt] Las constantes son “a”, “b” y […]