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 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 , è 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')