SAR雷达Yamaguchi四分量分解MATLAB代码（注释部分为极化角优化）

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