前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >[python从入门到放弃]手边没稿纸推公式sympy

[python从入门到放弃]手边没稿纸推公式sympy

作者头像
周星星9527
发布2019-08-14 16:30:17
1.1K0
发布2019-08-14 16:30:17
举报

sympy是一个非常好用的基于python的符号计算库,科技做微分、积分、极限等一系列高等数学运算。

SymPy is a Python library for symbolic mathematics. It aims to become a full-featured computer algebra system (CAS) while keeping the code as simple as possible in order to be comprehensible and easily extensible. SymPy is written entirely in Python. https://www.sympy.org/en/index.html

例如验证等熵过程的熵变为0。本文python环境安装了anaconda全家桶,并在jupyter notebook中运行,首先导入计算库sympy:

代码语言:javascript
复制
import sympy

为了让公式更好看,使其排版如同?????LATEX,我们如下设置:

代码语言:javascript
复制
from sympy import init_printing
init_printing(use_latex=True)

好了定义一些变量,比如熵变,压强等等:

代码语言:javascript
复制
Cp,Rg,T1,T2,P1,P2,kappa = sympy.symbols('Cp Rg T1 T2 P1 P2 kappa')
dS = Cp*sympy.ln(T2/T1)-Rg*sympy.ln(P2/P1)
dS

运行程序,结果如下图,可见公式排版杠杠的:

代码语言:javascript
复制
#又根据等熵过程有P2/P1为:
ratio=sympy.Pow(T2/T1,kappa/(kappa-1))
ratio
代码语言:javascript
复制
#Cp与kappa的关系为:
kappa*Rg/(kappa-1)

稍微整理,可知熵变为0

早上喝了一碗羊汤,要求开发票!老板问:"啥单位?"我说:"中华人民共和国内蒙古自治区巴彦淖尔市乌拉特右翼后旗乌兰哈达苏木那令格嘎查那令格村那令格二社村民委员会火山历史文化遗址申报联合国教育科学文化组织世界非物质文化遗产工作领导小组办公室主任乌力图乌力吉吉尔格拉.巴彦尔太!"老板说:"羊汤钱不要了!

本文参与 腾讯云自媒体分享计划,分享自微信公众号。
原始发表:2019-08-08,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 传输过程数值模拟学习笔记 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档