<< Chapter < Page Chapter >> Page >

Prony's method for time-domain design of iir filters

The problem of designing an IIR digital filter with a prescribed time-domain response is addressed in this section.Most formulations of time-domain design of IIR filters result in nonlinear equations for the same reasons as for frequency-domaindesign. Prony, in 1790, derived a special formulation for the analysis of elastic properties of gases, which resulted in linearequations. A more general form of Prony's method can be applied to the IIR filter design by use of a matrix description [link] , [link] .

The transfer function of an IIR filter is given by

H ( z ) = B ( z ) A ( z ) = b 0 + b 1 z - 1 + + b M z - M 1 + a 1 z - 1 + + a N z - N = h 0 + h 1 z - 1 + h 2 z - 2 +

and the impulse response h ( n ) is related to H ( z ) by the z transform.

H ( z ) = n = 0 h ( n ) z - n

[link] can be written

B ( z ) = H ( z ) A ( z )

which is the z-transform version of convolution. This convolution can be written as a matrix multiplication. Using the first K+1terms of the impulse response, this is written

b 0 b 1 b M 0 0 = h 0 0 0 0 h 1 h 0 0 h 2 h 1 h 0 h L h 0 1 a 1 a N 0 0

In order to uncouple the calculations of the a n and the b n , the matrices are partitioned to give

b 0 = H 1 h 1 H 2 1 a *

where b is the vector of the M + 1 numerator coefficients of [link] , a * is the vector of the N denominator coefficients ( a o = 1 ), h 1 is the vector of the last ( K - M ) terms of the impulse response, H 1 is the M + 1 by N + 1 partition of [link] , and H 2 is the ( K - M ) by N remaining part. The lower K - M equations are written

0 = h 1 + H 2 a *

or

h 1 = - H 2 a *

which must be solved for a * . The upper M + 1 equations of [link] are written

b = H 1 a

which allows the calculation of b .

If L = N + M , then H 2 is square. If H 2 is nonsingular, [link] can be solved exactly for the denominator coefficients in a * , which are augmented by the unity term to give a . From [link] , the numerator coefficients in b are found. If H 2 is singular [link] and there are multiple solutions, a lower order problem can be posed. If there are nosolutions, the methods of the next section must be used.

Note that any order numerator and denominator can be prescribed. If the filter is in fact an FIR filter, a is unity and a * does not exist. Under these conditions, [link] states that b n = h n , which is one of the cases of FIR frequency sampling covered in Section 3.1 of [link] . Also note that there is no control over the stability of the filter designed by thismethod.

Although Prony's method, applied to the time-domain design problem here, is similar to the solution of the frequency-samplingIIR design problem, there are important differences. The inverse DFT is used to obtain the matrix in the frequency domain problem, which iscyclic convolution. Equation [link] is noncyclic convolution and the K + 1 terms of h ( n ) , used to form H , result from a truncation of the infinitely long sequence.

An approximate solution or the least equation error problem

In order to obtain better practical filter designs, the interpolation scheme of the previous section is extended to give an approximation designmethod [link] . It should be noted at the outset that the method developed in this section minimizes an equation-error measure and not the usualfrequency-response error measure.

The number of samples specified, L + 1 , will be made larger than the number of filter coefficients, M + N + 1 . This means that H 2 is rectangular and, therefore, [link] cannot in general be satisfied. To formulate an approximation problem, a length- ( L + 1 ) error vector ε is introduced in [link] and [link] to give

b 0 = H 0 a + ε

[link] becomes

h 1 - ε = - H 2 a *

where now H 2 is rectangular with ( L - M ) > N . Using the same methods as used to derive [link] , the error ε is minimized in a least-squared error sense by the solution of the normal equations [link]

H 2 T h 1 = - H 2 T H 2 a *

If the equations are not singular, the solution is

a * = - [ H 2 T H 2 ] - 1 H 2 T h 1 .

