新浪博客

[转载]matlab小波变换

2018-01-28 15:13阅读:

迷糊了好几天突然顿悟了,博主写的很好

原文作者:夏天的味道

1. 离散傅立叶变换的 Matlab实现
Matlab 函数 fft、fft2 和 fftn 分别可以实现一维、二维和 N 维 DFT 算法;而函数 ifft、ifft2 和 ifftn 则用来计算反 DFT 。这些函数的调用格式如下:
A=fft(X,N,DIM)
其中,X 表示输入图像;N 表示采样间隔点,如果 X 小于该数值,那么 Matlab 将会对 X 进行零填充,否则将进行截取,使之长度为 N ;DIM 表示要进行离散傅立叶变换。
A=fft2(X,MROWS,NCOLS)
其中,MROWS 和 NCOLS 指定对 X 进行零填充后的 X 大小。
A=fftn(X,SIZE)
其中,SIZE 是一个向量,它们每一个元素都将指定 X 相应维进行零填充后的长度。
函数 ifft、ifft2 和 ifftn的调用格式于对应的离散傅立叶变换函数一致。
例子:图像的二维傅立叶频谱
% 读入原始图像
I=imread('lena.bmp');
imshow(I)
% 求离散傅立叶频谱
J=fftshift(fft2(I));
figure;
imshow(log(abs(J)),[8,10])
2. 离散余弦变换
的 Matlab 实现

2.1. dct2 函数
功能:二维 DCT 变换
格式:B=dct2(A)
B=dct2(A,m,n)
B=dct2(A,[m,n])
说明:B=dct2(A) 计算 A 的 DCT 变换 B ,A 与 B 的大小相同;B=dct2(A,m,n) 和 B=dct2(A,[m,n]) 通过对 A 补 0 或剪裁,使 B 的大小为 m×n。
2.2. dict2 函数
功能:DCT 反变换
格式:B=idct2(A)
B=idct2(A,m,n)
B=idct2(A,[m,n])
说明:B=idct2(A) 计算 A 的 DCT 反变换 B ,A 与 B 的大小相同;B=idct2(A,m,n) 和 B=idct2(A,[m,n]) 通过对 A 补 0 或剪裁,使 B 的大小为 m×n。
2.3. dctmtx函数
功能:计算 DCT 变换矩阵
格式:D=dctmtx(n)
说明:D=dctmtx(n) 返回一个 n×n 的 DCT 变换矩阵,输出矩阵 D 为 double 类型。
3. 图像小波变换的 Matlab 实现

