%freeman分解
clc
clear?all
close?all
%load?E:\SanFrancisco\T.mat
load（‘T‘）;
row?=900;?????
col?=1024;?
for?i=1:row
????for?j=1:col
????????t?=?[T（ij1）?T（ij2）?T（ij3）;T（ij4）?T（ij5）?T（ij6）;?T（ij7）?T（ij8）?T（ij9）];??
????????fv=（3/2）*t（22）;
????????pv=（8/3）*fv;?????%求pv
????????if?real（t（13））>0
????????????alph=-1;
????????????beta=1+（t（11）-t（33））/（t（13）+t（33）-2*t（22））;
????????????be=（t（11）+t（13）-2*t（22））/（t（13）+t（33）-2*t（22））;

%%求ps和pd

%????????????fs=（t（11）-t（33））/（（abs（beta））^2-1）;?不好
????????????fs=（t（13）+t（33）-2*t（22））/（1+beta）;
????????????ps=fs*（1+（abs（beta））^2）;
????????????fd=t（33）-fs-fv;
????????????pd=2*fd;
????????else
??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1424??2013-03-27?11:06??freeman.m