Main Content

このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。

大きくサンプリングされているウェーブレット パケット解析

この例では、1 次元信号のウェーブレット パケット変換を求める方法を示します。この例では、周波数順序が Paley 順序とは異なることも示します。

加法性ホワイト ガウス N(0,1/4) ノイズを伴う、周波数が 7π/8 ラジアン/サンプルの正弦波を作成します。正弦波は、信号のサンプル 128 ~ 512 で発生します。dwtmode を周期化に設定し、例の最後でこれを元の設定に戻します。

rng default
st = dwtmode('status','nodisplay');
dwtmode('per','nodisp');

n = 0:1023;
indices = (n>127 & n<=512);
x = cos(7*pi/8*n).*indices+0.5*randn(size(n));

4 つの消失モーメントを持つ Daubechies 最小非対称ウェーブレットを使用し、レベル 2 まで下げてウェーブレット パケット変換を求めます。ウェーブレット パケット ツリーをプロットします。

T = wpdec(x,2,'sym4');
plot(T)

Figure contains 2 axes objects and other objects of type uimenu. Axes object 1 with title Tree Decomposition contains 13 objects of type line, text. Axes object 2 with title data for node: 0 or (0,0). contains an object of type line.

終端ノードの Paley および周波数順序を検出します。

[tn_pal,tn_freq] = otnodes(T);

tn_freq にはベクトル [3 4 6 5] が含まれています。これは、最も高い周波数範囲 [3π/4,π) が実際には Paley 順序のウェーブレット パケット ツリーのノード 5 であることを示しています。

ウェーブレット パケット ツリーのノード (2,2) をクリックして、周波数順序が正弦波の存在を正しく予測していることを確認します。

2 次元イメージのウェーブレット パケット変換により、4 つの要素から成るウェーブレット パケット ツリーが得られます。イメージ例を読み込みます。再構成ウェーブレットに 3 つの消失モーメントと分解ウェーブレットに 5 つの消失モーメントを持つ双直交 B スプライン ウェーブレットを使用します。4 つの要素から成るウェーブレット パケット ツリーをプロットします。

load tartan
T = wpdec2(X,2,'bior3.5');
plot(T)

Figure contains 2 axes objects and other objects of type uimenu. Axes object 1 with title Tree Decomposition contains 41 objects of type line, text. Axes object 2 with title data for node: 0 or (0,0). contains an object of type image.

dwtmode(st,'nodisplay')

参考

アプリ

関数

関連するトピック