%?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;
?