脚本之家 服务器常用软件
快捷导航
您的位置：首页 → 脚本专栏 → python → Python Sympy计算梯度、散度和旋度

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计算梯度、散度和旋度的实例就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持脚本之家。

您可能感兴趣的文章:

相关文章

  • 浅谈django rest jwt vue 跨域问题

    浅谈django rest jwt vue 跨域问题

    这篇文章主要介绍了django rest jwt vue 跨域问题，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧
    2018-10-10
  • Python中变量的作用域的具体使用

    Python中变量的作用域的具体使用

    本文主要介绍了Python中变量的作用域的具体使用，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
    2023-03-03
  • python实现监控指定进程的cpu和内存使用率

    python实现监控指定进程的cpu和内存使用率

    这篇文章主要为大家详细介绍了python实现监控指定进程的cpu和内存使用率，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
    2022-01-01
  • Python matplotlib实现多重图的绘制

    Python matplotlib实现多重图的绘制

    Matplotlib作为Python的2D绘图库，它以各种硬拷贝格式和跨平台的交互式环境生成出版质量级别的图形。本文将利用Matplotlib库绘制多重图，感兴趣的可以了解一下
    2022-03-03
  • Django将默认的SQLite更换为MySQL的实现

    Django将默认的SQLite更换为MySQL的实现

    今天小编就为大家分享一篇Django将默认的SQLite更换为MySQL的实现，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
    2019-11-11
  • pyramid配置session的方法教程

    pyramid配置session的方法教程

    这篇文章主要介绍了pyramid如何配置session，大家可以参考使用
    2013-11-11
  • python之用Numpy和matplotlib画一个魔方

    python之用Numpy和matplotlib画一个魔方

    这篇文章主要介绍了如何用Numpy和matplotlib画一个魔方，本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
    2021-08-08
  • Python根据当前日期取去年同星期日期

    Python根据当前日期取去年同星期日期

    最近做项目，遇到这样的业务开发需求，需要对比当前时间段和去年同星期的时间段的数据，下面小编通过实例代码给大家分享Python根据当前日期取去年同星期日期，需要的朋友参考下
    2019-04-04
  • win10下Python3.6安装、配置以及pip安装包教程

    win10下Python3.6安装、配置以及pip安装包教程

    下面小编就为大家带来一篇win10下Python3.6安装、配置以及pip安装包教程。小编觉得挺不错的，现在就分享给大家，也给大家做个参考。一起跟随小编过来看看吧
    2017-10-10
  • Windows下将Python文件打包成.EXE可执行文件的方法

    Windows下将Python文件打包成.EXE可执行文件的方法

    这篇文章主要介绍了Windows下将Python文件打包成.EXE可执行文件的方法,需要的朋友可以参考下
    2018-08-08

最新评论