Pseudo-aleatorio ortogonal de la matriz de generación

Question

Un poco de tiempo, J. M. sugirió el uso de la Stewart algoritmo para la generación de $n$ $n$ pseudo-aleatorios ortogonal de matrices en $O(n^2)$ del tiempo. Además, señaló que esta metodología se implementa en Nick Higham la matriz de cálculo (caja de herramientasde matlab). Ahora este paquete contiene un montón de .m archivos. Puede que nadie me dirija a lo que uno se implementa la pseudoaleatoria ortogonal de la matriz generador? Alternativamente, podría alguien sabe de un R la implementación de este algoritmo?

Answer 1

Es en la Matriz de Prueba caja de herramientas, no de la Matriz de Cálculo de herramientas. El M-archivo qmult.m (premultiplication por un Haar-distribuido pseudoaleatoria ortogonal de la matriz) se puede encontrar aquí o aquí.