Estimate the order a Chebyshev type I IIR filter.
public int ChebyshevIOrder(double[] BEdges, double PassRipple, double StopRipple, TFilterType FilterType, ref double[] CutoffFreq, bool Analog);
Returns the order of the Chebyshev type I filter. Bedg array must contain the band edges of the transition region(s) sorted in ascending order. PassRipple defines the ripple of the passband and StopRipple defines the ripple of the stopband. The length of the CutoffFreq array must be equal to one half of the length of the BEdg array and must match the specified FilterType. The routine returns the estimated order as a result and fill's the CutoffFreq array. This array can then be passed to the ChebyshevIFilter routine.
Design an analog highpass filter with transition band between 1 and 4 rad/sec and with at least 50dB attenuation at the end of the transition band and. The passband should not have more then 0.1dB ripple.
using Dew.Math; using Dew.Math.Editors; using Dew.Math.Units; using Dew.Signal; using Dew.Signal.Units; using Dew.Math.Tee; using Dew.Signal.Tee; private void button1_Click(object sender, EventArgs e) { Vector z = new Vector(0); Vector p = new Vector(0); Vector num = new Vector(0); Vector den = new Vector(0); Vector Response = new Vector(0); Vector FreqFr = new Vector(0); double k; double[] Wc = new double[1]; int Order; //design a fifth order filter. Order = IIRFilters.ChebyshevIOrder(new double[2] { 1, 4 }, 0.2,40,TFilterType.ftHighPass,ref Wc,true); //design analog protype IIRFilters.ChebyshevIAnalog(Order, 0.1, z, p, out k); //design analog protype LinearSystems.LowpassToHighpass(z, p, ref k, Wc[0]); LinearSystems.ZeroPoleToTransferFun(num, den, z, p, k); FreqFr.Length = 1000; SignalUtils.LogRamp(FreqFr, -1, 1); SignalUtils.FrequencyResponses(num, den, FreqFr, Response, 0); TeeChart.DrawIt(Response, "Frequency response", false); }
|
What do you think about this topic? Send feedback!
|
|
Copyright (c) 1999-2010 by Dew Research. All rights reserved.
|