3.1 一维小波变换的 Matlab 实现
(1) dwt 函数
功能:一维离散小波变换
格式:[cA,cD]=dwt(X,'wname')
[cA,cD]=dwt(X,Lo_D,Hi_D)
说明:[cA,cD]=dwt(X,'wname') 使用指定的小波基函数 'wname' 对信号X 进行分解,cA、cD 分别为近似分量和细节分量;[cA,cD]=dwt(X,Lo_D,Hi_D) 使用指定的滤波器组 Lo_D、Hi_D 对信号进行分解。
(2) idwt 函数
功能:一维离散小波反变换
格式:X=idwt(cA,cD,'wname')
X=idwt(cA,cD,Lo_R,Hi_R)
X=idwt(cA,cD,'wname',L)
X=idwt(cA,cD,Lo_R,Hi_R,L)
说明:X=idwt(cA,cD,'wname') 由近似分量 cA 和细节分量 cD 经小波反变换重构原始信号 X 。
'wname' 为所选的小波函数
X=idwt(cA,cD,Lo_R,Hi_R) 用指定的重构滤波器 Lo_R 和 Hi_R 经小波反变换重构原始信号 X 。
X=idwt(cA,cD,'wname',L) 和 X=idwt(cA,cD,Lo_R,Hi_R,L) 指定返回信号 X 中心附近的 L 个点。
3.2 二维小波变换的 Matlab 实现
二维小波变换的函数
-------------------------------------------------
函数名 函数功能
---------------------------------------------------
dwt2 二维离散小波变换
wavedec2 二维信号的多层小波分解
idwt2 二维离散小波反变换
waverec2 二维信号的多层小波重构
wrcoef2 由多层小波分解重构某一层的分解信号
upcoef2 由多层小波分解重构近似分量或细节分量
detcoef2 提取二维信号小波分解的细节分量
appcoef2 提取二维信号小波分解的近似分量
upwlev2 二维小波分解的单层重构
dwtpet2 二维周期小波变换
idwtper2 二维周期小波反变换
-------------------------------------------------------------
(1) wcodemat 函数
功能:对数据矩阵进行伪彩色编码
格式:Y=wcodemat(X,NB,OPT,ABSOL)
Y=wcodemat(X,NB,OPT)
Y=wcodemat(X,NB)
Y=wcodemat(X)
说明:Y=wcodemat(X,NB,OPT,ABSOL) 返回数据矩阵 X 的编码矩阵 Y ;NB 伪编码的最大值,即编码范围为 0~NB,缺省值 NB=16;
OPT 指定了编码的方式(缺省值为 'mat'),即:
OPT='row' ,按行编码
OPT='col' ,按列编码
OPT='mat' ,按整个矩阵编码
ABSOL 是函数的控制参数(缺省值为 '1'),即:
ABSOL=0 时,返回编码矩阵
ABSOL=1 时,返回数据矩阵的绝对值 ABS(X)
(2) dwt2 函数
功能:二维离散小波变换
格式:[cA,cH,cV,cD]=dwt2(X,'wname')
[cA,cH,cV,cD]=dwt2(X,Lo_D,Hi_D)
说明:[cA,cH,cV,cD]=dwt2(X,'wname')使用指定的小波基函数 'wname' 对二维信号 X 进行二维离散小波变幻;cA,cH,cV,cD 分别为近似分量、水平细节分量、垂直细节分量和对角细节分量;[cA,cH,cV,cD]=dwt2(X,Lo_D,Hi_D) 使用指定的分解低通和高通滤波器 Lo_D 和 Hi_D 分解信号 X 。
(3) wavedec2 函数
功能:二维信号的多层小波分解
格式:[C,S]=wavedec2(X,N,'wname')
[C,S]=wavedec2(X,N,Lo_D,Hi_D)
说明:[C,S]=wavedec2(X,N,'wname') 使用小波基函数 'wname' 对二维信号 X 进行 N 层分解;[C,S]=wavedec2(X,N,Lo_D,Hi_D) 使用指定的分解低通和高通滤波器 Lo_D 和 Hi_D 分解信号 X 。
(4) idwt2 函数
功能:二维离散小波反变换
格式:X=idwt2(cA,cH,cV,cD,'wname')
X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R)
X=idwt2(cA,cH,cV,cD,'wname',S)
X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R,S)
说明:X=idwt2(cA,cH,cV,cD,'wname') 由信号小波分解的近似信号 cA 和细节信号 cH、cH、cV、cD 经小波反变换重构原信号 X ;X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R) 使用指定的重构低通和高通滤波器 Lo_R 和 Hi_R 重构原信号 X ;X=idwt2(cA,cH,cV,cD,'wname',S) 和 X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R,S) 返回中心附近的 S 个数据点。
(5) waverec2 函数
说明:二维信号的多层小波重构
格式:X=waverec2(C,S,'wname')
X=waverec2(C,S,Lo_R,Hi_R)
说明:X=waverec2(C,S,'wname') 由多层二维小波分解的结果 C、S 重构原始信号 X ,'wname' 为使用的小波基函数;X=waverec2(C,S,Lo_R,Hi_R) 使用重构低通和高通滤波器 Lo_R 和 Hi_R 重构原信号。
Allnodes 计算树结点
  appcoef 提取一维小波变换低频系数
  appcoef2 提取二维小波分解低频系数
  bestlevt 计算完整最佳小波包树
  besttree 计算最佳(优)树
*  biorfilt 双正交样条小波滤波器组
  biorwavf 双正交样条小波滤波器
*  centfrq 求小波中心频率
  cgauwavf Complex Gaussian小波
  cmorwavf coiflets小波滤波器
  cwt 一维连续小波变换
  dbaux Daubechies小波滤波器计算
  dbwavf Daubechies小波滤波器 dbwavf(W) W='dbN' N=1,2,3,...,50
  ddencmp 获取默认值阈值(软或硬)熵标准
  depo2ind 将深度-位置结点形式转化成索引结点形式
  detcoef 提取一维小波变换高频系数
  detcoef2 提取二维小波分解高频系数
  disp 显示文本或矩阵
  drawtree 画小波包分解树(GUI)
  dtree 构造DTREE类
  dwt 单尺度一维离散小波变换
  dwt2 单尺度二维离散小波变换
  dwtmode 离散小波变换拓展模式
