[信息与通信]3频率变换.ppt

频率变换 1 频率概述 2 频率变换概述 3 离散傅里叶变换 4 其他频率变换 5 滤波处理 1 频率概述 1 频率概述 1 频率概述 1 频率概述 2 频率变换概述 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 3 离散傅里叶变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 4 其他频率变换 5 滤波处理 5 滤波处理 5 滤波处理 5 滤波处理 5 滤波处理 5 滤波处理 5 滤波处理 频率变换 MATLAB Functions Terms Terms Terms 正弦型变换 DCT的核矩阵为： 与DFT不同的是，DCT是实值的，在图像编码领域有广泛的应用。可以用快速傅立叶变换算法实现快速离散余弦变换。MATLAB中进行二维离散余弦变换的函数为dct2。 余弦变换（DCT） 例：图像的二维余弦变换示例。 I1=imread(whitehouse.bmp); %读入一幅图像 subplot(2,2,1),imshow(I1) ,title(图像1);%显示图像1 C1=dct2(I1); %对图像1进行离散余弦变换 C2=log(abs(C1)); %对频谱的幅度取对数 subplot(2,2,2),imshow(C2,[]) ,title(图像1的频谱); I2=imread(appletree.bmp); %读入另一幅图像 subplot(2,2,3),imshow(I2) ,title(图像2); %显示图像2 C4=dct2(I2); %对图像2进行离散余弦变换 C5=log(abs(C4)); %对频谱的幅度取对数 subplot(2,2,4),imshow(C5,[]) ,title(图像2的频谱); Hartley变换 二维离散Hartley变换（DHT，Discrete Hartley Transform）为 基函数为 正弦型变换 Hartley变换 二维离散Hartley逆变换 二维DHT的核矩阵为 正弦型变换 DHT实际上是DFT的另一种计算方法，由于它避免了复数计算，从而可以大大地减少计算量。 沃尔什-哈达玛变换 为了减少计算量，还有一类图像变换采用方波型基函数，即沃尔什-哈达玛变换（WHT, Walsh-Hardamard Transform）。 沃尔什-哈达玛变换变换是对称的、可分离的酉变换，它的核矩阵中只有+1和-1两种元素。它的维数N=2n，其中n是整数。离散沃尔什-哈达玛变换变换为 方波型变换 逆变换为 沃尔什-哈达玛变换 核矩阵HN按如下方式构造： 方波型变换 最小阶（N＝2）的核矩阵为： 利用这个迭代性质求解N阶（N=2n）的哈达玛变换矩阵，便给出快速沃尔什-哈达玛变换算法WHT。 例：图像的沃尔什-哈达玛变换示例。 I=imread(cameraman.tif); subplot(121),imshow(I), title(原始图像);%显示原始图像 I=double(I); [m,n]=size(I); for j=1:n whtj(:,j)=hadamard(m)*I(:,j)/m; end for i=1:m whti(:,i)=hadamard(n)*whtj(i,:)/n; end wht=whti; subplot(122),imshow(wht) , title(频谱);%显示原始图像的频谱 滤波处理的目的： 使某些频率通过，使某些频率阻断。 让我们通过设定右图所示的参数a和b的值，使a以上、b以下的频率（斜线表示的频率）通过，其它的频率阻断来进行滤波处理。 例：小于40的频率通过，其余的频率阻断的MATLAB程序示例。 I=imread(whitehouse.bmp);%读入一幅BMP灰度图像 figure,subplot(2,2,1),imshow(I); title(原始图像);%把图形窗口划分为2×2显示区域，在左上角的区域显示图像I F=fft2(I); %傅立叶变换 F1=fftshift(F);%频谱中心化 F2=log(abs(F1)); %为了显示F1的幅值对其取对数 subplot(2,2,2),imshow(F2,[]); title(频谱);%在右上角的区域显示频谱F2 [M,N]=size(F1);%获取频谱F1的高度和宽度 m=round(M/2);n=round