matlab-线性代数 判断 合同矩阵

原创
2019/02/21 20:41
阅读数 841

     matlab : R2018a 64bit
      OS : Windows 10 x64
typesetting : Markdown
       blog : my.oschina.net/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,用于仿真和基于模型的设计,值得学习。
该博文仅可用于测试与参考。

展开阅读全文
加载中
点击引领话题📣 发布并加入讨论🔥
打赏
0 评论
0 收藏
0
分享
返回顶部
顶部