No self-introduction
Loading
Profile
Follow
codes (1)
Wavelet transform and inverse transform
no vote
%Examples% generate one-dimensional signal & nbsp; CLC; clearwname '; db3'; randn ('seed ', 531316785); s = 2 + Kron (ones (1,8), [1 - 1]) + ((1:16). ^ 2) / 32 + 0.2 * randn (1,16); figure; plot (s); title (' (a) original image ');% use db3 wavelet for two-scale decomposition% wavelet decomposition function DWT [CA1, CD1] = DWT (s, wname); figure; plot (CA1); title (' (b) approximation coefficient CA1 '); figure; plot (CD1);% use db3 wavelet for two-scale decomposition; Title ('(c) detail coefficient CD1');% decompose the approximate coefficients on a scale [Ca2, CD2] = DWT (CA1, wname); figure; plot (Ca2); title ('(d) approximation coefficient Ca2'); figure; plot (CD2); title ('(e) detail coefficient CD2'); [lo_ R,Hi_ R]=wfilters(wname,'r');ss=idwt(cA2,cD2,Lo_ R,Hi_ R);sss=idwt(ss,cD1,Lo_ R,Hi_ R) Error between reconstructed signal and original signal err = norm (s-sss) figure; plot (SSS); title ('(f) reconstructed original signal'); xlabel (['relative error = ', num2str (ERR)]);
why_dvdv 0 0
1
No more~
Follows0
Fans0