Similar documents
Linear System Fundamentals

One Class of Splitting Iterative Schemes

semigroups are consistent, we shall simply use the notation T (t). Assume further that U and Y are separable Hilbert spaces (the input and output spac

A THEOREM OF ROLEWICZ S TYPE FOR MEASURABLE EVOLUTION FAMILIES IN BANACH SPACES

Manprit Kaur and Arun Kumar

Bogoliubov Transformation in Classical Mechanics

Optimal Hankel norm approximation for infinite-dimensional systems

General System of Nonconvex Variational Inequalities and Parallel Projection Method

TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL

Approximate Analytical Solution for Quadratic Riccati Differential Equation

A CATEGORICAL CONSTRUCTION OF MINIMAL MODEL

Minimal state space realization of MIMO systems in the max algebra

A Note on the Sum of Correlated Gamma Random Variables

Factor Analysis with Poisson Output

On the Unit Groups of a Class of Total Quotient Rings of Characteristic p k with k 3

Practice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions

Unbounded solutions of second order discrete BVPs on infinite intervals

Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

Laplace Transformation

arxiv: v1 [math.ca] 23 Sep 2017

Lecture 8: Period Finding: Simon s Problem over Z N

7 Minimal realization and coprime fraction

Hyperbolic Partial Differential Equations

Computers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order

Lecture 10 Filtering: Applied Concepts

SOME RESULTS ON INFINITE POWER TOWERS

Research Article A New Kind of Weak Solution of Non-Newtonian Fluid Equation

MAE140 Linear Circuits Fall 2012 Final, December 13th

6. KALMAN-BUCY FILTER

IEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Chapter 4. The Laplace Transform Method

THE SPLITTING SUBSPACE CONJECTURE

696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of

Problem Set 8 Solutions

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho

1. The F-test for Equality of Two Variances

The Laplace Transform

REVERSE HÖLDER INEQUALITIES AND INTERPOLATION

arxiv: v2 [math.nt] 30 Apr 2015

RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

Robustness analysis for the boundary control of the string equation

The Hassenpflug Matrix Tensor Notation

The Laplace Transform , Haynes Miller and Jeremy Orloff

Digital Control System

ON THE SMOOTHNESS OF SOLUTIONS TO A SPECIAL NEUMANN PROBLEM ON NONSMOOTH DOMAINS

POINCARE INEQUALITY AND CAMPANATO ESTIMATES FOR WEAK SOLUTIONS OF PARABOLIC EQUATIONS

REPRESENTATION OF ALGEBRAIC STRUCTURES BY BOOLEAN FUNCTIONS. Logic and Applications 2015 (LAP 2015) September 21-25, 2015, Dubrovnik, Croatia

EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal

Exponential stabilization of a Rayleigh beam - actuator and feedback design

March 18, 2014 Academic Year 2013/14

A relationship between generalized Davenport-Schinzel sequences and interval chains

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems

Lecture 7: Testing Distributions

arxiv: v1 [math.mg] 25 Aug 2011

Introduction to Laplace Transform Techniques in Circuit Analysis

FOURIER SERIES AND PERIODIC SOLUTIONS OF DIFFERENTIAL EQUATIONS

Research Article Existence for Nonoscillatory Solutions of Higher-Order Nonlinear Differential Equations

Left invertible semigroups on Hilbert spaces.

1 Routh Array: 15 points

( 1) EE 313 Linear Signals & Systems (Fall 2018) Solution Set for Homework #10 on Laplace Transforms

A Limiting Property of the Matrix Exponential

A new scaled boundary finite element method using Fourier shape. functions

Codes Correcting Two Deletions

On mild solutions of a semilinear mixed Volterra-Fredholm functional integrodifferential evolution nonlocal problem in Banach spaces

THE HAUSDORFF MEASURE OF SIERPINSKI CARPETS BASING ON REGULAR PENTAGON

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem

Suggestions - Problem Set (a) Show the discriminant condition (1) takes the form. ln ln, # # R R

Fault Isolation Filter with Linear Matrix Inequality Solution to Optimal Decoupling

ON ASYMPTOTIC FORMULA OF THE PARTITION FUNCTION p A (n)

Spacelike Salkowski and anti-salkowski Curves With a Spacelike Principal Normal in Minkowski 3-Space

Some Applications of Spanning Trees in K s,t

Modeling the scalar wave equation with Nyström methods

Asymptotic behavior of solutions of mixed problem for linear thermo-elastic systems with microtemperatures

EFFECT ON PERSISTENCE OF INTRA-SPECIFIC COMPETITION IN COMPETITION MODELS

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

Flag-transitive non-symmetric 2-designs with (r, λ) = 1 and alternating socle

A PROOF OF TWO CONJECTURES RELATED TO THE ERDÖS-DEBRUNNER INEQUALITY

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Extension of Inagaki General Weighted Operators. and. A New Fusion Rule Class of Proportional Redistribution of Intersection Masses

Some Sets of GCF ϵ Expansions Whose Parameter ϵ Fetch the Marginal Value

Clustering Methods without Given Number of Clusters

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004

The continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.

INFINITESIMAL GENERATORS OF INVERTIBLE EVOLUTION FAMILIES

INTEGRAL CHARACTERIZATIONS FOR TIMELIKE AND SPACELIKE CURVES ON THE LORENTZIAN SPHERE S *

FILTERING OF NONLINEAR STOCHASTIC FEEDBACK SYSTEMS

Week 3 Statistics for bioinformatics and escience

Multi-dimensional Fuzzy Euler Approximation

Second order conditions for k means clustering: partitions vs. centroids

Research Article Triple Positive Solutions of a Nonlocal Boundary Value Problem for Singular Differential Equations with p-laplacian

EE Control Systems LECTURE 14

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis

where F (x) (called the Similarity Factor (SF)) denotes the function

TAYLOR POLYNOMIALS FOR NABLA DYNAMIC EQUATIONS ON TIME SCALES

Transcription:

Explicit formulae for J pectral factor for well-poed linear ytem Ruth F. Curtain Amol J. Saane Department of Mathematic Department of Mathematic Univerity of Groningen Univerity of Twente P.O. Box 800 P.O. Box 27 9700 AV Groningen 7500 AE Enchede The Netherland The Netherland R.F.Curtain@math.rug.nl A.J.Saane@math.utwente.nl Abtract The tandard way to obtain explicit formula for pectral factorization problem for rational tranfer function i to ue a minimal realization and then obtain formulae in term of the generator A B C and D. For well-poed linear ytem with unbounded generator thee formulae will not alway be well-defined. Intead we ugget another approach for the cla of well-poed linear ytem for which zero i in the reolvent et of A. Such a ytem i related to a reciprocal ytem having bounded generating operator depending on B C D and the invere of A. There are nice connection between well-poed linear ytem and their reciprocal ytem which allow u to tranlate a factorization problem for the well-poed linear ytem into one for it reciprocal ytem the latter having bounded generating operator. We illutrate thi general approach by giving explicit olution to the ub-optimal Nehari problem. Introduction The tandard way to obtain explicit formula for pectral factorization problem for rational tranfer function i to ue a minimal realization G) = D + CI A) B and then obtain formulae in term of the generator A B C and D ee Ball and Ran ). Thee formula typically depend on the controllability and obervability Gramian L B L C or on olution of variou Lyapunov equation. Such an approach ha been extended to certain clae of infinite-dimenional linear ytem ee Curtain and Ran 5 Curtain and Zwart 6 Kaahoek et al. 9) but the limiting factor i the difficulty in manipulating with the unbounded operator B and C. For example in Saane and Curtain and 0 where explicit olution to the ub-optimal Hankel norm approximation problem for exponentially table mooth Pritchard-Salamon ytem and exponentially table analytic ytem were obtained via the olution to the appropriate J pectral factorization problem uing the moothing propertie of thee clae. However it wa not poible to extend thi technique to more general well-poed linear ytem. In general it i not clear that the candidate pectral factor i even well-poed ee Staffan 3). While an obviou approach to get around the problem with unbounded operator would be to obtain factorization via the

