denoising an image

I have performed denoising of images using 2d-dwt,now how to plot each subbands such as

LL2,HL2,LH2,HH2,HL1,LH1,HH1,

my code is

I=imread('lena.bmp');

s=imresize(I,[256 256])

N = 256;

% Noise variance
sigma_n = 3;
n = sigma_n*randn(N);

% Add noise 
x = double(s) + double(n);

% Run local adaptive image denoising algorithm using dual-tree DWT. 
y = denoising_dtdwt(x);

please help