极化SAR目标分解中的cloude分解程序

% Function : Cloude decomposition
clear;clc;
fid_T11=fopen（‘c:\T\T11.bin‘‘r‘）;
fid_T22=fopen（‘c:\T\T22.bin‘‘r‘）;
fid_T33=fopen（‘c:\T\T33.bin‘‘r‘）;
fid_T12real=fopen（‘c:\T\T12_real.bin‘‘r‘）;
fid_T12imag=fopen（‘c:\T\T12_imag.bin‘‘r‘）;
fid_T13real=fopen（‘c:\T\T13_real.bin‘‘r‘）;
fid_T13imag=fopen（‘c:\T\T13_imag.bin‘‘r‘）;
fid_T23real=fopen（‘c:\T\T23_real.bin‘‘r‘）;
fid_T23imag=fopen（‘c:\T\T23_imag.bin‘‘r‘）;
Crow=1024;
Ccol=900;
T11=fread（fid_T11[Crow Ccol]‘float‘）;
T22=fread（fid_T22[Crow Ccol]‘float‘）;
T33=fread（fid_T33[Crow Ccol]‘float‘）;
T12_real=fread（fid_T12real[Crow Ccol]‘float‘）;
T12_imag=fread（fid_T12imag[Crow Ccol]‘float‘）;
T13_real=fread（fid_T13real[Crow Ccol]‘float‘）;
T13_imag=fread（fid_T13imag[Crow Ccol]‘float‘）;
T23_real=fread（fid_T23real[Crow Ccol]‘float‘）;
T23_imag=fread（fid_T23imag[Crow Ccol]‘float‘）;
T11 = T11‘;
T22=T22‘;
T33=T33‘;
T12_real=T12_real‘;
T12_imag=T12_imag‘;
T13_real=T13_real‘;
T13_imag=T13_imag‘;
T23_real=T23_real‘;
T23_imag=T23_imag‘;
fclose（fid_T11）;fclose（fid_T22）;fclose（fid_T33）;
fclose（fid_T12real）;
fclose（fid_T12imag）;
fclose（fid_T13real）;
fclose（fid_T13imag）;
fclose（fid_T23real）;
fclose（fid_T23imag）;
%----------------------------------------------
A=zeros（CcolCrow）; 
H=zeros（CcolCrow）;  
alpha=zeros（CcolCrow）;
for