test.m

matlab求小波变换模极大值的代码

clc
close all
clear all
%%% you have to generate the signal to be estimated first.

fm=120; %% max- freq you wanna process.

  [T1,F1,W_orig]=Wigner(D,fm);  %%  
  [T2,F2,W_sp] = STFT(D,64,fm);         %%% 

figure;
tfmax = max(max(W_sp));
tfmin = min(min(W_sp));
colormap(1-gray(256))
image(T2,F2,256*(W_sp-tfmin)/(tfmax-tfmin));
colorbar
%contour(T2,F2,256*(W_sp-tfmin)/(tfmax-tfmin));
axis('xy')
title('STFT');
xlabel('Time')
ylabel('Frequency')
% 
figure;
tfmax = max(max(W_orig));
tfmin = min(min(W_orig));
colormap(1-gray(256))
image(T1,F1,256*(W_orig-tfmin)/(tfmax-tfmin));
colorbar
%contour(T1,F1,256*(W_orig-tfmin)/(tfmax-tfmin));
axis('xy')
title('WVD');
xlabel('Time')
ylabel('Frequency')