You are here: Symbol Reference > Dew Namespace > Dew.Signal Namespace > Dew.Signal.Units Namespace > Classes > SignalUtils Class > SignalUtils Methods > SignalUtils.FractionalFirImpulse Method
Dew Signal for .NET
ContentsIndexHome
PreviousUpNext
SignalUtils.FractionalFirImpulse Method

Design an oversampled FIR filter with rectangular window.

Syntax
C#
Visual Basic
public static void FractionalFirImpulse(int Length, [In] TVec H, [In] double[] W, double Offset, double Step, TFilterType FilterType, double FS);

Compute a FIR impulse response filter (no window applied) and place the result in H. The transition regions are defined with the W array. There must be at least one (lowpass, highpass) and at most two (bandpass, bandstop) transition regions (2 or 4 elements). Filter type is defined with TFilterType. 20*Log10(Ripple) is also the required attenuation of the stop band in decibel. FS is the sampling frequency. Length specifies the length of the original filter and Step defines the oversampling factor. The actual length of the impulse response vector is computed like this: H.Length := Round(Length/Step) Step and offset must be bigger then 0. If the Offest is 0 and step is 1, the routine returns the same result as FirImpulse. The resulting H vector contains FIR type impulse response, which can be passed to an interpolation routine (linear, cubic, lagrange, etc..). Oversampled FIR filters are used for resampling with an arbitrary sampling frequency and, if Offset > 0 and Step = 1, fractional delay FIR filters can be implemented. If Step = 1 then the resulting impulse response can be passed directly to the FirInit and FirFilter routines. When setting Step bigger then 1, the filter designed must also work as an anti-aliasing filter (low-pass) or aliasing will occur.

Filter interpolation. 

 

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 h = new Vector(0); Vector x = new Vector(0); double FS = 2; double Step = 1; //interpolate by 1x Step = 1; SignalUtils.FractionalFirImpulse(30,h,new double[2] {0.5,0.6},0,Step,TFilterType.ftLowpass,FS); x = MtxExpr.Ramp(h,0,1.0/h.Length); MtxVecTee.DrawIt(x,h,"Interpolate by 1x",false); Step = 0.5; //interpolate by 2x SignalUtils.FractionalFirImpulse(30,h,new double[2] {0.5,0.6},0,Step,TFilterType.ftLowpass,FS); x = MtxExpr.Ramp(h,0,1.0/h.Length); MtxVecTee.DrawIt(x,h,"Interpolate by 2x",false); Step = 0.25; //interpolate by 4x SignalUtils.FractionalFirImpulse(30,h,new double[2] {0.5,0.6},0,Step,TFilterType.ftLowpass,FS); x = MtxExpr.Ramp(h,0,1.0/h.Length); MtxVecTee.DrawIt(x,h,"Interpolate by 4x",false); Step = 0.125; //interpolate by 8x SignalUtils.FractionalFirImpulse(30,h,new double[2] {0.5,0.6},0,Step,TFilterType.ftLowpass,FS); x = MtxExpr.Ramp(h,0,1.0/h.Length); MtxVecTee.DrawIt(x,h,"Interpolate by 8x",false); }
Copyright (c) 1999-2024 by Dew Research. All rights reserved.
What do you think about this topic? Send feedback!