01: /*
02: * Copyright (c) 2000 Silvere Martin-Michiellot All Rights Reserved.
03: *
04: * Silvere Martin-Michiellot grants you ("Licensee") a non-exclusive,
05: * royalty free, license to use, modify and redistribute this
06: * software in source and binary code form,
07: * provided that i) this copyright notice and license appear on all copies of
08: * the software; and ii) Licensee does not utilize the software in a manner
09: * which is disparaging to Silvere Martin-Michiellot.
10: *
11: * This software is provided "AS IS," without a warranty of any kind. ALL
12: * EXPRESS OR IMPLIED CONDITIONS, REPRESENTATIONS AND WARRANTIES, INCLUDING ANY
13: * IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE OR
14: * NON-INFRINGEMENT, ARE HEREBY EXCLUDED. Silvere Martin-Michiellot
15: * AND ITS LICENSORS SHALL NOT BE LIABLE FOR ANY DAMAGES
16: * SUFFERED BY LICENSEE AS A RESULT OF USING, MODIFYING
17: * OR DISTRIBUTING THE SOFTWARE OR ITS DERIVATIVES. IN NO EVENT WILL
18: * Silvere Martin-Michiellot OR ITS LICENSORS BE LIABLE
19: * FOR ANY LOST REVENUE, PROFIT OR DATA, OR FOR DIRECT,
20: * INDIRECT, SPECIAL, CONSEQUENTIAL, INCIDENTAL OR PUNITIVE DAMAGES, HOWEVER
21: * CAUSED AND REGARDLESS OF THE THEORY OF LIABILITY, ARISING OUT OF THE USE OF
22: * OR INABILITY TO USE SOFTWARE, EVEN IF Silvere Martin-Michiellot HAS BEEN
23: * ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
24: *
25: * This software is not designed or intended for use in on-line control of
26: * aircraft, air traffic, aircraft navigation or aircraft communications; or in
27: * the design, construction, operation or maintenance of any nuclear
28: * facility. Licensee represents and warrants that it will not use or
29: * redistribute the Software for such purposes.
30: *
31: */
32:
33: // This code is repackaged after the code from Craig A. Lindley, from Digital Audio with Java
34: // Site ftp://ftp.prenhall.com/pub/ptr/professional_computer_science.w-022/digital_audio/
35: // Email
36: package com.db.media.audio.dsp.filter;
37:
38: // The IIR highpass filter designed here has unity gain
39:
40: public class IIRHighpassFilterDesign extends IIRFilterDesignBase {
41:
42: public IIRHighpassFilterDesign(int cutoffFrequency, int sampleRate,
43: double dampingFactor) {
44:
45: super (cutoffFrequency, sampleRate, dampingFactor);
46:
47: }
48:
49: // Do the design of the filter. Filter response controlled by
50: // just three coefficients: alpha, beta and gamma.
51: public void doFilterDesign() {
52:
53: // Get radians per sample at cutoff frequency
54: double thetaZero = getThetaZero();
55:
56: double theSin = parameter / (2.0 * Math.sin(thetaZero));
57:
58: // Beta relates gain to cutoff freq
59: beta = 0.5 * ((1.0 - theSin) / (1.0 + theSin));
60:
61: // Final filter coefficient
62: gamma = (0.5 + beta) * Math.cos(thetaZero);
63:
64: // For unity gain
65: alpha = (0.5 + beta + gamma) / 4.0;
66:
67: }
68:
69: // Entry point for testing
70: public static void main(String[] args) {
71:
72: if (args.length != 3) {
73: System.out.println("\nIIR Highpass Filter Design Program");
74: System.out.println("\nUsage:");
75: System.out.println("\tIIRHighpassFilterDesign "
76: + "cutoffFrequency sampleRate dampingFactor");
77:
78: System.exit(1);
79: }
80: // Parse the command line arguments
81: int cutoffFrequency = new Integer(args[0]).intValue();
82: int sampleRate = new Integer(args[1]).intValue();
83: double dampingFactor = new Double(args[2]).doubleValue();
84:
85: // Instantiate highpass filter designer
86: IIRHighpassFilterDesign d = new IIRHighpassFilterDesign(
87: cutoffFrequency, sampleRate, dampingFactor);
88: // Do the design
89: d.doFilterDesign();
90:
91: // Print out the coefficients
92: d.printCoefficients();
93:
94: }
95:
96: }
|