【文章內(nèi)容簡介】 y 2 1/2 31 數(shù)學(xué)表達(dá)式 除了使用 SymPy中預(yù)先定義好的具有特殊運(yùn)算含義的數(shù)學(xué)函數(shù)之外，還可以使用Function()創(chuàng)建自定義的數(shù)學(xué)函數(shù)： 請注意 Function雖然是一個類，但是上面的語句所得到的 f并不是 Function類的實(shí)例。和預(yù)定義的數(shù)學(xué)函數(shù)一樣， f是一個類，它從Function類繼承： f = Function(f) isinstance(f, Function) False 32 數(shù)學(xué)表達(dá)式 當(dāng)我使用 f創(chuàng)建一個表達(dá)式時，就相當(dāng)于創(chuàng)建它的一個實(shí)例： f的實(shí)例 t可以參與表達(dá)式運(yùn)算： t = f(x,y) isinstance(t, Function) True type(t) f （ 其中 func和 args是 Basic類的兩個非常重要的屬性，分別表示對象的類和對象的參數(shù) ） f (x, y) t+t*t f(x, y)**2 + f(x, y) 33 符號運(yùn)算 ?表達(dá)式變換 和化簡 simplify()可以對數(shù)學(xué)表達(dá)式進(jìn)行化簡，例如 : simplify()調(diào)用 SymPy內(nèi)部的多種表達(dá)式變換函數(shù)對表達(dá)式進(jìn)行化簡運(yùn)算。但是數(shù)學(xué)表達(dá)式的化簡是 一件非常復(fù)雜的工作，并且對于同一個表達(dá)式，根據(jù)其使用目的可以有多種化簡方案。 simplify((x+2)**2 (x+1)**2) 2*x + 3 34 符號運(yùn)算 radsimp()對表達(dá)式的分母進(jìn)行有理化，它所得到的 表達(dá)式的分母部分將不含無理數(shù)。例如： 它也可以對帶符號的表達(dá)式進(jìn)行處理 : radsimp(1/(sqrt(5)+2*sqrt(2))) (sqrt(5) + 2*sqrt(2))/3 radsimp(1/(y*sqrt(x)+x*sqrt(y))) (sqrt(x)*y + x*sqrt(y))/(x*y*(x y)) 35 符號運(yùn)算 ratsimp()對表達(dá)式中的分母進(jìn)行通分運(yùn)算，即將表達(dá)式轉(zhuǎn)換為分 子除分母的形式： fraction()返回一個包含表達(dá)式的分子和分母的元組 ,用它可以獲得 ratsimp()通分之后的分 子或分母： 注意 fraction()不會自動對表達(dá)式進(jìn)行通分運(yùn)算，因此： ratsimp(x/(x+y)+y/(xy)) 2*y**2/(x**2 y**2) + 1 fraction(ratsimp(1/x+1/y)) (x + y, x*y) fraction(1/x+1/y) (1/y + 1/x, 1) 36 符號運(yùn)算 cancel()對分式表達(dá)式的分子分母進(jìn)行約分運(yùn)算，可以對純符號的分式表達(dá)式以及自定義函數(shù)表達(dá)式進(jìn)行約分，但是不能對內(nèi)部函數(shù)的表達(dá)式進(jìn)行約分。 cancel((x**21)/(1+x)) x1 cancel(sin((x**21)/(1+x))) cancel不能對函數(shù)內(nèi)部的表達(dá)式進(jìn)行約分 sin(x**2/(x + 1) 1/(x + 1)) cancel((f(x)**21)/(f(x)+1)) 能對自定義函數(shù)表達(dá)式進(jìn)行約分 f(x) 1 37 符號運(yùn)算 trigsimp()對表達(dá)式中的三角函數(shù)進(jìn)行化簡。它有兩個可選參數(shù) deep和 recursive，默認(rèn) 值都為 False。當(dāng) deep參數(shù)為 True時，將對表達(dá)式中的所有子表達(dá)式進(jìn)行簡化運(yùn)算；當(dāng) recursive 參數(shù)為 True時，將遞歸使用 trigsimp()進(jìn)行最大限度的化簡： trigsimp(sin(x)**2+2*sin(x)*cos(x)+cos(x)**2) sin(2*x) + 1 trigsimp(f(sin(x)**2+2*sin(x)*cos(x)+cos(x)**2)) 也能對自定義函數(shù)中的三角函數(shù)化簡，至今不知道 deep和 recursive是干嘛的 f(sin(2*x) + 1) 38 符號運(yùn)算 expand_trig()可以對三角函數(shù)的表達(dá)式進(jìn)行展開。它實(shí)際上是對 expand()的封裝，通過將 expand()的 trig參數(shù)設(shè)置為 True,實(shí)現(xiàn)三角函數(shù)的展開 計算。輸入 “expand_trig??” 來査看它調(diào)用 expand()時的參數(shù)。 expand()通用的展開運(yùn)算，根據(jù)用戶設(shè)置的標(biāo)志參數(shù)對表達(dá)式進(jìn)行展幵。默認(rèn)情況下，以下的標(biāo)志參數(shù)為 True。 mul： 展開乘法 expand_trig(sin(2*x+y)) (2*cos(x)**2 1)*sin(y) + 2*sin(x)*cos(x)*cos(y) 39 符號運(yùn)算 log:展開對數(shù)函數(shù)參數(shù)中的乘積和冪運(yùn)算 multinomial:展開加法式的整數(shù)次冪 power_base:展開冪函數(shù)的底數(shù)乘積 x,y=symbols(x,y,positive=True) expand(log(x*y**2)) log(x) + 2*log(y) expand((x+y)**3) x**3 + 3*x**2*y + 3*x*y**2 + y**3 expand(x**(y+z)) x**y*x**z 40 符號運(yùn)算 可以將默認(rèn)為 True的標(biāo)志參數(shù)設(shè)置為False,強(qiáng)制不展 開對應(yīng)的表達(dá)式。在下面的例子中， 將 mul設(shè)置為 False,因此不對乘法進(jìn)行展開： expand()的以下標(biāo)志參數(shù)默認(rèn)為 False。 plex:展開復(fù)數(shù)的實(shí)部和虛部，默認(rèn)不展開復(fù)數(shù)的實(shí)部和虛部 : x,y,z=symbols(x,y,z, positive=True) expand(x*log(y*z), mul=False) x*(log(y) + log(z)) x,y=symbols(x,y,plex=True) expand(x*y, plex=True) re(x)*re(y) + I*re(x)*im(y) + I*re(y)*im(x) im(x)*im(y) 41 符號運(yùn)算 func:對一些特殊函數(shù)進(jìn)行展 開 trig:展開三角函數(shù) expand_log()、 expand mul()、expand_plex()、 expand_trig()、expand_func()等函數(shù)則通過將相應(yīng)的標(biāo)志參數(shù)設(shè)置為 True,對 expand()進(jìn)行封裝。 expand (gamma (1+x),func=True) x*gamma(x) expand(sin(x+y), trig=True) sin(x)*cos(y) + sin(y)*cos(x) 42 符號運(yùn)算 factor()可以對多項(xiàng)式表達(dá)式進(jìn)行因式分解： collect()收集表達(dá)式中指定符號的有理指數(shù)次冪的系數(shù)。例如，希望獲得如下表達(dá)式中 x的各次冪的系數(shù)： factor(15*x**2+2*y3*x10*x*y) (3*x 2*y)*(5*x 1) factor(expand((x+y)**20)) (x + y)**20 a,b=symbols(39。a,b39。) eq = (1+a*x)**3 + (1+b*x)**2 43 符號運(yùn)算 首先需要對表達(dá)式 eq進(jìn)行展開，得到的表達(dá)式 eq2是一系列乘式的和 : 然后調(diào)用 collect(),對表達(dá)式 eq2中 X的冪的系數(shù)進(jìn)行收集 ： eq2 = expand (eq) eq2 a**3*x**3 + 3*a**2*x**2 + 3*a*x + b**2*x**2 + 2*b*x + 2 collect(eq2,x) a**3*x**3 + x**2*(3*a**2 + b**2) + x*(3*a + 2*b) + 2 44