小波时频变换matlab代码

以下是小波时频变换的MATLAB代码：

1. 离散小波变换（DWT）：

% 首先定义一个信号
t = 0:0.01:1;
x = sin(2*pi*20*t) + sin(2*pi*50*t);

% 定义小波参数
wname = 'db4';
level = 5;

% 进行离散小波变换
[c,l] = wavedec(x,level,wname);

% 重构信号
rec = waverec(c,l,wname);

% 绘制原始信号和重构信号
subplot(2,1,1);
plot(t,x);
title('Original Signal');
subplot(2,1,2);
plot(t,rec);
title('Reconstructed Signal');

2. 小波包变换（DWP）：

% 首先定义一个信号
t = 0:0.01:1;
x = sin(2*pi*20*t) + sin(2*pi*50*t);

% 定义小波包参数
wname = 'db4';
level = 5;
wp = wpdec(x,level,wname);

% 获取小波包节点
nodes = wpnodes(wp);

% 绘制小波包分解树
plot(wp);
title('Wavelet Packet Decomposition Tree');

% 获取小波包系数
c = wprcoef(wp,nodes);

% 重构信号
rec = wrcoef('p',wp,c);

% 绘制原始信号和重构信号
subplot(2,1,1);
plot(t,x);
title('Original Signal');
subplot(2,1,2);
plot(t,rec);
title('Reconstructed Signal');

3. 小波时频分析：

% 首先定义一个信号
t = 0:0.01:1;
x = sin(2*pi*20*t) + sin(2*pi*50*t);

% 定义小波参数
wname = 'db4';
level = 5;

% 进行小波分解
[c,l] = wavedec(x,level,wname);

% 定义时间窗口和频率窗口
twin = 0.05;
fwin = 5;

% 计算小波包系数的时频分布
[s,~,t,f] = spectrogram(c,twin*100,fwin,[],1/0.01);

% 绘制时频图
imagesc(t,f,abs(s));
set(gca,'YDir','normal');
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('Wavelet Time-Frequency Analysis');