matlab课设.docxVIP

《matlab》课程设计告专业电子信息工程班级 0901 设计题目图像的傅里叶变换及应用 设计人姓名 杨亚娟 学号2009001244同组人 杨鸿琼杨雪余道军 吴书侃 吴宗江 2011/6/25设计四：图像的傅立叶变换及其应用设计目的：通过该设计，掌握傅立叶变换的定义和含义二、设计内容：1.图像的离散傅立叶变换2.快速傅立叶变换的应用三、总体方案设计：图像的离散傅里叶变换：假设是一个离散空间的二维函数，则该函数的二维傅里叶变换定义为其中是频域变量，单位是弧度/采样单元。函数）为函数的频谱二维傅里叶反变换的定义为因此，函数可以用无数个不同频率的复指数信号的和表示，在频率处复指数信号的幅度和相位为MATLAB提供的快速傅里叶变换函数：fft2:用于计算机二维快速傅里叶变换，其语法格式为b=fft2(Ι),返回图像Ι的二维傅里叶变换矩阵，输入图像Ι和输出图像B大小相同；b=fft2(Ι,m,n),通过对图形Ι剪贴或补零，按用户指定的点数计算二维傅里叶变换，返回矩阵B的大小为m×n。很多MATLAB图像显示函数无法显示复数图像，为了观察图像傅里叶变换后的结果，应对变换后的结果求模，方法是对变换结果使用abs函数。fftn:用于计算n维快速傅里叶变换，其语法格式为fftn: b=fftn(I),计算图像的n维傅里叶变换，输出图像B和输入图像I大小相同；B=FFTN(I,size),通过对图像I剪贴过补零，按size指定的点数计算n维傅里叶变换，返回矩阵B的大小为 size。fftshift：用于将变换后的图像频谱中心从矩阵的原点移到矩阵的中心，其语法格式为：b=fftshift(Ι),将变换后的图像频谱中心从矩阵的原点移到矩阵的中心。2、快速傅里叶变换的应用1）滤波器频率响应滤波器冲激响应的傅里叶变换就是该滤波器的频率响应。MATLAB提供的freqz2函数可以同时计算和现示滤波器的频率响应，其语法格式为：freqz2（h），可同时计算和现示滤波器的频率响应。2）快速卷积傅里叶变换的另一个重要特性是能够实现快速卷积，由线性系统理论可知，两个系统卷积的傅里叶变换等于两个函数的傅里叶变换的乘积。该特性与快速傅里叶变换一起，可以快速计算函数的卷积。假设A为M×N矩阵，B为P×Q矩阵，则快速卷积方法如下（1）对A和B补0，使其大小都为(M+P-1)×(N+Q-1)（2）利用函数fft2对矩阵进行二维FFT变换；（3）将两个FFT结果相乘，利用函数ifft2对得到的乘积进行傅里叶反变换另外，利用MATLAB提供的卷积函数conv2可以直接对两个函数进行卷积，其语法格式为： C=conv2（A,B），计算A和B的二维卷积。四、设计报告要求(1)加载imdemos saturn2文件，显示文件中的图像saturn2，并对其进行傅里叶变换，给出源程序及结果，并显示其幅值的结果(2)对矩阵A=magic(3)和B=ones(3)按照快速计算卷积的方法计算两个矩阵的卷积，并用卷积函数conv2验证结果，给出源程序及结果。五、源程序(1) y=imread(saturn.png); image(y) b=fft2(y);Warning: FFT on values of class UINT8 is obsolete. Use FFT(DOUBLE(X)) or FFT(SINGLE(X)) instead. In uint8.fft at 10 In fft2 at 25 x=abs(b); c=uint8(x); imshow(c)Warning: Image is too big to fit on screen; displaying at 33% In imuitools\private\initSize at 86 In imshow at 196(2) A=magic(3)A = 8 1 6 3 5 7 4 9 2 B=ones(3)B = 1 11 1 11 1 11 A(5,5)=0A = 8 1 6 0 0 3 5 7 0 0 4 9 2 0 0 0 0000 0 0000 B(5,5)=0B = 1 11 0 0 1 11 0 0 1 11 0 0