Python-sympy计算问题入门-求极限

遇到复杂计算找python绝对不让你失望，Sympy是一个Python的科学计算库，用一套强大的符号计算体系完成诸如求极限、求导数、求微分、求积分等等计算问题。

虽然Matlab的类似科学计算能力也很强大，但是Python以其语法简单、易上手、异常丰富的三方库生态，个人认为可以更优雅地解决日常遇到的各种计算问题

正是因为有了丰富的 Python 库, Python 才发展得如此迅速👍。接下来我们将学习如何利用Sympy库求数列及函数极限。<br>

本片文章的目标:

1. 能够掌握Sympy库的安装与导入
2. 能够利用Sympy求数列极限与函数极限

01
sympy的安装与导入

1使用pip安装sympy

pip install sympy

2导入sympy库

from sympy import * #导入SymPy库

02
limit函数的使用

🎈 1. limit函数介绍

对于数列的计算，可以使用SymPy库中的limit函数实现，其语法格式如下。

✨ sympy.limit(e, z, z0, dir = ‘+’)

参数	说明
e	接收需要进行求极限的数列或函数表达式。无默认值。
z	接收需要进行求极限的数列的项或函数的自变量。无默认值。
z0	接收任何表达式，包括所有类型的数值、oo和-oo等，表示自变量趋于有限值或趋于无穷大，其中oo表示无穷大。无默认值
dir	接收+或-。取+时，表示求趋于有限值的右极限（z->z0)；取-时，表示求趋于有限值的左极限（z->zo-)。对于无穷大的zo（oo或-oo），dir参数无效。有默认值，默认为“+”。

🎈 2. 表达式符号写法

符号	写法
加号	+
减号	-
乘号	*
除号	/
指数	**
对数	log()
e的指数次幂	exp()
三角函数	sin(), cos(), tan(), cot()
反三角函数	asin(), acos(), atan(), acot()
Π	pi

sin(pi/2),cos(pi/2),tan(pi),cot(pi/2)
>>>(1, 0, 0, 0)

03
利用sympy求极限

✔ 例1：
数列1/2, 2/3, 3/4, ... n/(n+1),...的一般项为 *xn = n/(n+1)*，判断数列 {xn} 是否收敛。

n = Symbol('n') #将n定义为符号变量
xn = n/(n+1) #数列一般项

print('数列的极限为：',limit(xn,n,oo)) #数列的自变量都是趋于无穷大
>>>数列的极限为： 1

✔ 例2：
数列2, 4, 8, ... 2^n,...的一般项为 xn = 2^n，判断数列 {xn} 是否收敛。

n = Symbol('n') #将n定义为符号变量
xn = 2 ** n    # “**”符号表示乘方

print('数列的极限为：',limit(xn,n,oo))
>>>数列的极限为： oo

✔ 例3：
当x -> -1/2 时，计算函数f(x)=(1-4x^2)/(2x+1)的极限。

x = Symbol('x') #将x定义为自变量
y = (1-4*x**2)/(2*x+1)

print('函数的极限为：',limit(y,x,-1/2))
>>>函数的极限为： 2.00000000000000

✔ 例4：
当x -> 1 时，计算函数f(x)=(x^2-2x+1)/(x^2-1)的极限。

x = Symbol('x') #将x定义为自变量
y = (x**2-2*x+1)/(x**2-1)

print('函数的极限为：',limit(y,x,1))
>>>函数的极限为： 0

✔ 例5：
当x -> 0 时，计算函数f(x)= (tanx-sinx)/x的极限。

x = Symbol('x') #将x定义为自变量
y = (tan(x)-sin(x))/x

print('函数的极限为：',limit(y,x,0))
>>>函数的极限为： 0

04
总结

以上就是今天要讲的内容，本次主要介绍了Sympy计算库，利用该库进行求极限。最近该学一下数学了，数学对我们来说是走到哪都应该学的，我们应该从基础的开始学，了解极限的定义以及运算法则，进行一步一步推导。

对于Sympy能更方便的计算出结果，用于我们进行参考，我们学习数学，就应该懂得如何计算得出的。好啦，一起加油吧！

点个关注不迷路
觉得孔哥写的对你有帮助？请分享给更多的人
欢迎一起学习！博客平台同步发布，请搜索——和孔哥一起学

版权声明：

作者： 和孔哥一起学

导师：Qi Mengyao

本文版权归作者导师共有，欢迎转载，但未经作者同意必须在文章页面注明来源及原作者或原文链接，否则保留追究法律责任的权利。