hvsc.m

来自「混沌生成的很多形式」· M 代码 · 共 21 行

function y=hvsc(x)
[m,n]=size(x);
T=[x(m-4:m,:);x;x(1:4,:)];
I=[T(:,n-4:n) T T(:,1:4)];
gama=4;
for i=5:m+4,
    for j=5:n+4,
        It=[I(i-4:i+4,j-4:j+4)];
        %Mu=mean2(It);
        %a_9=(It-Mu).^gama;
        %a_9=std2(It).^2;%sample center moment
        %a_9=skewness(It);%third order moment
        Ckc=kurtscum(It);%forth order monent 
        %Ckr=kurtscum(It');
        a(i-4,j-4)=Ckc/9;
        %b(i-4,j-4)=sum(abs(Ckr))/9;
        %a(i-4,j-4)=sum(sum(abs(a_9)))/81;%9x9 array
    end
 end
%y=(a-min(min(a)))/(max(max(a))-min(min(a)))-0.5;%adding rule of generalization
y=a-(max(max(a))-min(min(a)))/2;%normalize only