马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
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); % 计算小波能谱 |