問題タブ [spectral]

スペクトル分析で時系列を分析しようとしています。1 週間 (24 * 7 = 168 測定) 記録された 1 時間ごとの測定値で構成されるデータの周期性を検出しようとしています。温度変化の日周成分を示すことを目指しています。これまでのところ、（たとえば）：

ここからどこへ行く？アドバイスをいただければ幸いです。

変更されました：

1時間あたりのサイクルではなく、x軸を時間に変換する方法を誰かが提案できますか? 私が試してみました

しかし、これはペリドグラムを台無しにします。

4 つの1xNサウンド信号があり、以下のリンクに示されているような平均スペクトル プロットを表示したいと考えています。

MATLAB 関数 PSD を使用してスペクトルをプロットしようとしましたが、代わりに別のプロットが得られます。両方のプロットに使用されるデータは同じであることに注意してください。

MATLAB で PSD 関数を使用して取得したプロット: http://i1233.photobucket.com/albums/ff396/sakurayen/Plot/PowerSpectralDensityofRJMF.png

ありがとう！

質問

バイナリ データをロードし、numpy 操作を実行し、結果を保存する小さな python バッチ プロセッサを作成しました。必要以上に多くのメモリを消費します。同様のスタックオーバーフローの議論を見て、さらに推奨事項を求めたいと思います。

バックグラウンド

スペクトル データを RGB に変換します。スペクトル データは、Band Interleaved by Line (BIL) イメージ ファイルに保存されます。そのため、データを 1 行ずつ読み取って処理します。numpy 配列を返すSpectral Python Libraryを使用してデータを読み取ります。hypは大きなスペクトル ファイルの記述子です: hyp.ncols=1600、hyp.nrows=3430、hyp.nbands=160

コード

メモリ消費量

すべてのビッグデータはループの外で初期化されます。私の素朴な解釈は、メモリ消費量が増えてはならないというものでした (Matlab を使いすぎたのでしょうか?) 増加係数 10 について説明してもらえますか? hyp.nrows = 3430 であるため、これは線形ではありません。メモリ管理を改善するための推奨事項はありますか?

ノート

rangeをxrangeに置き換えましたが、劇的な改善はありませんでした。三次補間が最速ではないことは承知していますが、これは CPU 消費に関するものではありません。

[-1 1] の範囲でホワイト ノイズを作成または生成したいのですが、正確な方法がわかりません。私が選んだプログラミング言語は matlab です。私の知る限り、 (ホワイト ガウス ノイズ)randnという名前の関数と、wgn(ホワイト ガウス ノイズ) という名前の関数が存在します。明確にするために、この問題について私を助けてください。たとえば、次の方程式を生成したいです。

ここでA3=amplitude 、ホワイト ノイズの範囲は [-1 1] です。私を助けて、それを行う方法を明確にしてください。私の混乱はホワイト ノイズに関連しており、他のノイズには関連していません。t が 1 から 100 に変化していると仮定しましょう。

次のモデルで次の正弦波モデルに適合させようとしました。これは、決定論的モデル + ホワイト ノイズを意味します。

しかし、次のコマンドを使用して実行すると:

次のエラーが発生しました：

では、何が問題なのですか？私が知っているように、rand関数に2つのパラメーターを使用できますか? ここで何が問題ですか？

更新：それは私がそれを修正するために行ったことです

と結果

高次元データに対して Mahout でスペクトル クラスタリング アルゴリズムを実行することに興味があります。私の質問は、高次元のデータ ベクトルのリストを取得し、最近傍グラフを作成する方法です。これはMahoutで行われますか、それともそのようなことを行うマップ削減可能な方法はありますか.

pwelchパワースペクトル密度をプロットするために使用しています。フォーマットを使いたい

ただし、x 軸に対数スケールがあります。

私も試してみました

しかし、上記と同じ結果のプロットは得られません。したがって、

良い解決策ではありません。