TypeError: Неверный тип операнда для abs (): 'ReferenceFrame', но у меня нет функции abs () - PullRequest
Новая
       59

TypeError: Неверный тип операнда для abs (): 'ReferenceFrame', но у меня нет функции abs ()

0 голосов
/ 21 апреля 2019

Я пытаюсь вычислить скручиваемость вектора с его оператором 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)

Я ожидаю, что оператор / набла пересечет вектор.

1 Ответ

0 голосов
/ 22 апреля 2019

Используйте 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
...