时频分析工具箱中的函数学习

AMGAUSS函数： Y=AMGAUSS（N,T0,T）
产生高斯幅值调制信号，N为信号点数，T0为中心，传播0.5的比例为T，即Y(T0)=1，Y(T0+T/2)=Y(T0-T/2)=0.5。T0默认为N/2，T默认为2*sqrt(N)。
例：z=AMGAUSS(200,90,40);plot(z);

FMCONST函数：[Y,IFLAW]=FMCONST(N,FNORM,T0)
产生一个固定频率的频率调制信号。N为产生信号的点数，FNORM为标椎化频率（默认为0.25），T0表示此时刻为正幅值（默认为round（N/2）），相当于规定了相位。
例：z=FMCONST(160,0.05,90);plot(real(z));

FMLIN函数：[Y,IFLAW]=FMLIN(N,FNORMI,FNORMF,T0)
产生一个线性调频信号。N为信号点数，FNORMI为开始标准频率（默认为0），FNORMF为结束标准频率（默认为0.5），相位由T0规定，即Y（T0）=1。IFLAW为设置的频率变化情况，这里就是一个线性变化的直线。
例：[z,f]=fmlin(128,0.05,0.3,50);plot(real(z));figure,plot(f);

短时傅里叶变换：
[tfr,t,f]=tfrstft(x,t,n,h,trace)
说明：
x–信号
t–时间（缺省值为1：length（x））
n–频率数（缺省值为length（x））
h–频率滑窗，h为归一化单位能量（缺省值为hamming（n/4））
trace–如果非零，显示算法的进程（缺省值为0）
tfr–时频分解（为复值），频率轴观察范围为-0.5~0.5
f–归一化频率
例：sig=[fmconst(128,0.2);fmconst(128,0.4)];
tfr=tfrstft(sig);
subplot(211);imagesc(abs(tfr));
subplot(212);imagesc(angle(tfr));

Wigner-Ville时频分布图：
[tfr,t,f]=tfrwv(x,t,ntrace)
例：
sig=fmlin(128,0.1,0.4);
tfrwv(sig);