Filter Design Results

Generated by:   http://www-users.cs.york.ac.uk/~fisher/mkfilter

Summary

You specified the following parameters:

Results

Command line: /www/usr/fisher/helpers/mkfilter -Ch -7.0000000000e-02 -Lp -o 4 -a 4.8828125000e-01 0.0000000000e+00
raw alpha1    =   0.4882812500
raw alpha2    =   0.4882812500
warped alpha1 =   8.6421677345
warped alpha2 =   8.6421677345
gain at dc    :   mag = 1.091256182e+00   phase =   0.0000000000 pi
gain at centre:   mag = 1.091256182e+00   phase =  -0.8699514648 pi
gain at hf    :   mag = 0.000000000e+00

S-plane zeros:

S-plane poles:
	-15.4841346398 + j  62.5630864614
	-37.3820078491 + j  25.9144789162
	-37.3820078491 + j -25.9144789162
	-15.4841346398 + j -62.5630864614

Z-plane zeros:
	 -1.0000000000 + j   0.0000000000	4 times

Z-plane poles:
	 -0.9834267115 + j   0.0593038260
	 -0.9291213337 + j   0.0466401741
	 -0.9291213337 + j  -0.0466401741
	 -0.9834267115 + j  -0.0593038260

Recurrence relation:
y[n] = (  1 * x[n- 4])
     + (  4 * x[n- 3])
     + (  6 * x[n- 2])
     + (  4 * x[n- 1])
     + (  1 * x[n- 0])

     + ( -0.8400367511 * y[n- 4])
     + ( -3.5058911151 * y[n- 3])
     + ( -5.4909777508 * y[n- 2])
     + ( -3.8250960906 * y[n- 1])

Ansi ``C'' Code

/* Digital filter designed by mkfilter/mkshape/gencode   A.J. Fisher
   Command line: /www/usr/fisher/helpers/mkfilter -Ch -7.0000000000e-02 -Lp -o 4 -a 4.8828125000e-01 0.0000000000e+00 -l */

#define NZEROS 4
#define NPOLES 4
#define GAIN   1.091256182e+00

static float xv[NZEROS+1], yv[NPOLES+1];

static void filterloop()
  { for (;;)
      { xv[0] = xv[1]; xv[1] = xv[2]; xv[2] = xv[3]; xv[3] = xv[4]; 
        xv[4] = next input value / GAIN;
        yv[0] = yv[1]; yv[1] = yv[2]; yv[2] = yv[3]; yv[3] = yv[4]; 
        yv[4] =   (xv[0] + xv[4]) + 4 * (xv[1] + xv[3]) + 6 * xv[2]
                     + ( -0.8400367511 * yv[0]) + ( -3.5058911151 * yv[1])
                     + ( -5.4909777508 * yv[2]) + ( -3.8250960906 * yv[3]);
        next output value = yv[4];
      }
  }

Download code and/or coefficients:

Magnitude (red) and phase (blue) vs. frequency

For an expanded view, enter frequency limits (as a fraction of the sampling rate) here:

Lower limit: Upper limit:

Impulse response

Step response

For a view on a different scale, enter upper time limit (integer number of samples) here:

Upper limit:

Tony Fisher fisher@minster.york.ac.uk