If the normal equations are singular, the pseudo-inverse [link] can be used to obtain a minimum norm or reduced order solution.

The numerator coefficients are found by the same techniques as before in [link]

b = H 1 a

which results in the upper M + 1 terms in ε being zero and the total squared equation error being minimum.

As is true for LS-error design of FIR filters, [link] is often numerically ill-conditioned and [link] should not be used to solve for a * . Special algorithms such as those used by Matlab and LINPACK [link] , [link] should be employed.

Various modifications can be made to the form of Prony's method presented. After the denominator is found by minimizingthe equation error, the numerator can be found by minimizing the solution error. It is possible to mix the exact and approximatemethods. The details can be found in [link] , [link] , [link] .

Several modifications to Prony's method have been made to use it to minimize the solution error. Most of these iterativelyminimize a weighted-equation error with Prony's method and update the weights from the previous determination of a [link] , [link] .

If an LS-error, time-domain approximation is the desired result, a minimization technique can be applied directly to thesolution error. The most successful method seems to be the Gauss- Newton algorithm with a step-size control. This combined withProny's method to find starting parameters is an effective design tool.

Questions & Answers

I'm interested in biological psychology and cognitive psychology
Tanya Reply
what does preconceived mean
sammie Reply
physiological Psychology
Nwosu Reply
How can I develope my cognitive domain
Amanyire Reply
why is communication effective
Dakolo Reply
Communication is effective because it allows individuals to share ideas, thoughts, and information with others.
effective communication can lead to improved outcomes in various settings, including personal relationships, business environments, and educational settings. By communicating effectively, individuals can negotiate effectively, solve problems collaboratively, and work towards common goals.
it starts up serve and return practice/assessments.it helps find voice talking therapy also assessments through relaxed conversation.
miss
Every time someone flushes a toilet in the apartment building, the person begins to jumb back automatically after hearing the flush, before the water temperature changes. Identify the types of learning, if it is classical conditioning identify the NS, UCS, CS and CR. If it is operant conditioning, identify the type of consequence positive reinforcement, negative reinforcement or punishment
Wekolamo Reply
please i need answer
Wekolamo
because it helps many people around the world to understand how to interact with other people and understand them well, for example at work (job).
Manix Reply
Agreed 👍 There are many parts of our brains and behaviors, we really need to get to know. Blessings for everyone and happy Sunday!
ARC
A child is a member of community not society elucidate ?
JESSY Reply
Isn't practices worldwide, be it psychology, be it science. isn't much just a false belief of control over something the mind cannot truly comprehend?
Simon Reply
compare and contrast skinner's perspective on personality development on freud
namakula Reply
Skinner skipped the whole unconscious phenomenon and rather emphasized on classical conditioning
war
explain how nature and nurture affect the development and later the productivity of an individual.
Amesalu Reply
nature is an hereditary factor while nurture is an environmental factor which constitute an individual personality. so if an individual's parent has a deviant behavior and was also brought up in an deviant environment, observation of the behavior and the inborn trait we make the individual deviant.
Samuel
I am taking this course because I am hoping that I could somehow learn more about my chosen field of interest and due to the fact that being a PsyD really ignites my passion as an individual the more I hope to learn about developing and literally explore the complexity of my critical thinking skills
Zyryn Reply
good👍
Jonathan
and having a good philosophy of the world is like a sandwich and a peanut butter 👍
Jonathan
generally amnesi how long yrs memory loss
Kelu Reply
interpersonal relationships
Abdulfatai Reply
What would be the best educational aid(s) for gifted kids/savants?
Heidi Reply
treat them normal, if they want help then give them. that will make everyone happy
Saurabh
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Digital signal processing and digital filter design (draft). OpenStax CNX. Nov 17, 2012 Download for free at http://cnx.org/content/col10598/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Digital signal processing and digital filter design (draft)' conversation and receive update notifications?

Ask