ウェーブレット変換を用いてsupectrumをプロットしたい

実現したいこと

ウェーブレット変換を使ってspectrumをプロットしたい。

前提

pythonで観測データ(dat形式)を読み込み、ウェーブレット変換を行い、supectrumをプロットしようとしています。
datファイルには、縦0列目に時間データ、縦1〜4列目までが観測データ(系4ch)あります。(要素数は110001)
各チャンネルごとにウェーブレット変換を行いたいです。
プログラミング初心者なため、調べながらやっていますが、どのサイトも解析対象の関数を作っているので観測データへの適用の仕方がわからず、どの部分を使えば良いかわからないです。
どなたかお教えしていただきたいです。

発生している問題・エラーメッセージ

該当のソースコード

試したこと

ファイルを読み込む
データを時間、観測データごとに読み込む
time = file[:,0]
data1 = file[:,1]
data2 = file[:,2] etc.
パラメータ
N = 110001 #要素数
dt = 0.02 #sampling-time
fs = 1/dt #sanpling-freqency
nq_f = fs/2 #ナイキスト周波数
wavelet_type = 'cmor1.5-1.0' # ウェーブレットのタイプ
追加すべきパラメータがあればご教示ください。

補足情報(FW/ツールのバージョンなど)

ここにより詳細な情報を記載してください。

コメントを投稿

0 コメント