hht3.m

来自「时频分布wigner-ville的仿真,可以通过改变信号,来看函数的分布情况」· M 代码 · 共 8 行

z1=amgauss(256,65).*fmconst(256,0.36,65);
z2=amgauss(256,65).*fmconst(256,0.24,65);
z3=amgauss(256,65).*fmconst(256,0.48,65);
z4=amgauss(256,195).*fmconst(256,0.36,195);
z=(z1+z2+z3+z4)/4;
n=1:256;
zz=z.*fmconst(256,0.07);
plot(real(zz))