20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

138 CHAPTER 8. H 2 MODEL-ORDER REDUCTION<br />

In the H 2 -approximation problem in continuous time, this enveloping space is<br />

taken <strong>to</strong> be the Hardy space H 2 (Π + ) of complex functions that are square integrable<br />

on the boundary iR of Π + = {s ∈ C | Re(s) > 0} and analytic on the closure of<br />

the right-half plane Π + . Such functions are in one-<strong>to</strong>-one relationship with square<br />

integrable functions on the imaginary axis ∂Π + = iR, which occur as their limits for<br />

s → iω with ω ∈ R. The associated H 2 -inner product of two functions G(s) and H(s)<br />

is given by:<br />

〈G, H〉 = 1 ∫<br />

G(s)H(−s ∗ ) ∗ ds = 1 ∫ ∞<br />

G(iω)H(iω) ∗ dω, (8.1)<br />

2πi i R<br />

2π −∞<br />

where the asterisk denotes complex conjugation. An important subset of H 2 (Π + )is<br />

formed by S, the class of strictly proper rational functions which have all their poles<br />

in the open left-half plane Π − = {s ∈ C | Re(s) < 0}. Such functions correspond <strong>to</strong><br />

transfer functions of continuous time stable and strictly causal LTI systems of finite<br />

order.<br />

When a state-space realization (A, B, C) is employed (thus D = 0) for H(s) (i.e.,<br />

it holds that H(s) =C(sI − A) −1 B), then it is well-known [102] that its squared<br />

H 2 -norm can be computed as:<br />

‖H‖ 2 H 2<br />

= CPC ∗ , (8.2)<br />

in which P satisfies the continuous time Lyapunov equation:<br />

AP + PA ∗ = −BB ∗ . (8.3)<br />

In discrete time, an analogous framework is constituted by the enveloping Hardy<br />

space H 2 (E) of complex functions which are square integrable on the boundary T =<br />

{z ∈ C |z| =1} of E = {z ∈ C |z| > 1} and analytic on the closure of E. These<br />

functions are in one-<strong>to</strong>-one correspondence with square integrable functions on the<br />

unit circle T = {z ∈ C ||z| =1}, which occur as the radial limits of these functions in<br />

H 2 (E). The associated H 2 -inner product is now given by:<br />

∮<br />

1<br />

〈Ĝ, Ĥ〉 =<br />

2πi T<br />

Ĝ(z)Ĥ( 1 dz<br />

)∗<br />

z∗ z = 1<br />

2π<br />

∫ 2π<br />

0<br />

Ĝ(e iΩ )Ĥ(eiΩ ) ∗ dΩ. (8.4)<br />

An important subset of H 2 (E) is the class of strictly proper rational functions which<br />

have all their poles inside the unit disk D = {z ∈ C ||z| < 1}. When represented as<br />

a Laurent series about infinity, a function Ĥ(z) in this class attains the form:<br />

of which the corresponding H 2 -norm satisfies:<br />

Ĥ(z) =Ĥ0 + Ĥ1z −1 + Ĥ2z −2 + ... (8.5)<br />

‖Ĥ‖2 H 2<br />

=<br />

∞∑<br />

|Ĥk| 2 < ∞. (8.6)<br />

k=0

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!