python特殊函数_SciPyTutorial-特殊函数
3. Scipy Tutorial- 特殊函数
在scipy里scipy的special子包和numpy的polynomial提供了一些数学上的著名的公式的生成函数,例如scipy.special里的log1p、exp10等函数;numpy.polynomial里的Chebyshev、Laguerre等模块。
3.1 scipy.special里的函数
在scipy.special里提供了很多数学公式的函数,例如:
求$10^x$
from scipy.special import exp10
a = exp10(-1)
print a
再比如赫尔维茨$ \zeta $函数公式:
scipy.special.zeta(x, q)
$$
\zeta(x, q) = \sum_{n=0}^{\infty} \frac{1} { (n+q)^x}\quad (x > 1 , q > 0)
$$
可以使用scipy.special.zeta来获得。
from scipy.special import zeta
c = zeta(2, 1)
print "c", c
黎曼ζ函数可以使用scipy.special.zetac函数来实现。
scipy.special.zetac(x)
$$
\zeta(x) = \sum_{k = 2}^{\infty}\frac{1}{k^x}
$$
from scipy.special import zeta, zetac
c = zeta(2, 1)
print "c", c
d = zetac(2)
print "d", d
程序执行结果:
c 1.6449340668482266
d 0.6449340668482264
3.2 numpy.polynomial里的函数
polynomial子包里也有很多的知名数学公式的实现函数,例如切比雪夫多项式。
import numpy.polynomial.chebyshev as chebyshev
import numpy as np
import numpy.linalg as linalg
x = np.array([1, 2, 3, 4])
y = np.array([1, 3, 5, 4])
deg = len(x) - 1
A = chebyshev.chebvander(x, deg)
print A, "# A"
c = linalg.solve(A, y)
print c,"# c"
执行结果:
[[ 1. 1. 1. 1.]
[ 1. 2. 7. 26.]
[ 1. 3. 17. 99.]
[ 1. 4. 31. 244.]] # A
[ 3.5 -3.875 1.5 -0.125] # c
3.3 Numpy.poly1d里的函数
poly1d可以构造一元多项式、求导、求积分以及求根等。
import numpy as np
p = np.poly1d([1, 2, 1])
print p
print np.polyder(p)
print np.polyint(p)
print np.roots([1, 2, 1])
print np.polyval([1, 2, 1], 3)
执行结果:
2
1 x + 2 x + 1
2 x + 2
3 2
0.3333 x + 1 x + 1 x
[-1. -1.]
16
python特殊函数_SciPyTutorial-特殊函数相关推荐
- Day10 Python基础之特殊函数(八)
一些特殊函数 1.递归函数(recursion) 递归函数的定义:在函数内部,可以调用其他函数.如果一个函数在内部调用自身本身,这个函数就是递归函数. 递归函数的优点:是定义简单,逻辑清晰.理论上,所 ...
- python中的特殊函数__call__
1 # _*_ coding:utf-8 _*_ 2 3 ''' 4 一个类实例可以变成一个可以调用的对象,只需要实现特殊方法__call__ 5 所以可以看出在python中函数也是对象 6 ''' ...
- python保存内容到文件(text、json、csv)
python保存内容到文件(text.json.csv) 在开发人员的日常中,将数据保存到文件是最常见的编程任务之一. 通常,程序需要一些输入并产生一些输出.在许多情况下,我们希望将这些结果持久化.我 ...
- python花式编码_Python——花式打印对象的若干种方法
今天是Python专题的第10篇文章,我们来聊聊Python当中的类. 打印实例 我们先从类和对象当中最简单的打印输出开始讲起,打印一个实例是一个非常不起眼的应用,但是在实际的编程当中却非常重要.原因 ...
- python格式化输出 xz_Python——详解__str__, __repr__和__format__
今天是Python专题的第10篇文章,我们来聊聊Python当中的类. 打印实例 我们先从类和对象当中最简单的打印输出开始讲起,打印一个实例是一个非常不起眼的应用,但是在实际的编程当中却非常重要.原因 ...
- 一篇文章带你了解Python运算符重载
回复"python"即可获赠从入门到进阶共10本电子书 今 日 鸡 汤 不堪玄鬓影,来对白头吟. 您可以根据所使用的操作数来更改Python中运算符的含义.这种做法称为运算符重载, ...
- python中len是什么函数_总结Python中len()函数的使用实例
这篇文章主要简单介绍了Python中的len()函数的使用,包括在四种情况下的使用小例子,是Python学习当中的基础知识,需要的朋友可以参考下 函数:len() 1:作用:返回字符串.列表.字典.元 ...
- python花式输出_Python——花式打印对象的若干种方法
今天是Python专题的第10篇文章,咱们来聊聊Python当中的类.编程 打印实例 咱们先从类和对象当中最简单的打印输出开始讲起,打印一个实例是一个很是不起眼的应用,可是在实际的编程当中却很是重要. ...
- python面试指南
1. a is b 和 a==b一样吗? 首先要了解对于Python对象而言,一般存在三个属性:type 类型, value 值 以及 地址id a==b是一个比较运算符,用于比较两个对象的value ...
- python中len的用法_总结Python中len()函数的使用实例
这篇文章主要简单介绍了Python中的len()函数的使用,包括在四种情况下的使用小例子,是Python学习当中的基础知识,需要的朋友可以参考下 函数:len() 1:作用:返回字符串.列表.字典.元 ...
最新文章
- 准备推出一款新的软件,提供给客户,请问需要做哪种签名?如何收费?
- 计算机主机电池,为什么电脑主机板上面有电池?很多人可能根本不知道
- numpy-自定义ufunc函数和广播
- 【渝粤教育】国家开放大学2018年春季 0674-21T财务管理 参考试题
- Java 面试之数据库
- ReactiveCocoa入门教程--第二部分
- 正在将错误信息上传至服务器,android - 将Base64编码的图像上传到服务器-400服务器错误 - 堆栈内存溢出...
- 【解决方案】SkeyeVSS石油油田钻井无线智能视频监控系统解决方案
- EdgeGallery — MEP — APIs
- nssa和stub_Stub区域和NSSA区域比较
- xp计算机怎样看ip地址,XP电脑ip地址怎么查看?XP系统怎么样查看电脑的IP地址?...
- unity 接入移动MM (3.1.10)
- 时间戳转换格林威治时间
- Maven项目管理工具
- Linux驱动开发: 杂项字符设备
- 六大基酒——朗姆酒的喝法
- python 听歌识曲_Shazam听歌识曲算法解析+python实现-3 检索歌曲
- hdu——4379 ——The More The Better
- 零基础怎么做好海报重叠设计
- rhel8安装libpcap-devel包