SciPy库（一）常数与特殊函数与插值

叁歲伎倆 2024-04-17 11:44 45阅读 0赞

## 一、概述 ##

> SciPy是一个开源的Python科学计算库，它建立在NumPy之上，提供了许多有用的科学计算功能。SciPy包括各种科学计算模块，如线性代数、优化、积分、插值、特殊函数、快速傅里叶变换、信号处理和图像处理等。SciPy库的主要特点是其高效性和易用性。

#### 常见模块 ####

<table> 
 <thead> 
  <tr> 
   <th>模块名</th> 
   <th>功能</th> 
   <th>参考文档</th> 
  </tr> 
 </thead> 
 <tbody> 
  <tr> 
   <td>scipy.cluster</td> 
   <td>向量量化</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/cluster.html" rel="nofollow" title="cluster API">cluster API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.constants</td> 
   <td>数学常量</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/constants.html" rel="nofollow" title="constants API">constants API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.fft</td> 
   <td>快速傅里叶变换</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/fft.html" rel="nofollow" title="fft API">fft API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.integrate</td> 
   <td>积分</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/integrate.html" rel="nofollow" title="integrate API">integrate API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.interpolate</td> 
   <td>插值</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/interpolate.html" rel="nofollow" title="interpolate API">interpolate API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.io</td> 
   <td>数据输入输出</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/io.html" rel="nofollow" title="io API">io API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.linalg</td> 
   <td>线性代数</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/linalg.html" rel="nofollow" title="linalg API">linalg API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.misc</td> 
   <td>图像处理</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/misc.html" rel="nofollow" title="misc API">misc API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.ndimage</td> 
   <td>N 维图像</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/ndimage.html" rel="nofollow" title="ndimage API">ndimage API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.odr</td> 
   <td>正交距离回归</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/odr.html" rel="nofollow" title="odr API">odr API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.optimize</td> 
   <td>优化算法</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/optimize.html" rel="nofollow" title="optimize API">optimize API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.signal</td> 
   <td>信号处理</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/signal.html" rel="nofollow" title="signal API">signal API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.sparse</td> 
   <td>稀疏矩阵</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/sparse.html" rel="nofollow" title="sparse API">sparse API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.spatial</td> 
   <td>空间数据结构和算法</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/spatial.html" rel="nofollow" title="spatial API">spatial API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.special</td> 
   <td>特殊数学函数</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/special.html" rel="nofollow" title="special API">special API</a></td> 
  </tr> 
  <tr> 
   <td>scipy.stats</td> 
   <td>统计函数</td> 
   <td><a href="https://docs.scipy.org/doc/scipy/reference/stats.mstats.html" rel="nofollow" title="stats.mstats API">stats.mstats API</a></td> 
  </tr> 
 </tbody> 
</table>

## 二、SciPy常数与特殊函数 ##

![在这里插入图片描述][watermark_type_d3F5LXplbmhlaQ_shadow_50_text_Q1NETiBAaHVzdGxlaQ_size_20_color_FFFFFF_t_70_g_se_x_16_pic_center]

### 2.1 数学常数 ###

import scipy.constants as C
    
    print("π = ", C.pi)
    print("e = ", C.e)
    print("自然对数的底数 = ", C.golden)

### 2.2物理常数 ###

> 常见的常数：
> 
>  *  constants.c：真空光速
>  *  constants.h：普朗克常数
>  *  constants.G：牛顿引力常数
>  *  constants.g：标准重力加速度
>  *  constants.e：基本电荷
>  *  constants.m\_e：电子质量
>  *  constants.m\_p：质子质量
>  *  constants.m\_n：中子质量

更多常量请看[SciPy 常量模块][SciPy]

### 2.3 特殊函数(scipy.special) ###

> scipy.specail模块包含了少量实用的函数，同时还包含很多基础的数学函数，比如贝塞尔基函数、椭圆函数、统计函数的基础函数、超几何函数、正交多项式、开尔文函数、伽玛相关函数等等。大多数是比较基础的函数，通常我们直接使用的较少。

#### 实用函数 ####

>  *  special.round(x)：四舍五入为最接近的整数
>  *  special.cbrt(x)：求立方根
>  *  special.sindg(x)：角度x(单位为度)的正弦
>  *  special.cosdg(x)：角度x(单位为度)的余弦
>  *  special.tandg(x)：角度x(单位为度)的正切
>  *  special.cotdg(x)：角度x(单位为度)的余切

