Python+SymPy實(shí)現(xiàn)秒解微積分詳解
之前我們分享過(guò)很多有用有趣的Python庫(kù)，今天繼續(xù)介紹一個(gè)：
SymPy 是一個(gè)Python庫(kù)，專注于符號(hào)數(shù)學(xué)，它的目標(biāo)是成為一個(gè)全功能的計(jì)算機(jī)代數(shù)系統(tǒng)，同時(shí)保持代碼簡(jiǎn)潔、易于理解和擴(kuò)展。
舉一個(gè)簡(jiǎn)單的例子，比如說(shuō)展開(kāi)二次方程：
from?sympy?import?*
x = Symbol('x')
y = Symbol('y')
d = ((x+y)**2).expand()
print(d)
# 結(jié)果：x**2 + 2*x*y + y**2你可以隨便輸入表達(dá)式，即便是十次方，它都能輕易的展開(kāi)，非常方便：
from?sympy?import?*
x = Symbol('x')
y = Symbol('y')
d = ((x+y)**10).expand()
print(d)
# 結(jié)果：x**10 + 10*x**9*y + 45*x**8*y**2 + 120*x**7*y**3 + 210*x**6*y**4 + 252*x**5*y**5 + 210*x**4*y**6 + 120*x**3*y**7 + 45*x**2*y**8 + 10*x*y**9 + y**10下面就來(lái)講講這個(gè)模塊的具體使用方法和例子。
1.準(zhǔn)備
請(qǐng)選擇以下任一種方式輸入命令安裝依賴：
1. Windows 環(huán)境 打開(kāi) Cmd (開(kāi)始-運(yùn)行-CMD)。
2. MacOS 環(huán)境 打開(kāi) Terminal (command+空格輸入Terminal)。
3. 如果你用的是 VSCode編輯器 或 Pycharm，可以直接使用界面下方的Terminal.
pip?install?Sympy
2.基本使用
簡(jiǎn)化表達(dá)式(化簡(jiǎn))
sympy支持三種化簡(jiǎn)方式，分別是普通化簡(jiǎn)、三角化簡(jiǎn)、指數(shù)化簡(jiǎn)。
普通化簡(jiǎn) simplify( )：
from?sympy?import?*
x = Symbol('x')
d = simplify((x**3?+ x**2?- x -?1)/(x**2?+?2*x +?1))
print(d)
# 結(jié)果：x - 1三角化簡(jiǎn) trigsimp( )：
from?sympy?import?*
x = Symbol('x')
d = trigsimp(sin(x)/cos(x))
print(d)
# 結(jié)果：tan(x)指數(shù)化簡(jiǎn) powsimp( )：
from?sympy?import?*
x = Symbol('x')
a = Symbol('a')
b = Symbol('b')
d = powsimp(x**a*x**b)
print(d)
# 結(jié)果：x**(a + b)解方程 solve()
第一個(gè)參數(shù)為要解的方程，要求右端等于0，第二個(gè)參數(shù)為要解的未知數(shù)。
如一元一次方程：
from?sympy?import?*
x = Symbol('x')
d = solve(x *?3?-?6, x)
print(d)
# 結(jié)果：[2]二元一次方程：
from?sympy?import?*
x = Symbol('x')
y = Symbol('y')
d = solve([2?* x - y -?3,?3?* x + y -?7],[x, y])
print(d)
# 結(jié)果：{x: 2, y: 1}求極限 limit()
dir=’+’表示求解右極限，dir=’-‘表示求解左極限：
from?sympy?import?*
x = Symbol('x')
d = limit(1/x,x,oo,dir='+')
print(d)
# 結(jié)果：0
d = limit(1/x,x,oo,dir='-')
print(d)
# 結(jié)果：0求積分 integrate( )
先試試求解不定積分：
from?sympy?import?*
x = Symbol('x')
d = integrate(sin(x),x)
print(d)
# 結(jié)果：-cos(x)再試試定積分：
from?sympy?import?*
x = Symbol('x')
d = integrate(sin(x),(x,0,pi/2))
print(d)
# 結(jié)果：1求導(dǎo) diff()
使用 diff 函數(shù)可以對(duì)方程進(jìn)行求導(dǎo)：
from?sympy?import?*
x = Symbol('x')
d = diff(x**3,x)
print(d)
# 結(jié)果：3*x**2
d = diff(x**3,x,2)
print(d)
# 結(jié)果：6*x解微分方程 dsolve( )
以 y′=2xy 為例:
from?sympy?import?*
x = Symbol('x')
f = Function('f')
d = dsolve(diff(f(x),x) -?2*f(x)*x,f(x))
print(d)
# 結(jié)果：Eq(f(x), C1*exp(x**2))3.實(shí)戰(zhàn)一下
今天群里有同學(xué)問(wèn)了這個(gè)問(wèn)題，“大佬們，我想問(wèn)問(wèn)，如果這個(gè)積分用Python應(yīng)該怎么寫(xiě)呢，謝謝大家”：
# Python 實(shí)用寶典
from?sympy?import?*
x = Symbol('x')
y = Symbol('y')
d = integrate(x-y, (y,?0,?1))
print(d)
# 結(jié)果：x - 1/2為了計(jì)算這個(gè)結(jié)果，integrate的第一個(gè)參數(shù)是公式，第二個(gè)參數(shù)是積分變量及積分范圍下標(biāo)和上標(biāo)。
運(yùn)行后得到的結(jié)果便是 x - 1/2 與預(yù)期一致。
如果大家也有求解微積分、復(fù)雜方程的需要，可以試試sympy，它幾乎是完美的存在。
到此這篇關(guān)于Python+SymPy實(shí)現(xiàn)秒解微積分詳解的文章就介紹到這了,更多相關(guān)Python SymPy微積分內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
利用Python如何實(shí)現(xiàn)數(shù)據(jù)驅(qū)動(dòng)的接口自動(dòng)化測(cè)試
這篇文章主要給大家介紹了關(guān)于利用Python如何實(shí)現(xiàn)數(shù)據(jù)驅(qū)動(dòng)的接口自動(dòng)化測(cè)試的相關(guān)資料，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家學(xué)習(xí)或者使用python具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們來(lái)一起看看吧2018-05-05
Python基礎(chǔ)語(yǔ)法(Python基礎(chǔ)知識(shí)點(diǎn))
這篇文章主要介紹了Python基礎(chǔ)語(yǔ)法(Python基礎(chǔ)知識(shí)點(diǎn)),需要的朋友可以參考下2016-02-02
這篇文章主要分享下在python中解析xml文件的模塊用法，以及對(duì)模塊封裝的方法，有需要的朋友參考下2014-02-02
Python實(shí)現(xiàn)將Sheet頁(yè)拆分成單獨(dú)的Excel文件
這篇文章主要為大家詳細(xì)介紹了如何使用 Python 將一個(gè) Excel 文件中的每個(gè)工作表（Sheet）保存成單獨(dú)的 Excel 文件,有需要的小伙伴可以了解下2025-02-02
當(dāng)前使用的是anaconda的3.8版本，無(wú)法正常下載tensorflow包，需要構(gòu)建虛擬環(huán)境使用3.7及以下的解釋器才可以，如何解決這個(gè)問(wèn)題呢，下面小編給大家?guī)?lái)了如何在pycharm中使用tensorflow，感興趣的朋友參考下吧2021-11-11
Python中執(zhí)行MySQL結(jié)果限制和分頁(yè)查詢示例詳解
這篇文章主要為大家介紹了Python中執(zhí)行MySQL結(jié)果限制和分頁(yè)查詢示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升職加薪2023-11-11
在本篇文章里小編給大家整理了一篇關(guān)于python unittest單元測(cè)試的步驟，對(duì)此有興趣的朋友們可以跟著學(xué)習(xí)下。2021-08-08
這篇文章主要介紹了Python 忽略文件名編碼的方法，文中講解非常細(xì)致，代碼幫助大家更好的理解和學(xué)習(xí)，感興趣的朋友可以了解下2020-08-08
使用Python識(shí)別和處理驗(yàn)證碼的代碼示例
驗(yàn)證碼作為一種常見(jiàn)的安全手段,廣泛應(yīng)用于各種網(wǎng)站和應(yīng)用中,以防止自動(dòng)化腳本的惡意攻擊,然而,在自動(dòng)化測(cè)試或數(shù)據(jù)抓取過(guò)程中,識(shí)別驗(yàn)證碼成為了一個(gè)不得不面對(duì)的問(wèn)題,本文將詳細(xì)介紹如何使用Python來(lái)識(shí)別和處理驗(yàn)證碼,通過(guò)實(shí)際案例和代碼,幫助讀者理解整個(gè)流程2025-01-01
python中k-means和k-means++原理及實(shí)現(xiàn)
本文主要介紹了python中k-means和k-means++原理及實(shí)現(xiàn)，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2022-05-05