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. nfp 网络共享服务器 搭建与配置
  2. 2019 课程设计个人报告
  3. 洛谷P2680 运输计划(倍增LCA + 树上差分 + 二分答案)
  4. 把老赵的页面缓存片断改一下,呵呵
  5. oracle lsnrctl status unknown,数据库监听卡住 Oracle lsnrctl status 卡在connecting ......
  6. 给你这张图,你能搜索到来历吗
  7. 基础学习笔记之opencv(14):随机数发生器绘制文字
  8. mysql join 索引 无效_ORACLE MYSQL中join 字段类型不同索引失效的情况-阿里云开发者社区...
  9. SSM项目-山东医院-可行性配置-1
  10. 遇到Python中文目录名问题,未解决
  11. 我新鲜出炉的 plist 图集拆图工具!github
  12. 求可以进行图片抽签的小程序
  13. python余弦函数_Python cos() 函数
  14. 云服务器文件传送,云服务器文件传送工具
  15. python爬取12306列车信息自动抢票并自动识别验证码(二)selenium登录验证篇
  16. 液化石油气采样钢瓶标准取样流程
  17. 【P2P租车】宝驾租车:学大创始人李如彬再创业
  18. Java基于JSP校园二手闲置商品交易系统
  19. Java工程师核心能力_java程序员的核心竞争力是什么?
  20. c盘里program files文件夹和program files(x86)文件夹

热门文章

  1. 命令方块召唤别墅指令_我的世界:如何在MC召唤实体303?告示牌的数字,才是关键!...
  2. 开源时代:刘韧对话任旭东崔宝秋章文嵩蒋涛
  3. 饿了么分布式服务治理及优化经验(含PPT)
  4. STM32多通道DMA—ADC采样
  5. linux系统etc什么意思,etc.是什么意思 linux下的etc是什么意思
  6. oracle设置查看格式,Oracle 设置显示格式
  7. 实验吧-web-天下武功唯快不破
  8. macOS Big Sur 11.7.1 (20G918) 正式版 ISO、PKG、DMG、IPSW 下载
  9. 想要快速、免费的翻译外文PDF文档?Get翻译神器
  10. java protobuf extend_protobuf中extension的使用