dicretized verion obtained by the uual Cayley tranform thi lead to horrible formulae. Intead we ugget tranlating the problem to the analogou one for reciprocal ytem which we now define. To motivate our definition we recall the definition of the tranfer function of a table well-poed linear ytem from Staffan 2: G) Gβ) = β )CβI A) I A) B for all β C + 0..) If 0 ρa) then we can ubtitute β = 0 to obtain the identity ) G) = G0) CA A A B.2) ) = G..3) We note that G ) i the tranfer function of the linear ytem with bounded generating operator A A B CA G0). We define thi linear ytem to be the reciprocal ytem of the well-poed linear ytem with generating operator A B C with A bounded. In addition to the relationhip.3) between their tranfer function the reciprocal ytem ha the ame controllability and obervability Gramian. Lemma.. Let A B C be generating operator of a well-poed linear ytem with tranfer function G. Suppoe that 0 ρa) and G i the tranfer function of it reciprocal ytem with generating operator A A B CA G0). Then the following hold:. C i an infinite-time admiible obervation operator for A iff CA i an infinitetime admiible obervation operator for A. If either C or CA i infinite-time admiible then the obervability Gramian are identical. 2. B i an infinite-time admiible control operator for A iff A B i an infinite-time admiible obervation operator for A. If either C or A B i infinite-time admiible then the obervability Gramian are identical. 3. G H LU Y )) iff G H LU Y )). Proof. From Hanen and Wei 8 ee alo Grabowki 7) we know that C i an infinite-time admiible obervation operator iff the Lyapunov equation Az L C z 2 + L C z Az 2 = Cz Cz 2.4) for all z and z 2 in DA) ha a nonnegative definite olution L C = L C 0. The equation.4) i clearly equivalent to the Lyapunov equation x L C A x 2 + LC A x x 2 = CA x CA x 2.5) for all x and x 2 in X which etablihe the equivalence. Moreover the obervability Gramian are the mallet poitive olution and o the Gramian are identical. 2

