function yamachuchi_matrix = yamaguchi_test(T,OA)
%对于给定的T矩阵做yamaguchi分解,返回结果是个各个散射的能量是一个row X col X 4的矩阵
%P_s(单词散射),P_d(二面角),P_v(体散射),P_c(旋状体)
row=size(T.T11,1);
col=size(T.T11,2);
i=sqrt(-1);
yamachuchi_matrix = zeros(row,col,4);
for p=1:col
for q=1:row
tmp(1,1) = T.T11(q,p);
tmp(1,2) = T.T12_real(q,p)+i*T.T12_imag(q,p);
tmp(1,3) = T.T13_real(q,p)+i*T.T13_imag(q,p);
tmp(2,1) = conj(tmp(1,2));
tmp(2,2) = T.T22(q,p);
tmp(2,3) = T.T23_real(q,p)+i*T.T23_imag(q,p);
tmp(3,1) = conj(tmp(1,3));
tmp(3,2) = conj(tmp(2,3));
tmp(3,3) = T.T33(q,p);
P_c = abs(2*imag(tmp(2,3))/i);
alfa = tm