2015-06-28 34 views
5

我無法找到一種方法讓SymPy將諸如cos(a)*cos(b)之類的產品擴展爲角度總和的trig函數的總和。Python3 - Sympy:展開trig函數的產品

from sympy import * 
init_printing() 
wrf,wlo,t = symbols('\omega_RF \omega_LO t') 
c = cos(wrf*t)*cos(wlo*t) 
expand_trig(c) 

保持產品完好無損。 simplify(c)trigsimp(c)也不提供任何替代形式。

我想要cos(a)*cos(b)擴展爲1/2*(cos(a+b) + cos(a-b)) ...任何提示?

回答

8

每文檔字符串,help(sympy.fu)

fu會盡量最小化目標函數measure。默認情況下,這個 首先最小化trig條件的數量,然後最小化總計操作的數量。

但是,如果你通過

measure=lambda x: -x.count_ops() 

然後fuwill try to maximize the op count


import sympy as sy 
sy.init_printing() 
wrf, wlo, t = sy.symbols('\omega_RF \omega_LO t') 
c = sy.cos(wrf*t)*sy.cos(wlo*t) 
print(sy.fu(c, measure=lambda x: -x.count_ops())) 

產生

cos(\omega_LO*t - \omega_RF*t)/2 + cos(\omega_LO*t + \omega_RF*t)/2 

或者,你可以call the Fu transformation TR8 directly

from sympy.simplify.fu import TR8 
print(TR8(c)) 

產生同樣的結果。文檔字符串help(sys.modules['sympy.simplify.fu'])解釋了可用的轉換。這是總結;檢查文檔字符串以獲取更多信息:

TR0 - simplify expression 
TR1 - sec-csc to cos-sin 
TR2 - tan-cot to sin-cos ratio 
TR2i - sin-cos ratio to tan 
TR3 - angle canonicalization 
TR4 - functions at special angles 
TR5 - powers of sin to powers of cos 
TR6 - powers of cos to powers of sin 
TR7 - reduce cos power (increase angle) 
TR8 - expand products of sin-cos to sums 
TR9 - contract sums of sin-cos to products 
TR10 - separate sin-cos arguments 
TR10i - collect sin-cos arguments 
TR11 - reduce double angles 
TR12 - separate tan arguments 
TR12i - collect tan arguments 
TR13 - expand product of tan-cot 
TRmorrie - prod(cos(x*2**i), (i, 0, k - 1)) -> sin(2**k*x)/(2**k*sin(x)) 
TR14 - factored powers of sin or cos to cos or sin power 
TR15 - negative powers of sin to cot power 
TR16 - negative powers of cos to tan power 
TR22 - tan-cot powers to negative powers of sec-csc functions 
TR111 - negative sin-cos-tan powers to csc-sec-cot