复杂chirp的瞬时频率

打开实时脚本

此示例显示了如何使用傅立叶同步转换来计算信号的瞬时频率。

产生带有正弦频率含量的呼叫。该信号嵌入白色高斯噪声中，并在3 kHz中采样1秒。

FS = 3000;t = 0：1/fs：1-1/fs;x = exp（2J*pi*100*cos（2*pi*2*t）） + randn（size（t））/100;

计算和绘制信号的傅立叶同步转换。在X- 轴和频率y-轴。

fsst（x，fs，'yaxis'）

图包含一个轴对象。带有标题傅里叶同步转换的轴对象包含类型图像的对象。

通过提取傅立叶同步转换的最大能源时频脊来找到信号的瞬时频率。

[sst，f，tfs] = fsst（x，fs）;冰箱= Tfridge（SST，F）;

覆盖转换图上的山脊。将时间转换为毫秒和频率为kHz。

抓住在情节（t*1000，冰箱/1000，'r'） 抓住离开

图包含一个轴对象。带有标题傅里叶同步转换的轴对象包含2个类型图像的对象。

对于真实信号，您可以使用Instfreq功能。例如，通过计算分析信号并区分其相位来显示复合chirp的实际部分的瞬时频率。

ax = real（x）;Instfreq（ax，fs，'方法'，，，，“希尔伯特”）

图包含一个轴对象。带有标题瞬时频率估计的轴对象包含类型线的对象。

也可以看看