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

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:

import sympy

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

from sympy import init_printing
init_printing(use_latex=True)

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

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

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

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

稍微整理,可知熵变为0

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

原文发布于微信公众号 - 传输过程数值模拟学习笔记(SongSimStudio)

原文发表时间:2019-08-08

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

扫码关注云+社区

领取腾讯云代金券