title: "SciPy-数值计算库" date: 2026-05-24T14:11:00Z 1import numpy as np 2import pylab as pl 1import matplotlib as mpl 2mpl.rcParams['font.sans-serif'] = ['SimHei'] SciPy-数值计算库 1import scipy 2scipy.version '1.0.0' 常数和特殊函数 1from scipy import constants as C 2print (C.c) # 真空中的光速 3print (C.h) # 普朗克常数 299792458.0 …

阅读更多

title: "10.scipy-spatial" date: 2026-05-24T14:11:03Z 1import numpy as np 2import pylab as pl 3from scipy import spatial 1import matplotlib as mpl 2mpl.rcParams['font.sans-serif'] = ['SimHei'] 空间算法库-spatial 计算最近旁点 1x = np.sort(np.random.rand(100)) 2idx = np.searchsorted(x, 0.5) 3print (x[idx], x[idx - 1]) #距离0.5最近的数是这两 …

阅读更多

title: "f计算方程组的误差，[1,1,1]是未知数的初始值" date: 2026-05-24T14:11:06Z 1import pylab as pl 2import numpy as np 1import matplotlib as mpl 2mpl.rcParams['font.sans-serif'] = ['SimHei'] 拟合与优化-optimize 非线性方程组求解 1from math import sin, cos 2from scipy import optimize 3 4def f(x): #❶ 5 x0, x1, x2 = x.tolist() #❷ 6 return [ 7 5*x1+3, …

阅读更多

title: "3.scipy-linalg" date: 2026-05-24T14:11:09Z 1import pylab as pl 2import numpy as np 3from scipy import linalg 1import matplotlib as mpl 2mpl.rcParams['font.sans-serif'] = ['SimHei'] 线性代数-linalg 解线性方程组 1import numpy as np 2from scipy import linalg 3 4m, n = 500, 50 5A = np.random.rand(m, m) 6B = …

阅读更多

title: "4.scipy-stats" date: 2026-05-24T14:11:12Z 1import numpy as np 2import pylab as pl 3from scipy import stats 1import matplotlib.pyplot as plt 2import matplotlib as mpl 3mpl.rcParams['font.sans-serif'] = ['SimHei'] 统计-stats 连续概率分布 1from scipy import stats 2[k for k, v in stats.dict.items() if isinstance(v, …

阅读更多

title: "调用ode对lorenz进行求解, 用两个不同的初始值" date: 2026-05-24T14:11:14Z 1import pylab as pl 2import numpy as np 3from scipy import integrate 4from scipy.integrate import odeint 5import matplotlib as mpl 6mpl.rcParams['font.sans-serif'] = ['SimHei'] 数值积分-integrate 球的体积 1def half_circle(x): 2 return (1-x**2)**0.5 1N = 10000 2x …

阅读更多

title: "设计一个带通滤波器：" date: 2026-05-24T14:11:17Z 1import pylab as pl 2import numpy as np 3from scipy import signal 1import matplotlib as mpl 2mpl.rcParams['font.sans-serif'] = ['SimHei'] 信号处理-signal 中值滤波 1#%fig=使用中值滤波剔除瞬间噪声 2t = np.arange(0, 20, 0.1) 3x = np.sin(t) 4x[np.random.randint(0, len(t), 20)] += …

阅读更多

title: "X-Y轴分为15*15的网格" date: 2026-05-24T14:11:20Z 1import numpy as np 2import pylab as pl 3from scipy import interpolate 1import matplotlib as mpl 2mpl.rcParams['font.sans-serif'] = ['SimHei'] 插值-interpolate 一维插值 WARNING 高次interp1d()插值的运算量很大，因此对于点数较多的数据，建议使用后面介绍的UnivariateSpline()。 1#%fig=interp1d的各阶插值 2from scipy …

阅读更多

title: "8.scipy-sparse" date: 2026-05-24T14:11:23Z 1%matplotlib inline 2import numpy as np 3import pylab as pl 4from scipy import sparse 5from scipy.sparse import csgraph 稀疏矩阵-sparse 稀疏矩阵的储存形式 1from scipy import sparse 2a = sparse.dok_matrix((10, 5)) 3 4a[2, 3] = 1.0 5a[3, 3] = 2.0 6a[4, 3] = 3.0 7print(a.keys()) …

阅读更多

title: "9.scipy-ndimage" date: 2026-05-24T14:11:26Z 1import numpy as np 2import pylab as pl 图像处理-ndimage 形态学图像处理 1import numpy as np 2 3def expand_image(img, value, out=None, size = 10): 4 if out is None: 5 w, h = img.shape 6 out = np.zeros((wsize, hsize),dtype=np.uint8) 7 8 tmp = …

阅读更多