SciPy教程 - 常数constants和特殊函数库special

http://blog.csdn.net/pipisorry/article/details/39615047

scipy.constants常数模块

数学常量Mathematical constants

pi Pi
golden Golden ratio
golden_ratio Golden ratio

物理常数Physical constants

constants模块包含了众多的物理常数。

>>> from scipy import constants as C
>>> C.c #真空中的光速
299792458.0
>>> C.h #普朗克常数
6.62606930800080626-34

在字典physical_constants中，以物理常量名为键，对应的值是一个含有三个元素的元组， 分别为常数值、单位及误差。
例如下面的程序可以查看电子质量：
>>> C.physical_constants["electron mass"]
 (9.1093825999999998e-31,’kg', 1.5999999999999999e-37)

单位信息

constants模块中还包括许多单位信息
>>> C.mile # 1英里等于多少米 
1609.3439999999998 
>>>C.inch     # 1英寸等于多少米 
0.025399999999999999 
>>> C.gram   # 1克等于多少千克 
0.001
>>> C.pound   # 1磅等于多少千克 

0.45359236999999997

[Constants (scipy.constants)]

皮皮blog

scipy.special特殊函数模块

special模块是一个非常完整的函数库，其中包含了基本数学函数、特殊数学函数以及NumPy中出现的所有函数。

常用特殊函数求值：

伽玛函数

gamma函数是阶乘函数在实数和复数范围上的扩展。

伽玛函数值的计算

>>> import scipy.special as S 
>>> S.gamma(0.5)
1.772458509055159
>>> S.gamma(1+1j)  # gamma 函数支持复数 (0.49801566811835629-0.15494982836181106j ) 
>>> S.gamma(1000) 
inf 

对数伽玛函数

gamma函数的增长速度非常快，因为1000的阶乘 已经超过了双精度浮点数的表示范围，因此结果是无穷大。为了计算更大的范围，可以使用 S.gammaln()。

S.gammaln(x)计算ln(|gamma(x)|)的值，它使用特殊的算法，直接计算gamma函数的对数值，因此可以表示更大的范围。

>>>S.gammaln(1000)
5905.2204232091817

log1p(x)

special模块中的某些函数并不是数学意义上的特殊函数，例如log1p(x)计算log(1+x)的值。 这是由于浮点数的精度有限，无法很精确地表示十分接近1的实数。例如无法用浮点数表示”1 + 1e-20”的值，因此“log(1+1e-20)”的值为0,而当使用log1p()时，则可以很精确地计算。实际上当x非常小时，log1p(x)约等于x 。

>>> 1 + 1e-20
1.0
>>> log(1+1e-20)
0.0
>>> S.log1p(1e-20)
9.9999999999999995e-21

[Special functions (scipy.special)]

from:http://blog.csdn.net/pipisorry/article/details/39615047

ref: