给大家分享我自己编的程序-小波变换谱

pengzk

function [SG, PS, WES, FreqBins, Scales] = WTSpectrum(Sig,WinLen,nLevel);
%============================================================%
%  Continuous Wavelet Transform using Morlet function               
%            Sig : 信号                                          
%    WinLen : 小波函数在尺度参数a=1时的长度   （默认为 10）                 
%      nLevel : 频率轴划区间数   （默认为1024）      
%
%            SG:   返回的小波尺度谱计算结果
%            PS:    返回的小波相位谱计算结果
%         WES:   返回的小波能谱计算结果
%  FreqBins :  返回频率轴划分结果（归一化频率，最高频率为0.5）
%       Scales:   返回与频率轴划分值相对应的尺度划分 （频率0.5对应的尺度为1）
%============================================================%

if (nargin == 0),
     error('At least 1 parameter required');
end;

[WT, FreqBins, Scales] = Wavelet_Morl(Sig,WinLen,nLevel); %计算信号的连续小波变换
   SG = abs(WT).^2;                             % 计算小波尺度谱
    PS = angle(WT);                               % 计算小波相位谱
WES = sum(SG,2);                               % 计算小波能谱

jerry0204103

怎么应用呢？

leaukai

看着不错，运行一下试试。楼主有关于小波谱的介绍性文章吗？

leaukai

要是把画图的功能加上就好了