Python Sympy计算梯度、散度和旋度的实例
更新时间：2019年12月06日 17:15:08 作者：落叶_小唱
今天小编就为大家分享一篇Python Sympy计算梯度、散度和旋度的实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
sympy有个vector 模块，里面提供了求解标量场、向量场的梯度、散度、旋度等计算,官方参考连接：
http://docs.sympy.org/latest/modules/vector/index.html
sympy中计算梯度、散度和旋度主要有两种方式：
一个是使用∇∇算子，sympy提供了类Del()，该类的方法有：cross、dot和gradient，cross就是叉乘，计算旋度的，dot是点乘，用于计算散度，gradient自然就是计算梯度的。
另一种方法就是直接调用相关的API：curl、divergence和gradient，这些函数都在模块sympy.vector 下面。
使用sympy计算梯度、散度和旋度之前，首先要确定坐标系，sympy.vector模块里提供了构建坐标系的类，常见的是笛卡尔坐标系， CoordSys3D，根据下面的例子可以了解到相应应用。
（1）计算梯度
## 1 gradient
C = CoordSys3D('C')
delop = Del() # nabla算子
# 标量场 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）计算散度
## divergence
C = CoordSys3D('C')
delop = Del() # nabla算子
# 向量场 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，向量场的散度是标量
（3）计算旋度
## curl
C = CoordSys3D('C')
delop = Del() # nabla算子
# 向量场 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，向量场的旋度是向量
以上这篇Python Sympy计算梯度、散度和旋度的实例就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持脚本之家。
相关文章
下面小编就为大家带来一篇spyder常用快捷键(分享)。小编觉得挺不错的，现在就分享给大家，也给大家做个参考。一起跟随小编过来看看吧2017-07-07
本文主要介绍了基于Python socket实现简易网络聊天室，本文将通过pyqt5作为桌面应用框架，socket作为网络编程的框架，从而实现包括客户端和服务端的网络聊天室的GUI应用，需要的可以参考一下2022-07-07
这篇文章主要为大家详细介绍了python实现word 2007文档转换为pdf文件，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2018-03-03
今天教大家如何利用Python爬虫批量下载喜马拉雅音频,文中有非常详细的代码示例,对正在学习python的小伙伴们很有帮助,需要的朋友可以参考下2021-05-05
在本篇内容里小编给大家分享的是关于python如何导入库的步骤和方法，有需要的朋友们可以学习下。2020-06-06
本文主要介绍了Python读取串口数据的实现方法,可以使用pySerial库来读取串口数据,具有一定的参考价值,感兴趣的可以了解一下2024-02-02
pyinstaller打包python3.6和PyQt5中各种错误的解决方案汇总
pyinstaller是打包python很方便的一个套件,我们可以很轻易地使用他,下面这篇文章主要给大家介绍了关于pyinstaller打包python3.6和PyQt5中各种错误解决的相关资料,需要的朋友可以参考下2022-08-08
这篇文章主要介绍了python使用Queue在多个子进程间交换数据的方法,实例分析了Queue实现进程间数据交互的技巧,需要的朋友可以参考下2015-04-04
本文是面向SEO人群的Python程序语言入门教程，也适用于其他没有程序基础但想学习些程序，以解决简单的实际应用需求的人群。在后面会尽量用最基础的角度来介绍这门语言2012-07-07
这篇文章主要介绍了Python多进程Process和管道Pipe的使用方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教2024-02-02
最新评论