一维离散小波算法

最新推荐文章于 2022-07-19 14:39:57 发布

原创最新推荐文章于 2022-07-19 14:39:57 发布 · 677 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #buffer #delete

STL/Algorithm 专栏收录该内容

32 篇文章

订阅专栏

本文介绍了一种实现一维离散小波变换(DWT)的方法，详细展示了如何通过循环遍历输入缓冲区并应用小波滤波器来完成变换过程。文中涉及了不同边界条件下的延拓方式，包括周期延拓、对称延拓和0填充。

void CDwt::Dwt1D(double *buffer, int buflen, int hBen, int hEnd, int h1Ben, int h1End, double *h, double *g) { double *p; int m,n,tempn; p = new double [buflen]; for (m=0; m<buflen; m++){="" *(p+m)="*(buffer+m);" *(buffer+m)="0;" }="" for="" (m="0;" m=""><(buflen>>1); m++){ for (n=hBen+2*m; n<=hEnd+2*m; n++){ tempn=n; /* 恒等延拓 */ /*if ( tempn<0 ) tempn=0; if ( tempn>= buflen ) tempn=buflen-1;*/ /* 周期延拓/ * if ( tempn<0 ) tempn=tempn+buflen; if ( tempn >buflen-1 ) tempn=tempn-buflen; /* 对称延拓 */ /**/if ( tempn<0 ) tempn=-tempn; if ( tempn >buflen-1 ) tempn=2*buflen-tempn-2; /* 0 填充 */ /*if ( tempn<0 ) tempn=0; if ( tempn >buflen-1 ) tempn=0;*/ *(buffer+m) = *(buffer+m)+*(p+tempn)*(*(h+n-2*m-hBen)); }// for n for (n=2*m+1-h1End; n<=2*m+1-h1Ben; n++){ tempn=n; /* 恒等延拓 */ /*if ( tempn<0 ) tempn=0; if ( tempn>= buflen ) tempn=buflen-1;*/ /* 周期延拓 /* if ( tempn<0 ) tempn=tempn+buflen; if ( tempn >buflen-1 ) tempn=tempn-buflen; /* 对称延拓 */ /**/if ( tempn<0 ) tempn=-tempn; if ( tempn >buflen-1 ) tempn=2*buflen-tempn-2; /* 0 填充 */ /*if ( tempn<0 ) tempn=0; if ( tempn >buflen-1 ) tempn=0;*/ *(buffer+m+(buflen>>1)) = *(buffer+m+(buflen>>1))+*(p+tempn)*(*(g+n-2*m-1+h1End)); }// for n }// for m delete []p; }