文档帮助中心
估计瞬时频率
不= 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,fs)
不
X
fs
不= instfreq(X,Ť)估计的瞬时频率X在存储在所述时间值采样Ť。
不= instfreq(X,Ť)
Ť
不= instfreq(xt)估计存储在MATLAB中的信号的瞬时频率®时间表xt。该功能将每一个独立变量里面的时间表所有变量和所有列。
不= instfreq(xt)
xt
不= instfreq(TFD,FD,TD)估计其时频分布为TFD,在存储在所述频率值被采样FD以及存储的时间值TD。
不= instfreq(TFD,FD,TD)
TFD
FD
TD
不= instfreq(___,名称,值)指定使用任何名字 - 值对参数的先前语法的附加选项。您可以指定用于估计瞬时频率或在计算中使用的频率范围的算法。
不= instfreq(___,名称,值)
名称,值
[不,Ť] = 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
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功能与默认设置估计信号的功率谱。利用估计值计算瞬时频率。
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 *π* (0:127)/ 16)+ randn (1128) / 100
例子:[2 1]。* sin(2 *π* (0:127)”。/ [64])指定一个双通道正弦波。
[2 1]。* sin(2 *π* (0:127)”。/ [64])
数据类型:单|双复数的支持:万博1manbetx是
单
双
采样率,指定为正的标量。采样率是单位时间内采样的数量。如果时间单位为秒,则采样率为Hz。
数据类型:单|双
持续时间
约会时间
采样时间,指定作为一个真正的载体,持续时间标量,一持续时间阵列,或一个约会时间数组中。
持续时间的连续样本之间的时间间隔X。
真正的向量,持续时间数组,或约会时间阵列 - 对应于每个元素的时刻X。
例子:秒(1)指定信号的连续测量之间的1秒延时。
秒(1)
例子:秒(0:8)指定一个信号是在1赫兹采样持续8秒。
秒(0:8)
数据类型:单|双|持续时间|约会时间
输入时间表。xt必须包含递增的、有限行的时间。
如果一个时间表丢失或重复的时间点,你可以使用提示修复清洁时间表与缺失,重复,或不均匀时报(MATLAB)。
例子:时间表(秒(0:4)”,randn(5,1))指定以1赫兹采样4秒的随机进程。
时间表(秒(0:4)”,randn(5,1))
例子:时间表(秒(0:4)”,randn(5,3),randn(5,4))包含一个三通道随机过程和一个四通道随机过程,均以1赫兹采样4秒。
时间表(秒(0:4)”,randn(5,3),randn(5,4))
时间频率分布,指定为在存储在所述频率采样的矩阵FD以及存储的时间值TD。此输入参数仅在以下情况下受支持万博1manbetx'方法'设定为'tfmoment'。
'方法'
'tfmoment'
例子:[P,F,T] = pspectrum(SIN(2 * PI *(0:511)/ 4),如图4所示, '频谱')指定一个1赫兹的时间 - 频率分布在4赫兹正弦波采样128秒,以及在其被计算出的频率和时间。
[P,F,T] = pspectrum(SIN(2 * PI *(0:511)/ 4),如图4所示, '频谱')
对于时间 - 频率分布的频率和时间的值,指定为矢量。这些输入参数,只有当支持万博1manbetx'方法'设定为'tfmoment'。
指定可选的逗号分隔的对名称,值参数。名称参数名和值是对应的值。名称必须出现在引号内。可以按任意顺序指定多个名称和值对参数名1,值1,...,NameN,值N。
名称
值
名1,值1,...,NameN,值N
'方法', 'tfmoment', 'FrequencyLimits',[25 50]
“FrequencyLimits”
(0 f / 2)
[fs / fs / 2)
频率范围内,指定为逗号分隔的一对组成的“FrequencyLimits”和以Hz一个两元素矢量。如果没有指定,“FrequencyLimits”默认为(0 f / 2)为实值信号,并[fs / fs / 2)为复值信号。只有在以下情况下才支持此论点万博1manbetx'方法'设定为'tfmoment'。
希尔伯特的
计算方法,指定为逗号分隔的一对组成的'方法',要么'tfmoment'要么希尔伯特的。
'tfmoment'-计算瞬时频率作为时间-频率分布的第一个条件谱矩X。如果X是不是非均匀采样instfreq内插信号到一个均匀网格,以计算瞬时频率。
希尔伯特的-计算瞬时频率,作为解析信号相位的导数X发现使用希尔伯特变换。该方法只接受均匀采样的实值信号,不支持时频分布输入。万博1manbetx
瞬时频率,返回作为载体,矩阵,或具有相同的尺寸与输入的时间表。
频率估计的时间,恢复作为一个真正的载体,持续时间阵列,或一个约会时间数组中。
非平稳信号的瞬时频率是一个时变参数,它与信号演化过程中出现的频率平均值有关[1],[2]。
如果'方法'设定为'tfmoment', 然后instfreq估计瞬时频率作为输入信号的时间 - 频率分布的第一条件的光谱的时刻。功能:
计算频谱功率谱P(Ť,F)中,输入用pspectrum函数,并使用频谱作为时间-频率分布。
使用估计的瞬时频率
F 研究所 ( Ť ) = ∫ 0 ∞ F P ( Ť , F ) d F ∫ 0 ∞ P ( Ť , F ) d F 。
如果'方法'设定为希尔伯特的, 然后instfreq将瞬时频率估计为输入解析信号相位的导数。功能:
计算解析信号,X一个中,输入用希尔伯特功能。
希尔伯特
F 研究所 ( Ť ) = 1 2 π d φ d Ť ,
哪里φ是输入的分析信号的相位。
[1] Boashash,布阿莱姆。“信号 - 第1部分的估计和解释瞬时频率:基本原理”。在IEEE论文集®。卷。80,1992年4月,第520-538。
[2] Boashash,布阿莱姆。“估计和信号 - 第2部分的解释瞬时频率:算法和应用。”在IEEE论文集。卷。80,1992年4月,第540-568。
希尔伯特|pspectrum
您的系统中存在此示例的修改版本。你想打开这个版本吗?
你点击了一个链接,对应于这个MATLAB命令:
通过在MATLAB命令窗口中输入该命令来运行它。Web浏览器不支持MATLAB命令。万博1manbetx
选择一个网站,以获得翻译的内容,其中可看到当地的活动和优惠。根据您的位置,我们建议您选择:。
您还可以选择从下面的列表中的网站:
选择最佳的网站性能的中国网站(在中国或英文)。其他MathWorks的国家网站都没有从您的位置访问进行了优化。
联系你当地的办公室