主要内容

instfreq

估计瞬时频率

折叠页面

语法

Ifq = instfreq(x,fs)
Ifq = instfreq(x,t)
Ifq = instfreq(xt)
Ifq = instfreq(tfd,fd,td)
Ifq = instfreq()＿＿＿、名称、值)
[ifq,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千赫采样4秒的信号。信号由一组持续时间递减的脉冲组成，这些脉冲由振荡幅度和波动频率呈上升趋势的区域隔开。绘制信号。

Fs = 5000;T = 0:1/fs:4-1/fs;S = besselj(0,1000*(sin(2* *t.^2/8).^4));%要听，输入声音(s,fs)情节(t, s)

图包含一个轴对象。axes对象包含line类型的对象。

估计信号的时间相关频率作为功率谱图的第一矩。绘制功率谱图，叠加瞬时频率。

instfreq(年代,fs)

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为image, line的对象。该对象表示瞬时频率。

打开实时脚本

生成一个复值信号，该信号由具有正弦变化频率内容的啁啾组成。信号以3khz采样1秒，并嵌入高斯白噪声中。

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

估计信号的时间相关频率作为功率谱图的第一矩。这是唯一的方法instfreq金宝app支持复值信号。绘制功率谱图，叠加瞬时频率。

instfreq (x, t)

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为image, line的对象。该对象表示瞬时频率。

打开实时脚本

创建一个双通道信号，以1khz采样2秒，由两个电压控制振荡器组成。

  • 在一个通道中，瞬时频率以锯齿波的形式随时间变化，锯齿波的最大值为周期的75%。

  • 在另一个信道中，瞬时频率随时间变化为占空比为30%的方波。

绘制两个通道的频谱图。为锯齿通道指定0.1秒的时间分辨率，为方形通道指定10hz的频率分辨率。

Fs = 1000;T = (0:1/fs:2)';X = vco(锯齿(2*pi*t,0.75)，[0.1 0.4]*fs,fs);Y = vco(square(2*pi*t,30)，[0.1 0.3]*fs,fs);次要情节(1、2、1)pspectrum (x, fs,的谱图，“TimeResolution”，0.1) subplot(1,2,2) pspectrum(y,fs，的谱图，“FrequencyResolution”, 10)

图中包含2个轴对象。标题为Fres = 25.6675 Hz, Tres = 100 ms的轴对象1包含一个类型为image的对象。标题为Fres = 10.0263 Hz, Tres = 256 ms的轴对象2包含一个类型为image的对象。

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

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

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为line的对象。这些对象表示x\_instfreq, y\_instfreq。

用解析信号重复计算。

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

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为line的对象。这些对象表示x\_instfreq, y\_instfreq。

打开实时脚本

产生一个二次啁啾调制的高斯。指定采样率为2khz，信号持续时间为4秒。

Fs = 2000;T = 0:1/fs:4-1/fs;Q = chirp(t-1,0,1/2,20，“二次”, 100,“凸”)。* exp (-1.7 * (2) ^ 2);情节(t, q)

图包含一个轴对象。axes对象包含line类型的对象。

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

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

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为image, line的对象。该对象表示瞬时频率。

使用同步压缩傅里叶变换重复计算。使用500个样本的汉恩窗口将信号分成段并对其进行窗口处理。

[s,sf,st] = fsst(q,fs,hann(500));instfreq (abs (s) ^ 2,科幻,圣)

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为image, line的对象。该对象表示瞬时频率。

比较使用两种不同方法得到的瞬时频率。

[psf,pst] = instfreq(p,f,t);[fsf,fst] = instfreq(abs(s).^2,sf,st);浮置板轨道,情节(fsf、pst psf)

图包含一个轴对象。axes对象包含2个line类型的对象。

打开实时脚本

生成一个以1 kHz采样0.3秒的正弦信号，并嵌入方差为1/16的高斯白噪声中。指定正弦波频率为200hz。估计和显示瞬时频率的信号。

Fs = 1000;T = (0:1/fs:0.3-1/fs)';X = sin(2*pi*200*t) + randn(size(t))/4;instfreq (x, t)

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为image, line的对象。该对象表示瞬时频率。

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

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

图包含一个轴对象。标题为瞬时频率估计的轴对象包含2个类型为image, line的对象。该对象表示瞬时频率。

打开实时脚本

产生一个由啁啾组成的信号，其频率在300 Hz和1200 Hz之间呈正弦变化。信号以3khz采样2秒。

Fs = 3e3;T = 0:1/fs:2;Y = chirp(t, 100,1200，“二次”);Y = vco(cos(2*pi*t)，[0.1 0.4]*fs,fs);

使用instfreq计算信号的瞬时频率和相应的采样次数。验证输出对应于信号时频分布的集中一阶条件谱矩，计算公式为tfsmoment(预测性维护工具箱)。

[z,tz] = instfreq(y,fs);[a,ta] = tfsmoment(y,fs,1,Centralize=false);情节(tz, z,助教,“。”)传说(“instfreq”，“tfsmoment”）

图包含一个轴对象。axes对象包含2个line类型的对象。这些对象代表instfreq, tfmoment。

使用instbw计算信号的瞬时带宽和相应的采样时间。指定比例因子为1。验证输出对应于信号时间分布的非集中二阶条件谱矩的平方根。换句话说，instbw产生标准差和tfsmoment生成方差。

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

图包含一个轴对象。axes对象包含2个line类型的对象。这些对象代表instfreq, tfmoment。

输入参数

全部折叠

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

例子:罪(2 *π* (0:127)/ 16)+ randn (1128) / 100指定一个带噪声的正弦波

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

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

采样率，指定为正标量。样本率是单位时间内的样本数量。如果时间的单位是秒，那么采样率的单位是Hz。

数据类型:单|双

采样时间，指定为实向量a持续时间标量,持续时间数组，或者datetime数组中。

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

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

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

例子:秒(0:8)指定信号以1hz的频率采样8秒。

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

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

如果时间表中缺少或重复的时间点，您可以使用提示来修复它清除时间表中缺少、重复或不一致的时间。

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

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

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

时频分布，指定为在存储的频率上采样的矩阵fd以及存储的时间道明。只有在以下情况下才支持此输入参数金宝app“方法”设为“tfmoment”。

例子:[p f t] = pspectrum(罪(2 *π*(0:511)/ 4),4,谱图)指定以4hz采样128秒的1hz正弦波的时频分布，以及计算它的频率和时间。

数据类型:单|双

时间-频率分布的频率和时间值，指定为向量。只有在以下情况下才支持这些输入参数金宝app“方法”设为“tfmoment”。

例子:[p f t] = pspectrum(罪(2 *π*(0:511)/ 4),4,谱图)指定以4hz采样128秒的1hz正弦波的时频分布，以及计算它的频率和时间。

数据类型:单|双

名称-值参数

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

在R2021a之前，使用逗号分隔每个名称和值，并将其括起来名字在报价。

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

频率范围，指定为由逗号分隔的对组成“FrequencyLimits”以及以Hz为单位的双元素矢量。如果未指定，“FrequencyLimits”默认为(0 f / 2)对于实值信号和[fs / fs / 2)对于复值信号。只有在以下情况下才支持此论点金宝app“方法”设为“tfmoment”。

数据类型:单|双

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

  • “tfmoment”-计算瞬时频率作为时频分布的第一个条件谱矩x。如果x是非均匀抽样的instfreq将信号内插到均匀网格中以计算瞬时频率。

  • 希尔伯特的-计算瞬时频率作为解析信号的相位的导数x用希尔伯特变换得到的。该方法只接受均匀采样的实值信号，不支持时频分布输入。金宝app

输出参数

全部折叠

瞬时频率，作为与输入相同维度的矢量、矩阵或时间表返回。

乘以频率估计，作为实向量返回，a持续时间数组，或者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 ，

      在哪里ϕ是输入的解析信号的相位。

参考文献

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

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

扩展功能

版本历史

在R2018a中引入

另请参阅

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

主题