伪+平滑wigner-ville的程序

clear all；clc；
N=1000;
fs=2000;
n=0:N;
t=n*1/2000;

x1=cos（2*pi*200*t）.*（t>=0 & t<=0.8）;
x2=cos（2*pi*50*t）.*（t>=0.2 & t<=1）;
x=x1+x2;
x=（x-mean（x））/std（x1）;  
figure（1）;
plot（tx）;
title（‘原始信号‘）;
xlabel（‘时间 t‘）;
ylabel（‘幅值 A‘）;
grid on;

x=（x-mean（x））/std（x1）;
x=hilbert（x）;
[tfrtf]=tfrwv（x‘1:NN）;
figure（2）;
contour（t/fsf（1:N）*fsabs（tfr））;
xlabel（‘时间 t/s‘）;
ylabel（‘频率 f/Hz‘）;
title（‘Wigner-Ville分布时频图‘）;
grid on
axis （[0 0.5 0 1000]）;
figure（3）;
mesh（t/fsf（1:N）*fsabs（tfr））;
xlabel（‘时间 t‘）;
ylabel（‘频率 f‘）;
zlabel（‘幅值 A‘）;
title（‘三维Wigner-Ville分布时频图‘）;

[tfr1t1f1]=tfrpwv（x‘1:NN）;
figure（4）;
contour（t1/fsf1（1:N）*fsabs（tfr1））;
xlabel（‘时间 t/s‘）;
ylabel（‘频率 f/Hz‘）;
title（‘伪Wigner-Ville分布时频图‘）;
grid on
axis （[0 0.