主要内容

instfreq

估计瞬时频率

所有的页面崩溃

语法

不= instfreq (x, fs)
不= instfreq (x, t)
不= instfreq (xt)
不= instfreq (tfd、fd td)
不= instfreq (___、名称、值)
[不,t] = instfreq (___)
instfreq (___)

描述

例子

不= instfreq (x,fs)估计信号的瞬时频率,x采样率fs。如果x是一个矩阵,然后独立函数估计瞬时频率为每一列对应的列和返回结果的不。

例子

不= instfreq (x,t)瞬时频率的估计x采样值存储在t。

例子

不= instfreq (xt)估计信号的瞬时频率存储在MATLAB中®时间表xt。函数对所有变量的时间表和所有列在每个独立变量。

例子

不= instfreq (tfd,fd,道明)估计信号的瞬时频率的时频分布,tfd,在频率采样值存储在fd中存储的时间值道明。

例子

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

例子

(不,t)= instfreq (___)同样的回报t样本的向量乘以对应不。

例子

instfreq (___)没有输出参数情节估计瞬时频率。

例子

全部折叠

打开生活的脚本

生成一个信号采样5 kHz 4秒。信号由一组脉冲减少持续时间由区域的振动振幅和频率波动增加的趋势。画出信号。

fs = 5000;t = 0:1 / fs: 4 - 1 / f;s = besselj(0, 1000 *(罪(2 *π* t。^ 2/8)。^ 4));%听,声音类型(年代,fs)情节(t, s)

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象。

估计的时间频率信号的一阶矩功率谱图。情节瞬时频率的功率谱图和叠加。

instfreq(年代,fs)

图包含一个坐标轴对象。坐标轴对象2标题瞬时频率估计包含对象类型的形象,线。这个对象表示瞬时频率。

打开生活的脚本

生成一个啁啾与呈现正弦的复值信号,由不同频率的内容。信号采样在1秒3千赫,嵌入在高斯白噪声。

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

估计的时间频率信号的一阶矩功率谱图。这是唯一的方法instfreq金宝app为复值信号。情节瞬时频率的功率谱图和叠加。

instfreq (x, t)

图包含一个坐标轴对象。坐标轴对象2标题瞬时频率估计包含对象类型的形象,线。这个对象表示瞬时频率。

打开生活的脚本

创建一个双通道信号,采样1 kHz 2秒钟,两个压控振荡器组成。

  • 在一个频道,瞬时频率随时间作为一个锯齿波的最大周期的75%。

  • 在其他频道,瞬时频率随时间一个方波占空比为30%。

图的谱图两个渠道。指定一个0.1秒的时间分辨率锯齿通道和10赫兹的频率分辨率广场通道。

fs = 1000;t = (0:1 / fs: 2) ';x = vco(锯齿(2 *π* t, 0.75), (0.1 - 0.4) * fs, fs);y = vco(广场(2 *π* 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)

图包含2轴对象。坐标轴对象1标题耐火的= 25.6675赫兹,女士非常= 100包含一个类型的对象的形象。坐标轴对象2标题耐火的= 10.0263赫兹,女士非常= 256包含一个类型的对象的形象。

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

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

图包含一个坐标轴对象。坐标轴对象与标题瞬时频率估计包含2线类型的对象。这些对象代表x \ _instfreq y \ _instfreq。

用解析信号重复计算。

instfreq (xt,“方法”,希尔伯特的)

图包含一个坐标轴对象。坐标轴对象与标题瞬时频率估计包含2线类型的对象。这些对象代表x \ _instfreq y \ _instfreq。

打开生活的脚本

生成一个二次啁啾高斯调制。指定2 kHz的采样率和信号持续时间4秒。

fs = 2000;t = 0:1 / fs: 4 - 1 / f;q =唧唧声(t - 1, 0, 1/2, 20日“二次”,100,“凸”)。* exp (-1.7 * (2) ^ 2);情节(t, q)

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象。

使用pspectrum与默认设置函数来估计信号的功率谱。使用计算瞬时频率的估计。

