20得票3回答
使用SciPy中的wavfile.write在Python中写入wav文件

我有这段代码:import numpy as np import scipy.io.wavfile import math rate, data = scipy.io.wavfile.read('xenencounter_23.wav') data2 = [] for i in rang...

9得票1回答
使用Scipy和Matplotlib将数据拟合到多模态分布

我有一个数据集,想要将其拟合到已知的概率分布中。目的是在数据生成器中使用拟合后的PDF,以便可以从已知(拟合)的PDF中抽取数据样本。数据将用于模拟目的。目前,我只是从正态分布中抽样,这与真实数据不一致,因此模拟结果不准确。 我最初想使用以下方法: 如何使用Scipy(Python)将经验分...

7得票2回答
Scipy无循环快速一维插值

我是一名有用的助手,可以为您进行文本翻译。 我有两个二维数组,x(ni,nj)和y(ni,nj),需要沿一个轴插值。我想对每个ni沿着最后一个轴进行插值。 我写了 import numpy as np from scipy.interpolate import interp1d z =...

7得票1回答
Python检验我的数据是否符合对数正态分布

您好,我有一组呈正偏态的结果分布,因此我想测试它是否适合于对数正态分布或Gumbell分布。 我以前使用过scipy,因此特别希望找到使用该工具的解决方案。 我的数据如下: listofdata =[2.6631285853098805, 1.948723030911822, 0.732...

9得票3回答
迭代拟合多项式曲线

我希望使用以下方法在Python中对数据进行迭代拟合曲线: 拟合多项式曲线(或任何非线性方法) 丢弃偏离曲线均值2个标准差以上的值 重复步骤1和2,直到所有值都在曲线的置信区间内。 我可以按照以下方式拟合多项式曲线: vals = array([0.00441025, 0.00490...

11得票4回答
如何在Python中应用Box-Cox转换?

我有以下形式的数据:X Y 3.53 0 4.93 50 5.53 60 6.21 70 7.37 80 9.98 90 16.56 100 我想找到n,以便将其适配为以下形式的函数: 我正在尝试通过Box-Cox变换确定n。如何在Python...

11得票2回答
Voronoi单元的体积(Python)

我正在使用Python 2.7中的Scipy 0.13.0计算3D Voronoi单元集。我需要获取每个单元的体积以便于加权输出专有模拟结果。是否有简单的方法来做到这一点?毕竟这是一个常见问题或Voronoi单元的常见用途,但我找不到任何信息。以下代码可以运行,并输出了scipy.spatia...

7得票1回答
使用未观测分量模型模拟时间序列

使用statsmodels的UnobservedComponents拟合本地级别模型后,我们正在尝试找到使用结果模拟新时间序列的方法。类似于: import numpy as np import statsmodels as sm from statsmodels.tsa.statespac...

7得票2回答
如何在使用外部C库时与scipy.weave.inline一起使用?

我正在尝试理解weave.inline,以在我的Python程序中包装C代码。下面的代码仅将Numpy数组的所有元素乘以2。 inl.py import numpy import scipy.weave a = numpy.array([1.0, 2.0, 3.0]) N = a.sh...

7得票4回答
在NumPy中从不均匀采样的数据生成均匀采样的数组

问题是我想要减少绘图和分析的数据量。我正在使用Python和Numpy。数据是不均匀采样的,因此有一个时间戳数组和一个对应值的数组。我希望在数据点之间至少有一定的时间间隔。我有一个简单的解决方案,用Python编写,在这里找到了至少在样本之间有一秒钟的指数: import numpy as ...