instfreq

估计瞬时频率

描述

= instfreq(Xfs估计信号的瞬时频率,X在的速率采样fs。如果X是一个矩阵,则该函数独立地估计瞬时频率为每列并返回结果中的相应的列

= instfreq(XŤ估计的瞬时频率X在存储在所述时间值采样Ť

= instfreq(xt估计存储在MATLAB中的信号的瞬时频率®时间表xt。该功能将每一个独立变量里面的时间表所有变量和所有列。

= instfreq(TFDFDTD估计其时频分布为TFD,在存储在所述频率值被采样FD以及存储的时间值TD

= instfreq(___名称,值指定使用任何名字 - 值对参数的先前语法的附加选项。您可以指定用于估计瞬时频率或在计算中使用的频率范围的算法。

[Ť] = instfreq(___同样的回报Ť的采样时间对应于矢量

instfreq(___没有输出参数绘制了估计的瞬时频率。

例子

全部折叠

产生一个信号采样在5千赫为4秒。信号由一组持续时间逐渐减小的脉冲组成,脉冲之间由振荡幅度和波动频率区域分开,并有增大的趋势。画出信号。

FS = 5000;t = 0时:1 / FS:4-1 / FS;S = BESSELJ(0,1000 *(SIN(2 * PI * T ^ 2/8)^ 4)。);图(T,S)

以功率谱图的第一阶矩估计信号的时变频率。绘制功率谱图,叠加瞬时频率。

instfreq(S,FS)

产生复值信号,即由具有正弦变化频率内容的线性调频脉冲的。信号3kHz处1秒进行采样,并嵌入在白高斯噪声。

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

以功率谱图的第一阶矩估计信号的时变频率。这是唯一的方法instfreq万博1manbetx支持复值信号。绘制功率谱图,叠加瞬时频率。

instfreq(X,T)

创建一个双声道信号,在1千赫取样2秒,由两个电压控制振荡器的。在一个通道中,瞬时频率以时间为锯齿波,其最大值是在时间段的75%而变化。在其他通道中,瞬时频率与时间的方波的30%的占空比变化。画出两个通道的频谱。指定0.1秒为锯齿波信道和10Hz的频率分辨率为方形信道的时间分辨率。

FS = 1000;T =(0:1 / fs的:2)';X = VCO(锯齿(2 * PI *吨,0.75),[0.1 0.4] * FS,FS);Y = VCO(正方形(2 * PI * T,30),[0.1 0.3] * FS,FS);副区(1,2,1)pspectrum(X,FS,“频谱”'TimeResolution'0.1)副区(1,2,2)pspectrum(Y,FS,“频谱”'FrequencyResolution',10)

将信号存储在时间表中。计算并显示瞬时频率。

xt =时间表(秒(t), x, y);clf instfreq (xt)

重复使用解析信号计算。

instfreq (xt,'方法'希尔伯特的

生成由高斯函数调制的二次啁啾。指定采样速率为2khz,信号持续时间为4秒。

FS = 2000;t = 0时:1 / FS:4-1 / FS;Q =啁啾(叔1,0,1 / 2,20,“二次”,100,“凸”)。* exp (-1.7 * (2) ^ 2);情节(t, q)

使用pspectrum功能与默认设置估计信号的功率谱。利用估计值计算瞬时频率。

[P,F,T] = pspectrum(Q,FS,“频谱”);instfreq(P,F,T)

使用同步压缩的傅里叶变换重复计算。使用一个500样本的Hann窗口将信号分成若干段并对它们进行窗口处理。

[S,SF,ST] = FSST(Q,FS,汉恩(500));instfreq(ABS(S)^ 2,SF,ST)

比较两种方法的瞬时频率。

(psf, pst) = instfreq (p f t);(fsf,置)= instfreq (abs (s) ^ 2、科幻、st);浮置板轨道,情节(fsf、pst psf)

产生在1kHz 0.3秒采样,并嵌入在方差1/16的高斯白噪声的正弦信号。指定200Hz的正弦波的频率。估计和显示信号的瞬时频率。

FS = 1000;T =(0:1 / FS:0.3-1 / FS)';X = SIN(2 * PI * 200 * T)+ randn(大小(T))/ 4;instfreq(X,T)

再次估计信号的瞬时频率,但现在使用粗频率分辨率为25hz的时频分布作为输入。

[P,FD,TD] = pspectrum(X,T,“频谱”'FrequencyResolution',25);instfreq(P,FD,TD)

输入参数

全部折叠

输入信号,指定为矢量或矩阵。如果X是矢量,然后instfreq把它当作一个单独的频道。如果X是一个矩阵,然后instfreq独立地计算瞬时频率为每列并返回结果中的相应的列

例子:罪(2 *π* (0:127)/ 16)+ randn (1128) / 100指定嘈杂的正弦波

例子:[2 1]。* sin(2 *π* (0:127)”。/ [64])指定一个双通道正弦波。

数据类型:|
复数的支持:万博1manbetx

采样率,指定为正的标量。采样率是单位时间内采样的数量。如果时间单位为秒,则采样率为Hz。

数据类型:|

采样时间,指定作为一个真正的载体,持续时间标量,一持续时间阵列,或一个约会时间数组中。

  • 持续时间的连续样本之间的时间间隔X

  • 真正的向量,持续时间数组,或约会时间阵列 - 对应于每个元素的时刻X

例子:秒(1)指定信号的连续测量之间的1秒延时。

例子:秒(0:8)指定一个信号是在1赫兹采样持续8秒。

数据类型:||持续时间|约会时间

输入时间表。xt必须包含递增的、有限行的时间。

如果一个时间表丢失或重复的时间点,你可以使用提示修复清洁时间表与缺失,重复,或不均匀时报(MATLAB)。

例子:时间表(秒(0:4)”,randn(5,1))指定以1赫兹采样4秒的随机进程。

例子:时间表(秒(0:4)”,randn(5,3),randn(5,4))包含一个三通道随机过程和一个四通道随机过程,均以1赫兹采样4秒。

数据类型:|
复数的支持:万博1manbetx

时间频率分布,指定为在存储在所述频率采样的矩阵FD以及存储的时间值TD。此输入参数仅在以下情况下受支持万博1manbetx'方法'设定为'tfmoment'

例子:[P,F,T] = pspectrum(SIN(2 * PI *(0:511)/ 4),如图4所示, '频谱')指定一个1赫兹的时间 - 频率分布在4赫兹正弦波采样128秒,以及在其被计算出的频率和时间。

数据类型:|

对于时间 - 频率分布的频率和时间的值,指定为矢量。这些输入参数,只有当支持万博1manbetx'方法'设定为'tfmoment'

例子:[P,F,T] = pspectrum(SIN(2 * PI *(0:511)/ 4),如图4所示, '频谱')指定一个1赫兹的时间 - 频率分布在4赫兹正弦波采样128秒,以及在其被计算出的频率和时间。

数据类型:|

名称-值对的观点

指定可选的逗号分隔的对名称,值参数。名称参数名和是对应的值。名称必须出现在引号内。可以按任意顺序指定多个名称和值对参数名1,值1,...,NameN,值N

例子:'方法', 'tfmoment', 'FrequencyLimits',[25 50]通过找到时间 - 频率分布的第一条件的光谱时刻计算从25赫兹到50赫兹的范围内的输入的瞬时频率。

频率范围内,指定为逗号分隔的一对组成的“FrequencyLimits”和以Hz一个两元素矢量。如果没有指定,“FrequencyLimits”默认为(0 f / 2)为实值信号,并[fs / fs / 2)为复值信号。只有在以下情况下才支持此论点万博1manbetx'方法'设定为'tfmoment'

数据类型:|

计算方法,指定为逗号分隔的一对组成的'方法',要么'tfmoment'要么希尔伯特的

  • 'tfmoment'-计算瞬时频率作为时间-频率分布的第一个条件谱矩X。如果X是不是非均匀采样instfreq内插信号到一个均匀网格,以计算瞬时频率。

  • 希尔伯特的-计算瞬时频率,作为解析信号相位的导数X发现使用希尔伯特变换。该方法只接受均匀采样的实值信号,不支持时频分布输入。万博1manbetx

输出参数

全部折叠

瞬时频率,返回作为载体,矩阵,或具有相同的尺寸与输入的时间表。

频率估计的时间,恢复作为一个真正的载体,持续时间阵列,或一个约会时间数组中。

更多关于

全部折叠

瞬时频率

非平稳信号的瞬时频率是一个时变参数,它与信号演化过程中出现的频率平均值有关[1][2]

  • 如果'方法'设定为'tfmoment', 然后instfreq估计瞬时频率作为输入信号的时间 - 频率分布的第一条件的光谱的时刻。功能:

    1. 计算频谱功率谱PŤF中,输入用pspectrum函数,并使用频谱作为时间-频率分布。

    2. 使用估计的瞬时频率

      F 研究所 Ť = 0 F P Ť F d F 0 P Ť F d F

  • 如果'方法'设定为希尔伯特的, 然后instfreq将瞬时频率估计为输入解析信号相位的导数。功能:

    1. 计算解析信号,X一个中,输入用希尔伯特功能。

    2. 使用估计的瞬时频率

      F 研究所 Ť = 1 2 π d φ d Ť

      哪里φ是输入的分析信号的相位。

参考

[1] Boashash,布阿莱姆。“信号 - 第1部分的估计和解释瞬时频率:基本原理”。在IEEE论文集®。卷。80,1992年4月,第520-538。

[2] Boashash,布阿莱姆。“估计和信号 - 第2部分的解释瞬时频率:算法和应用。”在IEEE论文集。卷。80,1992年4月,第540-568。

介绍了在R2018a