wv2.m

时频分布wigner-ville的仿真,可以通过改变信号,来看函数的分布情况

z1=amgauss(256,90).*fmconst(256,0.26,90);
z2=amgauss(256,90).*fmconst(256,0.13,90);
z3=amgauss(256,90).*fmconst(256,0.39,90);
z4=amgauss(256,225).*fmconst(256,0.26,225);
z=(z1+z2+z3+z4)/4;
n=1:256;
plot(real(z));
[trf,t,f]=tfrwv(z);
contour(t,f,abs(trf));