※ 引述《cy10911 (Ann)》之铭言：
: 如题
: 请问要怎么产生一个非方阵的正交矩阵
: 每个行向量跟别人内积为0
: 自己跟自己内积为1
http://elliottback.com/wp/modified-gram-schmidt-orthogonalization-in-matlab/
B = rand(50, 50);
A = QR(B);
all(all(abs(A * A' - eye(50)) < 1e-6)) % true
这样试试看，基本上不太会出现rank不足的问题 (就确定rank(B) == size(B, 1))
不过随机还是会发生rank(B) < size(B, 1)，所以要小心这件事情