维离散傅里叶变换及相关性质实验.docVIP

撰写人姓名： 撰写时间： 审查人姓名： 实 验 全 过 程 记 录 实验 名称 二维离散傅里叶变换及相关性质实验 时间 2 学时 地点 数学实验室 姓 名 学 号 班 组 同实验者 学 号 班 组 一、实验目的 1、掌握二维傅里叶变换的原理和方法； 2、掌握编程实现生成图像频谱图、中心化频谱及相互比较等方法。 二、实验内容： 编程实现图像频谱图、中心化频谱，并进行相互比较。 三、实验用仪器设备及材料 软件需求： 操作系统：Windows Xp或更新的版本 开发工具：MATLAB 7.01、Photoshop CS3 硬件需求： Pentium Ⅲ 1G以上的CPU处理器、256MB以上的内存、1.5G以上自由硬盘空间、CD-ROM驱动器、打印机、打印纸等。 四、实验步骤： 1、复习二维傅里叶变换的原理和方法； 2、启动MATLAB软件； 3、编程实现如下要求： ⑴产生右图所示亮块图像（128×128，暗处=0，明处=255），对其进行FFT（快速傅里叶变换）； ⑵同屏显示原图和FFT()的频谱图； ⑶令，重复以上过程，比较二者频谱的异同，并简述理由； ⑷将顺时针旋转得到，显示FFT()的频谱，并与FFT()的频谱进行比较。 4、撰写实验报告内容 ⑴简述实验过程； ⑵复制实验程序； ⑶复制实验结果（如结果数据、生成的图像等）； ⑷心得体会。 五、撰写实验报告： %实验一：二维离散傅里叶变换及相关性质实验 %构造原始图像 f1=zeros(128,128); f1(24:104,48:80)=1; subplot(3,2,1); imshow(f1); xlabel(构造原始图像); %求原始图像的傅里叶频谱 J=fft2(f1); F=abs(J); J1=fftshift(F); subplot(3,2,2); imshow(J1,[5 50]); xlabel(原始图像的傅里叶频谱); %构造f2的图像 f2=f1; for x=1:128 for y=1:128 f2(x,y)=(-1)^(x+y)*f1(x,y); end end subplot(3,2,3); imshow(f2); xlabel(构造f2的图像); %求f2的傅里叶频谱 I=fft2(f2); F=abs(I); J2=fftshift(F); subplot(3,2,4); imshow(J2,[5 50]); xlabel(f2的傅里叶频谱); %将f2顺时针旋转得到f3 f3=imrotate(f2,-45,bilinear,crop); subplot(3,2,5); imshow(f3); xlabel(f2图像顺时针旋转45度); %求f3的傅里叶频谱 K=fft2(f3); F=abs(K); J3=fftshift(F); subplot(3,2,6); imshow(J3,[5 50]); xlabel(f3的傅里叶频谱); 成绩评定： 指导教师： 年 月 日