草庐IT

matrix_size

全部标签

python - scipy.stats.multivariate_normal 提高 `LinAlgError: singular matrix` 即使我的协方差矩阵是可逆的

我在尝试使用scipy.stats.multivariate_normal时遇到问题,希望你们中的某个人能够提供帮助。我有一个2x2矩阵,可以找到使用numpy.linalg.inv()的逆矩阵,但是当我尝试将其用作multivariate_normal中的协方差矩阵时我收到LinAlgError声明它是一个奇异矩阵:In[89]:cov=np.array([[3.2e5**2,3.2e5*0.103*-0.459],[3.2e5*0.103*-0.459,0.103**2]])In[90]:np.linalg.inv(cov)Out[90]:array([[1.23722158e-1

python - urllib2.urlopen() : getting the size of the content

只要工作允许,我仍然会围绕python工作......我正在使用使用urllib2.urlopen的脚本查询大量内部webUI。我想知道如何从每个请求中获取页面内容的大小。我似乎无法弄清楚这一点。提前致谢MHibbin 最佳答案 printlen(urlopen(url).read())或>>>result=urllib2.urlopen('http://www.spiegel.de')>>>result.headers['content-length']'181291' 关于pytho

python - NumPy 错误 : Singular matrix

错误Numpyerror:Matrixissingular具体是什么意思(使用linalg.solve函数时)?我在Google上查看过,但找不到任何可以说明此错误发生时间的信息。 最佳答案 奇异矩阵是不可逆的。这意味着您要求解的方程组没有唯一解;linalg.solve无法处理这个问题。您可能会发现linalg.lstsq提供了一个可用的解决方案。 关于python-NumPy错误:Singularmatrix,我们在StackOverflow上找到一个类似的问题:

python - 索引错误 : index 1 is out of bounds for axis 0 with size 1/ForwardEuler

我正在对一阶微分方程组的x(t)进行数值求解。该系统是:dy/dt=(C)\*[(-K\*x)+M*A]我已经实现了正向欧拉方法来解决这个问题,如下所示:这是我的代码:importmatplotlibimportnumpyasnpfromnumpyimport*fromnumpyimportlinspacefrommatplotlibimportpyplotaspltC=3K=5M=2A=5#------------------------------------------------------------------------------defeuler(f,x0,t):n=l

python - 访问 coo_matrix 中的元素

这是一个非常简单的问题。对于像coo_matrix这样的SciPy稀疏矩阵,如何访问单个元素?类比本征线性代数库。可以使用coeffRef访问元素(i,j),如下所示:myMatrix.coeffRef(i,j) 最佳答案 来自coo_matrix的文档:|IntendedUsage|-COOisafastformatforconstructingsparsematrices|-Onceamatrixhasbeenconstructed,converttoCSRor|CSCformatforfastarithmeticandmatr

Python os.stat(file_name).st_size 与 os.path.getsize(file_name)

我有两段代码,它们都是为了做同样的事情——坐在一个循环中,直到一个文件被写入完成。它们都主要用于通过FTP/SCP传入的文件。代码的一个版本使用os.stat()[stat.ST_SIZE]:size1,size2=1,0whilesize1!=size2:size1=os.stat(file_name)[stat.ST_SIZE]time.sleep(300)size2=os.stat(file_name)[stat.ST_SIZE]另一个版本使用os.path.getsize():size1,size2=0,0whileTrue:size2=os.path.getsize(file

python - numpy matrix trickery - 逆时矩阵之和

我正在尝试执行以下操作,并重复直到收敛:其中每个Xi是nxp,还有r他们中的一个rxnxp名为samples的数组.U是nxn,V是pxp.(我得到了matrixnormaldistribution的MLE。)尺寸都可能很大;我期待的事情至少在r=200的顺序上,n=1000,p=1000.我当前的代码可以V=np.einsum('aji,jk,akl->il',samples,np.linalg.inv(U)/(r*n),samples)U=np.einsum('aij,jk,alk->il',samples,np.linalg.inv(V)/(r*p),samples)这没问题,但

python - Apache Spark : How to create a matrix from a DataFrame?

我在ApacheSpark中有一个带有整数数组的DataFrame,源是一组图像。我最终想对其进行PCA,但我无法从我的数组创建矩阵。如何从RDD创建矩阵?>imagerdd=traindf.map(lambdarow:map(float,row.image))>mat=DenseMatrix(numRows=206456,numCols=10,values=imagerdd)Traceback(mostrecentcalllast):File"",line2,inmat=DenseMatrix(numRows=206456,numCols=10,values=imagerdd)Fil

python ,Scipy : Building triplets using large adjacency matrix

我正在使用邻接矩阵来表示可以在视觉上解释为的friend网络Mary0111Joe1011Bob1101Susan1110MaryJoeBobSusan使用这个矩阵,我想编译所有可能的友谊三角列表,条件是用户1是用户2的friend,用户2是用户3的friend。对于我的列表,用户1不需要是用户3的friend。(joe,mary,bob)(joe,mary,susan)(bob,mary,susan)(bob,joe,susan)我有一些代码可以很好地处理小三角形,但我需要它来缩放非常大的稀疏矩阵。fromnumpyimport*fromscipyimport*defbuildTri

python - 将稀疏矩阵 (csc_matrix) 转换为 pandas 数据帧

我想将这个矩阵转换为Pandas数据框。csc_matrix括号中的第一个数字应该是索引,第二个数字是列和最后的数字是数据。我想这样做是为了在文本分析中进行特征选择,第一个数字代表文档,第二个数字是单词的特征,最后一个数字是TFIDF分数。获取数据框帮助我将文本分析问题转化为数据分析。 最佳答案 fromscipy.sparseimportcsc_matrixcsc=csc_matrix(np.array([[0,0,4,0,0,0],[1,0,0,0,2,0],[2,0,0,1,0,0],[0,0,0,0,0,1],[4,0,3,