scipy

    0热度

    1回答

    我目前正在使用scipy绘制数据矩阵X中的一组特征变量的正常可靠性图。但是,我使用的模块没有返回我的r^2值。这是我的简单代码: Data_Matrix=pd.read_csv('My csv') My_datum=My_Data.as_matrix() #loop through all feature variables for i in range(My_datum.shape[1])

    0热度

    1回答

    获取python中的sparse.csr_matrix矩阵的中值(沿着行和列)的最佳方法是什么? PS:本webpage可是没有中位数

    0热度

    1回答

    我确实找到了计算点群集的中心坐标的方法。然而,当初始坐标的数量增加时(我有大约100 000个坐标),我的方法非常慢。 瓶颈是代码中的for循环。我试图通过使用np.apply_along_axis来删除它,但发现这只不过是一个隐藏的Python循环。 是否有可能以矢量化的方式检测并平均出各种大小的过于接近点的聚类? import numpy as np from scipy.spatial i

    1热度

    1回答

    这是我第一次使用scipy.signal库,并且正在用方法filtfilt()试验一个错误。 这是我试图执行代码: Fs = 1000 # s is an array of numbers a=signal.firwin(10, cutoff=0.5/(Fs/2)) ss = s - np.mean(s) se = signal.filtfilt(a, 1, ss, method="gu

    0热度

    1回答

    我有一个非常具体的任务,我需要找到我的指数函数的斜率。 我有两个阵列,一个表示波长范围在400和750 nm之间,另一个表示吸收光谱。 x =波长,y =吸收。 我的拟合函数应该是这个样子的是: y_mod = np.float(a_440) * np.exp(-S*(x - 440.)) 其中S是斜率和图像等于0.016,这应该是在价值观我应该得到的范围(+/- 0.003)。 a_440是在4

    1热度

    1回答

    我有以下脚本来使用odeint来计算dRho。 P_r = 10e5 rho_r = 900 L = 750 H = 10 W = 150 A = H * W V = A * L fi = 0.17 k = 1.2e-13 c = 12.8e-9 mu = 2e-3 N = 50 dV = V/N dx = L/N P_in = P_r rho_in = rho

    1热度

    2回答

    在Matlab中/八度音功能不是,spdiags([-8037.500 50.000 -12.500], 0:2, 1, 51)给出以下输出: (1, 1) -> -8037.5 (1, 2) -> 50 (1, 3) -> -12.500 然而,当我使用在Python下面,它不会产生类似的结果作为在Matlab /八度: ​​ Python的spdiags()产生下面的输出,这是在缺少一

    1热度

    1回答

    如何调整N-D numpy图像的大小? 我不只是想对它进行二次采样,而是对像素进行插值/平均。 例如,如果我开始与 array([[[3, 1, 3, 1], [3, 1, 3, 1], [3, 1, 3, 1], [3, 1, 3, 1]], [[3, 1, 3, 1], [3, 1, 3, 1], [3, 1, 3, 1

    1热度

    1回答

    我想使用scipy和python 3.x三态高斯函数。我认为我真的差不多在那里,但我在这里摸不着头脑,因为我无法弄清楚它出了什么问题。 data =np.loadtxt('mock.txt') my_x=data[:,0] my_y=data[:,1] def gauss(x,mu,sigma,A): return A*np.exp(-(x-m

    1热度

    1回答

    我不能够使用PIP安装SciPy的库并且得到错误 Collecting scipy Using cached scipy-0.19.1.tar.gz Installing collected packages: scipy Running setup.py install for scipy ... error Complete output from comman