xieyafeng.m

这么代码是用Matlab编写的

%图像置乱程序
x0=imread('e:\work\lena1.jpg');%读入图像
x=double(x0);   %将读入图像由unit8型转化为double型
e=eye(256);   %生成单位矩阵
X=fft2(x);
E=fft2(e);
y=ifft2(X.*E);   %对X，E在频域中点乘后再进行快速傅立叶反变换
subplot(1,2,1),imshow(x0),title('图1')
subplot(1,2,2),imshow(abs(y),[]),title('图2')
%图象恢复程序
% x0=imread('d:\lena.jpg');   %读入图像
% x=double(x0);    %图像由unit8型转化为double型
% e=eye(256);     %生成单位矩阵
% X=fft2(x);
% E=fft2(e);
% M=X.*E;       %恢复开始
% y=ifft2(M);
% Y=fft2(y);     
% Y=M./E;       %点除，结果为M对应元素除以H对应元素
% y0=ifft2(Y);
% subplot(1,2,1),imshow(abs(y),[]),title('图3')
% subplot(1,2,2),imshow(abs(y0),[]),title('图4')  % ？？？