|
给你一个例子,你可以看一下,这是CWD分布,要是用平滑伪可以换成tfrspwv(sig);
clc
clear
N=1024;
fs=1024;
t=1/fs*(0:N/4-1);
x1=sin(2*pi*150*t);
t=1/fs*(N/4:N/2-1);
x2=sin(2*pi*350*t);
t=1/fs*(N/2:3*N/4-1);
x3=sin(2*pi*50*t);
t=1/fs*(3*N/4:N-1);
x4=sin(2*pi*20*t);
x=[x1 x2 x3 x4];
t=1/fs*(0:(N-1));
figure(1)
plot(t,x);xlabel('时间 t/s');ylabel('幅值 A/um');
%计算Choi-Williams分布
sig=x';
sig1=hilbert(sig);
[tfr,t,f]=tfrcw(sig1);
figure(2)
contour(t/fs,f(1:N)*fs,abs(tfr));
xlabel('Choi-Williams分布时频图 时间 t/s');
ylabel('频率 f/Hz');
grid on |
|