问题标签 [spectral-density]

让我们考虑以下代码

正如我所知，傅里叶域中的峰值代表这个频率，它们存在于信号中，例如让我们绘制这个信号的傅里叶变换图让我们运行这个信号

和情节

但很明显，它没有显示信号中给定频率的峰值，有什么问题？请帮助我，通常它是固定信号，为什么这张图应该显示准确的图片但它没有，为什么？

编辑：

y1=generate1(3,500,0);

我有两个信号，我希望其中一个对另一个做出响应，但有一定的相移。

现在我想计算相干性或归一化交叉谱密度，以估计输入和输出之间是否存在任何因果关系，以找出这种相干性出现在哪些频率上。

例如，请参阅此图像（来自此处），它似乎在频率 10 处具有高相干性：

现在我知道我可以使用互相关来计算两个信号的相移，但是如何使用相干性（频率为 10）来计算相移？

图像代码：

.
.
编辑：

对于它的价值，我已经添加了一个答案，也许是对的，也许是错的。我不确定..

我正在尝试生成 100 个随机位（1，-1）的数据包（每个大小为 8192 的行向量），使用巴特沃斯滤波器对它们进行过滤，然后绘制它们的平均功率谱密度。我必须使用 MATLAB 来执行此操作。我的 o/p 应该是具有非常尖锐峰值的滤波 sinc。当我将此代码用于较小尺寸的数据包时，它说 100 它可以工作。但对于 8192 则不然。我希望有人检查我的代码是否有错误。

我试图找到在不均匀时间测量的信号的功率谱密度。数据看起来像这样：

其中第一列是自第一次测量以来的时间（以秒为单位），第二列是测量值。

目前，使用 Matlab 中的周期图函数，我已经能够使用以下方法估计功率谱密度：

现在，为了绘制这个我一直在使用

这显然没有在功率谱密度上放置正确的 x 轴（并产生类似于下图的东西），它需要在频率空间中。鉴于数据的采样不均匀，我对如何解决这个问题感到困惑。

在估计不均匀采样的数据的功率谱密度时，转换到（然后绘制）频率空间的正确方法是什么？我也有兴趣从 python/numpy/scipy 的角度解决这个问题，但到目前为止只查看了 Matlab 函数。

我有一段时间内的脑电图信号数据集，我需要计算它的 PSD 观察结果。

在 matlab 中，我可以计算信号的 PSD，因为我有一个内置函数。

R 是否具有计算数据集的 PSD 观测值的功能？

我在 matlab 中使用 pwelch 方法来计算一些风速测量的功率谱。因此，到目前为止，我已经编写了以下代码作为示例：

我无法包含该情节，因为我没有足够的声望点

这有意义吗？我正在为在情节右侧产生噪音的想法而苦苦挣扎。分解后的信号是一个没有噪声的正弦波，这个噪声是从哪里来的呢？y轴上的值为负的事实是否表明这些频率可以忽略不计？此外，如果风速以 m/s 为单位测量，那么在 y 轴上写单位的最佳方法是什么，这可以转换为对环境科学家更有意义的东西吗？

matplotlib 中的 mlab.csd：http://matplotlib.org/api/mlab_api.html#matplotlib.mlab.csd可用于获得实值交叉谱密度。如果我想从光谱密度中获取相位信息，我需要一个返回复数值的 csd 计算。有吗？

我正在使用Lomb-Scargle代码来估计短时间序列中的一些频率，时间序列显示在第一张图像中。Lomb-Scargle 分析的结果显示在第二个中，我放大了每天大约 2 个周期的显着峰值。然而，该峰值被涂抹，因此证明难以解析该分量的实际频率。是否有任何其他方法或改进我正在使用的方法来准确解析这个短时间序列内的重要频率分量？

这里有一些关于短时间序列方法使用的信息，但不清楚它们是否需要定期采样。理想情况下，我正在寻找一种适用于不规则采样数据的方法，从一些研究看来，最大熵方法是答案，但我不确定这些方法是否已在 MATLAB 中实现？虽然从这个链接，似乎有一个等效的方法，'Yule-Walker AR 方法产生与最大熵估计器相同的结果。但是，是否需要对数据进行均匀采样，这又是一个不清楚的地方？

我一直在使用psd()来计算 .wav 文件上的功率谱密度。我已经把它展示给我的主管，他不希望它平均来计算 Pxx：

每个段的|FFT(i)|^2 被平均以计算Pxx

他建议我使用 PSD，但每次手动将其重叠一帧，而不是传递整个数据数组。我已经尝试过了，它看起来像这样：

imag_array波形文件中的数据数组在哪里。我已经把它发给他了，他对 Python 不是很了解，因为他不能运行它，所以他不能调试它。这看起来正确吗？

我正在尝试将模型拟合到您在下面看到的曲线：

我需要为这个问题找到一个例程，因此它为这些参数提供了优化值（即使曲线如此嘈杂）。我可以找到参数的一些初步猜测，但我需要精确拟合。这是我用来将模型拟合到数据的脚本：

但是，我得到这样的东西：

我不知道出了什么问题，或者还有什么其他方法可以拟合自定义模型。任何人都可以提出一些想法吗？归根结底，我希望为我的数据获得类似的东西：

谢谢。