hht2.m

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

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));