tfw2.m

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

z1=amgauss(256,65).*fmconst(256,0.26,65);
plot(real(z1));
z2=amgauss(256,65).*fmconst(256,0.13,65);
z3=amgauss(256,65).*fmconst(256,0.39,65);
z4=amgauss(256,195).*fmconst(256,0.26,195);
z=(z1+z2+z3+z4)/4;
figure
plot(real(z));
n=1:256;
zz=z.*fmconst(256,0.07);
plot(real(zz))
[y,iflaw]=fmlin(256,0.01);
zzz=z.*y;

[trf,t,f]=tfrwv(zzz);
contour(t,f,abs(trf));
tfrwv(zz)