matlab 如何实现哈达玛

clc;
cr=05;
I=imread('8tif');
I1=double(I)/255;%图像为256级灰度图像，对图像进行归一化 *** 作
subplot(1,2,1);imshow(I1);
xlabel('(a) 原始图像');%显示原始图像
disp('原始图像的大小为：');
whos('I1')
%对图像进行哈达玛变换
T=hadamard(8);%产生88的哈达玛矩阵
htcoe=blkproc(I1,[8 8],'P1xP2',T,T);%将图像分割为88的子图像进行FFT
coevar=im2col(htcoe,[8 8],'distinct');%降变换系数矩阵重新排列
coe=coevar;
[y,ind]=sort(coevar);
[m,n]=size(coevar);%根据压缩比确定要变0的系数个数

clc;
cr=05;
I=imread('8tif');
I1=double(I)/255;%图像为256级灰度图像，对图像进行归一化 *** 作
subplot(1,2,1);imshow(I1);
xlabel('(a) 原始图像');%显示原始图像
disp('原始图像的大小为：');
whos('I1')
%对图像进行哈达玛变换
T=hadamard(8);%产生88的哈达玛矩阵
htcoe=blkproc(I1,[8 8],'P1xP2',T,T);%将图像分割为88的子图像进行FFT
coevar=im2col(htcoe,[8 8],'distinct');%降变换系数矩阵重新排列
coe=coevar;
[y,ind]=sort(coevar);
[m,n]=size(coevar);%根据压缩比确定要变0的系数个数

---