tfridge
时频脊
句法
描述
例子
查找噪声信号的岭
创建一个类似于用锋利的脊时频矩阵的矩阵。可视化的矩阵在三个维度。
t = 0时：0.05：10;F = 0：0.2：8;RV = 1;[F，T] = ndgrid（F，T）;S =零（大小（T））;S（ABS（（F-4）-cos（（T-6）^ 2））<0.1。）= RV;目（T，F，S）视图（-30,60）
噪声添加到矩阵，并重新显示情节。
S = S +兰特（大小（S））/ 10;目（T，F，S）视图（-30,60）xlabel（'时间'）ylabel（'频率'）
提取脊和绘制的结果。
[冰箱，〜，lridge] = tfridge（S，F）;rvals = S（lridge）;保持上plot3（吨，电冰箱，rvals，数k，'行宽'4）保持离
从多组信号中提取高能量脊
生成被3kHz处取样的一秒钟的信号。所述信号由两个色调和一个二次啁啾。
第一音具有1000Hz与单元的幅度的频率。
第二音调具有1200赫兹和单位幅度的频率。
该线性调频脉冲具有500赫兹的初始频率并达到在采样结束750赫兹。它有6个的幅度。
FS = 3000;t = 0时：1 / FS：1-1 / FS;X1 = 6 *啁啾（T，FS / 6，T（端部），FS / 4，“二次”）;X2 = SIN（2 * PI * FS / 3 * T）;X3 = SIN（2 * PI * FS / 2.5 * T）;X = X1 + X2 + X3;
计算并显示傅立叶变换synchrosqueezed的信号。
[SST中，f] = FSST（X，FS）;MX = MAX（ABS（SST（：）））*酮（大小（T））;目（T，F，ABS（SST））的视图（2）
提取并绘制两个最高能量的信号分量。设置没有惩罚改变频率。
penval = 0;冰箱= tfridge（SST，F，penval，'NumRidges'，2）;保持上plot3（吨，冰箱，MX，'W'，'行宽'5）保持离
两个音调具有相同的幅度，以及它们之间的算法跳转。设置改变频率为1的惩罚。
penval = 1;冰箱= tfridge（SST，F，penval，'NumRidges'，2）;目（T，F，ABS（SST））的视图（2）xlabel（'时间（s）'）ylabel（'频率（Hz）'）保持上plot3（吨，冰箱，MX，'W'，'行宽'5）保持离
刑罚设置为比较高的值。啁啾受到处罚，因为它的频率不恒定。
penval = 1000;冰箱= tfridge（SST，F，penval，'NumRidges'，2）;目（T，F，ABS（SST））的视图（2）xlabel（'时间（s）'）ylabel（'频率（Hz）'）保持上plot3（吨，冰箱，MX，'W'，'行宽'5）保持离
邻居单元去除的影响。
生成两个二次啁啾组成的信号。该信号在1 kHz采样为3秒。的线性调频脉冲是这样的瞬时频率是对称采样间隔的中间点。一个线性调频脉冲是凹的，而其他的啁啾是凸的。凹啁啾具有凸线性调频脉冲的振幅的两倍。
FS = 1E3;t = 0时：1 / fs的：3;X =啁啾（叔1.5,100,1.1,200，“二次”，[]“凸”）;Y = 2 *啁啾（叔1.5,300,1.1,400，“二次”，[]'凹'）;％听到，类型soundsc（X + Y，FS）
计算并显示傅立叶变换synchrosqueezed的信号。
SIG = X + Y;[SST，F，T] = FSST（SIG，FS）;FSST（SIG，FS，“Y轴”）
提取具有最高的能量两个时频脊。指定频率的变化为1的损失。提取第二脊前取出与周围的最高能量的脊1的频点。画出脊。
NFB = 1;[FR，IR] = tfridge（SST中，f，1，'NumRidges'，2，'NumFrequencyBins'，NFB）;情节（T，FR）
一个二进制位是不够的：该函数查找第二凸脊是部分地在所述第一的斜率。提高到50仓的数量，以除去和重复计算。
NFB = 50;[FR，IR] = tfridge（SST中，f，1，'NumRidges'，2，'NumFrequencyBins'，NFB）;情节（T，FR）
卸下太多箱扭曲降低能源脊。减少的数量为15，重复计算。
NFB = 15;[FR，IR] = tfridge（SST中，f，1，'NumRidges'，2，'NumFrequencyBins'，NFB）;情节（T，FR）
反变换对应于两个脊。添加脊重建信号。画出重构信号与啁啾之间的差异。
ITR = ifsst（SST，[]，IR，'NumFrequencyBins'，NFB）;xrec =总和（ITR'）;情节（T，xrec-（X + Y））ylim（[ - 1 0.1]）
％听到，类型soundsc（xrec，FS）
该协议是很好的大部分在两端，因为后者的频率变化最迅速的时间，但恶化。
输入参数
输出参数
算法
该函数使用一个惩罚前后贪婪算法来提取一个时间 - 频率矩阵的最大能量脊。该算法通过最小化发现的最大的时频脊-ln一个在每个时间点，在那里一个是矩阵的绝对值。最小化-ln一个等同于最大化的值一个。该算法可选地约束在频率跳跃带罚成比例的频率窗口之间的距离。
下面的示例示出了使用惩罚是两倍频率仓之间的距离的时间 - 频率脊算法。具体而言，元件之间的距离（J，K）和（M，N）被定义为（J-M）2。时间 - 频率矩阵具有三个频率仓和三个时间步长。矩阵的列对应于时间步骤，并且矩阵的行对应于频率仓。在第二行中的值表示的正弦波。
假设你有矩阵：
1 4 4 2 2 2 5 5 4
更新值用于（1,2）元素如下。
在第一时间点不变保留值。与基质，其中介绍在第二时间点在第一频率上的（1,2）元素开始的算法。的bin值是4惩罚在基于从（1,2）元素的距离的第一列中的值。运用刑罚的第一列产生
原始值+罚×距离1 + 2×0 = 1 2 + 2×1 = 4 5 + 2×4 = 13
第一列的最小值为1，这是在仓1。1 4 4 2 13 5
在第1列添加最小值到当前二进制值，4.（1,2）的更新后的值变为5，其从箱1来了。
更新用于在列中的其余元件2如下的值。
重新计算原始列1点的值使用相同的处理步骤2a中的惩罚因子。获得使用相同的过程剩余的第二列的值作为在步骤2b。例如，在更新（2,2）元件，其具有单元条值时2中，施加惩罚到列产量
添加最小值，2，当前二进制值。为（2,2）的更新值变为4更新（3,2）元素之后，矩阵是原始值+罚×距离1 + 2×1 = 3 2 + 2×0 = 2 5 + 2×1 = 7
只有第二栏已更新。该标表示斌在从一个值来上一列的索引。1 5（1）4 2 4（2）2 5 9（2）4
将第三列重复执行步骤2。但是，现在的处罚适用于更新的第二列。例如，在更新（1,3）元素时，该惩罚是
最小值，如图5所示，这是在第一仓，被添加到（1,3）的合并值。在第三列中更新所有的值后，最终的基体是5 + 2×0 = 5 4 + 2×1 = 6 9 + 2×4 = 17
1 5（1）9（1）2 4（2）6（2）5月9日（2）10（2）
在矩阵的最后一列开始，找到最小值。通过在前面的时间点，从目前的仓去那个箱子的起源穿行矩阵时间回来。跟踪箱索引，其形成通道组成的山脊。该算法通过使用起源箱，而不是与最小值的bin平滑的过渡。在这个例子中，脊指数
2，2，2，其在步骤1中所示的矩阵的第2行的正弦波的能量路径相匹配。
如果要提取的多个脊，所述算法删除从时间 - 频率矩阵的第一个凸脊和重复该过程。
扩展功能
介绍了在R2016b
您还可以选择从下面的列表中的网站：
