scipy相关内容

从顶点坐标创建三角网格

给定一组坐标为x和y的二维数据点(左图),有没有一种简单的方法可以在它的顶部(右图)构建一个三角形网格?即返回指示连接了哪些顶点的元组列表。解决方案不是唯一的,但任何合理的网格都可以。 推荐答案 您可以尝试scipy.spatial.Delaunay。从该链接: points = np.array([[0, 0], [0, 1.1], [1, 0], [1, 1]]) from ..
发布时间:2022-09-06 09:37:22 Python

两类稀疏CSR矩阵的差值

我有2个scipy.sparse.csr_matrix如下: A = [ 1 0 1 0 0 1 1 0 0 1 0 0 0 1 0 0 0 0 ] B = [ 1 0 1 0 1 1 1 1 0 1 0 0 1 1 1 0 0 0 ] 我愿意购买在B中出现但在A中没有的“新的”。 C = [ 0 0 0 0 1 0 0 ..
发布时间:2022-08-12 18:55:10 Python

将KDE与来自scipy.Integrate.quad的奇怪行为和设置的带宽进行集成

我正在寻找一种方法来从绘制的分布中获得平均值(期望值),我使用该分布来拟合scipy.stats.Gauss_kde的核密度估计。我记得在我的统计课程中,期望值就是pdf(X)*x从-无穷大到无穷大的积分: 我在代码中使用了scipy.Integrate.quad函数来执行此任务,但我遇到了这种明显奇怪的行为(这可能与KDE中的band with参数有关)。 问题 import ..
发布时间:2022-08-09 15:26:03 Python

一组坐标中的边界点

我有一组很长的点,我想从这些点中提取形成边界的点,我使用了凸包,但对于我的包来说,这是不够的,因为Transporthull只返回形成多边形的最远的点,在所有点都适合的地方,我需要形成透射器的所有点,类似于我附加的图像。我能做什么呢?是否有某种包可以随时使用,而不是实现任何空间算法? 谢谢 推荐答案 您必须使用用于凸面的包。下面是一个例子: import alphashape i ..
发布时间:2022-07-14 16:49:42 Python

将加载的MAT文件转换回NumPy数组

我使用scipy.io.savemat()将大小为5000,96,96的NumPy数组中的图像保存到.mat文件中。 当我想要将这些图像加载回Python时,我使用scipy.io.loadmat(),但这次它们被放入词典中。 如何才能将它们从字典整齐地放到NumPy数组中? 我正在使用scipy.io.loadmat加载MATLAB文件,并希望将其放入NumPy数组中。图像的亮 ..
发布时间:2022-07-11 22:30:52 Python

如何为差分进化增加几个约束?

我的问题与this question中的问题相同,但我不想只向优化问题添加一个约束,而是希望添加多个约束。 所以,例如,我想在x1和x2的和小于5和x2小于3的约束下最大化x1 + 5 * x2(不用说,实际问题要复杂得多,不能像这个问题那样抛到scipy.optimize.minimize中;它只是用来说明问题...)。 我可以对付这样的丑陋黑客: from scipy.opt ..
发布时间:2022-07-11 22:19:12 Python

Interp1d给出了NaN的外推

我正在尝试使用interp1d带外推的interp1d插补具有NaN值的数据: import numpy as np import matplotlib.pyplot as plt x = np.arange(11) y = np.array([np.nan, np.nan, np.nan, 1, np.nan, np.nan, 9, 7, 6, np.nan, np.nan]) f = ..
发布时间:2022-07-05 13:58:46 Python

如何在.wav文件中找到特定声音的时间戳?

我有一个.wav文件,我录制了自己的声音并讲了几分钟。比方说,我想找到我在音频中提到的确切时间。我研究了语音识别,并用Google Speech API进行了一些测试,但我得到的时间戳远远不准确。 作为另一种选择,我录制了一个非常短的.wav文件,我刚才说了&Mike";。我正在尝试比较这两个.wav文件,找出在较长的.wav文件中提到的每个时间戳。我遇到了令人惊叹的SleuthEy ..
发布时间:2022-06-07 21:28:36 Python

Monod生长/降解方程对实验数据的曲线拟合

因此,这里面临的问题是Monod方程与实验数据的曲线拟合。细菌生长和降解有机碳的模型如下: dx/dt=(u*S*X)/(K+S) ds/dt=((-1/Y)*u*S*X)/(K+S) 这些方程使用scipy odeint函数求解。积分后的结果存储在两个向量中,一个用于生长,另一个用于降解。下一步是将该模型与实验观测数据进行曲线拟合,并估计模型参数:U、K和Y。一旦运行代码,就会 ..
发布时间:2022-05-23 18:09:04 Python

Scipy Curve_Fit和局部极小值:尽可能快地获得全局极小值

我手头的问题:我正在使用scipycurve_fit来拟合一条曲线(https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.curve_fit.html),但在许多情况下,为此类曲线估计的参数引用了许多局部“局部”最小值之一,而不是“全局”最小值。考虑到curve_fit的设计方式,这是意料之中的。尽管如此,我真的需要 ..
发布时间:2022-05-23 17:27:29 Python