signal reconstruction using wavelets
in [Matlab]
|
From: FAISAL PEER MOAHMED on 9 May 2010 15:46 Hi I am using wavelets to denoise a noisy signal. Perfect reconstruction depends upon the wavelet basis functions used. cross correlation of the signal with wavelts basis can help us to select the optimum wavelet. Comeback to wavelet decomposition , signal is convoluted with filter coefficients. Can any one please help me the relation between cross corelation and convolution in wavelet decomposition ? Regards Faisal
From: Wayne King on 9 May 2010 17:27 "FAISAL PEER MOAHMED" <pfaisalbe(a)gmail.com> wrote in message <hs73dr$6ku$1(a)fred.mathworks.com>... > Hi > > I am using wavelets to denoise a noisy signal. Perfect reconstruction depends upon the > wavelet basis functions used. cross correlation of the signal with wavelts basis can help us to select the optimum wavelet. Comeback to wavelet decomposition , signal is convoluted with filter coefficients. > > Can any one please help me the relation between cross corelation and convolution in wavelet decomposition ? > > Regards > > Faisal Hi Faisal, the cross correlation of a signal with an FIR filter is the same as convolving the signal with the complex-conjugated and time-reversed filter. At the decimation, the wavelet coefficients are obtaining by cross-correlating the signal with the wavelet filter at even lags, but that is equivalent to downsampling the convolution of the signal with the time-reverse and complex-conjugate of the wavelet filter. Note for [LoD,HiD,LoR,HiR] = wfilters('db1'); that the difference between HiD and HiR is just that they are time-reversed versions of each other. Since the filters are real-valued, you obviously don't see the effect of taking the conjugate. Next note that: dwtmode('per','nodisp'); reset(RandStream.getDefaultStream); [ca,cd] = dwt(x,'db1'); produces wavelet coefficients, cd, that are equivalent to: downsample(conv(x,HiD,'valid'),2) Hope that helps, Wayne
From: FAISAL PEER MOAHMED on 10 May 2010 08:52 Many Thanks. As you mentioned "At the decimation, the wavelet coefficients are obtaining by cross-correlating the signal with the wavelet filter at even lags". So in which domain this cross correlation is taking place ? i mean in time domain or frequency domain Regards Faisal "Wayne King" <wmkingty(a)gmail.com> wrote in message <hs79b7$k5i$1(a)fred.mathworks.com>... > "FAISAL PEER MOAHMED" <pfaisalbe(a)gmail.com> wrote in message <hs73dr$6ku$1(a)fred.mathworks.com>... > > Hi > > > > I am using wavelets to denoise a noisy signal. Perfect reconstruction depends upon the > > wavelet basis functions used. cross correlation of the signal with wavelts basis can help us to select the optimum wavelet. Comeback to wavelet decomposition , signal is convoluted with filter coefficients. > > > > Can any one please help me the relation between cross corelation and convolution in wavelet decomposition ? > > > > Regards > > > > Faisal > > Hi Faisal, the cross correlation of a signal with an FIR filter is the same as convolving the signal with the complex-conjugated and time-reversed filter. At the decimation, the wavelet coefficients are obtaining by cross-correlating the signal with the wavelet filter at even lags, but that is equivalent to downsampling the convolution of the signal with the time-reverse and complex-conjugate of the wavelet filter. > > Note for > > [LoD,HiD,LoR,HiR] = wfilters('db1'); > > that the difference between HiD and HiR is just that they are time-reversed versions of each other. Since the filters are real-valued, you obviously don't see the effect of taking the conjugate. Next note that: > > dwtmode('per','nodisp'); > reset(RandStream.getDefaultStream); > [ca,cd] = dwt(x,'db1'); > > produces wavelet coefficients, cd, that are equivalent to: > > downsample(conv(x,HiD,'valid'),2) > > Hope that helps, > Wayne
From: Wayne King on 10 May 2010 09:00 "FAISAL PEER MOAHMED" <pfaisalbe(a)gmail.com> wrote in message <hs8vi7$6jl$1(a)fred.mathworks.com>... > Many Thanks. > > As you mentioned "At the decimation, the wavelet coefficients are obtaining by cross-correlating the signal with the wavelet filter at even lags". > > So in which domain this cross correlation is taking place ? i mean in time domain or frequency domain > > Regards > > Faisal > > > > > > "Wayne King" <wmkingty(a)gmail.com> wrote in message <hs79b7$k5i$1(a)fred.mathworks.com>... > > "FAISAL PEER MOAHMED" <pfaisalbe(a)gmail.com> wrote in message <hs73dr$6ku$1(a)fred.mathworks.com>... > > > Hi > > > > > > I am using wavelets to denoise a noisy signal. Perfect reconstruction depends upon the > > > wavelet basis functions used. cross correlation of the signal with wavelts basis can help us to select the optimum wavelet. Comeback to wavelet decomposition , signal is convoluted with filter coefficients. > > > > > > Can any one please help me the relation between cross corelation and convolution in wavelet decomposition ? > > > > > > Regards > > > > > > Faisal > > > > Hi Faisal, the cross correlation of a signal with an FIR filter is the same as convolving the signal with the complex-conjugated and time-reversed filter. At the decimation, the wavelet coefficients are obtaining by cross-correlating the signal with the wavelet filter at even lags, but that is equivalent to downsampling the convolution of the signal with the time-reverse and complex-conjugate of the wavelet filter. > > > > Note for > > > > [LoD,HiD,LoR,HiR] = wfilters('db1'); > > > > that the difference between HiD and HiR is just that they are time-reversed versions of each other. Since the filters are real-valued, you obviously don't see the effect of taking the conjugate. Next note that: > > > > dwtmode('per','nodisp'); > > reset(RandStream.getDefaultStream); > > [ca,cd] = dwt(x,'db1'); > > > > produces wavelet coefficients, cd, that are equivalent to: > > > > downsample(conv(x,HiD,'valid'),2) > > > > Hope that helps, > > Wayne Time domain, but it is equivalent a sum of multiplications in the frequency domain ( a sum of multiplications and not just one because of the decimation). Wayne
|
Pages: 1 Prev: Subtracting from cell arrays Next: two variables output checker |