开发者社区

文档建议反馈控制台

文章/答案/技术大牛

发布

社区首页 >问答首页 >SciPy与Numpy伪逆的区别

问SciPy与Numpy伪逆的区别
EN

Stack Overflow用户

提问于 2012-11-07 07:48:06

回答 1查看 10.4K关注 0票数 13

我发现有两个版本的pinv()函数，它在Scipy和numpy中计算矩阵的伪逆，文档可以在以下位置查看：

http://docs.scipy.org/doc/numpy/reference/generated/numpy.linalg.pinv.html

http://docs.scipy.org/doc/scipy/reference/generated/scipy.linalg.pinv.html

问题是我有一个50000*5000的矩阵，当使用scipy.linalg.pinv时，它花费了我超过20 it的内存。但是当我使用numpy.linalg.pinv时，只使用了不到1 1GB的内存。

我想知道为什么numpy和scipy在不同的实现下都有一个pinv。以及为什么他们的表现如此不同。

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2012-11-07 08:39:16

我无法解释为什么在scipy和numpy中都有实现，但我可以解释为什么行为是不同的。

numpy.linalg.pinv使用奇异值分解(精确地说是lapack方法dgesdd )近似摩尔-彭罗斯伪逆，而scipy.linalg.pinv在最小二乘意义下求解模型线性系统以近似伪逆(使用dgelss)。这就是为什么他们的表现是不同的。我希望得到的伪逆估计的总体精度也会有所不同。

您可能会发现scipy.linalg.pinv2的性能更类似于numpy.linalg.pinv，因为它也使用奇异值分解方法，而不是最小平方近似。

票数 17

EN

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/13265299

复制

相关文章

MoorePenrose伪逆

Moore-Penrose 伪逆常用于求解或简化非一致线性方程组的最小范数最小二乘解。其在实数域和复数域上都是唯一的，并且可以通过奇异值分解求得。

hotarugali

2022/03/03

5260

Moore-Penrose伪逆

对于非方阵矩阵而言，其逆矩阵没有定义。假设在下面的问题中。我们希望通过矩阵A的左逆B来求解线性方程：

狼啸风云

2019/09/18

1.8K0

Moore-Penrose伪逆

野生码农的逆袭之路：在跨界中找到自我

大数据数据分析

本文由CDA作者库成员HarryZhu原创，并授权发布。 CDA作者库凝聚原创力量，只做更有价值的分享。 Day Job and Night Job 我非常认同《黑客与画家》里的 Paul Graham 说的一句话：码农需要一个 day job for food，也需要一个 night job for fun。和格雷厄姆不同的是，我的night job不是一个画家，而是一个作家，是的，一个技术专栏的撰稿人。通常，晚餐之后，刷一遍自己的 Feedly 和 GitHub，搞搞黑科技，这就是一种

CDA数据分析师

2018/02/24

1.2K0

野生码农的逆袭之路：在跨界中找到自我

爬虫如何正确从网页中提取伪元素？

css html selenium javascript xslt & xpath

可以看到，在源代码里面没有请抓取我！这段文字。难道这个网页是异步加载？我们现在来看一下网页的请求：

崔庆才

2020/12/15

2.9K0

爬虫如何正确从网页中提取伪元素？

从Java视角理解系统结构（三）伪共享

从我的前一篇博文中, 我们知道了CPU缓存及缓存行的概念, 同时用一个例子说明了编写单线程Java代码时应该注意的问题. 下面我们讨论更为复杂, 而且更符合现实情况的多核编程时将会碰到的问题. 这些问

用户1263954

2018/01/30

6720

从Java视角理解系统结构（三）伪共享

运用伪逆矩阵求最小二乘解

线性回归 scipy 工具函数模型

已经有工具可以解很多最小二乘的模型参数了，但是几个专用的最小二乘方法最多支持一元函数的求解，难以计算多元函数最小二乘解，此时就可以用伪逆矩阵求解了。

为为为什么

2023/04/09

1.8K0

【答疑点评必看】如何从「数据范围」中找到解题「突破口」...

这是 LeetCode 上的「395. 至少有K个重复字符的最长子串」，难度为 Medium。

宫水三叶的刷题日记

2021/03/12

7290

如何从菜鸟程序员成长为(伪)高手

1.摘要最近有一些毕业不久的同事问我：“你工作的时候有没有什么窍门？怎么才能快速成为高手？” 想起当初刚入职，新人培训的时候，也跟其他同事讨论过这个问题：如何才能成为业界大牛？当时自己只是觉得兴趣是最好的老师，思路方法什么的没有多想。加入微博平台架构部的时间也不短了，趁着快过春节总结了一下自己入职微博以来的工作情况，从互联网开发的半个门外汉，到如今能设计一些架构、排查一些问题、分享一些经验，收获颇多，感想颇多,也逐渐意识到思路和方法的重要性，在此跟大家分享一下。主要分为学、做、想三方面。 2.学会学习

小小明童鞋

2018/06/13

8470

LeetCode 85 | 如何从矩阵当中找到数字围成的最大矩形的面积？

今天是LeetCode专题53篇文章，我们一起来看看LeetCode中的85题，Maximal Rectangle（最大面积矩形）。