from scipy import special
    
    x = 5.2
    
    # 四舍五入为最接近的整数
    print(special.round(x))
    
    # 求立方根
    print(special.cbrt(x))
    
    # 角度x(单位为度)的正弦
    print(special.sindg(x))
    
    # 角度x(单位为度)的余弦
    print(special.cosdg(x))
    
    # 角度x(单位为度)的正切
    print(special.tandg(x))
    
    # 角度x(单位为度)的余切
    print(special.cotdg(x))

#### 特别函数 ####

>  *  special.gamma(z)：伽玛函数。伽玛函数是阶乘函数在实数和复数系上的扩展
>  *  special.gammaln(z)：伽玛函数的对数
>  *  special.beta(a,b)：贝塔函数
>  *  special.betaln(a,b)：贝塔函数
>  *  special.ellipj(u,m)：雅可比椭圆函数
>  *  special.euler(n,x)：欧拉数计算
>  *  special.binom(n,k)：二项式系数
>  *  special.sici(x)：正弦和余弦的积分

### 2.4 输入输出(scipy.io) ###

#### 2.4.1 matlab数据文件 ####

mditc = io.loadmat(file)：加载.mat类型的matlab文件。返回字典mdict，key为变量名  
io.savemat(file, mdict)：把mdict字典数据保存到matlab格式的.mat文件  
io.whosmat(file)：列出matlab格式的.mat文件中的变量。

#### 2.2 wav声音文件 ####

rate,data = io.wavfile.read(file)：打开wav声音文件。rate为采样频率，data为数据  
io.wavfile.write(file, rate, data)：写入一个简单的未压缩的wav文件

### 2.5 Matlab 数组 ###

#### 转换 ####

import numpy as np
    mat = np.matrix('1,2;3,4')
    mat

![39bc260070a74377b6c3b6ade7eafac0.png][]

#### 单位矩阵 ####

import numpy as np
    import numpy.matlib
    mat1 = np.matlib.identity(5)
    mat1

![8a78eb20783b4a83a2598181c96ac654.png][]

#### 转置 ####

import numpy as np
    import numpy.matlib
    mat = np.matrix('1,2;3,4')
    print(mat)
    print(mat.T)

#### ![19ca907deb374e5abc039653f56b7604.png][] ####

#### 数乘 ####

> `np.multiply()` 函数用于让两个矩阵同样位置的元素相乘。因此，`mat3 = np.multiply(ma, ma1)` 表示将矩阵 `ma` 和 `ma1` 中同样位置的元素相乘，得到一个新的矩阵 `mat3`。

#mat2 = mat1*2
    mat2 = np.multiply(4,mat1)
    mat3 = np.multiply(mat2,mat1)#矩阵位乘
    mat3

![bdf41e63ca5e43148f6dabf1a1c162c8.png][]

####  ####

#### 矩阵加法 ####

import numpy as np
    import numpy.matlib
    ma = np.matrix('1,2;3,4')
    ma1 = np.matrix('3,4;1,2')
    #ma+ma1
    print(np.add(ma, ma1))

![a285ffdf44894dc9b8e8d300ae5abfb8.png][]

#### 矩阵乘法 ####

import numpy as np
    import numpy.matlib
    ma = np.matrix('1,2;3,4')
    print(ma)
    ma1 = np.matrix('3,4;1,2')
    print(ma1)
    # ma*ma1
    mat3 = np.dot(ma, ma1)  # 矩阵相乘
    mat3

#### ![cb1bc4ae69c54561b15ed9d7ae5744eb.png][] ####

#### 矩阵的逆 ####

![f4c8708f2f5e413ea4e17b016c9af481.png][]

import numpy as np
    import numpy.matlib
    ma = np.matrix('1,2;3,4')
    ma.I

![08961988d2cb4c8fbf6fcdf9c6d0e11a.png][]

## 三、插值与拟合 ##

> 插值和拟合是数学中的两个重要概念。插值是指通过已知数据点求出一个函数，使得这个函数在已知数据点上的函数值与已知数据点上的函数值相等。而拟合是指通过已知数据点求出一个函数，使得这个函数在已知数据点上的函数值与已知数据点上的函数值的差的平方和最小

