Minimum phase

Minimum phase

In control theory and signal processing, a linear, time-invariant system is said to be minimum-phase if the system and its inverse are causal and stable.[1][2]

For example, a discrete-time system with rational transfer function H(z) can only satisfy causality and stability requirements if all of its poles are inside the unit circle. However, we are free to choose whether the zeros of the system are inside or outside the unit circle. A system is minimum-phase if all its zeros are also inside the unit circle. Insight is given below as to why this system is called minimum-phase.

Contents

Inverse system

A system \mathbb{H} is invertible if we can uniquely determine its input from its output. I.e., we can find a system \mathbb{H}_{inv} such that if we apply \mathbb{H} followed by \mathbb{H}_{inv}, we obtain the identity system \mathbb{I}. (See Inverse matrix for a finite-dimensional analog). I.e.,

\mathbb{H} \, \mathbb{H}_{inv} = \mathbb{I}

Suppose that \tilde{x} is input to system \mathbb{H} and gives output \tilde{y}.

\mathbb{H} \, \tilde{x} = \tilde{y}

Applying the inverse system \mathbb{H}_{inv} to \tilde{y} gives the following.

\mathbb{H}_{inv} \, \tilde{y} = \mathbb{H}_{inv} \, \mathbb{H} \, \tilde{x} = \mathbb{I} \, \tilde{x} = \tilde{x}

So we see that the inverse system \mathbb{H}_{inv} allows us to determine uniquely the input \tilde{x} from the output \tilde{y}.

Discrete-time example

Suppose that the system \mathbb{H} is a discrete-time, linear, time-invariant (LTI) system described by the impulse response h(n) \, \forall \, n \, \in \mathbb{Z}. Additionally, \mathbb{H}_{inv} has impulse response h_{inv}(n) \, \forall \, n \, \in \mathbb{Z}. The cascade of two LTI systems is a convolution. In this case, the above relation is the following:

(h * h_{inv}) (n) = \sum_{k=-\infty}^{\infty} h(k) \, h_{inv} (n-k) =  \delta (n)

where δ(n) is the Kronecker delta or the identity system in the discrete-time case. Note that this inverse system \mathbb{H}_{inv} is not unique.

Minimum phase system

