%圖像置亂程序
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')  % ？？？