1回答
作者多次强调，用正交基计算最小二乘解非常容易，快速，因为Qᵀ*Q = I，其中Q是一个具有正交基的设计矩阵。所以你的方程变成了x̂ = Qᵀb。我得到的印象是，在应用最小二乘之前，每次计算QR分解都是个好主意。但后来我发现了QR分解的，结果是计算QR分解，然后应用最小二乘比普通的x̂ = inv(AᵀA)Aᵀb更昂贵。
说的对吗
浏览 3提问于2021-11-05得票数 1
回答已采纳
我一直在查看函数numpy.linalg.lstsq，它包含了LAPACK提供的所有类型的求解器，除了基于QR分解的GLES例程(实际上，持家反射器因为用Q计算QR是浪费计算时间)。而且，在我在网络上看到的基于QR的最小二乘的所有实现中，它们计算QR分解中的Q，然后使用它，但是它需要低效率，因为计算Q是一种浪费，我们只需要反射。LAPCAK的网站解释了最小二乘求解器的类型：，如你所
浏览 3提问于2020-08-13得票数 0
回答已采纳
4回答
它是在做最小二乘计算，还是有迭代的方式？编辑:我认为梯度下降是可行的。最小二乘计算可能会有点繁琐(必须做qr分解或其他事情来保持稳定)。
浏览 2提问于2010-10-28得票数 6
回答已采纳
因此，我试图找出如何使用R中的qr.solve命令来根据特定的数据集求解最佳匹配线(斜率和截距)的组件。虽然我已经知道如何使用lm函数来完成这个任务，但据我所知，矩阵的qr分解可能产生同样的结果。qr.solve函数的输入可以在R中实现，还是必须先自己求解，然后再插入R中呢？
我尝试将我的数据集(每行2列表示图上的点)输入到一个矩阵中，并在函数qr.solve中使用它作为x参数。
浏览 6提问于2022-11-26得票数 1
回答已采纳
1回答
QR分解
、、、
有没有像那样实现QR分解的方法？特别是，我对以下命令感兴趣：根据它的描述“返回一个排列矩阵P，它被选择来减少填充R，你可以使用C，R，和P来计算稀疏线性系统SX =B的最小二乘解，其中X= P(R我需要运行相当于MATLAB代码[Q, R, E] = qr(X,0);的Python。
浏览 9提问于2022-11-11得票数 0
nnmf仅允许用于r< min(size(A))
有没有办法把矩阵分解成两个矩阵？
浏览 7提问于2012-03-07得票数 1
我在usconsumption数据集上使用R中的forecast包拟合Arima(2,0,0)模型。然而，当我使用lm模拟相同的拟合时，我得到了不同的系数。我的理解是它们应该是一样的。下面是我的代码。> library(fpp)> #load data> > lagpad <- function(x, k=1) {
+
浏览 0提问于2016-07-12得票数 3
4回答
我正在寻找一个普通最小二乘的numpy-based实现，可以用更多的观测值更新拟合。类似于或R的的东西。有什么建议吗？
浏览 8提问于2011-05-10得票数 3
回答已采纳
我正在研究对一个或多个自变量进行线性回归的算法。C + D*am + E* bm = ym我将使用矩阵表示法其中Beta是向量C，D，E，其中这些值将是最佳拟合线。我是否应该计算出
或者我应该使用矩阵的LU分解/分解。在大量的数据上，每个数据的性能如何(即m的一个大值，按10^8 .)编辑
如果答案是使用Cholesky分解或QR分解，那么是否有实现提示/简
浏览 0提问于2013-10-21得票数 4
回答已采纳
1回答
我正在研究一种算法来查找列表对象中的峰值。我想出了一个很好(或者足够好)的算法，通过查看一个点和它的邻居，如果它是一个峰值，将它添加到结果列表中。我对LabView如何做到这一点做了一些研究，发现如下：
“这个峰值检测器VI基于一种算法，该算法将二次多项式拟合成数据点的顺序组，其中使用的数据点的数目由宽度指定。对于每个峰或谷，二次拟合是根据阈值进行测试的。峰高低于阈值，谷谷高于阈值则被忽略。只有在VI处理超过峰值或山谷位置的大约宽度/2个数据点之后才能检测到峰和谷。这种延迟只
浏览 5提问于2011-03-14得票数 1
例如，一种策略是计算和存储一个适当的矩阵分解(LUP或QR)，否则每次线性代数库都会调用它(内部)。但是，我几乎找不到在中间求解器(特别是在使用LAPACK时)使用哪种方法的指导，或者分解在速度或遇到坏边缘情况下的相对优点。
注意:由于速度和精度问题，“坏”策略是计算逆矩阵inv(A)。
浏览 0提问于2016-06-01得票数 2
回答已采纳
Apache mahout在没有Hadoop的情况下工作吗?如果不是这样的话，那么mahout的哪些部分依赖于Hadoop。我正在尝试Mahout集群实现。
浏览 4提问于2014-06-11得票数 0
回答已采纳
1回答
我用Matlab来估计一个具有普通最小二乘(OLS)的回归模型。[C,R] = qr(x,y,0)b = R\C
我的问题是，我是否需要担心这里的数字错误。我还需要做一些额外的迭代吗？我应该检查R或R'R的条件号吗？
浏览 3提问于2014-04-22得票数 2
1回答
我想问一下$coefficients和$effects在aov输出方面的区别。让我们考虑以下简单的数据集。f2 <- c(1,0,0,1,1,0,0,1)av <- aov(r ~ f1 * f2)av$coefficientsplot(f1, r)
由于f1 mean(r[f1=
浏览 1提问于2016-10-24得票数 5
回答已采纳
2回答
我需要在一个列表中solve超过上千个矩阵。但是，我得到了错误Lapack routine dgesv: system is exactly singular。我的问题是，我的输入数据是非奇异矩阵，但是在对这些矩阵进行计算之后，其中一些得到了奇异。但是，一个包含数据子集的可重复的示例是不可能的，因为它会很长(我已经尝试过了)。这里是我问题的一个基本例子(A是经过一些计算后的矩阵，R是我需要做的下一个计算)：R = solve(diag(4)-A)
你有什么想法如何“解决”这个问题，也许其
浏览 3提问于2015-11-06得票数 0
1回答
我正在编写一个用ruby编写的，并且我遇到了一个算法需要使用的地方。我正在寻找一个宝石来转换一个矩阵，以减少行梯形形式。
浏览 0提问于2013-02-14得票数 1
有人知道Armadillo库是如何找到具有非可逆A矩阵(Ax=b)的线性系统的近似解的吗？考虑的函数是解(x，A，b)。
浏览 0提问于2018-05-01得票数 2
回答已采纳
我有一个有点复杂的算法，需要将二次曲线拟合到一组点上。这个二次曲面是由它的参数化(u, v, f(u,v))给出的，其中f(u,v) = au^2+bv^2+cuv+du+ev+f。除了计算奇异值分解外，这个库中是否有更直接的算法可以在最小二乘意义上解决系统(同样，除了奇异值分解或使用朴素的inv(transpose(A)*A)*transpose(A)*b公式解决一般的最小二乘问
浏览 3提问于2013-05-18得票数 0
回答已采纳
1回答
我已经试验了两种实现最小二乘拟合算法的方法，如所示。V = [x1,1; x2,1; x3,1; ... xn,1] % this
浏览 6提问于2010-12-02得票数 2
回答已采纳
2回答
Lapack via方案
、、、、
这至少包括所有的BLAS操作，主要的分解(SVD，LU，QR，Cholesky)，为了方便，还包括最小二乘。我可以看到Chicken和Chicken有一些东西，但它没有覆盖上面的内容。
浏览 2提问于2011-10-28得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
云直播
腾讯会议活动推荐
腾讯云开发者
