给大家分享我自己编的程序-小波变换谱
function = 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;
= Wavelet_Morl(Sig,WinLen,nLevel); %计算信号的连续小波变换
SG = abs(WT).^2; % 计算小波尺度谱
PS = angle(WT); % 计算小波相位谱
WES = sum(SG,2); % 计算小波能谱 怎么应用呢? 看着不错,运行一下试试。楼主有关于小波谱的介绍性文章吗? 要是把画图的功能加上就好了
页:
[1]