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-特殊函数相关推荐

  1. Day10 Python基础之特殊函数(八)

    一些特殊函数 1.递归函数(recursion) 递归函数的定义:在函数内部,可以调用其他函数.如果一个函数在内部调用自身本身,这个函数就是递归函数. 递归函数的优点:是定义简单,逻辑清晰.理论上,所 ...

  2. python中的特殊函数__call__

    1 # _*_ coding:utf-8 _*_ 2 3 ''' 4 一个类实例可以变成一个可以调用的对象,只需要实现特殊方法__call__ 5 所以可以看出在python中函数也是对象 6 ''' ...

  3. python保存内容到文件(text、json、csv)

    python保存内容到文件(text.json.csv) 在开发人员的日常中,将数据保存到文件是最常见的编程任务之一. 通常,程序需要一些输入并产生一些输出.在许多情况下,我们希望将这些结果持久化.我 ...

  4. python花式编码_Python——花式打印对象的若干种方法

    今天是Python专题的第10篇文章,我们来聊聊Python当中的类. 打印实例 我们先从类和对象当中最简单的打印输出开始讲起,打印一个实例是一个非常不起眼的应用,但是在实际的编程当中却非常重要.原因 ...

  5. python格式化输出 xz_Python——详解__str__, __repr__和__format__

    今天是Python专题的第10篇文章,我们来聊聊Python当中的类. 打印实例 我们先从类和对象当中最简单的打印输出开始讲起,打印一个实例是一个非常不起眼的应用,但是在实际的编程当中却非常重要.原因 ...

  6. 一篇文章带你了解Python运算符重载

    回复"python"即可获赠从入门到进阶共10本电子书 今 日 鸡 汤 不堪玄鬓影,来对白头吟. 您可以根据所使用的操作数来更改Python中运算符的含义.这种做法称为运算符重载, ...

  7. python中len是什么函数_总结Python中len()函数的使用实例

    这篇文章主要简单介绍了Python中的len()函数的使用,包括在四种情况下的使用小例子,是Python学习当中的基础知识,需要的朋友可以参考下 函数:len() 1:作用:返回字符串.列表.字典.元 ...

  8. python花式输出_Python——花式打印对象的若干种方法

    今天是Python专题的第10篇文章,咱们来聊聊Python当中的类.编程 打印实例 咱们先从类和对象当中最简单的打印输出开始讲起,打印一个实例是一个很是不起眼的应用,可是在实际的编程当中却很是重要. ...

  9. python面试指南

    1. a is b 和 a==b一样吗? 首先要了解对于Python对象而言,一般存在三个属性:type 类型, value 值 以及 地址id a==b是一个比较运算符,用于比较两个对象的value ...

  10. python中len的用法_总结Python中len()函数的使用实例

    这篇文章主要简单介绍了Python中的len()函数的使用,包括在四种情况下的使用小例子,是Python学习当中的基础知识,需要的朋友可以参考下 函数:len() 1:作用:返回字符串.列表.字典.元 ...

最新文章

  1. 准备推出一款新的软件,提供给客户,请问需要做哪种签名?如何收费?
  2. 计算机主机电池,为什么电脑主机板上面有电池?很多人可能根本不知道
  3. numpy-自定义ufunc函数和广播
  4. 【渝粤教育】国家开放大学2018年春季 0674-21T财务管理 参考试题
  5. Java 面试之数据库
  6. ReactiveCocoa入门教程--第二部分
  7. 正在将错误信息上传至服务器,android - 将Base64编码的图像上传到服务器-400服务器错误 - 堆栈内存溢出...
  8. 【解决方案】SkeyeVSS石油油田钻井无线智能视频监控系统解决方案
  9. EdgeGallery — MEP — APIs
  10. nssa和stub_Stub区域和NSSA区域比较
  11. xp计算机怎样看ip地址,XP电脑ip地址怎么查看?XP系统怎么样查看电脑的IP地址?...
  12. unity 接入移动MM (3.1.10)
  13. 时间戳转换格林威治时间
  14. Maven项目管理工具
  15. Linux驱动开发: 杂项字符设备
  16. 六大基酒——朗姆酒的喝法
  17. python 听歌识曲_Shazam听歌识曲算法解析+python实现-3 检索歌曲
  18. hdu——4379 ——The More The Better
  19. 零基础怎么做好海报重叠设计
  20. rhel8安装libpcap-devel包

热门文章

  1. Koo叔说Shader—最基本的Shader
  2. 【NOIP2017提高A组冲刺11.2】救赎(数学期望)
  3. debian7系统设置固定IP
  4. 蓝桥杯2021年PYTHON 真题,跳房子
  5. Python3,一行代码实现文件夹共享,看到结果我酸了~
  6. 算法与数据结构 03 匹配问题
  7. Deep Learning论文翻译(Nature Deep Review)
  8. 【问链财经-区块链基础知识系列】 第二十三课 一张纸搞定B2B供应链金融产品设计思路
  9. 区块链代采供应链金融平台方案设计手稿
  10. matlab画图不显示中文_如何解决 Matlab 画图时中文显示乱码的问题?