0

好的,这是一个更新。您将拥有与我相同的输入矩阵。

我从以下矩阵 X 开始:

X <- matrix(0,10,25)
X[1,6] <- 1
X[2,11] <- 1
X[3,16] <- 1
X[4,21] <- 1
X[5,12] <- 1
X[6,17] <- 1
X[7,22] <- 1
X[8,18] <- 1
X[9,23] <- 1
X[10,21] <- -1.46114
X[10,22] <- 0.005024144
X[10,23] <- 0.005329883
X[10,24] <- 1.050356
X[10,25] <- 0.09473591

这个矩阵应该被转换成一个矩阵 Y,即(25x15)。由于这个很大,我只给 ua 几个值(右下角应该是这样的):

Y[24,14] = 0.0080693
Y[25,14] = -0.0894663
Y[24,15] = -0.0894663
Y[25,15] = 0.9919307

这种转换描述如下:

它使用奇异值分解来获取对应于零奇异值的列。

我试图找到解决方案svd()但没有成功。在 RATS 中,这种精确的转换由过程管理%perp。一个 Estima 的人向我解释说,这个过程“使用奇异值分解来获取对应于零奇异值的列”。

希望提供的数据有助于找到解决方案。

马丁

4

0 回答 0