时频特性分析（Matlab）

最近瞅了一些关于时频分析工具箱的matlab函数使用方法,总结一下吧.

我使用的是2011a的matlab,貌似没有自带的时频分析工具箱,可以到网上下载,google一搜就能搜到,安装后就可以使用了(所谓安装就是把工具箱的目录包含到matlab工作目录中即可).

下面说一些时频工具箱函数的用法(由于我下的工具箱没有html版的帮助,所以只能使用helpfunname的方法查看帮助信息).

AMGAUSS

Generate gaussian amplitude modulation.

Y=AMGAUSS(N,T0,T) generates a gaussian amplitude modulation

centered on a time T0, and with a spread proportional to T.

This modulation is scaled such that Y(T0)=1

and Y(T0+T/2) and Y(T0-T/2) are approximately equal to 0.5 .

产生高斯幅值调制信号，其以信号点数为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(160,90,40); plot(z);

 

FMCONST

Signal with constant frequency modulation.

[Y,IFLAW] = FMCONST(N,FNORM,T0) generates a frequencymodulation

with a constant frequency fnorm.

The phase of this modulation is such that y(t0)=1.

产生一个固定频率的频率调制信号。N为产生信号的点数，FNORM为标准化频率（默认为0.25），T0表示此时刻为正幅值（默认为round(N/2)），相当于规定了相位。IFLAW为设置的频率变化情况，这里就是一个常数。

如z=fmconst(128,0.05,50); plot(real(z));

 

FMLIN

Signal with linear frequency modulation.

[Y,IFLAW]=FMLIN(N,FNORMI,FNORMF,T0) generates a linearfrequency

modulation.

The phase of this modulation is such that Y(T0)=1.

产生一个线性调频信号。N为信号的点数，FNORMI为开始标准频率（默认为0.0），FNORMF为结束标准频率（默认为0.5），相位由T0规定，即Y(T0)=1。IFLAW为设置的频率变化情况，这里就是一个线性变化的直线。

如[z, f]=fmlin(128,0.05,0.3,50); plot(real(z));figure,plot(f)

以上是比较常用的产生信号的函数。

 

下面说一些时频分析常用的分析方法，主要有短时傅里叶变换STFT、Wigner-Ville分布WVD、伪Wigner-Ville分布PWVD等，其他还有很多分析方法，都有相关函数，具体可以查看帮助。

短时傅里叶变换STFT

功能：计算时间序列的短时离散傅里叶变换，得到瞬时频率。

格式：

[tfr, t, f] = tfrstft(x) %计算时间序列x的短时傅里叶变换，参数tfr为短时傅里叶变换系数，t为系数tfr对应的时刻，f为归一化频率向量

[tfr, t, f] = tfrstft(x, t) % 计算对应时刻t的短时傅里叶变换

[tfr, t, f] = tfrstft(x, t, n) % 计算n点对应时刻t的短时傅里叶变换

[tfr, t, f] = tfrstft(x, t, n, h) % 参数h为归一化频率平滑窗

[tfr, t, f] = tfrstft(x t, n, h, trace) % 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时频分布图

功能：计算时间序列的Wigner-Ville时频分布图，得到瞬时频率

格式：

[tfr, t, f] = tfrwv(x)

[tfr, t, f] = tfrwv(x, t)

[tfr, t, f] = tfrwv(x, t, n)

[tfr, t, f] = tfrwv(x, t, n, trace)

说明：

x--信号

t--时间（缺省值为1：length(x))

n--频率数（缺省值为length(x))

trace--如果非零，显示算法的进程（缺省值为0）

tfr--时频分解（为复值），频率轴观察范围为-0.5~0.5

f--归一化频率

如

sig=fmlin(128,0.1,0.4);

tfrwv(sig);

 

伪Wigner-Ville时频分布图

功能：计算时间序列的伪Wigner-Ville时频分布图，得到瞬时频率

格式：

[tfr, t, f] = tfrpwv(x)

[tfr, t, f] = tfrpwv(x, t)

[tfr, t, f] = tfrpwv(x, t, n)

[tfr, t, f] = tfrpwv(x, t, n, trace)

说明：

x--信号

t--时间（缺省值为1：length(x))

n--频率数（缺省值为length(x))

trace--如果非零，显示算法的进程（缺省值为0）

tfr--时频分解（为复值），频率轴观察范围为-0.5~0.5

f--归一化频率

如

sig=fmlin(128,0.1,0.4);

tfrpwv(sig);

还可以设置成三维显示，如下

其他时频分析

fmt

计算时间序列的快速梅林变换

ifmt

计算快速梅林逆变换。

tfrbj

生成时间序列的Bord-Jondan时频分布图，得到瞬时频率。

tfrbud

生成时间序列的Butterworth时频分布图，得到瞬时频率。

tfrcw

生成时间序列的Choi-Williams时频分布图，得到瞬时频率。

tfrgrd

生成时间序列的广义矩形时频分布图，得到瞬时频率。

tfrmh

生成时间序列的Margenau-Hill时频分布图，得到瞬时频率。

tfrmhs

生成时间序列的Margenau-HillP谱图时频分布图，得到瞬时频率。

tfrsp

计算时间序列的频谱图分布，得到瞬时频率。

tfrscalo

计算时间序列的尺度图时频表示，即连续小波变换的幅值的平方。

tfrrsp

计算时间序列的时频分布和它的重排形式，得到瞬时频率。

friedman

计算时间序列的瞬时频率密度。

ridges

从时间序列的重排时频分布中提取瘠。

tfrideal

计算给定时间序列的瞬时频率规律的理想时频表示。

以上给了很多函数，但常用的基本就是短时傅里叶变换STFT、Wigner-Ville分布WVD、伪Wigner-Ville分布PWVD和提取瘠ridges、得到理想时频表示tfrideal。