[求助]耦合矩阵变换问题

根据S曲线提取耦合矩阵后要经过相似变换为对应拓扑的耦合矩阵，Matlab里实现老是变不对，大佬们帮忙给看看吧 @BW~A@8  
拓扑结构如下： 'Ur$jW  
1    1    0    0    0    0 5p=T*Y  
0    1    1    1    0    0 z4{|?0=C  
0    0    1    1    0    0 !\?? [1_e  
0    0    0    1    1    0 G'{4ec0<{  
0    0    0    0    1    1 s)ZL`S?</  
0    0    0    0    0    1 J,+| Fb  
原始M矩阵如下： }ZvL%4jT  
0.0315   -0.6036   -0.2204   -0.2691    0.5697   -0.0000 i'`Z$3EF)  
   -0.6036   -0.4937   -0.0419    0.6875   -0.0616   -0.5517 kCN9`9XI{  
   -0.2204   -0.0419    0.1524    0.0100    0.7249    0.1892 9.1%T06$  
   -0.2691    0.6875    0.0100   -0.0887    0.0253   -0.2951 +uv]d
D*i  
    0.5697   -0.0616    0.7249    0.0253    0.4040   -0.6507 -] J V  
   -0.0000   -0.5517    0.1892   -0.2951   -0.6507    0.0728 ,pf<"^li  
变换后的拓扑和设置的不一样，如下： K[T?--H  
0.0315   -0.8999         0         0         0         0 "MK:y[+*  
   -0.8999    0.0225    0.6291   -0.0494         0         0 ;znIY&Z  
         0    0.6291    0.0959    0.5725         0   -0.1154 V/762&2X  
         0   -0.0494    0.5725    0.0853   -0.6802    0.199 #W!@j"8eK  
         0         0         0   -0.6802   -0.2296   -0.8932 .D)'ZY  
         0         0   -0.1154    0.1990   -0.8932    0.0728 |3"'>* J  
变换程序如下： Ej6vGC.,  
Ru(1,:)=[1 1 0 0 0 0]; 7c+TS--  
Ru(2,:)=[0 1 1 1 0 0]; kS=OX5  
Ru(3,:)=[0 0 1 1 0 0]; wCn W]<+  
Ru(4,:)=[0 0 0 1 1 0]; ">CjnF2>R  
Ru(5,:)=[0 0 0 0 1 1]; HRPTP+  
Ru(6,:)=[0 0 0 0 0 1]; 4Hk
Og)a  
nuKjp Ap!  
Rd=tril(ones(N));%取下三角矩阵，其余为0； Z4E:Z}~''  
RU=Ru+Rd FS1<f:  
Ro=RU; 10}\7p8  
row=1; Bv!j.$0d{  
colum=0; &#;UKk~)Of  
Mr=M; }@a_x,O/x}  
      for a=1:N bnUd !/;  
        if row>colum%先行开始消 m./lrz  
            for i=N:-1:1 ;'R{b$B;|  
                if Ro(row,i)==0;%则要消去 hQ8/-#LO_  
              pivot=[i-1,i]; ^U:pv0Qz  
              seta=-atan(Mr(row,i)/Mr(row,i-1)); MH!'g7iK8  
              Rt=eye(N); {T m-X`  
              Rt(i-1,i)=-sin(seta); S G]e^%i  
              Rt(i,i-1)=sin(seta); k\x>kJ}0  
              Rt(i-1,i-1)=cos(seta); rf]]I#C7  
              Rt(i,i)=cos(seta); $Wb"X=}tl  
              Mr=Rt*Mr*Rt'; z#bOFVg#  
              Ro(row,i)=1;%完成行消元 3Yg/-=U(  
                else ]hL 1qS  
                end d t0?4 d  
            end oDM}h +  
            colum=colum+1; XP6R$0yN  
        else ]}KmT"vA  
            for i=1:N %再列开始消 JBi<TDm/  
                if Ro(i,N+1-colum)==0;%则要消去 H.w
p{m{  
                    pivot=[i,i+1]; GaM#a[p  
                    seta=atan(Mr(i,N+1-colum)/Mr(i+1,N+1-colum)); 'iUfr@  
                    Rt=eye(N); hD # Yz<  
              Rt(i+1,i)=sin(seta); K~v"%sG{`  
              Rt(i,i+1)=-sin(seta); M<g>z6   
              Rt(i+1,i+1)=cos(seta); >8>`-  
              Rt(i,i)=cos(seta); 9s5CqB  
             .. VJZ   
ldi'@
^  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

有关于s参数曲线提取耦合矩阵的资料吗共享出来一起看看啊