Articles of sympy

Evaluando jacobianos en puntos específicos usando sympy

Estoy tratando de evaluar al jacobiano en (x, y) = (0,0) pero no puedo hacerlo. import sympy as sp from sympy import * import numpy as np x,y=sp.symbols(‘x,y’, real=True) J = Function(‘J’)(x,y) f1=-y f2=x – 3*y*(1-x**2) f1x=diff(f1,x) f1y=diff(f1,y) f2x=diff(f2,x) f2y=diff(f2,y) J=np.array([[f1x,f1y],[f2x,f2y]]) J1=J(0,0) print J1 El error correspondiente a —> 16 J1=J(0,0) es TypeError: ‘numpy.ndarray’ object […]

Simpy y raíces cuadradas de números complejos.

Cuando se utiliza solve para calcular las raíces de una ecuación cuadrática, SymPy devuelve expresiones que podrían simplificarse pero no puedo lograr que las simplifique. Un ejemplo mínimo se ve así: from sympy import * sqrt(-24-70*I) Aquí, SymPy simplemente devuelve sqrt(-24-70*I) mientras que Mathematica o Maple responderán con el equivalente de 5-7*I Soy consciente de […]

SymPy: ¿Cómo devolver una expresión en términos de otra (s) expresión (s)?

Soy bastante nuevo en SymPy y tengo lo que podría ser una pregunta básica. O simplemente podría estar malinterpretando cómo se debe usar SymPy. ¿Hay una manera de crear una expresión que no esté representada por átomos, sino por una combinación de otras expresiones? Ejemplo: >>> from sympy.physics.units import * >>> expr1 = m/s >>> […]

Cómo lambdificar una expresión SymPy que contiene la función erf para usar con NumPy

Me gustaría lambdificar una expresión simbólica que contiene la función erf con SymPy. Esto se puede hacer para los argumentos escalares de la siguiente manera: log_normal = 0.5 + 0.5 * sym.erf((sym.log(x) – mu) / sym.sqrt(2 * sigma**2)) F = sym.lambdify([x, mu, sigma], log_normal) F(1.0, 0.0, 1.0) Me gustaría vectorizar lo anterior. Normalmente haría lo […]

Sympy: resolviendo matrices en un campo finito

Para mi proyecto, necesito resolver para una matriz X dadas matrices Y y K. (XY = K) Los elementos de cada matriz deben ser números enteros módulo una prima aleatoria de 256 bits. Mi primer bash de resolver este problema usó la función mod_inv(n) . El problema con esto es que me estoy quedando sin […]

¿Está rota la impresión bonita de Sympy en el nuevo portátil de Jupyter?

Anteriormente he usado la impresión bonita de las matemáticas en el cuaderno de ipython. Después de actualizar a jupyter (también actualiza muchos otros paquetes relacionados con ipython), la impresión bastante ya no funciona como antes. Utilizo este código en la parte superior de mis cuadernos para configurarlo: import sympy as sp sp.init_printing() También he intentado […]

Añadir signos de multiplicación (*) entre coeficientes

Tengo un progtwig en el que un usuario ingresa una función, como sin(x)+1 . Estoy usando ast para tratar de determinar si la cadena es ‘segura’ al incluir en la lista blanca los componentes como se muestra en esta respuesta . Ahora me gustaría analizar la cadena para agregar signos de multiplicación ( * ) […]

Diagtwig de coordenadas esféricas en matplotlib

R (teta, phi) = cos (phi ^ 2), teta [0, 2 * pi], phi [0, pi] ¿Cómo dibujar una gráfica de esta función (R (teta, phi)) en coordenadas esféricas con la ayuda de matplotlib? La documentación no he encontrado coordenadas esféricas.

Convierta una fórmula LaTex a un tipo que pueda usarse dentro de SymPy

Quiero analizar las fórmulas de LaTeX y usarlas directamente como expresiones de SymPy. En otras palabras, lo que necesito es algo similar a sympify: from sympy import sympify f = sympify(‘x^2 + sin(y) + 1/2’) print f pero eso puede tomar expresiones LaTeX (cadenas) como entrada, por ejemplo: f = latex_sympify(‘\frac{x}{1+x}’) Dado que sympify es […]

¿Cómo realizar una optimización no lineal con scipy / numpy o sympy?

Estoy tratando de encontrar la solución óptima para el siguiente sistema de ecuaciones en Python: (x-x1)^2 + (y-y1)^2 – r1^2 = 0 (x-x2)^2 + (y-y2)^2 – r2^2 = 0 (x-x3)^2 + (y-y3)^2 – r3^2 = 0 Dados los valores, un punto (x, y) y un radio (r): x1, y1, r1 = (0, 0, 0.88) x2, […]