连续小波变换matlab代码

下面是一段连续小波变换的MATLAB代码：

% 连续小波变换函数
function CWT = ContinuousWaveletTransform(sig, scale, fs, wname)
 
% sig：信号，scale：尺度向量，fs：采样率，wname：小波函数名称
 
% 计算CWT的函数主体 
numOfScales = length(scale); % 尺度个数 
sigLen = length(sig); % 信号长度 
dt = 1/fs; % 采样时间间隔 
CWT = zeros(numOfScales,sigLen); % 初始化存储CWT矩阵 
for i = 1:numOfScales 
    s = scale(i); % pick one scale
    t = (1:sigLen)*dt; % 构建时间坐标轴
    coef = Morlet(s,t,wname); % 计算Morlet小波函数系数
    temp = conv(sig, coef,'same'); % 对信号进行小波变换 
    CWT(i,:) = temp; % 存储结果
end
end
 
 
% Morlet小波函数
function [w] = Morlet(s,t,wname)
 
% s:尺度
% t:时间坐标轴
% wname:小波函数名称
% w:小波函数
 
switch wname
    case 'morlet'
        f0=6;
        w=exp(2*1i*pi*f0*t).*exp(-t.^2./(2*(s^2)))./sqrt(s); % Morlet小波函数
    case 'mexican_hat'
        w=(1-2*(t.^2)./s^2).*exp(-t.^2./(2*s^2))./sqrt(s/(3*pi)); % Mexican hat小波函数
end
end

希望能够帮助到您！