Calcular la transformada wavelet unidimensional en MATLAB - tomando la señal Doppler como ejemplo

La función para ejecutar IDWT es: x=waverec(C,L,'wname')

Donde wname es el nombre de la wavelet; x es la señal de secuencia en el dominio del tiempo; N es la serie de la transformada de wavelet.

Utilice coeficientes wavelet parciales para reconstruir la señal: x=wrcoef('type',C,L,'wname',N)

donde type='a' es un N -aproximación de nivel La señal de reconstrucción de componentes

type='d' es la señal de reconstrucción de componentes de detalle de nivel N

Doppler= wnoise(4,10);

N = 5;

nombrew = 'db4';

[C,L] = wavedec(doppler,N,wname); waverec(C ,L,wname);

A5 = wrcoef('a',C,L, wname,5);

A4 = wrcoef('a',C, L,wname ,4);

A3 = wrcoef('a',C,L,wname,3);

A2 = wrcoef('a',C,L, wname,2 );

A1 = wrcoef('a ',C,L,wname,1);

D5 = wrcoef('d',C,L,wname, 5);

D4 = wrcoef('d',C,L,wname,4);

D3 = wrcoef('d',C,L,wname,3) ;

p>

D2 = wrcoef('d',C,L,wname,2);

D1 = wrcoef('d',C,L,wname ,1);

figure(1);

plot(doppler,'k','linewidth',2);

espera;

plot(A0 ,'-.r','linewidth',2);

espera;

plot(A5,'--b',' linewidth',2);

legend('señal Doppler original','señal IDWT----A0','señal de aproximación de nivel 5----A5');

eje([0 400 -1 1]);

set(gcf,'unidad','centímetros','posición',[10 10 30 10]) )

figura(2);

subplot(3,2,1);

plot(doppler,'k','linewidth',2);

título('señal Doppler');

eje([0 400 -1 1]

subplot(3,2,2);