La traduzione di questa pagina non è aggiornata. Fai clic qui per vedere l'ultima versione in inglese.

Analisi di pacchetti wavelet a campionamento critico

Questo esempio mostra come ottenere la trasformata di pacchetti wavelet di un segnale monodimensionale. L'esempio dimostra anche che l'ordinamento in frequenza è diverso dall'ordinamento di Paley.

Creare un segnale costituito da un’onda sinusoidale con una frequenza di $$7\pi /8$$ radianti/campione in rumore gaussiano bianco additivo N(0,1/4). L’onda sinusoidale si verifica tra i campioni 128 e 512 del segnale. Impostare dwtmode su periodizzazione e riportarlo sull'impostazione originale alla fine dell'esempio.

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));

Ottenere la trasformata dei pacchetti wavelet fino al livello 2 utilizzando la wavelet asimmetrica minima di Daubechies con 4 momenti di fuga. Tracciare l'albero dei pacchetti wavelet.

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

Trovare l’ordinamento di Paley e di frequenza dei nodi terminali.

[tn_pal,tn_freq] = otnodes(T);

tn_freq contiene il vettore [3 4 6 5], che mostra che l’intervallo di frequenza più alto $$[3\pi /4,\pi )$$, è in effetti il nodo 5 nell’albero dei pacchetti wavelet nell’ordinamento di Paley.

Facendo clic sul nodo (2,2) dell'albero dei pacchetti wavelet, è possibile notare che l'ordinamento delle frequenze prevede correttamente la presenza dell'onda sinusoidale.

La trasformata di pacchetti wavelet di un'immagine bidimensionale produce un albero dei pacchetti wavelet quaternario. Caricare un'immagine di esempio. Utilizzare la wavelet B-spline biortogonale con 3 momenti di fuga nella wavelet di ricostruzione e 5 momenti di fuga nella wavelet di scomposizione. Tracciare l'albero dei pacchetti wavelet quaternari risultante.

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

dwtmode(st,'nodisplay')

Vedi anche

App

• Wavelet Signal Analyzer

Funzioni

• dwpt | wavedec | wavedec2

Argomenti complementari

• Practical Introduction to Multiresolution Analysis
• Wavelet Packets: Decomposing the Details