DFT Matlab代码

Matlab代码开发DFT,里面有详细的过程和照片，是课程作业的一部分

function DFT_2Image（）
Dog=imread（‘Dog.jpg‘）;
Owl=imread（‘Owl.jpg‘）;

subplot（321）;
imshow（Dog）;
title（‘The picture of Dog‘）;
subplot（322）;
imshow（Owl）;
title（‘The picture of Owl‘）;

grayDog=rgb2gray（Dog）;
subplot（323）;
imshow（grayDog）;
title（‘The gray picture of Dog‘）;
grayOwl=rgb2gray（Owl）;
subplot（324）;
imshow（grayOwl）;
title（‘The gray picture of Owl‘）;

dogAmp=abs（fft2（double（grayDog）））;
dogPhase=angle（fft2（double（grayDog）））;
owlAmp=abs（fft2（double（grayOwl）））;
owlPhase=angle（fft2（double（grayOwl）））;

exchangeDog=ifft2（dogAmp.*exp（sqrt（-1）*owlPhase））;
exchangeOwl=ifft2（owlAmp.*exp（sqrt（-1）*dogPhase））;
subplot（325）;
imshow（abs（exchangeDog）[]）;
title（‘The picture of Dog‘）;
subplot（326）;
imshow（abs（exchangeOwl）[]）;
title（‘The picture of Owl‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      26104  2014-04-01 16:01  DFT\DFT_2Image.jpg

     文件        799  2014-04-01 15:59  DFT\DFT_2Image.m

     文件      21210  2014-04-01 16:02  DFT\DFT_IDFT.jpg

     文件        467  2014-04-01 16:20  DFT\DFT_IDFT.m

     文件      19859  2014-04-01 14:22  DFT\Dog.jpg

     文件      14078  2014-02-24 20:37  DFT\MyPicture.jpg

     文件      41836  2014-04-01 12:44  DFT\Owl.jpg

     目录          0  2014-04-01 16:20  DFT

----------- ---------  ---------- -----  ----

               124353                    8