scipy模块的简介与使用

阿原~

已于 2024-09-21 17:50:20 修改

阅读量793

点赞数 8

文章标签： scipy

于 2024-09-21 17:47:40 首次发布

本文链接：https://blog.csdn.net/m0_73901987/article/details/142370477

版权

一. scipy模块的简介

scipy作为世界上著名的，开源的高级科学计算库，从20世纪90年代发展到现在已经每年有千万的下载量。 Scipy是基于Numpy构建的一个集成了多种数学算法和方便的函数的Python模块。通过给用户提供一些高层的命令和类来操作和可视化数据，SciPy在python交互式会话中，大大增加了操作和可视化数据的能力。通过SciPy，Python的交互式会话变成了一个数据处理和一个系统原型system-prototyping环境，可以与MATLAB、IDL、Octave、R-Lab和SciLab等系统相匹敌。
更重要的是，在Python中使用SciPy，还可以同时用一门强大的语言—Python来开发复杂和专业的程序。用SciPy写科学应用，还能获得世界各地的开发者开发的模块的帮助，受益于世界各地的开发人员在软件领域的许多小众领域中开发的附加模块。从并行编程到web到数据库子例程到各种类，Python程序员都可以使用。这些强大的功能，SciPy都有，特别是它的数学库。

二. Python如何下载scipy模块

步骤一：安装pip

Python安装库的驱动是pip，在安装scipy库之前要确保自己的Python有着pip驱动。

# 下载并安装Python的pip
python get-pip.py

步骤二：使用pip安装scipy库

使用pip安装驱动来安装scipy库

# 使用pip来安装scipy库
pip install scipy

Python3.以来使用pip有时会产生报错，当产生报错的时候可以尝试使用pip3来进行安装，pip和pip3一般Python系统会自带安装

# 使用pip3来进行scipy库安装
pip3 install scipy

有的开发人员喜欢使用anconda来进行Python环境的搭建，在anconda中进行scipy库的安装

# anconda中安装scipy库
conda install scipy

三. scipy库的常见使用

主要模块与功能：

Integration（数值积分）:

提供了一套全面的数值积分方法，包括单变量和多变量函数的定积分、不定积分以及自适应辛普森法则等算法。

scipy库中存放了大量的数值积分的计算公式函数，如矩形法、梯形法、求和法等。

下面使用梯形法（一重积分）来演示数值积分

import numpy as np
from scipy.integrate import trapz
 
# 定义被积函数
def f(x):
    return x**2
 
# 定义积分区间
a = 0
b = 1
 
# 生成一系列点用于计算数值积分
x = np.linspace(a, b, 1000)
 
# 计算数值积分
integral = trapz(y=f(x), x=x)
 
print(f"数值积分结果为: {integral}")

scipy数值积分的方法包括

1.quad：用于一维积分的通用积分函数。

2.dblquad：用于二维积分的通用积分函数。

3.tplquad：用于三维积分的通用积分函数。

4.nquad：用于n维积分的通用积分函数。

5.fixed_quad：用于固定积分的数值积分方法。

6.quadrature：用于自适应积分的数值积分方法。

7.romberg：用于Romberg积分的数值积分方法。

8.simps：用于Simpson积分的数值积分方法。

9.trapz：用于梯形积分的数值积分方法。

10.cumtrapz：用于累积梯形积分的数值积分方法。

Optimization（优化）:

包括非线性最小化、线性规划、二次规划、约束优化等功能，适用于寻找函数的最小值、最大值或者满足条件下的最优解。

scipy优化器常常于机器学习，可以定义自己的优化函数来通过scipy优化器来进行封装使用。

如下面的约束优化案例

from scipy.optimize import minimize
 
# 定义目标函数
def objective(x):
    return x[0]**2 + x[1]**2
 
# 定义约束条件
constraints = ({'type': 'eq', 'fun': lambda x: x[0] - x[1]})
 
# 进行优化
result = minimize(objective, [1, 1], constraints=constraints)
print(f"Optimal solution with constraints: {result.x}")  # 输出符合约束条件的最小值点

constraints 参数用于定义等式约束（'type': 'eq'）或不等式约束（'type': 'ineq'）。
常见的scipy库优化

# minimize
# minimize 函数用于求解多维优化问题，可以处理带有约束条件的优化问题。
scipy.optimize.minimize(fun, x0, args=(), method='BFGS', jac=None, constraints=(), bounds=None, tol=None, options=None)
fun：目标函数，需要最小化的函数。
x0：初始猜测值。
method：优化算法，如 'BFGS'、'L-BFGS-B'、'SLSQP' 等。
constraints：优化问题的约束条件。
bounds：变量的上下界。
jac：目标函数的雅可比矩阵，用于梯度信息