TechFlow-承志

2020/07/14

1.5K0

如何从40亿个整数中找到不存在的一个

给定一个最多包含40亿个随机排列的32位的顺序整数的顺序文件，找出一个不在文件中的32位整数。（在文件中至少确实一个这样的数-为什么？）。在具有足够内存的情况下，如何解决该问题？如果有几个外部的“临时”文件可用，但是仅有几百字节的内存，又该如何解决该问题？

编程珠玑

2019/09/02

1.6K0

如何在以太坊合并中找到机会？

比特币数字货币区块链

以太坊合并可能导致一场经济上的混乱局面，而在这一系列不确定中，仍存在确定性的机会。

小将

2022/08/11

5430

我是如何从3亿IP中找到CISCO后门路由器的

安全安全漏洞

接到某单位通知让查找中国具有SYNful Knock后门的CISCO路由器，按照曼迪安特分析的报告称中国已经发现3台具有SYNful Knock后门的路由器，如何快速从全国3亿IP地址中快速查找出3个IP地址难度还是十分的大啊，而我经过查找发现中国已经有4个IP被植入了后门，现将检测过程分享给大家。一、获取IP地址为保证中国IP的全面性，从apnic重新获取亚洲区域所分配到的IP，过滤出CN的IP，结果如下。apnic文件中每行为一个IP地址段，以"|"作为分隔，第四个字段为IP起始地址，第五个字

FB客服

2018/02/06

1.7K0

我是如何从3亿IP中找到CISCO后门路由器的

Java 给定字符串中找到数字

public myTest{ public static StringBuilder hhhh(String value) { StringBuilder sb = new StringBuilder(); String strIndex = ""; // String regex = "\\d*"; 可以用String 但最终结果只有后面的数字例子中的只能显示222222222，因为会被替代 Pattern pattern = Pa

桑鱼

2020/04/24

3860

Java 给定字符串中找到数字

从图像中找到文字 -- 开源 OCR 工具 tesseract 使用简介

文字识别 jar java ubuntu https

OCR 是 Optical Character Recognition （光学字符识别）的缩写，指的是通过检测图像，从而识别出文字的技术。

用户3147702

2022/06/27

2.1K0

从图像中找到文字 -- 开源 OCR 工具 tesseract 使用简介

如何在Shopify中找到大量XSS漏洞

电商漏洞扫描服务安全漏洞数据安全安全运营平台

XSS攻击通常指的是通过利用网页开发时留下的漏洞，通过巧妙的方法注入恶意指令代码到网页，使用户加载并执行攻击者恶意制造的网页程序。

小生观察室

2021/08/29

8740

如何在Shopify中找到大量XSS漏洞

如何从有序数组中找到和为指定值的两个元素下标

例如:{2, 7, 17, 26, 27, 31, 41, 42, 55, 80} target=72.求得值为17和55,对应下标为:2,8

一个架构师

2022/06/20

2.3K0

如何从有序数组中找到和为指定值的两个元素下标

一日一技：如何从多个Jupyter Notebook中找到需要代码段

开源编程算法

某一天，我突然想找很久之前写的一段代码片段，但是我不记得是写在哪个Notebook里面了。应该怎么办呢？

青南

2019/08/05

1.9K0

java-逆波兰算法

四则运算是栈的重要应用之一中缀表达式转后缀表达式（逆波兰算法）过程从左到右遍历中缀表达式数字直接输出为后缀表达式一部分如果是符号，则判断与栈顶元素的优先级高于栈顶元素优先级直接入栈低于或等于栈顶优先级栈顶元素出栈并输出为后缀表达式一部分（注意这里是递归比较栈顶元素的优先级并出栈），最后将当前元素入栈直到遍历完中缀表达式，最终输出后缀表达式下面是自己的实现源码 package com.yhq.demospringboot; import org.apache.commons.lang3.St

IT架构圈

2018/06/01

7340

从微盟“删库”事件来谈谈企业的信息安全伪壁垒

运维网站机器人数据库 sql

背景：2月23日19点，我们收到系统监控报警，服务出现故障，随后我们立刻召集相关技术人员进行定位，发现大面积服务集群无法响应，生产环境及数据遭受严重破坏。我们立刻启动紧急响应机制，并与腾讯云技术团队一起研究制定生产环境和数据修复方案。

顾黄亮

2020/06/25

9510

从微盟“删库”事件来谈谈企业的信息安全伪壁垒

[教程]Typecho 如何开启伪静态？

如果在保存上述配置的时候，typecho 无法自动配置，那么你可能需要手动配置服务器的 rewrite 规则。 nginx 配置. conf：

字节星球Henry

2021/08/09

5.7K0

相似问题

阴影问题浮动操作按钮(api = 19)

362

用于扩展的Xamarin.Android浮动操作按钮的阴影填充错误与AndroidSDK API 19及更少

18

安卓浮动动作按钮API19 (KitKat)

20

Android浮动操作按钮为空

10

Android带有API级别19的RecyclerView (API 19)

25

添加站长进交流群

领取专属 10元无门槛券

AI混元助手 在线答疑

关注 腾讯云开发者公众号

洞察 腾讯核心技术

剖析业界实践案例