scipy 相关问题

SciPy是一个用于Python编程语言的算法和数学工具的开源库。

为什么 scipy 的快速傅立叶变换操作会给数据添加如此多的噪声?

如果我创建两个非常简单的信号,一个是指数衰减,另一个是高斯信号,我可以使用 scipy 的 fft 和 ifft 函数非常轻松地对这两个信号进行卷积: 卷积 = ifft(fft(衰变...

回答 1 投票 0

如何在指数(时间序列)和网格数据之间进行 Spearman 相关?

我需要执行类似于此处可用的 Spearman 相关性(皮尔逊相关性),其中我可以在网格数据和气候指数(或自定义时间序列)之间进行关联。

回答 1 投票 0

用户警告 scipy.stats.shapiro

我使用夏皮罗威尔克斯正态性检验得到此错误 UserWarning:scipy.stats.shapiro:对于 N > 5000,计算的 p 值可能不准确。当前 N 为 1048575。 stat, p = st.shapiro(

回答 1 投票 0

子类化 scipy.stats.rv_continuous 以实现截断分布

我想从截断的互补误差函数创建一个随机变量。为此,我想对 scipy.stats.rv_continuous 进行子类化。我已经检查了 Subclassing scipy's Continuous

回答 1 投票 0

无法上传包含稀疏矩阵列表的pickle对象。我收到 ModuleNotFoundError: No module named 'scipy.sparse._arrays'

进口泡菜 从 scipy.sparse 导入 csr_matrix 当我尝试打开包含我之前多次使用过的稀疏矩阵列表的 Pickle 文件时,没有任何问题,我得到以下信息...

回答 1 投票 0

scipy.ndimage 导入测量。如何导入“标签”

来自 pylab 导入 * 从 scipy.ndimage 导入测量值 长=100 pv = [0.2,0.3,0.4,0.5,0.6,0.7] z = 兰特(L,L) 对于范围内的 i(len(pv)): p = pv[i] m = z from pylab import * from scipy.ndimage import measurements L = 100 pv = [0.2,0.3,0.4,0.5,0.6,0.7] z = rand(L,L) for i in range(len(pv)): p = pv[i] m = z<p lw, num = measurements.label(m) area = measurements.sum(m, lw, index=arange(lw.max() + 1)) areaImg = area[lw] subplot(2,3,i+1) tit = 'p='+str(p) imshow(areaImg, origin='lower') title(tit) axis() 我正在尝试在 vscode 中运行此代码,但它在终端上显示此内容 DeprecationWarning:请从 label 命名空间导入 scipy.ndimage; scipy.ndimage.measurements 命名空间已弃用,并将在 SciPy 2.0.0 中删除。 lw, num = 测量值.label(m) DeprecationWarning:请从 sum 命名空间导入 scipy.ndimage; scipy.ndimage.measurements 命名空间已弃用,并将在 SciPy 2.0.0 中删除。 面积=测量值.sum(m,lw,index=arange(lw.max() + 1)) 应该显示这个 我实际上是新人,不知道如何导入标签和总和。 在第二行,您可以使用 导入 from scipy.ndimage import measurements 和 label,而不是 sum from scipy.ndimage import label, sum 然后您需要调整以从代码的其余部分中删除 measurements,如下所示: lw, num = label(m) # this works because you imported label area = sum(m, lw, index=arange(lw.max() + 1)) # this works because you imported sum

回答 1 投票 0

如何在python中找到原始wav文件中带有噪声的子波文件的位置?

原始wav文件的长度为20分钟。 sub wav 文件的长度为 20-40 秒。 1分钟内找到位置。 我尝试使用 scipy.io.wavfile 将文件读取为 numpy 数组并使用...

回答 1 投票 0

Python:ValueError:对象对于所需数组来说太深。我该如何解决这个问题?

VScode 在第 12 行和第 22 行显示错误。这是代码。 从 scipy.io 导入 wavfile 将 numpy 导入为 np def change_pitch(输入文件, 输出文件, 半音, nsemitones): # 读取 .wav 文件...

回答 1 投票 0

创建 numpy 数组的缩减长度重采样

我有一个一定长度的 numpy 数组,并且想要生成一个更短的数组,其中较短数组中的每个值都是基于原始数组中的值的 bin 聚合的。在我的具体说明中...

回答 1 投票 0

使用少量控制点(~200)的 RBF 对大量点(~150k)进行高效重复插值

我正在尝试解决最小二乘问题,其中我的自变量是数百个控制点中每个控制点的少数值,然后将其插值到大量其他点...

回答 1 投票 0

如何在Python中使用SVD求矩阵的一般逆?

所以我试图使用 svd 找到矩阵 A 的广义逆 A+ 这是我的 svd 代码 A4 = np.array([[3, 1,2,3,4], [5,2,3,4,5], [7,4,5,6,7]]) 打印(A4) U4,...

回答 1 投票 0

python / scipy 中的牛顿法

我正在尝试在Python上使用牛顿算法来获取函数的根。即使我更改精度级别,我也会出现运行时错误。您能帮我了解如何改进它吗? 最好的,...

回答 2 投票 0

如何使用 scipy 从频谱图中回归多个高斯峰值?

我有这段代码,旨在适应此处的数据 DriveFilelink 函数 read_file 用于以结构化方式提取信息。然而,我在高斯方面遇到了挑战......

回答 1 投票 0

如何打开和读取.nc文件?

我在打开 .nc 文件并将其转换为 .csv 文件时遇到问题,但我仍然无法读取它们(指第一部分)。我看到了这个链接也这个链接,但我找不到如何打开......

回答 1 投票 0

错误:无法安装 1.18.0,因为这些软件包版本具有冲突的依赖项

ubuntu:20.04.6 蟒蛇:3.10.14 还有其他问题,但我可以先问一下吗? 我在阅读论文和安装 scDALI 时陷入困境。我通过搜索找不到方法。什么该...

回答 1 投票 0

为什么即使数据相同,stats.chisquare 和 stats.chi2_contingency 的结果也不同

我有一个名为“data_count”的数组,存储9个数字组的计数。我使用本福德定律生成一个名为“expected_counts”的预期计数数组。我想测试两个数组是否有...

回答 1 投票 0

为什么这个约束函数比类似的函数慢得多?如何提高 Scipy Optimize 中的速度?

我有一个优化,大约有 22 个约束函数。我尝试将其转换为单个约束函数,但优化器花费了 10 倍的时间。有什么办法可以降低速度吗? 是...

回答 1 投票 0

尝试整合一个作为变量函数的矩阵,但遇到了一个我不知道如何解决的问题

我想对矩阵进行积分,作为求解具有时间相关势的薛定谔方程的一部分。 该矩阵由 4 个 8x8 矩阵和标量函数组成。第一个矩阵描述...

回答 1 投票 0

如何在使用 Hypergeometric Confluence 函数时调试此错误?

我正在尝试使用 scipy.special 库中的超几何汇合函数定义此函数: def Fl(x,l,p): F=sc.hyp1f1(l+1-(1j*p),2*(l+1),-2*1j*x) s=sc.gamma(l+1+(1j*p))/math.gamma...

回答 1 投票 0

使用 SciPy 进行曲线拟合未能给出正确的拟合

我有两个 NumPy 数组 x_data 和 y_data。当我尝试使用二阶阶跃响应模型和 scipy.optimize.curve_fit 与此代码来拟合我的数据时: 将 numpy 导入为 np 来自 scipy.optimize imp...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.