annotate src/samer/units/IIRFilter.java @ 8:5e3cbbf173aa tip

Reorganise some more
author samer
date Fri, 05 Apr 2019 22:41:58 +0100
parents bf79fb79ee13
children
rev   line source
samer@0 1 /*
samer@0 2 * Copyright (c) 2000, Samer Abdallah, King's College London.
samer@0 3 * All rights reserved.
samer@0 4 *
samer@0 5 * This software is provided AS iS and WITHOUT ANY WARRANTY;
samer@0 6 * without even the implied warranty of MERCHANTABILITY or
samer@0 7 * FITNESS FOR A PARTICULAR PURPOSE.
samer@0 8 */
samer@0 9
samer@0 10 package samer.units;
samer@0 11 import samer.core.*;
samer@0 12 import samer.maths.*;
samer@0 13 import samer.tools.*;
samer@0 14 import java.util.*;
samer@0 15
samer@0 16 public class IIRFilter implements Filter
samer@0 17 {
samer@0 18 int N; // order of filter
samer@0 19 double A[]; // coefficients
samer@0 20 double u[]; // circular buffer of previous values
samer@0 21 int k; // current position in u buffer
samer@0 22
samer@0 23 public IIRFilter( int n)
samer@0 24 {
samer@0 25 N=n;
samer@0 26 u=new double[n];
samer@0 27 A=new double[n+1];
samer@0 28 reset();
samer@0 29 }
samer@0 30
samer@0 31 public IIRFilter( double [] A)
samer@0 32 {
samer@0 33 N=A.length-1;
samer@0 34 u=new double[N];
samer@0 35 this.A=A;
samer@0 36 reset();
samer@0 37 }
samer@0 38
samer@0 39 public void dispose() {}
samer@0 40 public void reset() { Mathx.zero(u); k=0; }
samer@0 41
samer@0 42 /** FIR filter using array of filter coefficients */
samer@0 43 public final double filter( double z)
samer@0 44 {
samer@0 45 double y=z*A[0];
samer@0 46
samer@0 47 int i=1, j;
samer@0 48 for (j=k; j<N; j++) y += u[j]*A[i++];
samer@0 49 for (j=0; j<k; j++) y += u[j]*A[i++];
samer@0 50
samer@0 51 if (k==0) k=N-1; else k--;
samer@0 52 u[k]=y; // store output value in buffer
samer@0 53 return y;
samer@0 54 }
samer@0 55 }
samer@0 56