我在評(píng)估 Sympy 光學(xué)模塊的 TWave 實(shí)例時(shí)遇到問(wèn)題。起初，我嘗試lambdify()使用此代碼：from sympy import symbolsfrom sympy.physics.optics import TWaveimport numpy as npA, f, phi = symbols('A f phi')wave = TWave(A, f, phi)# the evaluationxpoint = tpoint = np.linspace(0, 10, 100)fb = lambdify([A, f, phi, x, t, n], wave.rewrite('cos'))fb(1, 1000, 0.5*np.pi, xpoint, tpoint, 1)但是，我收到了這個(gè)錯(cuò)誤：AttributeError: 'Float' object has no attribute 'cos'相反。然后我嘗試了一個(gè)簡(jiǎn)單evalf()的一點(diǎn)，只是為了確定。# evaluation with evalf()data = { 'A0': 1, 'f' : 100, 'phi' : np.pi/2,}ff = wave.subs(data).rewrite('cos')ff# result: TWave instance, rewritten in cosine term, with variable substituted現(xiàn)在，后來(lái)發(fā)生的事情讓我感到奇怪：我一眼就看到需要替換的變量有 3 個(gè)，即x, t, n. 但是，當(dāng)我嘗試這樣做時(shí)：ff.subs({'x':1, 't':10, 'n':1})我得到了ValueError: summation of quantities of incompatible dimensions錯(cuò)誤。我試圖n從 subs 字典中刪除：ff.subs({'x':1, 't':10})它起作用了，但n仍然是一個(gè)象征性的例子。因此，我無(wú)法從中獲得任何數(shù)值ff.subs({'x':1, 't':10}).evalf()任何想法？我只想能夠評(píng)估任何給定點(diǎn)的橫波值。謝謝，PS：可能不相關(guān)，但我在 Anaconda 上使用 Sympy 1.2
1 回答
一只甜甜圈
TA貢獻(xiàn)1836條經(jīng)驗(yàn) 獲得超5個(gè)贊
看起來(lái)問(wèn)題與TWave在不告訴您的情況下向表達(dá)式添加單位有關(guān)。
wave.rewrite('cos')給A*cos(-2*pi*f*t + phi + pi*second*f*n*x/(149896229*meter)). 如果您消除meter和second符號(hào)，您的代碼就可以工作。
from sympy import symbols
from sympy.physics.optics import TWave
from sympy.physics.units import meter, second
import numpy as np
A, f, phi = symbols('A f phi')
wave = TWave(A, f, phi)
# the evaluation
xpoint = tpoint = np.linspace(0, 10, 100)
fb = sympy.lambdify(
[A, f, phi, x, t, n],
wave.rewrite(
'cos'
).subs(
[(meter, 1), (second, 1)]
),
modules='numpy'
)
fb(1, 1000, 0.5*np.pi, xpoint, tpoint, 1)
添加回答
舉報(bào)
0/150
提交
取消