1次元連続ウェーブレット変換で、ウェーブレット係数を表示するとき、縦軸を擬似周波数として表示させることはできますか?
15 views (last 30 days)
Show older comments
MathWorks Support Team
on 25 Oct 2013
Answered: MathWorks Support Team
on 25 Oct 2013
1次元連続ウェーブレット変換で、ウェーブレット係数を表示するとき、縦軸を擬似周波数として表示させる方法を教えてください。
Accepted Answer
MathWorks Support Team
on 25 Oct 2013
スケールファクタから、擬似周波数への変換は、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) % ウェーブレット係数表示
0 Comments
More Answers (0)
See Also
Categories
Find more on 連続ウェーブレット変換 in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!