When we impose the constraints of causality and stability, the inverse system is unique; and the system \mathbb{H} and its inverse \mathbb{H}_{inv} are called minimum-phase. The causality and stability constraints in the discrete-time case are the following (for time-invariant systems where h is the system's impulse response):

Causality

h(n) = 0 \,\, \forall \, n < 0

and

h_{inv} (n) = 0 \,\, \forall \, n < 0

Stability

\sum_{n = -\infty}^{\infty}{\left|h(n)\right|} = \| h \|_{1} < \infty

and

\sum_{n = -\infty}^{\infty}{\left|h_{inv}(n)\right|} = \| h_{inv} \|_{1} < \infty

See the article on stability for the analogous conditions for the continuous-time case.

Frequency analysis

Discrete-time frequency analysis

Performing frequency analysis for the discrete-time case will provide some insight. The time-domain equation is the following.

(h * h_{inv}) (n) = \,\! \delta (n)

Applying the Z-transform gives the following relation in the z-domain.

H(z) \, H_{inv}(z) = 1

From this relation, we realize that

H_{inv}(z) = \frac{1}{H(z)}

For simplicity, we consider only the case of a rational transfer function H (z). Causality and stability imply that all poles of H (z) must be strictly inside the unit circle (See stability). Suppose

H(z) = \frac{A(z)}{D(z)}

where A (z) and D (z) are polynomial in z. Causality and stability imply that the poles – the roots of D (z) – must be strictly inside the unit circle. We also know that

H_{inv}(z) = \frac{D(z)}{A(z)}

So, causality and stability for Hinv(z) imply that its poles – the roots of A (z) – must be inside the unit circle. These two constraints imply that both the zeros and the poles of a minimum phase system must be strictly inside the unit circle.

Continuous-time frequency analysis

Analysis for the continuous-time case proceeds in a similar manner except that we use the Laplace transform for frequency analysis. The time-domain equation is the following.

(h * h_{inv}) (t) = \,\! \delta (t)

where δ(t) is the Dirac delta function. The Dirac delta function is the identity operator in the continuous-time case because of the sifting property with any signal x (t).

\delta(t) * x(t) = \int_{-\infty}^{\infty} \delta(t - \tau) x(\tau) d \tau = x(t)

Applying the Laplace transform gives the following relation in the s-plane.

H(s) \, H_{inv}(s) = 1

From this relation, we realize that

H_{inv}(s) = \frac{1}{H(s)}

Again, for simplicity, we consider only the case of a rational transfer function H(s). Causality and stability imply that all poles of H (s) must be strictly inside the left-half s-plane (See stability). Suppose

H(s) = \frac{A(s)}{D(s)}

where A (s) and D (s) are polynomial in s. Causality and stability imply that the poles – the roots of D (s) – must be inside the left-half s-plane. We also know that

H_{inv}(s) = \frac{D(s)}{A(s)}

So, causality and stability for Hinv(s) imply that its poles – the roots of A (s) – must be strictly inside the left-half s-plane. These two constraints imply that both the zeros and the poles of a minimum phase system must be strictly inside the left-half s-plane.

Relationship of magnitude response to phase response

A minimum-phase system, whether discrete-time or continuous-time, has an additional useful property that the natural logarithm of the magnitude of the frequency response (the "gain" measured in nepers which is proportional to dB) is related to the phase angle of the frequency response (measured in radians) by the Hilbert transform. That is, in the continuous-time case, let

H(j \omega) \ \stackrel{\mathrm{def}}{=}\  H(s) \Big|_{s = j \omega}  \

be the complex frequency response of system H(s). Then, only for a minimum-phase system, the phase response of H(s) is related to the gain by

 \arg \left[ H(j \omega) \right] = -\mathcal{H} \lbrace \log \left( |H(j \omega)| \right) \rbrace \

and, inversely,

 \log \left( |H(j \omega)| \right) = \log \left( |H(j \infty)| \right) + \mathcal{H} \lbrace \arg \left[H(j \omega) \right] \rbrace \ .

Stated more compactly, let

H(j \omega) = |H(j \omega)| e^{j \arg \left[H(j \omega) \right]} \ \stackrel{\mathrm{def}}{=}\  e^{\alpha(\omega)}  e^{j \phi(\omega)} = e^{\alpha(\omega) + j \phi(\omega)} \

where α(ω) and ϕ(ω) are real functions of a real variable. Then

 \phi(\omega) = -\mathcal{H} \lbrace \alpha(\omega) \rbrace \

and

 \alpha(\omega) = \alpha(\infty) + \mathcal{H} \lbrace \phi(\omega) \rbrace \ .

The Hilbert transform operator is defined to be

\mathcal{H} \lbrace x(t) \rbrace \ \stackrel{\mathrm{def}}{=}\  \widehat{x}(t) = \frac{1}{\pi}\int_{-\infty}^{\infty}\frac{x(\tau)}{t-\tau}\, d\tau \ .

An equivalent corresponding relationship is also true for discrete-time minimum-phase systems.

Minimum phase in the time domain

For all causal and stable systems that have the same magnitude response, the minimum phase system has its energy concentrated near the start of the impulse response. i.e., it minimizes the following function which we can think of as the delay of energy in the impulse response.

 \sum_{n = m}^{\infty} \left| h(n) \right|^2 \,\,\,\,\,\,\, \forall \, m \in \mathbb{Z}^{+}

Minimum phase as minimum group delay

For all causal and stable systems that have the same magnitude response, the minimum phase system has the minimum group delay. The following proof illustrates this idea of minimum group delay.

Suppose we consider one zero a of the transfer function H(z). Let's place this zero a inside the unit circle (\left| a \right| < 1) and see how the group delay is affected.

a = \left| a \right| e^{i \theta_a} \, \mbox{ where } \, \theta_a = \mbox{Arg}(a)

Since the zero a contributes the factor 1 − az − 1 to the transfer function, the phase contributed by this term is the following.

\phi_a \left(\omega \right) = \mbox{Arg} \left(1 - a e^{-i \omega} \right)
= \mbox{Arg} \left(1 - \left| a \right| e^{i \theta_a} e^{-i \omega} \right)
= \mbox{Arg} \left(1 - \left| a \right| e^{-i (\omega - \theta_a)} \right)
= \mbox{Arg} \left( \left\{ 1 - \left| a \right| cos( \omega - \theta_a ) \right\} + i \left\{ \left| a \right| sin( \omega - \theta_a ) \right\}\right)
= \mbox{Arg} \left( \left\{ \left| a \right|^{-1} - \cos( \omega - \theta_a ) \right\} + i \left\{ \sin( \omega - \theta_a ) \right\} \right)

ϕa(ω) contributes the following to the group delay.

-\frac{d \phi_a (\omega)}{d \omega} =

\frac{ \sin^2( \omega - \theta_a ) + \cos^2( \omega - \theta_a ) - \left| a \right|^{-1} \cos( \omega - \theta_a )

}{ 

\sin^2( \omega - \theta_a ) + \cos^2( \omega - \theta_a ) + \left| a \right|^{-2} - 2 \left| a \right|^{-1} \cos( \omega - \theta_a )

 }
 -\frac{d \phi_a (\omega)}{d \omega} =

\frac{ \left| a \right| - \cos( \omega - \theta_a )

}{ 

\left| a \right| + \left| a \right|^{-1} - 2 \cos( \omega - \theta_a )

 }

The denominator and θa are invariant to reflecting the zero a outside of the unit circle, i.e., replacing a with (a − 1) * . However, by reflecting a outside of the unit circle, we increase the magnitude of \left| a \right| in the numerator. Thus, having a inside the unit circle minimizes the group delay contributed by the factor 1 − az − 1. We can extend this result to the general case of more than one zero since the phase of the multiplicative factors of the form 1 − aiz − 1 is additive. I.e., for a transfer function with N zeros,

\mbox{Arg}\left( \prod_{i = 1}^N  \left( 1 - a_i z^{-1} \right) \right) = \sum_{i = 1}^N \mbox{Arg}\left( 1 - a_i z^{-1} \right)

So, a minimum phase system with all zeros inside the unit circle minimizes the group delay since the group delay of each individual zero is minimized.

Non-minimum phase

Systems that are causal and stable whose inverses are causal and unstable are known as non-minimum-phase systems. A given non-minimum phase system will have a greater phase contribution than the minimum-phase system with the equivalent magnitude response.

Maximum phase

A maximum-phase system is the opposite of a minimum phase system. A causal and stable LTI system is a maximum-phase system if its inverse is causal and unstable. That is,

  • The zeros of the discrete-time system are outside the unit circle.
  • The zeros of the continuous-time system are in the right-hand side of the complex plane.

Such a system is called a maximum-phase system because it has the maximum group delay of the set of systems that have the same magnitude response. In this set of equal-magnitude-response systems, the maximum phase system will have maximum energy delay.

For example, the two continuous-time LTI systems described by the transfer functions

\frac{s + 10}{s + 5} \qquad \text{and} \qquad \frac{s - 10}{s + 5}

have equivalent magnitude responses; however, the first system has a much larger contribution to the phase shift. Hence, in this set, the second system is the minimum-phase system and the first system is the maximum-phase system.

Mixed phase

A mixed-phase system has some of its zeros inside the unit circle and has others outside the unit circle. Thus, its group delay is neither minimum or maximum but somewhere between the group delay of the minimum and maximum phase equivalent system.

For example, the continuous-time LTI system described by transfer function

\frac{ (s + 1)(s - 5)(s + 10) }{ (s+2)(s+4)(s+6) }

is stable and causal; however, it has zeros on both the left- and right-hand sides of the complex plane. Hence, it is a mixed-phase system.

Linear phase

A linear-phase system has constant group delay. Non-trivial linear phase or nearly linear phase systems are also mixed phase.

See also

References

  1. ^ Hassibi, Babak; Kailath, Thomas; Sayed, Ali H. (2000). Linear estimation. Englewood Cliffs, N.J: Prentice Hall. pp. 193. ISBN 0-13-022464-2. 
  2. ^ J. O. Smith III, Introduction to Digital Filters with Audio Applications (September 2007 Edition).

Further reading

  • Dimitris G. Manolakis, Vinay K. Ingle, Stephen M. Kogon : Statistical and Adaptive Signal Processing, pp. 54-56, McGraw-Hill, ISBN 0-07-040051-2
  • Boaz Porat : A Course in Digital Signal Processing, pp. 261-263, John Wiley and Sons, ISBN 0-471-14961-6

Wikimedia Foundation. 2010.

Игры ⚽ Нужно сделать НИР?

Look at other dictionaries:

  • minimum phase shift system — minimaliosios fazės sistema statusas T sritis automatika atitikmenys: angl. minimum phase shift system vok. Phasenminimum System, n rus. минимально фазовая система, f pranc. système à déphasage minimal, m ryšiai: sinonimas – mažiausiosios fazės… …   Automatikos terminų žodynas

  • Phase margin — In electronic amplifiers, phase margin is the difference, measured in degrees, between the phase of the amplifier s output signal and 360°. In feedback amplifiers, the phase margin is measured at the frequency at which the open loop voltage gain… …   Wikipedia

  • phase — [ faz ] n. f. • 1661 astron.; h. 1544 fig.; répandu XIXe; gr. phasis « lever d une étoile » 1 ♦ Astron. Chacun des aspects que présentent la Lune et les planètes à un observateur terrestre, selon leur éclairement par le Soleil. ⇒ apparence. Les… …   Encyclopédie Universelle

  • Phase de la Lune — Phase lunaire En astronomie, une phase lunaire désigne une portion de Lune illuminée par le Soleil et vue à partir de la Terre. La Lune tournant en orbite autour de la Terre, les positions relatives du Soleil, de la Terre et de la Lune changent… …   Wikipédia en Français

  • Phase de la lune — Phase lunaire En astronomie, une phase lunaire désigne une portion de Lune illuminée par le Soleil et vue à partir de la Terre. La Lune tournant en orbite autour de la Terre, les positions relatives du Soleil, de la Terre et de la Lune changent… …   Wikipédia en Français

  • Phase distortion synthesis — is a synthesis method introduced in 1984 by Casio in its CZ range of synths, and similar to Frequency modulation synthesis in the sense that they are both built on phase modulation. Basically a sine wave is played, but by modifying the phase… …   Wikipedia

  • phase — phaseless, adj. phasic, phaseal, adj. /fayz/, n., v., phased, phasing. n. 1. any of the major appearances or aspects in which a thing of varying modes or conditions manifests itself to the eye or mind. 2. a stage in a process of change or… …   Universalium

  • Minimum-shift keying — Passband modulation v · d · e Analog modulation AM · …   Wikipedia

  • Minimum-shift keying — Pour les articles homonymes, voir MSK (homonymie). Diagramme de constellation CPM La modulation Minimum shift keying (MSK) est un type de modulation numéri …   Wikipédia en Français

  • Minimum energy performance standard — MEPS (Minimum Energy Performance Standard) is a specification containing a number of performance requirements for an energy using device, and that effectively limits the maximum amount of energy that may be consumed by a product in performing a… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”