1次元連続ウェーブレ​ット変換で、ウェーブ​レット係数を表示する​とき、縦軸を擬似周波​数として表示させるこ​とはできますか？

1次元連続ウェーブレット変換で、ウェーブレット係数を表示するとき、縦軸を擬似周波数として表示させる方法を教えてください。

Answer 1

[MathWorks Support Team]

スケールファクタから、擬似周波数への変換は、SCAL2FRQ 関数で求めることができます。SCAL2FRQ 関数は、下記の様な書式で使用します。

 f = scal2frq(scalef,wname,delta)

 f : 擬似周波数ベクトル (Hz) 
 scalef : スケールファクタ 
 wname : マザーウェーブレット 
 delta : サンプリング周波数 (sec)

なお、擬似周波数(Fa)は、マザーウェーブレットの中心周波数(Fc)をもとにに次式で計算されています。

 Fa = Fc ./ (scalef * delta)

適用するマザーウェーブレット(wname) の 中心周波数(Fc)は、CENTFRQ 関数より求めることができます。

Fc = centfrq(wname)

縦軸を擬似周波数として表示するサンプルプログラム wavefreq.mを、関連ドキュメントよりダウンロードすることができます。

(サンプルプログラム使用例)

ダウンロード後、MATLAB Command Window上で下記コマンドを実行してください。

Fs = 100; % サンプリング周波数  
time = 0:1/Fs:2; % 時間ベクトル 
sig = sin(2*pi*6*time)+sin(2*pi*15*time);  % 信号 (6Hz / 15 Hz の正弦波) 
wname = 'gaus6'; % マザーウェーブレット 
scalef = 1:20; % スケールファクタ 
scalinterv = 4:4:20; % スケールファクタのラベル表示間隔 

wavefrq(Fs,sig,wname,scalef,scalinterv) % ウェーブレット係数表示