极化目标分解中的四分量FOURCOMPONENT分解

% Function : Four-component decomposition
% Four-component分解把C矩阵分解为C=fv*Cv+fd*Cd+fs*Cs+fh*Ch;
clear;clc;
fid_c11=fopen（‘c:\C3\C11.bin‘‘r‘）;
fid_c22=fopen（‘c:\C3\C22.bin‘‘r‘）;
fid_c33=fopen（‘c:\C3\C33.bin‘‘r‘）;
fid_c13real=fopen（‘c:\C3\C13_real.bin‘‘r‘）;
fid_c13img=fopen（‘c:\C3\C13_imag.bin‘‘r‘）;
fid_c12img=fopen（‘c:\C3\C12_imag.bin‘‘r‘）;
fid_c23img=fopen（‘c:\C3\C23_imag.bin‘‘r‘）;
Crow=1024;
Ccol=900;
c11=fread（fid_c11[Crow Ccol]‘float‘）;
c22=fread（fid_c22[Crow Ccol]‘float‘）;
c33=fread（fid_c33[Crow Ccol]‘float‘）;
c13_real=fread（fid_c13real[Crow Ccol]‘float‘）;
c13_img=fread（fid_c13img[Crow Ccol]‘float‘）;
c12_img=fread（fid_c12img[Crow Ccol]‘float‘）;
c23_img=fread（fid_c23img[Crow Ccol]‘float‘）;
c11 = c11‘;
c22 = c22‘;
c33 = c33‘;
c13 = （c13_real + 1i*c13_img）‘;
c13_real=c13_real‘;
c12_img=c12_img‘;
c23_img=c23_img‘;
fclose（fid_c11）;
fclose（fid_c22）;
fclose（fid_c33）;
fclose（fid_c13real）;
fclose（fid_c13img）;
fclose（fid_c12img）;
fclose（fid_c23img）;
%-----------------------------------------------------
pv=zeros（CcolCrow）; 
ps=zeros（CcolCrow）;  
pd=zeros（CcolCrow）;
ph=zeros（CcolCrow）;
span1=zeros（CcolCrow）;
span2=zeros（CcolCrow）;
for ii=1 : 900
    for jj=1:1024
        
        R=10*log（c33（iijj）/c11（iijj））;
        
        if R<-2
            
            fh=2*（c12_img（iijj）+c23_img（iijj））;
            fv=（15/2）*（c22（iijj）/2-fh/4）;
            a=c11（iijj）-3*fv/15-fh/4;
            b=c33（iijj）-8*fv/15-fh/4;
            c=c13（iijj）-2*fv/15+fh/4;
            
        elseif R>2
            
            fh=2*（c12_img（iijj）+c23_img（iijj））;
            fv=（15/2）*（c22（iijj）/2-fh/4）;
            a=c11（iijj）-8*fv/15-fh/4;
            b=c33（iijj）-3*fv/15-fh/4;