Code covered by the BSD License  

Highlights from
DFiltFIR

  • DFiltFIR
  • Example1
  • Example2
  • Example3
  • IntMCPiecewise monotonic cubic interpolation.
  • IntMCIPiecewise monotonic cubic interpolation.
  • MMBSpecCheck and fill in missing band specifications
  • MMcosCofGiven the sample values for a polynomial, this routine returns the
  • MMdSpecConvert band-by-band specifications into a specification on a dense
  • MMdesLPFIRThis routine designs multiple band bandpass filters, differentiators,
  • MMdesRemezThis routine finds a minimum maximum deviation approximation to a
  • MMexchExtThe routine finds the locations of the extrema of the error function for
  • MMfiltCofThis routine finds the coefficents of a linear phase FIR filter given the
  • MMgridGiven the input band specifications, this routine generates the
  • MMintCofThis subroutine calculates a set of Lagrange interpolation coefficients
  • MMintValThis function calculates the interpolated value of a polynomial using the
  • MMmmDevThis routine calculates the deviation for a minimax approximation. The
  • MMmodSpecThis routine modifies the specifications of certain types of filters.
  • MMprintSpecThis routine prints the design parameters for an FIR filter.
  • MMtrChebThe inputs to this routine are a set of coefficients for an expansion
  • PrMatrixPretty print a matrix. For the most compact printout, transpose column
  • PrMatrixPretty print a matrix. For the most compact printout, transpose column
  • SlopeMCSlopes at each point for use with polynomial interpolation using a
  • StartDiaryFileThis routine is used for test scripts. All output from the test script
  • StartDiaryFileThis routine is used for test scripts. All output from the test script
  • tDFiltFIR
  • tMMgridSet the filter type
  • tinterp1Script illustrating an undesired behaviour in interp1
  • View all files
Be the first to rate this file! 13 Downloads (last 30 days) File Size: 433 KB File ID: #24662 Version: 1.0
image thumbnail

DFiltFIR

by

Peter Kabal (view profile)

 

Design linear-phase FIR filters - with upper/lower constraints and flexible specifications

| Watch this File

File Information
Description

DFiltFIR designs linear phase Finite (length) Impulse Response (FIR) filters. It can be considered to extend the functionality of the Matlab routine firpm.
DFiltFIR is a fully Matlab implementation. It started as a Fortran program based on the original McClellan-Parks code in [1]. The Fortran code was later modified to include constraints as described by Grenez [2]. At each stage of modification, more comments were added to help in understanding the code. The Fortran code was converted to the C-language and modularized allowing it to be called from programs which designed Nyquist filter and Minimum Phase filters. Finally, when desktop computation power had crossed an appropriate threshold, the program was converted into Matlab. The Matlab code is much easier to maintain and much easier to understand.

Features:
1. Each of the bands in the filter specifications is defined as a sequence of frequencies. At a minimum, bands of non-zero extent are defined by the end points of a band. Additional points within the band can be specified. The desired values, weights, and limits can be defined at those points. The design procedure interpolates these values, weights and limits onto a dense grid of frequency points.
2. The interpolation uses a monotonic cubic interpolation (similar to the Matlab routine pchip). Such an interpolation is applied segment by segment within a band and ensures that the interpolated values between the given frequency points are monotonic.
3. The dense grid of points has a density which varies within a band – the density is higher at the edges of a band than in the middle. The dense grid is applied separately between points within a band, i.e. the points in the band will appear as part of the dense grid.
4. Several versions of the interpolation are available. A function is used to transform the values to be interpolated, interpolation is applied, and then the inverse function is applied. The functions available are the identity function, log, and square root. The first is the default. The second allows for cubic interpolation of dB values of the desired values. The square root option is useful for designing minimum-phase filters (see the DFiltMPFIR package which calls the present package).
5. Relative to firpm, this program features the following.
  a. More flexible definition of values, weights.
  b. Band values are subject to limits.

Example Design:
The following design is a low pass filter (127 coefficients) with two bands: a passband from 0 to 3450 Hz, and a stopband from 4000 to 24000 Hz. The sampling frequency is 48000 Hz. The desired value in the passband is 1 and the desired value in the stopband is 0. The weight in the passband varies from 5 down to 1. This means that the passband ripple amplitude at near zero frequency is 1/5 of the ripple at the passband edge. The weights for the stopband are given in two segments. The limits on the passband are such that at zero frequency, the response is forced to be unity (the upper and lower limits are both 1 at that point). By setting the upper limit to infinity and the lower limit to minus infinity, there are effectively no limits at other frequencies.
  Fs = 24000;
  [B(1:2).Freq] = deal([0 3450], [4000, 6000, 12000]);
  [B(1:2).Value] = deal(1, 0);
  [B(1:2).Weight] = deal([5 1], [10 100 100]);
  B(1).LLimit = [1 -Inf];
  B(1).ULimit = [1 Inf];
  B(2).WeightInt = 'log';
  h = DFiltFIR(63, B, 'bpf', Fs);

1. J. H. McClellan, T. W. Parks and L. R. Rabiner, “A Computer Program for Designing Optimum FIR Linear Phase Digital Filters”, IEEE Trans. Audio and Electroacoustics, vol. 21, pp. 506–526, December 1973.
2. F. Grenez, “Design of Linear or Minimum-Phase FIR Filters by Constrained Chebyshev Approximation”, Signal Processing, vol. 5, pp. 325–332, July 1983.

MATLAB release MATLAB 7.8 (R2009a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.

Contact us