hht4.m

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

z1=amgauss(256,65).*fmconst(256,0.26,65);
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;
n=1:256;
[y,iflaw]=fmlin(256,0.01);
zzz=z.*y;
plot(real(zzz))