函数的频谱分析

信息处理 傅里叶变换 频谱
2022-02-19 20:03:21

我在 Math SE 上问过这个问题,但没有收到任何回复。我希望这是一个相关的论坛来提问。

我想分析以下函数的频谱:

f(t)=cos(ta(1+bcos(ct)))
a,b,c一些实际参数。由于该函数是非周期性的(也许除非a,b,c是微调),最自然的方法是对其进行傅里叶变换。但是手工做是相当可怕的,Wolfram Alpha aso 很挣扎,使用 Mathematica 我得到了
iacos(2π)δ[ω]abccos2(2π)δ[ω]
我只对结果的真实部分感兴趣。我现在正试图弄清楚我得到了什么。

1) Mathematica 结果对您来说是否合理?

2)从生成器函数的角度来看(如here),似乎有一个以零为中心的尖峰,在一些处伴随着更多的小峰。我发现这个结果令人惊讶,因为我预计频谱将以为中心,作为“主”频率,然后是一些较小的谐波或某种类似的函数。±ΩΩ>0asinc

3) 遵循比完整傅里叶变换简单得多的参数,我尝试将其视为具有在“时间”( ) 中振荡的频率的函数。然后也许作为一个类比,我可以问自己哪些频率最有可能。在这种情况下,我可以将项视为均匀随机变量的余弦,因此概率密度函数实际上是 U 形的以的一般形式的函数给出(归一化等)tcos(ct)ct:=YUni(π,π)a

fY(x)=1π1x2

但这两种分析相距甚远,我想了解发生了什么。

对我的概率解释的一个反驳是,决定频率分量强度的项是傅立叶分解的系数(如果它是周期性的),这里这不是具有不同概率的频率的线性组合,而是一些别的。不过,我需要一些帮助来解释结果。谢谢!

1个回答

最外层余弦函数的相位项是随时间变化的:瞬时频率由特别是,在附近的一个很小的邻域中,我们可以认为只有一个频率分量在较大的值下,项获得控制。ϕ(t):=ta(1+bcos(ct))ϕt=a+abcos(ct)abctsin(ct).t=0f(t)a+abtabctsin(ct)

请注意,频率随时间而变化。这使得很难直观地解释假设所有频率始终存在的傅立叶变换。对于此类信号,最好通过查看的窗口段在不同值周围的频谱来使用时频分析。f(t)t