Cómo utilizar matlab para encontrar la distribución pseudo-Wigner-Weil

Borrar todo;clc;

%Señal experimental

sig=load('D:\data \ zheng1234.txt');

N=length(sig);

%Normalizar

sig=(sig-mean(sig))/std( sig ,1);

%Frecuencia de muestreo

fs=10000;

% Calcular la distribución pseudo-Wigner-Veile

sig= (sig-media(sig))/std(sig,1);

sig=hilbert(sig);

[tfr,t,f]=tfrpwv(sig);

figura(1);

contorno(t/fs,f(1：N)*fs,abs(tfr));

xlabel('tiempo t/s ');

ylabel('frecuencia f/Hz');

cuadrícula en

eje ([0 0.02 0 5000]);

figura(2);

mesh(t/fs,f(1:N)*fs,abs(tfr));

xlabel(' tiempo t') ;

ylabel('frecuencia f');

zlabel('amplitud A');

title('gráfico 3D');