问题标签 [spectral]

我正在尝试使用频谱分析来分析时间序列。我试图检测我的数据中的任何周期性，这些周期性由记录一周的每小时测量值组成（24 * 7 = 168 次测量值），我的目标是显示温度变化的昼夜分量。到目前为止，我有（例如）：

我会从这里去哪里？任何建议将不胜感激。

更改：

谁能建议我如何将 x 轴转换为小时而不是每小时的周期？我试过了

但这弄乱了周期图。

我有四个 1xN 声音信号，我想查看以下链接中给出的平均频谱图：http: //i1233.photobucket.com/albums/ff396/sakurayen/Plot/AMaximumLikelihoodApproachtoSinglechannelSourceseparationpdf-AdobeReader.jpg

我尝试使用 MATLAB 函数 PSD 来绘制光谱，但我得到了不同的图。请注意，用于两个图的数据是相同的。

在 MATLAB 中使用 PSD 函数获得的绘图：http: //i1233.photobucket.com/albums/ff396/sakurayen/Plot/PowerSpectralDensityofRJMF.png

谢谢！

问题

我编写了一个小型 python 批处理器，它加载二进制数据、执行 numpy 操作并存储结果。它消耗的内存比它应该消耗的多得多。我查看了类似的堆栈溢出讨论，并想寻求进一步的建议。

背景

我将光谱数据转换为 rgb。光谱数据存储在带线交错 (BIL) 图像文件中。这就是我逐行读取和处理数据的原因。我使用Spectral Python Library读取数据，它返回一个 numpy 数组。hyp是大型光谱文件的描述符：hyp.ncols=1600, hyp.nrows=3430, hyp.nbands=160

代码

内存消耗

所有大数据都在循环外初始化。我天真的解释是内存消耗不应该增加（我使用了太多的 Matlab 吗？）有人可以解释一下增加因子 10 吗？这不是线性的，因为 hyp.nrows = 3430。有什么建议可以改善内存管理吗？

笔记

我用 xrange 替换了range ，没有大幅改进。我知道三次插值并不是最快的，但这与 CPU 消耗无关。

我想创建或生成 [-1 1] 范围内的白噪声，但我不知道具体该怎么做。我选择的编程语言是matlab。据我所知，存在一个名为的函数randn和一个名为wgn（高斯白噪声）的函数。所以请帮我解决这个问题，澄清一下，例如我想生成以下等式：

其中A3=amplitude ，白噪声在 [-1 1] 范围内。请帮助我并澄清如何做。我的困惑与白噪声有关，与其他人无关，让我们假设 t 从 1 变为 100。

我试图通过以下模型拟合以下正弦模型，这意味着确定性模型+白噪声：

但是当我使用以下命令运行时：

我收到以下错误：

那么有什么问题呢？据我所知，我可以为rand函数使用两个参数，对吗？这里有什么问题？

更新：这就是我为修复它所做的

和结果

我有兴趣在 Mahout 中对高维数据运行谱聚类算法。我的问题是如何获取高维数据向量列表并创建最近邻图？这是在 Mahout 中完成的，还是有地图可简化的方法来做这样的事情。

我pwelch用来绘制功率谱密度。我想使用格式

但在 x 轴上有一个对数刻度。

我也试过

但它没有给出与上面相同的结果图。所以，

不是一个好的解决方案。

现在我正在用声音进行频谱分析，奇怪的是用 绘制图形与 Periodogram[]使用PeridogramArray[]然后使用的结果不同ListLinePlot[]。有什么我理解错了吗？

源代码：

我在 Matlab 中阅读 Trefethen 的光谱方法。

在创建微分矩阵时，

列= [任何东西]

D=toeplitz(列,列([1 N:-1:2]))

有人可以解释一下上面行中的 [ ... ] 内到底发生了什么。

我了解您正在移动列，但我不了解该语法。

我将我的 MacBook 更新为 Mavericks，重新安装了 Macports 和我通常使用的所有 Python 2.7 模块。在运行 Python 时，我收到以下消息：

导入 mlab 时：

from mayavi import lab

(process:1146): Gtk-WARNING **: Locale not supported by C library. Using the fallback 'C' locale.

运行 mlab 命令（例如 mlab.mesh()）时，显示窗口打开，不显示任何内容并冻结。

导入光谱时我没有收到此消息，但在运行view_cube()显示图像立方体的显示窗口时收到此消息，冻结但显示数据立方体。Xterm似乎有问题，但我无法弄清楚。如何防止显示窗口冻结并摆脱 Gtk-WARNING？

我检查了 locale 和locale -a，但没有发现任何异常： locale:

locale
LANG=
LC_COLLATE="C"
LC_CTYPE="C"
LC_MESSAGES="C"
LC_MONETARY="C"
LC_NUMERIC="C"
LC_TIME="C"
LC_ALL=