I already solve this problem. This is my solution, in case it helps someone:
[c,l]=wavedec(bandpass_data,3,'db2'); %Descomposición de la señal en 3 niveles
Lev=1;
ca1=appcoef(c,l,'db2',Lev); %Coeficientes de aproximación (Frecuencias bajas)
Lev2=2;
ca2=appcoef(c,l,'db2',Lev2); %Coeficientes de aproximación (Frecuencias bajas)
Lev3=3;
ca3=appcoef(c,l,'db2',Lev3); %Coeficientes de aproximación (Frecuencias bajas)
[cd1,cd2,cd3]=detcoef(c,l,[1,2,3]); %Coeficientes de detalle (Frecuencias altas)
subplot(7,1,1)
plot(bandpass_data)
title('Señal Original')
subplot(7,1,2)
plot(ca1)
title('Coeficientes de Aproximación - Nivel 1')
subplot(7,1,3)
plot(ca2)
title('Coeficientes de Aproximación - Nivel 2')
subplot(7,1,4)
plot(ca3)
title('Coeficientes de Aproximación - Nivel 3')
subplot(7,1,5)
plot(cd1)
title('Coeficientes de Detalle - Nivel 1')
subplot(7,1,6)
plot(cd2)
title('Coeficientes de Detalle - Nivel 2')
subplot(7,1,7)
plot(cd3)
title('Coeficientes de Detalle - Nivel 3')