www.pudn.com > hdjm.rar > hdjm.m

function output=hdjm(inputx) 
yuantu=inputx;             %读入原图像 
[m,n]=size(yuantu);                         %计算图象大小 
 
% 1 从原图中抽取一个128*128的子图像 
 
           %% 将原图分解为1024个8*8小块 
point1=zeros(1024,2);                       %point1存储各个小块的序号及首个坐标值(x,y) 
s=1;                                        %s用来标记1024个小块的序号 
for i=1:8:m-7                                                  
    for j=1:8:n-7                                     
        point1(s,:)=[i j];                  %将每个小块的首个坐标值存入数组 
        s=s+1; 
    end                                      
end     
 
           %% 从1024个小块中选取256个                 
x=zeros(1,1024);                            %创建1*1024的数组x 
y=zeros(1,256);                             %创建1*256的数组y，存放x中的前256个不同的数 
x(1)=0.7;                                   %赋初值，满足0