matlab-線性代數 判斷 合同矩陣

     matlab : R2018a 64bit
      OS : Windows 10 x64
typesetting : Markdown
       blog : my.oschina.net/zhichengjiu
      gitee : gitee.com/zhichengjiu

 

code

clear
clc

% 輸入須要斷定的矩陣
 
a=[1 2 3;2 1 0;3 0 1]
b=[1 2 3;8 9 8;8 7 6]
 
c=size(a);
d=size(b);
 
if(c==d)
    flag=1;
    fprintf('兩矩陣同階\n');
else
    flag=0;
    fprintf('兩矩陣不合同\n');
end
 
% 矩陣中有複數，進行轉置，取共軛.
% 因此不管實數,複數不同。  只有實數去共軛纔是它自己。
if(flag==1)
    if((a==a')&(b==b')) 
        fprintf('該矩陣是實對稱矩陣\n');
        flagg=1;
    else
        fprintf('兩矩陣不合同\n');
        flagg=0;
    end
end

if(flagg==1)
% 計算兩個矩陣的特徵值並排序
e=eig(a);
ee=sort(e);
 
f=eig(b);
ff=sort(f);
 
ee
ff
    if(ee==ff)

        fprintf('兩矩陣合同.\n');
    else
        fprintf('兩矩陣不合同.\n');
    end
end

 

result

a =

     1     2     3
     2     1     0
     3     0     1


b =

     1     2     3
     8     9     8
     8     7     6

兩矩陣同階
兩矩陣不合同
>>

 

resource

• [文檔] ww2.mathworks.cn/help/matlab
• [文檔] ww2.mathworks.cn/help/simulink
• [平臺] www.oschina.net
• [平臺] gitee.com

 

---

感謝幫助 志成就 的人們。
matlab優秀，值得學習。基礎知識 + 專業知識 + matlab = ?
Simulink，用於仿真和基於模型的設計，值得學習。
該博文僅可用於測試與參考。