Re: [問題] 如何產生非方陣的正交矩陣
※ 引述《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)，所以要小心這件事情
--
R資料整理套件系列文：
magrittr #1LhSWhpH (R_Language) http://tinyurl.com/j3ql84c
data.table #1LhW7Tvj (R_Language) http://tinyurl.com/hr77hrn
dplyr(上) #1LhpJCfB (R_Language) http://tinyurl.com/jtg4hau
dplyr(下) #1Lhw8b-s (R_Language)
tidyr #1Liqls1R (R_Language) http://tinyurl.com/jq3o2g3
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.74.87
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1466151781.A.319.html
※ 編輯: celestialgod (140.109.74.87), 06/17/2016 16:28:12
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 2 篇)：
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
