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

镇场诗：慈心积善融学习，技术誓为大家学。善心速造多好事，前人栽树后乘凉。我今于此写经验，愿见文者得启发。
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

%输入需要判定的矩阵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;endendif(flagg==1)%计算两个矩阵的特征值并排序e=eig(a);ee=sort(e);f=eig(b);ff=sort(f);eeffif(ee==ff)        fprintf('两矩阵合同.\n');else    fprintf('两矩阵不合同.\n');            endend

           
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
感恩帮助过我的人。博客的精髓在技术部分，更在镇场一诗。
我是一个新手，代码还有许多不完善的地方，请您看代码的时候多多思考。
MATLAB是一个优秀的软件，值得学习。如果您有一些不会的知识，咱们可以相互讨论。
如果您认为代码可以有改进的地方，有错误的地方，请留下评论，我会处理的。
注：如果我的博文无意中侵犯了您的权益，请告知。看到您的告知后，我将及时作出处理。