matlab求傅里叶变换
matlab⼆维离散傅⽴叶变换,⼆维离散傅⾥叶变换.PPT
⼆维离散傅⾥叶变换
3.1 ⼆维离散傅⾥叶变换(DFT) 3.1.1 ⼆维连续傅⾥叶变换 ⼆维连续函数 f (x, y)的傅⾥叶变换定义如下: 设 是独⽴变量 的函数,且在 上绝对可积,则定义积分 为⼆维连续函数 的付⾥叶变换,并定义 为 的反变换。 和 为傅⾥叶变换对。 3.1.2 ⼆维离散傅⾥叶变换 尺⼨为
M×N的离散图像函数的DFT 反变换可以通过对F(u,v) 求IDFT获得 DFT变换进⾏图像处理时有如下特点: (1)直流成分为F(0,0)。 (2)幅度谱|F(u,v)|对称于原点。 (3)图像f (x, y)平移后,幅度谱不发⽣变化,仅有相位发⽣了变化。 3.1.3 ⼆维离散傅⾥叶变换的性质 1.周期性和共轭对称性 周期性和共轭对称性来了许多⽅便。 我们⾸先来看⼀维的情况。 设有⼀矩形函数为,求出它的傅⾥叶变换: 幅度谱: DFT取的区间是[0,N-1],在这个区间内频谱是由两个背靠背的半周期组成的 ,要显⽰⼀个完整的周期,必须将变换的原点移⾄u=N/2点。 根据定义,有 在进⾏DFT之前⽤(-1)x 乘以输⼊的信号 f (x) ,可以在⼀个周期的变换中(u=0,1,2,…,N-1),求得⼀个完整的频谱。 推⼴到⼆维情况。在进⾏傅⾥叶变换之前⽤(-1)x+y 乘以输⼊的图像函数,则有: DFT的原点,即F(0,0)被设置在u=M/2和v=N/2上。(0,0)点的变换值为: 即 f (x,y) 的平均值。 如果是⼀幅图像,在原点的傅⾥叶变换F(0,0)等于图像的平均灰度级,也称作频率谱的直流成分。 2.可分性 离散傅⾥叶变换可以⽤可分离的形式表⽰ 这⾥ 对于每个x值,当
v=0,1,2,…,N-1时,该等式是完整的⼀维傅⾥叶变换。 3.离散卷积定理 设f(x,y)和g(x,y) 是⼤⼩分别为A×B和C×D的两个数组,则它们的离散卷积定义为 卷积定理 【例3.2】⽤MATLAB实现图像的傅⾥叶变换。 解:MATLAB程序如下: A=imread('pout.tif'); %读⼊图像 imshow(A); %显⽰图像 A2=fft2(A); %计算⼆维傅⾥叶变换 A2=fftshift(A2); %将直流分量移到频谱图的中⼼ figure, imshow(log(abs(A2)+1),[0 10]); %显⽰变换后的频谱图3.2 ⼆维离散余弦变换(DCT) 任何实对称函数的傅⾥叶变换中只含余弦项,余弦变换是傅⾥叶变换的特例,余弦变换是简化DFT的重要⽅法。 3.2.1 ⼀维离散余弦变换 将⼀个信号通过对折延拓成实偶函数,然后进⾏傅⾥叶变换,我们就可⽤2N点的DFT来产⽣N点的DCT。1.以x=-1/2为对称轴折叠原来的实序列f(n) 得: 3.对0到2N-1的2N个点的离散周期序列 作DFT,得 令i=2N-m-1,则上式为 为了保证变换基的规范正交性,引⼊常量,定义: DCT逆变换为 【例3.3】应⽤MATLAB实现图像的DCT变换。 解:MATLAB程序如下:A=imread('pout.tif'); %读⼊图像 I=dct2(A); %对图像作DCT变换 subplot(1,2,1),imshow(A); %显⽰原图像
subplot(1,2,2),imshow(log(abs(I)),[0 5]); 3.3 ⼆维离散沃尔什-哈达玛变换(DHT) 前⾯的变换都是余弦型变换,基底函数选⽤的都是余弦型。 图像处理中还有许多变换常常选⽤⽅波信号或者它的变形。 沃尔什(Walsh)变换。 沃尔什函数是⼀组矩形波,其取值为1和-1,⾮常便于计算机运算。 沃尔什函数有三种排列或编号⽅式,以哈达玛排列最便于快速计算。 采⽤哈达玛排列的沃尔什函数进⾏的变换称为
沃尔什-哈达玛变换,简称WHT或直称哈达玛变换。 3.3.1 哈达玛变换 哈达玛矩阵:元素仅由+1和-1组成的正交⽅阵。 正交⽅阵:指它的任意两⾏(或两列)都彼此正交,或者说它们对应元素之和为零。 哈达玛变换要求图像的⼤⼩为N=2n 。 ⼀维哈达玛变换

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。