有時候做實驗需要使用對稱矩陣,這里介紹如何使用Matlab產(chǎn)生隨機的對稱矩陣。
用例子說明一下:我要產(chǎn)生4X4的隨機矩陣,要求是對稱矩陣。
產(chǎn)生對稱矩陣
A?=?rand(4); B?=?tril(A,-1)+triu(A',0);
方法一:B是A的下三角矩陣改造的對稱矩陣
A?=?rand(4); B?=?triu(A,0)?+?tril(A',-1);
方法二:B是A的上三角矩陣改造的對稱矩陣
驗證對稱矩陣
只需要這樣驗證:
sum(sum(B==B'))==numel(B)
只要上面的代碼輸出結(jié)果是1就對了
附帶一個5X5的實對稱正定矩陣
????2.1012,0.1300,-1.6081,-1.1935,0.3851, ????0.1300,0.6209,-0.2666,-0.3431,-0.5251, ????-1.6081,-0.2666,2.6402,1.0969,-0.5009, ????-1.1935,-0.3431,1.0969,3.3753,-0.6894, ????0.3851,-0.5251,-0.5009,-0.6894,1.5310