![588a2bece7d9426d8842f7df9f87ee22.png][]

### 3.1插值法 ###

![46ae54484fbf415fa66b8a52094b5b7f.png][]

原理

![72614b0af1bc4438a0f0b2ed5aed3960.png][]

#### ![bcdf63cb44d2453892718298dc955856.png][]3.1.1单变量插值， 一维插值interpld() ####

> interp1d函数是SciPy库中interpolate模块中最常用的插值函数之一。它可以根据给定的数据点和插值方法，返回一个可调用的函数对象，该对象可以用于估计新点的值。interp1d函数的详细参数如下
> 
>     scipy.interpolate.interp1d(x, y, kind='linear', axis=-1, copy=True,
>     bounds_error=None, fill_value=nan, assume_sorted=False)
> 
> 其中，x和y是一维数组，表示要进行插值的数据点。kind是指定插值方法的字符串参数，可以取值为’linear’、‘nearest’、‘zero’、‘slinear’、‘quadratic’、‘cubic’。bounds\_error是一个布尔型参数，表示是否在超出插值范围时引发异常。fill\_value是一个浮点型参数，表示在超出插值范围时要使用的填充值，assume\_sorted参数是一个布尔型参数，表示输入的x数组是否已经按升序排列。

>  *  kind：指定插值的样条曲线的阶数，可以是’linear’（线性），‘nearest’（最近邻），‘zero’（零阶），‘slinear’（一阶样条），‘quadratic’（二阶样条），‘cubic’（三阶样条）等。

import numpy as np
    from scipy import interpolate as intp
    import matplotlib.pyplot as plt
    
    # 生成等差数列，从0到4，共12个数
    x = np.linspace(0,4,12)
    # 计算cos(x^2/3+4)的值
    y = np.cos(x**2/3+4)
    
    print (x)
    print (y)
    
    # 画出散点图
    plt.plot(x,y,"o")
    # 画出平滑曲线
    f = intp.interp1d(x,y,kind="cubic")
    xnew = np.linspace(0,4,30)
    ynew = f(xnew)
    plt.plot(xnew,ynew,"-")
    
    plt.show()

#### ![76b7b1272fdd4b0887f1fdd3075e5fb2.png][] ####

[watermark_type_d3F5LXplbmhlaQ_shadow_50_text_Q1NETiBAaHVzdGxlaQ_size_20_color_FFFFFF_t_70_g_se_x_16_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/b937f3f3d53841d2b8e64d40c0f6b4e6.png
[SciPy]: https://www.runoob.com/scipy/scipy-constants.html
[39bc260070a74377b6c3b6ade7eafac0.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/ee78d785a1fe42f083d69ad7b3c3b258.png
[8a78eb20783b4a83a2598181c96ac654.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/2ed2402cef9a46e886430bc266c8ba4f.png
[19ca907deb374e5abc039653f56b7604.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/974594884a6b4aff8fdbe38cefcbdaef.png
[bdf41e63ca5e43148f6dabf1a1c162c8.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/419f9619c4a8480ca11411baa90ab8c8.png
[a285ffdf44894dc9b8e8d300ae5abfb8.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/fce581c45db4493c9d2e198eb7d21cb6.png
[cb1bc4ae69c54561b15ed9d7ae5744eb.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/3eb8646f0bb24004bb93d89019e43843.png
[f4c8708f2f5e413ea4e17b016c9af481.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/9222386cef8f4bc1b9a8918b11278ea0.png
[08961988d2cb4c8fbf6fcdf9c6d0e11a.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/91a3acd03c1a49fba7ffb15b909c2711.png
[588a2bece7d9426d8842f7df9f87ee22.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/94d371dfab2c4a798c62156e0a050957.png
[46ae54484fbf415fa66b8a52094b5b7f.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/8d29afdf8b78479c994248c9c575bad8.png
[72614b0af1bc4438a0f0b2ed5aed3960.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/74c80e0a82f744e2a282f26d46aff464.png
[bcdf63cb44d2453892718298dc955856.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/079e4756229949d09860c4607fa47f55.png
[76b7b1272fdd4b0887f1fdd3075e5fb2.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/12/370197acf1d845938b89f9f95f4712f5.png