A=imread（‘5.jpg‘）;%讀入原始圖像
figure;
subplot（141）;
imshow（A）title（‘原始圖像‘）;%顯示原始圖像
B=imresize（A[512?512]）;
subplot（142）;
imshow（B）title（‘轉化圖像‘）;%顯示轉化圖像
C=rgb2gray（B）;%獲取灰度圖
subplot（143）;?
imshow（C）title（‘灰度圖像‘）;%顯示灰度圖像
imwrite（C‘lotus_grayscale.bmp‘‘bmp‘）;?%?存儲為灰度圖像
R=rand（512512）;
CC=double（C）;
M=CC.*exp（i*2*pi.*R）;%數字化圖像乘以隨機相位因子平滑傅里葉變換譜
Am=abs（M）;