用四阶统计量方法对图像特征进行提取的MATLAB仿真程序

用四阶统计量方法对图像特征进行提取，进行分析，所用的是MATLAB仿真，仿真的程序和结果都有！

%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