allpass.c

来自「包含有各种各样的数字信号处理经典算法源代码,很有用的.」· C语言 代码 · 共 20 行

/* allpass.c - allpass reverberator with circular delay line */

double tap();
void cdelay();

double allpass(D, w, p, a, x)                   /* usage: y=allpass(D,w,&p,a,x); */
double *w, **p, a, x;                           /* \(p\) is passed by address */
int D;
{
       double y, s0, sD;

       sD = tap(D, w, *p, D);                   /* \(D\)th tap delay output */
       s0 = x + a * sD;
       y  = -a * s0 + sD;                       /* filter output */
       **p = s0;                                /* delay input */
       cdelay(D, w, p);                         /* update delay line */

       return y;
}