2. Thi i dual to part above. 3. Thi follow from.3). The idea i then to tranlate a factorization problem for the well-poed linear ytem with tranfer function G into one for the ytem with tranfer function G the latter having bounded generating operator. We illutrate thi general approach by giving explicit olution to the ub-optimal Nehari problem for the cla of well-poed linear ytem with B C finite rank A bounded which are input output and input-output table. Thi cla include exponentially table well-poed linear ytem with finite rank B C. The aumption that A i not eential and can be removed. A mentioned in Curtain 2 for any iω in the reolvent et of A it i alway poible to define a reciprocal ytem baed on A ω = A iωi. For then the new reciprocal ytem with generating operator A ω A ω B CA ω Giω) and tranfer function G ω atifie G ω and an analogou verion of Lemma. hold. ) = G + iω) 2 The Nehari problem We olve thi problem for the well-poed linear ytem on a Hilbert pace X with generating operator A B C) under the following aumption: A0. The input and output pace are finite-dimenional that i U = C m and Y = C p. A. 0 ρa). A2. B i an infinite-time admiible control operator for {T t)} t 0. A3. C i an infinite-time admiible obervation operator for {T t)} t 0. A4. G ) H C p m ). The ub-optimal Nehari problem i the following: If σ > H G the Hankel norm of the ytem then find all K ) H C p m ) uch that Gi ) + Ki ) σ. K i then called a olution of the ub-optimal Nehari problem. 3

