Daubechies小波（Daubechies wavelet）.docVIP

Daubechies小波（Daubechies wavelet） Construction of wavelet bases of Daubechies This program implements the construction of wavelet bases % periodic_wavelet.m Function ss=periodic_wavelet; Clear; clc; % global MOMENT;% vanishing order of moments % global LEFT_SCALET;% scale function; left support interval % global RIGHT_SCALET;% scaling function; right support interval % global LEFT_BASIS;% wavelet base function; left support interval % global RIGHT_BASIS;% wavelet function; right support section % global MIN_STEP;% minimum discrete step size % global LEVEL;% calculated the number of layers required (discrete precision) % global MAX_LEVEL;% periodic wavelet maximum number of layers [s2, h]=scale_integer; [test, h]=scalet_stretch (S2, H); Wave_base=wavelet (test, H); Ss=periodic_waveletbasis (wave_base); Function, [s2, h]=scale_integer; % this function to solve the wavelet scaling function, discrete integer points % sacle_integer.m MOMENT=10;% vanishing moment order LEFT_SCALET=0;% scale function; left support interval RIGHT_SCALET=2*MOMENT-1;% scaling function; right support interval LEFT_BASIS=1-MOMENT;% wavelet base function; left support interval RIGHT_BASIS=MOMENT; the right support section of the wavelet base function MIN_STEP=1/512;% minimum discrete step size LEVEL=-log2 (MIN_STEP);% calculated the number of layers required (discrete precision) MAX_LEVEL=8;% periodic wavelet; maximum number of layers H=wfilters (db10,r) -% filter factor H=h*sqrt (2);% FI (T) =SQRT (2) *SUM (H (N) *FI (2T-N)) N=0:2*MOMENT-1; For i=LEFT_SCALET+1:RIGHT_SCALET-1 For j=LEFT_SCALET+1:RIGHT_SCALET-1 K=2*i-j+1; If (k=1k=RIGHT_SCALET+1) A (I, J) =h (k);% matrix coefficient matrix Else A (I, J) =0; End End End [s, w]=eig (a);% for characteristic vectors, bases of Solutions S1=s ((:: 1); S2=[0; s1/sum (S1); 0];%, according to the conditional SUM (FI (T)) =1; This function extends the discrete values of the scaling function % scalet_stretch.m Function, [s2, h]=scalet_stretch (S2, H); MOMENT