TypeError: неверный тип операнда для abs(): 'ReferenceFrame', но у меня нет функции abs()
Я пытаюсь вычислить завиток вектора с помощью его оператора del, поэтому я использую завиток из sympy и просто пересекаю вектор с системой отсчета, но получаю эту ошибку: что именно я делаю неправильно?
import sympy
import numpy as np
import math
from sympy import Symbol, diff, Array, sin, cos, curl
from sympy import init_printing
from sympy.physics.vector import curl, ReferenceFrame
init_printing()
# Variables being used - B0 is the initial mag field, alpha is a constant. x/y/z are for the direction
alpha = Symbol('\u03B1')
B0 = Symbol('B0')
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
print ('Symbols: ', alpha, B0, x, y, z)
# This has the reference frame and the vector
R = ReferenceFrame('R')
V = B0*sin(alpha*x)*R.y + B0*cos(alpha*x)*R.z
print(V)
C = curl(V,R)
print(C)
Я ожидаю, что оператор /nabla пересечет вектор.
Что такое B0, альфа, R? СИМВОЛЫ? Покажите, как вы инициализируете их. Дайте нам некоторые, которые работают с копированием и вставкой.
Я отредактировал свой вопрос с константой и для чего они используются
Согласно документам и коду, curl расширяет расчет, с dot и так далее. Он работает для меня и печатает 0.
Я буквально должен вернуть alpha* вектор как завиток. Но я не понимаю, как я получаю ошибку, которая у меня есть
Ответы 1
Используйте R[0] вместо x. Я сделал вывод из примеров curl (и пытался воспроизвести некоторые примеры из Вики)
Забегаю isympy
(Невозможно импортировать curl из sympy)
In [2]: from sympy.physics.vector import curl, ReferenceFrame
...: alpha = Symbol('\u03B1')
...: B0 = Symbol('B0')
...: R = ReferenceFrame('R')
...: V = B0*sin(alpha*R[0])*R.y + B0*cos(alpha*R[0])*R.z
...: print(V)
...: C = curl(V,R)
...: print(C)
...:
...:
B0*sin(R_x*α)*R.y + B0*cos(R_x*α)*R.z
B0*α*sin(R_x*α)*R.y + B0*α*cos(R_x*α)*R.z
Может быть способ разложить на множители alpha или иным образом показать, что C есть alpha*V. Я просто работаю над учебниками sympy.
In [4]: expand(C-alpha*V)
Out[4]: 0
In [8]: V
Out[8]: B₀⋅sin(Rₓ⋅α) r_y + B₀⋅cos(Rₓ⋅α) r_z
In [9]: C
Out[9]: B₀⋅α⋅sin(Rₓ⋅α) r_y + B₀⋅α⋅cos(Rₓ⋅α) r_z
Показать полную трассировку (как правильно отформатированный текст в вопросе).