Mercurial > hg > jslab
annotate src/samer/units/FIRFilter.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 FIRFilter 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 FIRFilter( int n) |
samer@0 | 24 { |
samer@0 | 25 N=n; |
samer@0 | 26 u=new double[n]; |
samer@0 | 27 A=new double[n]; |
samer@0 | 28 reset(); |
samer@0 | 29 } |
samer@0 | 30 |
samer@0 | 31 public FIRFilter( double [] A) |
samer@0 | 32 { |
samer@0 | 33 N=A.length; |
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 if (k==0) k=N-1; else k--; |
samer@0 | 46 double y=0; |
samer@0 | 47 u[k]=z; |
samer@0 | 48 |
samer@0 | 49 int i=0, j; |
samer@0 | 50 for (j=k; j<N; j++) y += u[j]*A[i++]; |
samer@0 | 51 for (j=0; j<k; j++) y += u[j]*A[i++]; |
samer@0 | 52 |
samer@0 | 53 return y; |
samer@0 | 54 } |
samer@0 | 55 } |
samer@0 | 56 |