×

Efficient digital filter design tool for approximating an FIR filter with a low-order linear-phase IIR filter

  • US 7,373,367 B2
  • Filed: 04/19/2004
  • Issued: 05/13/2008
  • Est. Priority Date: 04/19/2004
  • Status: Expired due to Fees
First Claim
Patent Images

1. A method of approximating an FIR filter with low-order linear-phase IIR filters by the rational Arnoldi algorithm with adaptive orders containing the following steps:

  • a) initialize the first vector of the Krylov sequence for each expansion point;

    b) in the jth iteration of the algorithm, choosing an expansion frequency wherein the heuristics of selecting expansion frequencies in advance for the proposed rational Arnoldi method we given by(a) low-pass filters;

    the proposed method with the expansion point ω

    1=0;

    (b) high-pass filters;

    the special structures of state-space matrices used to present the duality between low-pass and light pass filters;

    let state matrices become Ā

    =−

    A, b=b, c=c, and h0=−

    h0, the expansion point ω

    1=0 chosen to perform the Arnoldi algorithm;

    when the corresponding orthonormal matrix Vq is obtained and then the high-pass IIR filter, which satisfies the same specifications as the original FIR filter; and

    (c) band-pass/band-stop filters;

    the passband edge and the stopband edge frequencies being the appropriate candidate expansion points in meeting the specifications of the design, and other expansion points with uniform spacing recommended to be selectedsuch that the frequency gives the greatest difference between the (j+1)st-order output moment of the original FIR filter H(z) and that of the lower-order IIR filter Ĥ

    (z) wherein the expression of output moment errors between the ĵ

    ith-order moments H

    i)(zi) and Ĥ



    i)(zi) at each expansion point zi are expressed as follows;



    H ( j ^ i )

    ( z i )
    - H ^ ( j ^ i )

    ( z i )


    =

    h π



    c T

    r ( j ^ i - 1 )

    ( z i )


    ,

    where




    h π



    ( z i )
    =

    j








    r ( j - 1 )

    ( z i )




    is the normalization coefficient when an expansion frequency zi is selected in the jth iteration;

    vector c contains the last n impulse response coefficients of a FIR filter with length n+1; and

    r(j−

    1)
    (zi) is the residual vector in the (j−

    1)st iteration of the disclosed adaptive rational Arnoldi algorithm at the expansion frequency zi;

    c) after the choosing the expansion point in jth iteration being determined, the single-point Arnoldi method applied at the expansion point to generate the new orthnormal vector; and

    d) determine a new residual at each expansion point for next iteration;

    whereby, after the giving total iteration number of the algorithm, outputting the resulting orthogonal projection matrix.

View all claims
  • 1 Assignment
Timeline View
Assignment View
    ×
    ×