我有 129 个信号样本。我必须发现它是什么信号(它有什么组件)并尝试重建它。所以在我完成 FFT 之后,我想找到信号的相位,所以我尝试:
phs=unwrap(angle(y));
pha_wek=(0:length(y)-1)'/length(y)*129;
plot(pha_wek,phs)
结果我得到了它不是我所期望的。你能指出我有什么问题吗?
这个情节对我来说看起来完全正确,这取决于原始信号是什么。相位角的不连续性是由于幅度过零造成的。但是,相位跳变小于 pi(如果输入信号最初是真实的并且在应用 FFT 时已加窗,则这是您可能期望的)。所以我的猜测是,这个信号是某种形式的复杂调制,带有一些离散的频率。
通常,查看相位比信号幅度更能说明信号。它对于重建信号至关重要,但更难以解释(特别是因为它会随着时域的简单变化而发生剧烈变化)。
在不使用 unwrap 命令的情况下绘制相位可能会更好。在这种情况下,问题在于 unwrap 命令 - 它的目的是将相位跳跃限制为小于 pi,这就是为什么最终相位约为 -4pi 的原因。
我不确定你的挑战是什么,但不清楚你重建样本的意思。你已经有了这些,那有什么问题呢?(是否要求您重建插值样本)?