tzz.m

klt经典程序

%%%%%%%%%%%%%%%%%%%%%%求一个8X8的矩阵的特征矩阵，被调用函数%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%求一个8X8的矩阵的特征矩阵，被调用函数%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%求一个8X8的矩阵的特征矩阵，被调用函数%%%%%%%%%%%%%%%%%%%%%%%%%%
 k=1;
 for n1=1:8      %%%%%%%%%%%%%%%%%%%%%%%%扫描8X8矩阵
     for n2=1:8
         fz(k)=A(n1,n2);    
         k=k+1;
     end
 end             %%%%%%%%%%%%%%%%%%%%%%%%扫描8X8矩阵
 k=k-1;
 aa=64;
 for k=1:64      %%%%%%%%%%%%%%%%%%%%将元素变换成一维
     X(k,1)=fz(k); 
 end             %%%%%%%%%%%%%%%%%%%%将元素变换成一维
 for k=1:64      %%%%%%%%%%%%%%%%%%减去均值后的元素值
     fz(k)=fz(k)-jz;
 end             %%%%%%%%%%%%%%%%%%减去均值后的元素值           
 B=zeros(64);
 for n1=1:64     %%%%%%%%%%%%%%%%%%%%%%%%求协方差矩阵
     for n2=1:64
        B(n1,n2)=fz(n1)*fz(n2); 
     end
 end             %%%%%%%%%%%%%%%%%%%%%%%%求协方差矩阵