Now ) G) + K) = G + K) = CA I A ) A B + K) + G0) = CA I A ) A B + K where K ) := G0) + K ). Clearly G H C p m ) iff G H C p m ) and K ) H C p m ) iff K ) H C p m ) and in the L -norm G + K = G + K. Thi mean that the Hankel norm of G i equal to that of G. So intead of olving the uboptimal Nehari problem for G we olve the uboptimal Nehari problem for the reciprocal ytem with the bounded generating operator A A B CA. Thi ytem atifie all the condition in Curtain and Ootveen 4: it tranfer function i in H C p m ) and the operator A B and CA are infinite-time admiible with the ame Gramian L B and L C a the original ytem. So we tranlate the reult to thi ytem. Let N σ = ) I L σ 2 B L C LX). Then all olution K to the ub-optimal Nehari problem CA i I A ) A B + K i ) σ are given by where K ) = R )R 2 ) R ) R 2 ) Q ) H C p m ) atifie Q and Λ ) = + CA L B 0 σi m σ 2 σ A B) = Λ ) Q ) I m ) N σ I + A ) CA ) L C A B. Appealing to Theorem. in Wei and Wei 4 we ee that Λ ) H 2 C p+m) p+m) ) 0 σi m but not all component need be in H. Λ ) i invertible for every C + 0 and Λ ) I p H 2 H C p ) and Λ ) i invertible for every C + 0. Furthermore Λ ) atifie the following J pectral factorization problem Ip G iω) I p 0 G iω) I m 0 σ 2 I m 0 I m 4 = Λiω) 0 I m Λiω)

for ω R and ince G ) = G ) we obtain a J pectral factorization over H2 for G: I p 0 Ip Giω) Giω) I m 0 σ 2 = Λ I m 0 I m iω ) Λ 0 I m iω ) for ω R. The final olution to our original uboptimal Nehari problem i K ) = G0) + K ) = G0) + R ) R2 where ) R ) R 2 = Λ ) Q0 ) I m and Q 0 ) H C p m ) atifie Q 0. From Curtain and Ootveen 4 we have the following formula. Λ = ) 0 I σ m + σ 2 CA L B σ A B) I A ) ) N σ ) CA ) σ L CA B for all C + 0. While it i tempting to try to write thi in term of it reciprocal we know that thi will not in general be well-defined ee Staffan 3). So we leave the explicit olution a it tand. A imilar approach to the optimal Hankel norm problem for well-poed linear ytem ha been taken in Curtain and Saane 3. For other application of reciprocal ytem ee Curtain 2. Reference J.A. Ball and A.C.M. Ran. Optimal Hankel norm model reduction and Wiener-Hopf factorization I: The canonical cae. SIAM Journal on Control and Optimization 252):362 382 987. 2 R.F. Curtain. Reciprocal of Regular linear ytem: a urvey. Proceeding MTNS 2002. 3 R.F. Curtain and A.J. Saane. Hankel norm approximation for well-poed linear ytem. manucript 4 R.F. Curtain and J.C. Ootveen. The Nehari problem for nonexponentially table ytem. Integral Equation and Operator Theory 3:307 320 998. 5 R.F. Curtain and A. Ran. Explicit formula for Hankel norm approximation of infinitedimenional ytem. Integral Equation and Operator Theory 3:455 469 989. 5

6 R.F. Curtain and H.J. Zwart. The Nehari problem for the Pritchard-Salamon cla of infinite-dimenional linear ytem: a direct approach. Integral Equation and Operator Theory 8:30 53 994. 7 P. Grabowki. On the pectral-lyapunov approach to parametric optimization of ditributed parameter ytem. IMA Journal of Mathematical Control and Information 7:37 338 990. 8 S. Hanen and G. Wei. New reult on the operator Carleon meaure criterion. IMA Journal of Mathematical Control and Information 4:3 32 997. 9 M.A. Kaahoek C.V.M van der Mee and A.C.M. Ran. Wiener-Hopf factorization of tranfer function of extended Pritchard-Salamon realization. Math. Nachrichten 96:7 02 998. 0 A.J. Saane and R.F. Curtain. Sub-optimal Hankel norm approximation for the analytic cla of infinite-dimenional ytem. Integral Equation and Operator Theory to appear. A.J. Saane and R.F.Curtain. Optimal Hankel norm approximation for the Pritchard- Salamon cla of infinite-dimenional ytem. Integral Equation and Operator Theory 39:98-26 200. 2 O.J. Staffan. Well-Poed Linear Sytem I: General Theory. To be publihed 2002. draft available at http://www.abo.fi/ taffan/publ.htm. 3 O.J. Staffan. Admiible factorization of Hankel operator induced by well-poed linear ytem. Sytem and Control Letter 37:30 307 999. 4 M. Wei and G. Wei. Optimal control of table weakly regular linear ytem. Mathematic of Control Signal and Sytem 0:287 330 997. 6