python 如何使用Symy计算向量的范数

bfnvny8b  于 2023-02-02  发布在  Python
关注(0)|答案(4)|浏览(203)

我需要计算一个向量的范数使用symy和符号从symy。摘录的代码如下:

from sympy import *
x = Symbol('x')
sb = [2,1]
func = sympy.exp(-(sympy.sqrt((x.norm() + (x-sb).norm())**2 - sb**2)/(2)))
func_prime = func.diff(x)
ff = lambdify(x, func_prime, 'numpy')
f = -1*ff(np.array(r))

上面的实现给了我错误:AttributeError: 'Symbol' object has no attribute 'norm'.
如果我使用func = V0 * sympy.exp(-(sympy.sqrt((sympy.sqrt(sum(x**2)) + sympy.sqrt(sum((x-sb*e)**2)))**2 - sb**2)/(2*sig))),我会得到以下错误:TypeError: 'Pow' object is not iterable.
sympy不接受numpy. linalg. norm()。
在这种情况下,我怎样才能求出向量的范数呢?
谢谢

tez616oj

tez616oj1#

Sympy有一个vector模块可以使用。

from sympy.vector import CoordSys3D
C = CoordSys3D('C')
v = 3*C.i + 4*C.j + 5*C.k
v.dot(v)

Output: 50

您还可以将其用于符号:

from sympy import Symbol
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
v = x*C.i + y*C.j + z*C.k
sympy.sqrt(v.dot(v))

Output: sqrt(x**2 + y**2 + z**2)
tvmytwxo

tvmytwxo2#

你好:)我建议你如下:

import sympy as sp
import sympy.physics.vector as spv

M = spv.ReferenceFrame("M")

vector=V1*M.x+V2*M.y+V3*M.z

其中M.x,M.y,M.z是物理学中的酉向量i,j,k,我们知道范数是向量与自身的点积的平方根,所以

norm=sp.sqrt(spv.dot(vector, vector))
print(norm)

如果要以LaTeX格式打印结果

print(sp.latex(norm))

如果你想简化表达式,

print(norm.simplify())
2g32fytz

2g32fytz3#

x, y, z = symbols('x y z')
vec = Matrix([x, y, z])
vec_norm = sqrt(sum(sympy.matrices.dense.matrix_multiply_elementwise(vec, vec)))
muk1a3rh

muk1a3rh4#

这个怎么样?

from sympy import *
x = Symbol('x')
y = Symbol('y')
a = Matrix([x,y])
anorm = sqrt(a.T @ a)

相关问题