欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

Python Sympy計(jì)算梯度、散度和旋度的實(shí)例

 更新時(shí)間:2019年12月06日 17:15:08   作者:落葉_小唱  
今天小編就為大家分享一篇Python Sympy計(jì)算梯度、散度和旋度的實(shí)例,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧

sympy有個(gè)vector 模塊,里面提供了求解標(biāo)量場(chǎng)、向量場(chǎng)的梯度、散度、旋度等計(jì)算,官方參考連接:

http://docs.sympy.org/latest/modules/vector/index.html

sympy中計(jì)算梯度、散度和旋度主要有兩種方式:

一個(gè)是使用∇∇算子,sympy提供了類Del(),該類的方法有:cross、dot和gradient,cross就是叉乘,計(jì)算旋度的,dot是點(diǎn)乘,用于計(jì)算散度,gradient自然就是計(jì)算梯度的。

另一種方法就是直接調(diào)用相關(guān)的API:curl、divergence和gradient,這些函數(shù)都在模塊sympy.vector 下面。

使用sympy計(jì)算梯度、散度和旋度之前,首先要確定坐標(biāo)系,sympy.vector模塊里提供了構(gòu)建坐標(biāo)系的類,常見的是笛卡爾坐標(biāo)系, CoordSys3D,根據(jù)下面的例子可以了解到相應(yīng)應(yīng)用。

(1)計(jì)算梯度

## 1 gradient

C = CoordSys3D('C')
delop = Del() # nabla算子

# 標(biāo)量場(chǎng) f = x**2*y-xy
f = C.x**2*C.y - C.x*C.y

res = delop.gradient(f, doit=True) # 使用nabla算子
# res = delop(f).doit()
res = gradient(f) # 直接使用gradient

print(res) # (2*C.x*C.y - C.y)*C.i + (C.x**2 - C.x)*C.j

(2)計(jì)算散度

## divergence

C = CoordSys3D('C')
delop = Del() # nabla算子

# 向量場(chǎng) f = x**2*y*i-xy*j
f = C.x**2*C.y*C.i - C.x*C.y*C.j

res = delop.dot(f, doit=True)

# res = divergence(f)

print(res) # 2*C.x*C.y - C.x,即2xy-x,向量場(chǎng)的散度是標(biāo)量

(3)計(jì)算旋度

## curl

C = CoordSys3D('C')
delop = Del() # nabla算子

# 向量場(chǎng) f = x**2*y*i-xy*j
f = C.x**2*C.y*C.i - C.x*C.y*C.j

res = delop.cross(f, doit=True)

# res = curl(f)

print(res) # (-C.x**2 - C.y)*C.k,即(-x**2-y)*k,向量場(chǎng)的旋度是向量

以上這篇Python Sympy計(jì)算梯度、散度和旋度的實(shí)例就是小編分享給大家的全部?jī)?nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持腳本之家。

相關(guān)文章

  • Django對(duì)接支付寶實(shí)現(xiàn)支付寶充值金幣功能示例

    Django對(duì)接支付寶實(shí)現(xiàn)支付寶充值金幣功能示例

    今天小編就為大家分享一篇Django對(duì)接支付寶實(shí)現(xiàn)支付寶充值金幣功能示例,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧
    2019-12-12
  • Python實(shí)現(xiàn)switch/case語(yǔ)句

    Python實(shí)現(xiàn)switch/case語(yǔ)句

    與Java、C\C++等語(yǔ)言不同,Python中是不提供switch/case語(yǔ)句的,這一點(diǎn)讓我感覺到很奇怪。我們可以通過如下幾種方法來實(shí)現(xiàn)switch/case語(yǔ)句
    2021-08-08
  • pip install過程中出現(xiàn)error: subprocess-exited-with-error錯(cuò)誤的解決辦法

    pip install過程中出現(xiàn)error: subprocess-exited-with-erro

    在Python環(huán)境配置過程中,經(jīng)常會(huì)遇到各種安裝錯(cuò)誤,其中一個(gè)常見的錯(cuò)誤是 pip install 過程中出現(xiàn)的 error: subprocess-exited-with-error,在本文中,我將詳細(xì)解析這個(gè)錯(cuò)誤的原因,并提供具體的解決方案和步驟,確保大家能夠順利完成安裝,避免類似問題的再次發(fā)生
    2024-07-07
  • Pandas實(shí)現(xiàn)數(shù)據(jù)類型轉(zhuǎn)換的一些小技巧匯總

    Pandas實(shí)現(xiàn)數(shù)據(jù)類型轉(zhuǎn)換的一些小技巧匯總

    這篇文章主要給大家匯總介紹了關(guān)于Pandas實(shí)現(xiàn)數(shù)據(jù)類型轉(zhuǎn)換的一些小技巧,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧
    2018-05-05
  • 在Python下利用OpenCV來旋轉(zhuǎn)圖像的教程

    在Python下利用OpenCV來旋轉(zhuǎn)圖像的教程

    這篇文章主要介紹了在Python下利用OpenCV來旋轉(zhuǎn)圖像的教程,代碼和核心的算法都非常簡(jiǎn)單,需要的朋友可以參考下
    2015-04-04
  • python+matplotlib實(shí)現(xiàn)鼠標(biāo)移動(dòng)三角形高亮及索引顯示

    python+matplotlib實(shí)現(xiàn)鼠標(biāo)移動(dòng)三角形高亮及索引顯示

    這篇文章主要介紹了Python+matplotlib實(shí)現(xiàn)鼠標(biāo)移動(dòng)三角形高亮及索引顯示,具有一定借鑒價(jià)值,需要的朋友可以參考下
    2018-01-01
  • python矩陣列的實(shí)現(xiàn)示例

    python矩陣列的實(shí)現(xiàn)示例

    在Python和NumPy庫(kù)的幫助下,矩陣列可以很容易地進(jìn)行各種操作,本文主要介紹了python矩陣列的實(shí)現(xiàn)示例,具有一定的參考價(jià)值,感興趣的可以了解一下
    2024-02-02
  • 詳解Python中for循環(huán)是如何工作的

    詳解Python中for循環(huán)是如何工作的

    如果你對(duì)python中的for循環(huán)不是很清楚,那么建議你看看這篇文章,本文主要給大家介紹了關(guān)于Python中for循環(huán)是如何工作的相關(guān)資料,介紹的非常詳細(xì),對(duì)大家具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面來一起看看吧。
    2017-06-06
  • Python基礎(chǔ)之Spyder的使用

    Python基礎(chǔ)之Spyder的使用

    Spyder是一個(gè)用于科學(xué)計(jì)算的使用Python編程語(yǔ)言的集成開發(fā)環(huán)境(IDE),它結(jié)合了綜合開發(fā)工具的高級(jí)編輯、分析、調(diào)試等功能,需要的朋友可以參考下
    2023-05-05
  • Python利用Matplotlib庫(kù)實(shí)現(xiàn)繪制餅形圖

    Python利用Matplotlib庫(kù)實(shí)現(xiàn)繪制餅形圖

    這篇文章主要為大家分享了基于python+matplotlib庫(kù)的餅形圖繪制,具體內(nèi)容涉及一般的餅圖、分裂餅圖、以及環(huán)形圖,感興趣的小伙伴可以了解一下
    2022-04-04

最新評(píng)論