matlab对矩阵谱分解

syms p real; %定义实数p
A = [ ]; %输入一个矩阵A
[V,D] = eig(A) %A的特征值为D的对角线元素，V的列向量为其对应的特征向量
simplify(VDinv(V)); %通过特征值和特征向量求出矩阵
O = orth(V); %将特征向量正交规范化
[la,lb,lc,ld] = deal(D(1,1),D(2,2),D(3,3),D(4,4)); %将特征值赋值给la，lab，lc，ld
[a,b,c,d] = deal(O(:,1),O(:,2),O(:,3),O(:,4)); %将正交规范化后的特征向量赋值给a，b，c，d
simplify(laaa’ + lbbb’ + lccc’ + lddd’) %a’为a的共轭转置，|a>’=<a|