function[]=pauli_cov（c11c22c33c12c13c23str）
%%??使用協(xié)方差矩陣數(shù)據(jù)的?PAULI?分解??（R:HH+VVG:HVB:HH-VV）
if?nargin==6
????str=‘pauli‘;
end

B=[1?0?1;1?0?-1;0?1.414?0]/1.414;%%??原協(xié)方差矩陣中的HV相已經(jīng)有了1.414，所以B中不是2
[s1s2]=size（c11）;
pauli_out=zeros（s1s23）;
for?ii=1:s1
????for?jj=1:s2
????????cov=[c11（iijj）c12（iijj）c13（iijj）;c12（iijj）‘c22（iijj）c23（iijj）;c13（iijj）‘c23（iijj）‘c33（iijj）];
????????t=abs（B*cov*B‘）;
????????if?t（11）>1
????????????pauli_out（iijj1）=0.9999;
????????else
????????????pauli_out（iijj1）=t（11）;
????????end
????????if?t（22）>1
????????????pauli_out（iijj3）=0.9999;
????????else
????????????pauli_out（iijj3）=t（22）;
????????end
????????if?t（33）>1
????????????pauli_out（iijj2）=0.9999;
????????else
????????????pauli_out（iijj2）=t（33）;
????????end
????end
end
figure;imshow（?（pauli_out））;title（str）;
????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????907??2006-01-11?14:37??pauli_cov.m

-----------?---------??----------?-----??----

?????????????????1125????????????????????2