(R) 生成半正交随机矩阵答案

【问题标题】：(R) Generate a semi-orthogonal random matrix(R) 生成半正交随机矩阵
【发布时间】：2020-05-07 04:30:48
【问题描述】：

我需要生成一个随机 NxK 矩阵（其中 N > K），其中 K 列是正交随机向量。

我尝试的一个选项是生成一个大小为 NxN 的平方正交矩阵，然后选择前 K 列，但我想知道是否有更有效的方法。

当前代码（在 R 中）：

library(pracma)
Z <- rortho(N)[,1:K]

【问题讨论】：

请阅读有关how to ask a good question 的信息以及如何提供reproducible example。这将使其他人更容易帮助您。另外，添加您现有的代码。

【解决方案1】：

您可以使用 rstiefel 包做到这一点：

library(rstiefel)
rustiefel(N, K)

【讨论】：