If you are sure you just want the 1st-level (and no further), you can use dwt2.
If you need multiple levels, use wavedec2().
However, your statement " mean and standard deviation of each coefficient and form one dimensional feature vector." does not make sense.
Do mean the mean and standard deviation of each subband image?
You can easily obtain the subband images
load woman;
dwtmode('per');
[CA,CH,CV,CD] = dwt2(X,'haar');
CH, CV, and CD are the wavelet subband images. Now can you point me to some text which describes how to find your feature vector from those matrices?
