首页 > 编程语言> 文章详细

【Python】SymPy符号计算库学习

2020-11-26 20:32:20 阅读：325 来源： 互联网

标签：tmp 函数 Python 积分 sqrt SymPy integrate 符号计算表达式

SymPy符号整理

定义变量（符号）：symbols
定义函数：Function

SymPy函数整理

积分与泰勒展开

表达式展开：expand()

expand(,complex=True):表达式分为实数、虚数两部分

泰勒展开：series(函数表达式,自变量,0,余项次数)
不定积分运算：integrate(表达式，自变量)

定积分运算：integrate（表达式，（自变量，积分下界，积分上届））

算式中x换成y：subs(x,y)

IN：

import numpy as np
from sympy import *
#将x定义为符号
x=Symbol("x",real=True)
#创建多个符号：x,y,r=symbols('x,y,r')
#参数：positive=True，表示符号为正
#var():快速创建变量和Symbol对象
#泰勒展开
tmp=series(exp(I*x),x,0,10)
tmp
#获得tmp实部
re(tmp)
#获得tmp虚部
im(tmp)

OUT:

1 + I*x - x**2/2 - I*x**3/6 + x**4/24 + I*x**5/120 - x**6/720 - I*x**7/5040 + x**8/40320 + I*x**9/362880 + O(x**10)
x**8/40320 - x**6/720 + x**4/24 - x**2/2 + re(O(x**10)) + 1

表达式化简

化简数学表达式：simplify()
分母有理化：radsimp()
通分：fraction()
约分：trim()
三角函数化简：trigsimp()

deep=True：对所有子表达式化简

recursive=True：递归进行最大程度化简

乘法展开：mul()
整数次幂展开：multinomial()
log展开：expand()
合并同类项：collect()

IN:

var("a,b,x")
#展开eq得到eq2
eq=(1+a*x)**3+(1+b*x)**2
eq2=expand(eq)
eq2
#合并同类项
collect(eq2,x)
#得到x的各次幂系数，如获得x的二次项系数
p=collect(eq2,x,evaluate=False)
p[x**2]

OUT:

a**3*x**3 + 3*a**2*x**2 + 3*a*x + b**2*x**2 + 2*b*x + 2
a**3*x**3 + x**2*(3*a**2 + b**2) + x*(3*a + 2*b) + 2
3*a**2 + b**2

解方程

解方程：solve()

IN:

var("a,b,c")
#解一元二次方程
solve(a*x**2+b*x+c,x)
#解二元二次方程，结果每个元组表示方程的一组解
solve((x**2+x*y+1,y**2+x*y+2),x,y)

OUT:

[(-b + sqrt(-4*a*c + b**2))/(2*a), -(b + sqrt(-4*a*c + b**2))/(2*a)]
[(-sqrt(3)*I/3, -2*sqrt(3)*I/3), (sqrt(3)*I/3, 2*sqrt(3)*I/3)]

微分

得到导函数：t=Derivative(f(x),x)
计算出导函数：t.doit()
高阶导函数：Derivative(f(x),x,3)
多个变量导函数：Derivative(f(x,y),x,2,y,3)
微分方程符号求解：dsolve((Derivative(f(x),x)-f(x),f(x))

hint=best，返回最简单的解，得到最简单的显函数表达式

积分

不定积分运算：integrate(表达式，自变量)
定积分运算：integrate（表达式，（自变量，积分下界，积分上届））
二重不定积分运算：integrate（表达式，x,y）
二重定积分运算：integrate（表达式，(x,a,b),(y,c,d)）
计算出积分函数：doit()
无法表示为初等函数的积分——数值计算：evalf()

（不够精确，且不适合计算无穷积分）

其他

一些基础语句补充：

E:自然底数
I:虚数单位
pi:π
**:^
exp():e^x

参考资料

《Python科学计算》

标签：tmp,函数,Python,积分,sqrt,SymPy,integrate,符号计算,表达式
来源： https://www.cnblogs.com/fighterkaka22/p/14044367.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9