%4階統計量方法
clear?all;
filename1=[‘E:\電影\課堂作業\神經網絡\NN\風景圖\3.bmp‘];
???imgdata=double（imread（filename1‘bmp‘））;????????????????%%%讀取圖像3
????????????????
[hw]=size（imgdata）;??????????????????????%%%圖像的大小為[hw]

?len=100000;??????????????????????????????????%%%%選取10萬個樣本
?X=zeros（256len）;
for?i=1:len???????????????????????????????%%%%
????x=floor（rand（）*（w-16））+1;
????y=floor（rand（）*（h-16））+1;
????TI=imgdata（y:（y+15）x:（x+15））;?%%%從圖中隨機選取，小圖的大小?16*16???
????X（:i）=TI（:）;????????????%%將每個小圖整形成[2561]的數據；
end?????????????????????????????????????
?

X=X-mean（X2）*ones（1len）;
?Cx=X*X‘/len;????????????????????????????%%樣本去均值??????
[U?V]=svd（Cx）;??????????????????????????%%%
?epsilon=1e-10;
n=length（find（diag（V）>epsilon））;????%%%將大于>epsilon的奇異值取出

U=U（:1:n）;???????????????????????????%%%[256n]
V=V（1:n1:n）;
V=diag（（diag（V）.^（-0.5）））;????????????%%%[nn]
M=V*U‘;????????????????????????????????%%%[nn]*[n256]->[n256]?
Z=M*X;?????????????????????????????????%%[nn]*[256len]=[nlen]
sum=0;
for?i=1:len
????sum=sum+（norm（Z（:i）））^2*Z（:i）*Z（:i）‘;
end????????????????????????????????????%%%[nn]
Rz4=sum/len;???????????????????????????%%%?四階矩大小為[nn]??
[B?S]=svd（Rz4）;????????????????????????%%B:[nn]
A=pinv（M）*B;???????????????????????????%%[256n]A的每一列極為一個基；
for?i=1:n
????A（:i）=255*（A（:i）-min（A（:i）））/（max（A（:i））-min（A（:i）））;?
end
%imgdata=uint8（imgdata）;
A=uint8（A）;
figure;
for?i=1:n
????base=reshape（A（:i）1616）;?????%%%%將基轉換為[1616]的圖
????subplot（ceil（n/floor（sqrt（n）））fix（sqrt（n））i）imshow（base）;?
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????618746??2010-01-03?14:41??4階統計法\100000.fig

?????文件?????619005??2010-01-03?14:45??4階統計法\100000_e-10.fig

?????文件???????1649??2010-01-03?13:44??4階統計法\j_four_1.asv

?????文件???????1668??2010-01-03?14:41??4階統計法\j_four_1.m

?????文件??????57924??2009-12-30?17:52??4階統計法\ref1000_e2.jpg

?????文件??????51744??2009-12-30?17:45??4階統計法\ref100_e15.jpg

?????文件??????15120??2009-12-30?17:40??4階統計法\ref100_e2.jpg

?????目錄??????????0??2010-06-26?14:32??4階統計法

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

??????????????1365856????????????????????8