wavelets

DWT フィルターバンクの時間領域ウェーブレット

ページ内をすべて折りたたむ

構文

psi = wavelets(fb)

[psi,t] = wavelets(fb)

説明

psi = wavelets(fb) は、離散ウェーブレット変換 (DWT) フィルターバンク fb のウェーブレット通過帯域フィルターに対応する時間領域と中央揃えされたウェーブレットを返します。

例

[psi,t] = wavelets(fb) は、サンプリング時点 t を返します。

例

すべて折りたたむ

DWT フィルターバンクウェーブレット

ライブスクリプトを開く

Daubechies db2 ウェーブレットと 1 kHz のサンプリング周波数を使用して、信号長が 1,000 サンプルの 7 レベル DWT フィルターバンクを作成します。

wv = "db4";
len = 1000;
lev = 7;
Fs = 1e3;
fb = dwtfilterbank('Wavelet',wv,'SignalLength',len,'Level',lev,'SamplingFrequency',Fs);

ウェーブレットバンドパスフィルターに対応する時間領域と中央揃えされたウェーブレットをプロットします。

[psi,t] = wavelets(fb);
plot(t,psi')
grid on
title('Time-domain Wavelets')

Figure contains an axes object. The axes object with title Time-domain Wavelets contains 7 objects of type line.

最も細かいスケールの時間領域ウェーブレットと、対応するウェーブレットバンドパスフィルターの片側振幅周波数応答をプロットします。

sc = 1;
[psidft,f] = freqz(fb);
subplot(2,1,1)
plot(t,psi(sc,:))
grid on
xlabel('Time (sec)')
ylabel('Magnitude')
title(['Level ',num2str(sc),' Time-Domain Wavelet'])
subplot(2,1,2)
plot(f(len/2:end),abs(psidft(sc,len/2:end)))
grid on
xlabel('Hz')
ylabel('Magnitude')
title('Magnitude Frequency Response')

Figure contains 2 axes objects. Axes object 1 with title Level 1 Time-Domain Wavelet contains an object of type line. Axes object 2 with title Magnitude Frequency Response contains an object of type line.

入力引数

すべて折りたたむ

`fb` — 離散ウェーブレット変換フィルターバンク
`dwtfilterbank` オブジェクト

離散ウェーブレット変換 (DWT) フィルターバンク。dwtfilterbank オブジェクトとして指定します。

出力引数

すべて折りたたむ

`psi` — 時間中心のウェーブレット
実数値の行列

ウェーブレット通過帯域フィルターに対応する時間中心のウェーブレット。L 行 N 列の行列として返されます。ここで、L はフィルターバンクの Level、N は SignalLength です。psi のウェーブレットは、最も細かいスケール分解能から最も粗いスケール分解能の順に並べられます。

`t` — サンプリング時点
実数値のベクトル

サンプリング時点。長さ N の実数値ベクトル t として返されます。ここで、N はフィルターバンクの SignalLength です。サンプリング時点は、間隔 $[- ½ N D T, ½ N D T)$ にあります。ここで、DT はフィルターバンクのサンプリング周期 (フィルターバンクのサンプリング周波数の逆数) です。