有時(shí)候做實(shí)驗(yàn)需要使用對(duì)稱(chēng)矩陣，這里介紹如何使用Matlab產(chǎn)生隨機(jī)的對(duì)稱(chēng)矩陣。

用例子說(shuō)明一下：我要產(chǎn)生4X4的隨機(jī)矩陣，要求是對(duì)稱(chēng)矩陣。

產(chǎn)生對(duì)稱(chēng)矩陣

A?=?rand(4);
B?=?tril(A,-1)+triu(A',0);

方法一：B是A的下三角矩陣改造的對(duì)稱(chēng)矩陣

A?=?rand(4);
B?=?triu(A,0)?+?tril(A',-1);

方法二：B是A的上三角矩陣改造的對(duì)稱(chēng)矩陣

驗(yàn)證對(duì)稱(chēng)矩陣

只需要這樣驗(yàn)證：

sum(sum(B==B'))==numel(B)

只要上面的代碼輸出結(jié)果是1就對(duì)了

附帶一個(gè)5X5的實(shí)對(duì)稱(chēng)正定矩陣

????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