[p f t] = pspectrum (q, fs,的谱图);instfreq (p f t)

图包含一个坐标轴对象。坐标轴对象2标题瞬时频率估计包含对象类型的形象,线。这个对象表示瞬时频率。

重复使用synchrosqueezed傅里叶变换的计算。使用500 -样本损害窗口把信号分成段和窗口。

[年代,科幻,st] = fsst (q, fs,损害(500));instfreq (abs (s) ^ 2,科幻,圣)

图包含一个坐标轴对象。坐标轴对象2标题瞬时频率估计包含对象类型的形象,线。这个对象表示瞬时频率。

瞬时频率进行比较发现使用两种不同的方法。

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

图包含一个坐标轴对象。坐标轴对象包含2线类型的对象。

打开生活的脚本

生成一个正弦信号采样1 kHz 0.3秒和嵌入在高斯白噪声方差的1/16。指定一个正弦信号的频率200赫兹。估计和显示信号的瞬时频率。

fs = 1000;t = 0.3 (0:1 / fs: 1 / fs) ';x =罪(2 *π* 200 * t) + randn(大小(t)) / 4;instfreq (x, t)

图包含一个坐标轴对象。坐标轴对象2标题瞬时频率估计包含对象类型的形象,线。这个对象表示瞬时频率。

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

[p fd td] = pspectrum (x, t,的谱图,“FrequencyResolution”25);instfreq (p fd td)

图包含一个坐标轴对象。坐标轴对象2标题瞬时频率估计包含对象类型的形象,线。这个对象表示瞬时频率。

打开生活的脚本

产生一个信号,由频率啁啾的变化呈现正弦300赫兹到1200赫兹。信号在3千赫采样2秒。

fs = 3 e3;t = 0:1 / fs: 2;y =唧唧声(t, 100, 1200,“二次”);y = vco (cos(2 *π* t) (0.1 - 0.4) * fs, fs);

使用instfreq计算信号的瞬时频率和相应的样品时间。验证输出对应于集中的一阶条件光谱信号的时频分布的计算tfsmoment(预测维护工具箱)。

(z, tz) = instfreq (y, fs);[一个,ta] = tfsmoment (y, fs,集中= false);情节(tz, z,助教,“。”)传说(“instfreq”,“tfsmoment”)

图包含一个坐标轴对象。坐标轴对象包含2线类型的对象。这些对象代表instfreq tfsmoment。

使用instbw计算信号的瞬时带宽和相应的样品时间。指定一个比例因子为1。验证输出对应的平方根noncentralized二阶条件光谱信号的时间分布的时刻。换句话说,instbw生成一个标准差和tfsmoment生成一个方差。

[w, tw] = instbw (y, fs, ScaleFactor = 1);[m, tm] = tfsmoment (y, fs, 2);情节(tw, w, tm, sqrt (m),“。”)传说(“instfreq”,“tfsmoment”)

图包含一个坐标轴对象。坐标轴对象包含2线类型的对象。这些对象代表instfreq tfsmoment。

输入参数

全部折叠

输入信号,指定为一个向量或矩阵。如果x是一个矢量,然后呢instfreq将它视为一个单一的通道。如果x是一个矩阵,然后呢instfreq独立计算瞬时频率为每一列,并返回对应的列的结果不。

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

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

数据类型:单|双
复数的支持:金宝app是的

采样率,指定为一个积极的标量。采样率是单位时间内样品的数量。如果时间的单位是秒,然后采样率在赫兹。

数据类型:单|双

样品时间,指定为一个真正的向量,持续时间标量,持续时间数组,或datetime数组中。

  • 持续时间标量-连续采样之间的时间间隔x。

  • 真正的向量,持续时间数组,或datetime数组——即时的每个元素对应的时间x。

例子:秒(1)指定了一个秒的连续测量信号之间的失误。

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

数据类型:单|双|持续时间|datetime

输入时间表。xt必须包含增加,有限的行乘以。