# minimize_scalar
# minimize_scalar 用于一维函数的最小化，适用于只有一个变量的优化问题。
scipy.optimize.minimize_scalar(fun, args=(), method='brent', bounds=None, bracket=None, tol=None, options=None)
method：优化方法，如 'brent'、'bounded'。
bounds 和 bracket 用于限制优化范围。

Interpolation（插值）:

提供多种插值方法，如一维、二维和多维插值，包括线性插值、多项式插值（包括样条插值）、不规则网格插值等。

在SciPy库中，插值函数主要位于scipy.interpolate模块中。这个模块提供了多种插值方法，适用于不同的数据和应用场景。

interp1d函数（一维插值）

interp1d函数用于对一维数据点进行插值。它可以创建一个插值函数，该函数可以用来预测新的数据点的值。

from scipy.interpolate import interp1d

# 给定的数据点
x = np.linspace(0, 10, 11)
y = np.sin(x)

# 创建一个线性插值函数
interp_func = interp1d(x, y, kind='linear')

# 使用插值函数预测新的数据点
x_new = np.linspace(0, 10, 101)
y_new = interp_func(x_new)

kind参数可以设置插值的类型如'linear', 'nearest', 'zero', 'slinear', 'quadratic', 'cubic'等。

常见的scipy库的插入方式

# interp1d函数（一维插值）
# interp1d函数用于对一维数据点进行插值。它可以创建一个插值函数，该函数可以用来预测新的数据点的值。
interp1d(x, y, kind='linear')

# UnivariateSpline类（一维样条插值）
# UnivariateSpline类用于对一维数据点进行样条插值。它构造一个使用B样条的插值函数
UnivariateSpline(x, y, s=0)
s参数控制样条的平滑度，值越大，插值曲线越平滑。

# griddata函数（多元插值）
# griddata函数用于在多维空间中进行插值。它适用于不规则分布的数据点，可以在规则或不规则的网格上进行插值。
griddata(points, values, (xi, yi), method='cubic')
method参数可以设置插值算法，如'linear', 'nearest', 'cubic'等

# Rbf类（径向基函数插值）
# Rbf类用于使用径向基函数进行插值，特别适合于稀疏数据集

Signal Processing（信号处理）:

实现了滤波、FFT变换、窗口函数、信号恢复、频谱分析等相关功能。

SciPy库的scipy.signal模块提供了一系列用于信号处理的函数和类。

常见的信号处理方法有

# filtfilt()函数：零相位滤波
# filtfilt()函数用于对信号进行零相位滤波，即前后滤波相结合，以消除滤波引入的相位延迟。这个函数接收滤波器的系数和要滤波的信号，并返回滤波后的信号。
from scipy.signal import filtfilt

# 使用butter()设计的滤波器系数进行滤波
y = filtfilt(b, a, x)

# fft()函数：快速傅里叶变换
fft()函数用于计算信号的快速傅里叶变换（FFT），它将信号从时域转换到频域，便于分析信号的频率成分。
from scipy.fftpack import fft

# 计算信号的FFT
Y = fft(x)

# lfilter()函数：线性滤波
# lfilter()函数用于对信号进行线性滤波，它根据滤波器的分子和分母系数对数据进行滤波处理。
from scipy.signal import lfilter

# 使用lfilter()应用滤波器
y = lfilter(b, a, x)

# periodogram()函数：计算功率谱估计
# periodogram()函数用于计算信号的周期图，它是功率谱密度的一种估计方法，通常用于分析信号的频率特性。
from scipy.signal import periodogram

# 计算信号的周期图
Pxx, f = periodogram(x, fs=1000)

Linear Algebra（线性代数）:

提供了大量线性代数操作，如解线性方程组、特征值与特征向量计算、奇异值分解、LU分解等。

scipy的线性代数的方法

基础线性代数运算
inv()：计算矩阵的逆。
det()：计算方阵的行列式。
solve()：求解线性方程组Ax = b。
norm()：计算向量或矩阵的范数。
矩阵分解
lu()：计算矩阵的LU分解。
svd()：执行奇异值分解（SVD）。
eig()：计算矩阵的特征值和特征向量。
eigvals()：仅计算矩阵的特征值。
cholesky()：计算对称正定矩阵的Cholesky分解。
qr()：计算矩阵的QR分解。
矩阵函数
expm()：计算矩阵的指数。
logm()：计算矩阵的对数。
cosm()、sinm()、tanm()：计算矩阵的三角函数。
特殊矩阵运算
block_diag()：创建块对角矩阵。
circulant()：创建循环矩阵。
线性方程组求解
lstsq()：计算线性最小二乘问题的解。
其他功能
linspace()：生成线性间隔的数组，类似于NumPy的numpy.linspace。
logspace()：生成对数间隔的数组，类似于NumPy的numpy.logspace。