复啁啾的瞬时频率

打开实时脚本

这个例子展示了如何使用傅里叶同步压缩变换计算信号的瞬时频率。

产生具有正弦变化频率内容的啁啾。信号嵌入高斯白噪声中，以3千赫采样1秒。

Fs = 3000;T = 0:1/fs:1-1/fs;x = exp (2 j *π* 100 * cos(2 *π* 2 * t)) + randn(大小(t)) / 100;

计算并绘制信号的傅里叶同步压缩变换。上显示时间x-轴和频率在y设在。

fsst (x, fs,“桠溪”）

图中包含一个轴对象。标题为傅里叶同步压缩变换的坐标轴对象包含一个图像类型的对象。

通过提取傅里叶同步压缩变换的最大能量时频脊来找到信号的瞬时频率。

[sst,f,tfs] = fsst(x,fs);冰箱= tfridge(sst,f);

在变换图上叠加山脊。将时间转换为毫秒，频率转换为千赫。

持有在情节(t * 1000冰箱/ 1000,“r”)举行从

图中包含一个轴对象。标题为傅里叶同步压缩变换的坐标轴对象包含两个类型为image、line的对象。

对于一个真实的信号，你可以更容易地找到瞬时频率使用instfreq函数。例如，通过计算解析信号并求导其相位，显示复啁啾实部的瞬时频率。

Ax = real(x);instfreq (ax, fs,“方法”，希尔伯特的）

图中包含一个轴对象。标题为“瞬时频率估计”的axis对象包含一个类型为line的对象。

另请参阅