如果时间表有缺失或重复的时间点,你可以修理它使用的技巧清洁时间表与失踪、重复或不均匀。

例子:时间表(秒(0:4)”,randn (1))指定了一个随机过程样本1赫兹为4秒。

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

数据类型:单|双
复数的支持:金宝app是的

时频分布,指定为一个矩阵在存储在频率采样fd中存储的时间值道明。只有当这个输入参数支持金宝app“方法”被设置为“tfmoment”。

例子:[p f t] = pspectrum(罪(2 *π*(0:511)/ 4),4,谱图)指定1 Hz正弦信号的时频分布采样4赫兹为128秒,而且计算的频率和时间。

数据类型:单|双

时频分布、频率和时间值指定为向量。只有当这些输入参数支持金宝app“方法”被设置为“tfmoment”。

例子:[p f t] = pspectrum(罪(2 *π*(0:511)/ 4),4,谱图)指定1 Hz正弦信号的时频分布采样4赫兹为128秒,而且计算的频率和时间。

数据类型:单|双

名称-值参数

指定可选的双参数作为Name1 = Value1,…,以=家,在那里的名字参数名称和吗价值相应的价值。名称-值参数必须出现在其他参数,但对的顺序无关紧要。

R2021a之前,用逗号来分隔每一个名称和值,并附上的名字在报价。

例子:“方法”、“tfmoment’,‘FrequencyLimits’, 50 [25]计算瞬时频率的输入范围从25 Hz至50 Hz找到第一个条件时频分布谱的时刻。

频率范围,指定为逗号分隔组成的“FrequencyLimits”双元素向量和赫兹。如果不指定,“FrequencyLimits”默认为(0 f / 2)并为实值信号[fs / fs / 2)为复值信号。只有当这个论点是支持金宝app“方法”被设置为“tfmoment”。

数据类型:单|双

计算方法,指定为逗号分隔组成的“方法”,要么“tfmoment”或希尔伯特的。

  • “tfmoment”——计算瞬时频率作为第一个条件的时频分布谱的时刻x。如果x非均匀采样,然后instfreq篡改的信号均匀网格计算瞬时频率。

  • 希尔伯特的——计算瞬时频率的导数的阶段分析的信号x发现利用希尔伯特变换。这个方法只接受均匀采样,实值信号和输入不支持时频分布。金宝app

输出参数

全部折叠

瞬时频率,作为一个矢量,返回一个矩阵,或相同的维数作为输入的时间表。

频率估计的时候,作为一个真正的返回向量,一个持续时间数组,或datetime数组中。

更多关于

全部折叠

瞬时频率

的瞬时频率的非平稳信号是时变的参数与频率的平均值信号中,因为它的发展[1],[2]。

  • 如果“方法”被设置为“tfmoment”,然后instfreq估计瞬时频率作为第一个条件谱输入信号的时频分布的时刻。功能:

    1. 计算光谱图功率谱P(t,f)的输入使用pspectrum功能和使用频谱作为时频分布。

    2. 估计瞬时频率使用

      f 本月 ( t ) = ∫ 0 ∞ f P ( t , f ) d f ∫ 0 ∞ P ( t , f ) d f 。

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

    1. 计算分析信号,x一个输入使用希尔伯特函数。

    2. 估计瞬时频率使用

      f 本月 ( t ) = 1 2 π d ϕ d t ,

      在哪里ϕ的阶段分析信号的输入。

引用

[1]Boashash Boualem。“估计和解释一个信号的瞬时频率。即基本面因素。”IEEE学报》®80年,没有。4(1992年4月):520 - 538。https://doi.org/10.1109/5.135376。

[2]Boashash Boualem。“估计和解释一个信号的瞬时频率。二世。算法和应用程序”。IEEE学报》80年,没有。4(1992年5月):540 - 568。https://doi.org/10.1109/5.135378。

扩展功能

版本历史

介绍了R2018a

另请参阅

|||(预测维护工具箱)|(预测维护工具箱)|(预测维护工具箱)

主题