matlab 小波包分解信号获取时频图

Posted 2023-03-04 studyer_domi

1、内容简介

略

643-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

略

clc
close all
clear
%% 小波包分解
load sep_wave1.mat;
x2=sep_wave;
fs = 1e6; % 采样频率
wavename='cmor3-3'; % 小波名称cmor3-3
totalscal=512*4; %totalscal是对信号进行小波变换时所用尺度序列的长度(通常需要预先设定好)
Fc=centfrq(wavename); % 小波的中心频率
c=2*Fc*totalscal;    % 测得
scals=c./(1:totalscal);
f=scal2frq(scals,wavename,1/fs); % 将尺度转换为频率
coefs = cwt(x2,scals,wavename); % 求连续小波系数

clear all  
clc
fs=1024;  %采样频率
f1=100;   %信号的第一个频率
f2=300;   %信号第二个频率
t=0:1/fs:1;
s=sin(2*pi*f1*t)+sin(2*pi*f2*t);  %生成混合信号
[tt]=wpdec(s,3,'dmey');  %小波包分解，3代表分解3层，'dmey'使用meyr小波
plot(tt)               %画小波包树图
wpviewcf(tt,1);        %画出时间频率图

4、参考论文

略