Electrical computers: arithmetic processing and calculating – Electrical digital calculating computer – Particular function performed
Reexamination Certificate
2000-04-14
2003-02-25
Malzahn, David H. (Department: 2124)
Electrical computers: arithmetic processing and calculating
Electrical digital calculating computer
Particular function performed
C708S323000
Reexamination Certificate
active
06526429
ABSTRACT:
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates to filter design, and more particularly to a method for determining coefficients for an all pass filter.
2. Description of the Related Art
All pass filters are used to realize a given group delay response in linear systems. Various methods have been proposed for designing all pass filters to satisfy a given group delay specification. For example, in “Design of Recursive Group Delay Filters by Autoregressive Modeling,” IEEE Trans. Acoust., Speech and Signal Processing, Vol. ASP-30, pp. 632-637, August 1982, by B. Yegnanarayana, autoregressive modeling is used to derive the all pass filter coefficients from the cepstral coefficients. G. R. Reddy and M. N. S. Swamy, in “Digital All-pass Filter Design Through Discrete Hilbert Transform,” pp. 646-649, ISCAS 1990, used the discrete Hilbert transform relationship between the log magnitude function and the phase of the Fourier transform to obtain the filter coefficients of an all pass filter. P. A. Bernhardt, in “Simplified Design of Higher-order Recursive Group Delay Filters,” IEEE Trans. Acoust. Speech and Signal Processing, vol. ASSP-28, pp. 498-503, October 1980, used an iterative computer algorithm for solving a set of non-linear equations to get the parameters of the filter.
Although the above method provide alternatives to determining coefficients for all pass filter designs, the methods and algorithms are expensive in terms of time for processing. Therefore, a need exists for a method, which reduces the computation time and the complexity for determining all pass-filter coefficients.
SUMMARY OF THE INVENTION
A method for designing an all pass filter, in accordance with the present invention, includes providing a plurality of electrical or optical signals having a group delay, and normalizing the group delay to determine a normalized group delay function. Cepstral coefficients of the normalized group delay function are determined for to a denominator polynomial of a transfer function of the all pass filter wherein a magnitude function and the normalized group delay function of the denominator polynomial are related through the cepstral coefficients. The denominator polynomial coefficients are determined through a non-linear recursive difference equation by employing the cepstral coefficients.
In other methods, the step of determining cepstral coefficients may include the steps of determining the cepstral coefficients c(k) using the relation:
τ
(
w
)
=
∑
k
=
1
∞
kc
(
k
)
cos
kw
where &tgr;(w) is a group delay function k is an index and w is 2&pgr;f and f is a normalized sampling frequency. The relation
τ
(
w
)
=
∑
k
=
1
∞
kc
(
k
)
cos
kw
may be solved by an Inverse Fast Fourier Transform (IFFT). The step of determining the denominator polynomial coefficients through a non-linear recursive difference equation by employing the cepstral coefficients may further include employing the non-linear recursive difference equation:
a
n
=
∑
k
=
0
n
(
k
n
)
c
(
k
)
a
n
-
k
,
n
>
0
to determine the denominator polynomial coefficients, wherein a
n
represents the denominator polynomial coefficients, c(k) represents the cepstral coefficients, k is an index, and n is the order of the all pass filter. The transfer function may given by H(z)=N(z)/D(z) where N(z) is a numerator polynomial and D(z) is the denominator polynomial. The coefficients of the denominator polynomial are preferably the same as the numerator polynomial with mirror-image symmetry such that determining the denominator polynomial coefficients determines H(z). The method steps of the present invention may be implemented on a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform method steps for is designing an all pass filter given a plurality of electrical or optical signals having a group delay.
A method for designing an all pass filter includes the steps of providing a plurality of electrical or optical signals having a group delay, normalizing the group delay to determine a normalized group delay function, and determining coefficients for the all pass filter. The step of determining coefficients for the all pass filter includes determining cepstral coefficients of the normalized group delay function by employing only one inverse Fourier transform operation and determining polynomial coefficients of the all pass filter through a non-linear recursive difference equation by employing the cepstral coefficients.
In other methods, the step of determining cepstral coefficients may include the steps of determining the cepstral coefficients c(k) using the relation:
τ
(
w
)
=
∑
k
=
1
∞
kc
(
k
)
cos
kw
where &tgr;(w) is a group delay function k is an index and w is 2&pgr;f and f is a normalized sampling frequency. The non-linear recursive difference equation may include:
a
n
=
∑
k
=
0
n
(
k
n
)
c
(
k
)
a
n
-
k
,
n
>
0
,
wherein a
n
represents the polynomial coefficients, c(k) represents the cepstral coefficients, k is an index, and n is the order of the all pass filter. The transfer function may be given by H(z)=N(z)/D(z) where N(z) is a numerator polynomial and D(z) is a denominator polynomial, wherein the polynomial coefficients of the denominator polynomial are the same as the numerator polynomial with mirror-image symmetry such that determining the denominator polynomial coefficients determines H(z). The method steps of the present invention may be implemented on a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform method steps for designing an all pass filter given a plurality of signals having a group delay.
REFERENCES:
patent: 3997772 (1976-12-01), Crochiere et al.
patent: 5144569 (1992-09-01), Kunold
patent: 5654909 (1997-08-01), Sun
patent: 5687101 (1997-11-01), Lee
Reddy et al., “Digital All-Pass Filter Design Through Discrete Hilbert Transform,” 1990 IEEE International Symposium on Circuits and Systems, vol. 1, pp. 646-649, May 1990.
Paul A. Bernhardt, “Simplified Design of High-Order Recursive Group-Delay Filters,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-28, No. 5, pp. 498-503, Oct. 1980.
B. Yegnanarayana, “Design of Recursive Group-Delay Filters by Autoregressive Modeling,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-30, No. 4, pp. 632-637, Aug. 1982.
Lai Yhean-Sen
Rajamani Kannan
Agere Systems Inc.
Malzahn David H.
LandOfFree
Method for designing all pass digital filters does not yet have a rating. At this time, there are no reviews or comments for this patent.
If you have personal experience with Method for designing all pass digital filters, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Method for designing all pass digital filters will most certainly appreciate the feedback.
Profile ID: LFUS-PAI-O-3156185