*  dyaddown 二元取样
*  dyadup 二元插值
  entrupd 更新小波包的熵值
  fbspwavf B样条小波
  gauswavf Gaussian小波
  get 获取对象属性值
  idwt 单尺度一维离散小波逆变换
  idwt2 单尺度二维离散小波逆变换
  ind2depo 将索引结点形式转化成深度—位置结点形式
*  intwave 积分小波数
  isnode 判断结点是否存在
  istnode 判断结点是否是终结点并返回排列值
  iswt 一维逆SWT(Stationary Wavelet Transform)变换
  iswt2 二维逆SWT变换
  leaves   Determine terminal nodes
  mexihat 墨西哥帽小波
  meyer Meyer小波
  meyeraux Meyer小波辅助函数
  morlet Morlet小波
  nodease 计算上溯结点
  nodedesc 计算下溯结点(子结点)
  nodejoin 重组结点
  nodepar 寻找父结点
  nodesplt 分割(分解)结点
  noleaves   Determine nonterminal nodes
  ntnode   Number of terminal nodes
  ntree   Constructor for the class NTREE
*  orthfilt 正交小波滤波器组
  plot 绘制向量或矩阵的图形
*  qmf 镜像二次滤波器
  rbiowavf   Reverse biorthogonal spline wavelet filters
  read 读取二进制数据
  readtree 读取小波包分解树
*  scal2frq   Scale to frequency
  set  
  shanwavf   Shannon wavelets
  swt 一维SWT(Stationary Wavelet Transform)变换
  swt2 二维SWT变换
  symaux   Symlet wavelet filter computation.
  symwavf Symlets小波滤波器
  thselect 信号消噪的阈值选择
  thodes   References
  treedpth 求树的深度
  treeord 求树结构的叉数
  upcoef 一维小波分解系数的直接重构
  upcoef2 二维小波分解系数的直接重构
  upwlev 单尺度一维小波分解的重构
  upwlev2 单尺度二维小波分解的重构
  wavedec 单尺度一维小波分解
  wavedec2 多尺度二维小波分解
  wavedemo 小波工具箱函数demo
*  wavefun 小波函数和尺度函数
*  wavefun2 二维小波函数和尺度函数
  wavemenu 小波工具箱函数menu图形界面调用函数
*  wavemngr 小波管理函数
  waverec 多尺度一维小波重构
  waverec2 多尺度二维小波重构
  wbmpen   Penalized threshold for wavelet 1-D or 2-D de-noising
  wcodemat 对矩阵进行量化编码
  wdcbm   Thresholds for wavelet 1-D using Birge-Massart strategy
  wdcbm2  Thresholds for wavelet 2-D using Birge-Massart strategy
  wden 用小波进行一维信号的消噪或压缩
  wdencmp  De-noising or compression using wavelets
  wentropy 计算小波包的熵
  wextend  Extend a vector or a matrix
*  wfilters 小波滤波器
  wkeep 提取向量或矩阵中的一部分
*  wmaxlev 计算小波分解的最大尺度
  wnoise 产生含噪声的测试函数数据
  wnoisest 估计一维小波的系数的标准偏差
  wp2wtree 从小波包树中提取小波树  
  wpcoef 计算小波包系数
  wpcutree 剪切小波包分解树
  wpdec 一维小波包的分解
  wpdec2 二维小波包的分解
  wpdencmp 用小波包进行信号的消噪或压缩
  wpfun 小波包函数
wpjoin  重组小波包
  wprcoef 小波包分解系数的重构
  wprec 一维小波包分解的重构
  wprec2 二维小波包分解的重构
  wpsplt 分割(分解)小波包
  wpthcoef 进行小波包分解系数的阈值处理
  wptree   显示小波包树结构
  wpviewcf   Plot the colored wavelet packet coefficients.
  wrcoef 对一维小波系数进行单支重构
  wrcoef2 对二维小波系数进行单支重构
  wrev 向量逆序
  write 向缓冲区内存写进数据
  wtbo   Constructor for the class WTBO
  wthcoef 一维信号的小波系数阈值处理
  wthcoef2 二维信号的小波系数阈值处理
  wthresh 进行软阈值或硬阈值处理
  wthrmngr 阈值设置管理
  wtreemgr 管理树结构

我的更多文章

下载客户端阅读体验更佳

APP专享