arxiv: v1 [math.oc] 22 Jan 2019

Transcription

1 SAMPLED-DATA OUTPUT REGULATION OF WELL-POSED INFINITE-DIMENSIONAL SYSTEMS MASASHI WAKAIKI AND HIDEKI SANO arxiv: v1 [math.oc] 22 Jan 219 Abstract. We solve the sample-data control problem of output tracking and disturbance rejection for unstable well-posed linear infinite-dimensional systems. We obtain a sufficient condition for the existence of finite-dimensional sampled-data controllers that are solutions of this control problem with constant disturbance and reference signals. To obtain this result, we study the problem of output tracking and disturbance rejection in the discrete-time setup and propose a design method of finite-dimensional controller by using a solution of a Nevanlinna-Pick interpolation problem with both interior and boundary conditions. To illustrate our main result, we describe an application to a delay system of retarded type with output delays. Key words. Frequency-domain methods, output regulation, sampled-data control, state-space methods, well-posed infinite-dimensional systems AMS subject classifications. 93B52, 93C5, 93C25, 93C35, 93C57, 93D15 1. Introduction. Due to the development of computer technology, we implement digital controllers for continuous-time plants. We call such closed-loop systems sampled-data systems. In addition to their practical motivation, sampled-data systems yield theoretically interesting problems arising from the combination of both continuous-time and discrete-time dynamics. Sampled-data control theory for infinitedimensional systems has been developed, e.g., in [11 13, 18 21, 29 31, 35]. Several specifically relevant studies will be cited below again. In this paper, we study the problem of sampled-data output regulation for unstable well-posed systems. The main objective in our control problem is to find a finite-dimensional digital controller achieving convergence of the output to a given constant reference signal in the presence of external constant disturbances. A theory for well-posed systems has been extensively developed; see, e.g., the survey [37] and the book [36]. Well-posed systems allow unbounded control and observation operators and provide a framework to formulate control problems for systems governed by partial differential equations with point control and observation and by delay differential equations with delayed control action and measurements. Our output regulation method is based on the internal model principle, which was originally developed for finite-dimensional systems in [7] and was later generalized for infinite-dimensional systems with finite-dimensional and infinite-dimensional exosystems in [25 28] and references therein. In particular, output regulation of nonsmooth periodic signals has been extensively studied as repetitive control [1]. For regular systems, which is a subclass of well-posed systems, the authors of [2, 24, 25, 42] have provided design methods of continuous-time controllers for robust output regulation. For stable well-posed systems with finite-dimensional exosystems, low-gain controllers suggested by the internal model principle have been constructed for the continuoustime setup in [22,32] and for the sampled-data setup in [11 13,21]. The difficulty of the problem we consider arises from the instability of well-posed systems. If the system is unstable, then low-gain controllers cannot achieve closed-loop stability. Ukai and Submitted to the editors DATE. Funding: This work was supported by JSPS KAKENHI Grant Numbers JP17K Graduate School of System Informatics, Kobe University, Nada, Kobe, Hyogo , Japan wakaiki@ruby.kobe-u.ac.jp, sano@crystal.kobe-u.ac.jp). 1

2 2 MASASHI WAKAIKI AND HIDEKI SANO v u <latexit sha1_base64="qxfejlz+ao3iuapy7fgab7zpre8=">aaac6xichvfllwrbec7jtczrczg4iarku2nc3eslo7rszbsymzmtq3o9vrmzmy3wrn/wil4yq4id/giou7opgj4kji4qb2dkqq1ks7qr+q7+ueksmr3pmrhxquxqbmltzym9rednvhe/pxfhssmuytgkl211di8jllna575gq47ldmsqlgou5mr5tiw5hrflsl912lqlfyxf4kbue5spj+cmvuqyauvyrftnwfyoubvlxkbwlejyllwkc3he5jeaz/bloujccylb2/gxwuwaytymabawk+xqj8ojlggyidmhrebdmustdpintuilbpipgftqhjdqldmpgqkrztdml2sbdimi5xbyeybzhc3zgw7zar3z7vssinwpvqzi36lzm5lt3+pde/2vz5h3y/gt9wfiodsm7qev3/wvohdan5tpynvnnl3zsmwp7wkkntojuumpwnspb+89l4vdwqie4bp15rgf8iy6iysv5tkcwzwclqarfr/jz2bliqlhuluysmzmrioowqamwsjncrjmyb5skkz79+asrubaksm7yofywc9vgijoh3wx5fqdkugumw==</latexit> <latexit <latexit sha1_base64="vf3qrdyv628yl8sjtzzare31a=">aaac6xichvfllwrbec7jtczrczg4iarku2nc3eslo7rszbsymzmtq3o9vrmzmy3wrn/wil4yq4id/giou7opgj4kji4qb2dkqq1ks7qr+q7+ueksmr3pmrhxquxqbmltzym9rednvhe/pxfhssmuytgkl211di8jllna575gq47ldmsqlgou5mr5tiw5hrflsl912lqlfyxf4kbue5spj+cmvuqyauvyrftnwfyoubvlxkbwlyjyllwkc3he5jeaz/bloujccylb2/gxwuwaytymabawk+xqj8ojlggyidmhrebdmustdpintuilbpipgftqhjdqldmpgqkrztdml2sbdimi5xbyeybzhc3zgw7zar3z7vssinwpvqzi36lzm5lt3+pde/2vz5h3y/gt9wfiodsm7qev3/wvohdan5tpynvnnl3zsmwp7wkkntojuumpwnspb+89l4vdwqie4bp15rgf8iy6iysv5tkcwzwclqarfr/jz2bliqlhuluysmzmrioowqamwsjncrjmyb5skkz79+asrubaksm7yofywc9vgijoh3wx5fqdjtqumg==</latexit> <latexit y <latexit sha1_base64="7qdr5hkurub/sc3lsuomhjpngii=">aaac6xichvfllwrbec7jtczrczg4iarku2nc3eslo7rszbsymzmtq3o9vrmzmy3wrn/wil4yq4id/giou7opgj4kji4qb2dkqq1ks7qr+q7+ueksmr3pmrhxquxqbmltzym9rednvhe/pxfhssmuytgkl211di8jllna575gq47ldmsqlgou5mr5tiw5hrflsl912lqlfyxf4kbue5spj+cmvuqyauvyrftnwfyoubvlxkbwlujyllwkc3he5jeaz/bloujccylb2/gxwuwaytymabawk+xqj8ojlggyidmhrebdmustdpintuilbpipgftqhjdqldmpgqkrztdml2sbdimi5xbyeybzhc3zgw7zar3z7vssinwpvqzi36lzm5lt3+pde/2vz5h3y/gt9wfiodsm7qev3/wvohdan5tpynvnnl3zsmwp7wkkntojuumpwnspb+89l4vdwqie4bp15rgf8iy6iysv5tkcwzwclqarfr/jz2bliqlhuluysmzmrioowqamwsjncrjmyb5skkz79+asrubaksm7yofywc9vgijoh3wx5fqdmhaung==</latexit> <latexit sha1_base64="55ujujwtfcvvrteupuu7npi1ho=">aaac7xichvhpl8rqeb71u34tlhixsshisjm6ecfh4oi1luzkpu+9bkvr3b3wq1/gf3hiieret8c25c3mxbnycojc4ozrsvqtdnezpvm/nmvx5juil7aejdg9ly1nzs2taudnr2dfckevtwfkfsmsxjoslxvg3dz4jllgl4iniq6zhdngtlgqw5wj5byz7phbkcvf22buuw5bvc1aoccgnmg8zimtqft+t4tpo3hfhqza5um2nuwxrbdjz/kgojjkywsigfgzahcqhtguncqd5kiidjptbbgysaoof6odtlwmnefzc1iekzkoir3kg26astxvjcpqkuw3tkxaikc62nh7fnukxq8og5bcn4j+f4jld4gy/49muvmopreuvvfhnmrfqszoqfv2xzzmpypot9qfjo9aj7aat3/8vphfynz3xxnpzp7uiqn9upaentdwiqalj1ntutvaf9nli+eonuat6xkmd3hfysjki3m6yjyoqkxbat/h+dnymuhpmniwj5izs/gi22aqhmgm5jgjmzapc5che/fgeq7hrnguxevqoaqxkgxpx++mhl2dqvssja=</latexit> <latexit S H <latexit sha1_base64="98blja2zuctayfyuc3odsp++kg=">aaac93ichvhnsuraec7j+hf/rveysbdxumtdupgiejl1skd1heezyndjtgnjpxosnoexzavsc3jyw/7aloonn2fbvkghnwe8ajgxym1myioqbw6q/qr+r7uvdmbfjfgvogwoj9dff9pn9a4ndw5mrbxir4uul7i+9mn1hvccsulwmjj14oqm8+rvojsllbyxtopi+grvdi+kbhqkpsczdpgsqzedvhvafijkvvttdn2/fljwp45glt9pjedpmm882yrvmnlxvhmrlmszmmlgxl4gvbllibcnpnienffdbhrp4wegbplgcg4i+eliaebc2ctfhiuuiyxnogkncglvxqmce7tbepvmprrwdw5prwnbpfkkrjoyytoa5hua1nuahxuldq1pxotf6s4o8+byodz87vp+9l2wr17d9iprdczdbujzlvqv/19m7/aentems8u3e6fjby7pgacebans6o7b1qjv7l3n51cm4kn8jvful194gf+om6p+4/5z5ivfwatbws/h+djym8lzmlowz7ilx9ir98jngicpmumslmbxwiic3fstjue/nbon44fx19hvlxodkecjpdhj6b7habut</latexit> <latexit sha1_base64="1eo9ctrecawni11/nelmtzwopho=">aaac93ichvhnsuraec7juqvrxwf1ingrbu8dbuvu3is9elrv3guiqydtds2djoh6rmydfmcvoahd/6aonjygqtxtgulbx9bpcp48wbnjiiquhw6q/qr+r7uvdmbfjfgvgoz2j91fp7s2wv293z91pv53rcu+dxq5xnxl3647lcis6f4xgst+xiqcuy5khecjelmvldjysr8tajrav/1wewjneeytvapm2k7vcjuzksoqpggatu+led1j1xs2h7t6y6t8uyjzgtwptrx5afauialouxs6g1gpuewupv1m//ahjl44eivpocgqfmsgufexxesqagiw4wysjaikeq5nmakbpwqofuwqjdor9cpmkkkzk3nkgg7diukfrjzcebweo/xfi/wbk/x4v2tonfovqvo3mlxevdq3rxyup8vyyovyf2z9qhjqtak7bqt9/8vpnd4z5or48xch73qppcz6yggngqjuyo29ko/d66xzich4lhcr9vqc97eivn1blvu3mp5vj8npgwov1gn8gsxm5c3pw3er26lc64k4yhgeyozn+gcmygvni727ca5/4k9rn3amq+oovwqpzx+eghg6spgxbui</latexit> yd <latexit sha1_base64="izhb1jgzxutupvezshzj/hhrxf=">aaac8xichvhnssnaeb7jf/yrehg8fisihspebmwt6mvjq9yqrzqkxevszsykaaggvoaviopfbq/iswhe1kt48bheo4ixd7tikioe3zn9pv5vt3mgi7gno/4kktrv3dhz1qz29ffdscghdc+uucblmlaw3q1d95jgkmv87gu24bhmtwzbskz5qzhpvpnrcvuu+twhbvl6sfjtbuo+qyxytn5gjs4dxfcsnkqrecmwra9mkqvuwihiu1a8wfcpwwtxo6oqs2asq4v/dlqosebkktt2b3koggmvmacbhj8igxo4ngxawqhmk2icdmpyiheqz1uilbospgftqhzdpldmpfqkrzq9ml2sbdimi5xizdon7jot7jnv7gi779qhwego231mgbts5zcgn7i6uv/7is8j7sfll+yhzuuptdpvx7/wxduttc3mql/2zfz7pzyu94nqtjqa3tgbgtxd/efv+zxxyajp8yn6coipeewdkdux8yznvo5apcfq38f4m1iftmqy1nltiyxfamrdmapjmelzniufwiyuzojeq7ieg7hvpovaovzom6vks8qzhi+mxlwdmlmyba==</latexit> d <latexit sha1_base64="shwvl11nxv2ee7vhy/wiblxmq=">aaac9nichvhntttaeb5cogd+qntb4hirgrchajwlhrnqlzcauiqgfn+wcwq67vloxgplrfoc/rqqvuriys48a6oeqf2gdjwckjhipxsaxphcaecxtqd2w/m+3y9ywwcrzhixydyynbo+oxiqdo2pje5lzt+trx5zdbmfdsxflizzigjllkl5rfgtsbkpmcjvrua73r5aoufefflh7gdsg3pdcxf4byze2tkvnsluvwmpucuxoqouuulj2p75bjvl3pxm41ed7281epy6rzu2rkclje1pipaylcpdzup87axc8mggjnjaqejmsqativrqoafcqng2jisffpez6adkngbvmwowisqbtlp3qgsjr3nkoubdmtglzizdzm4zkeyhdp8qgv8f+jwkmqxtlm7zv57lampo8u/73lmsjh8pulesjxk1tsnkdwo//xlv8g41py7vq/2iia9n2kpopuksjfed+y+ruvtl255dxm+wccf+jf68gmv8iq6i1tx9v4g2/wkkg1wuz/gh8fwqahhudsofdbezimegvmygwa4zkswxtyhwrd+x1+ww/4o+wp35r95abfqgxknndwx5tjawzctao=</latexit> ud - yref <latexit sha1_base64="pzccwqtgpp1ssio3xfmndmuudb8=">aaac8xichvhnssnaeb7jf/yrehg8fisihspebmwt6mvjq9yqrzqkxevszsykaaggvoaviopfbq/iswhe1kt48bheo4ixd7tikioe3zn9pv5vt3mgi7gno/4kktrv3dhz1qz29ffdscghdc+uucblmlaw3q1d95jgkmv87gu24bhmtwzbskz5qzhpvpnrcvuu+twhbvl6sfjtbuo+qyxytn5gjs4dxfcsnkqrecmwra9mkqvusihiu1a8wfcpwwtxo6oqs2asq4v/dlqosebkktt2b3koggmvmacbhj8igxo4ngxawqhmk2icdmpyiheqz1uilbospgftqhzdpldmpfqkrzq9ml2sbdimi5xizdon7jot7jnv7gi779qhwego231mgbts5zcgn7i6uv/7is8j7sfll+yhzuuptdpvx7/wxduttc3mql/2zfz7pzyu94nqtjqa3tgbgtxd/efv+zxxyajp8yn6coipeewdkdux8yznvo5apcfq38f4m1iftmqy1nltiyxfamrdmapjmelzniufwiyuzojeq7ieg7hvpovaovzom6vks8qzhi+mxlwdjv2yaa==</latexit> <latexit sha1_base64="ivpwjimwnjua88/lnpsvv6vews=">aaac83ichvgxtunaejy7gjmoalxneqtxqrcfngyahavguzkklwikegr7wzcivxawjurgpuf4acoqgpicot4cwiqju5bz9wogsjhokjy4tueddw7sy+mfd2pwn5gvsb4nvcgxj88heooaypfbr9pjyan1j33aaywdf2has2lnnngktwdhgg2ianmolygpwsndvevtriyueu/bgpbblma3j69w2a4kqqfmkxrpchqbgdtnbsuwk2p+2yex6u1qwffowrf6r6cuk7wnumoqg1mmlpymoiga7gtuanfuieaugbdexxgicggwiajpn1lmadbi2wlqsiurtzkm+iattwmvtgqmandob1bp3kmsjr3np2ibdmtgpyizhqm8a+e4c1e4in+xydxtcjio/ewnnmrz2vedwzva+h+xzzdpodtz9ybjkdardk6rdf/l6r3om+am7plpu9cehvmeobp554edlrjt3xao3u3xaw8lphnb7hdfxlek/xgjojw3f2zxzldgnwrr/j/flsd6xntbr5oyyyyvxijmwcd9ghua4amvwhdagspcewdn8giutqxw1i+24x6olys4x+mesclj7lx</latexit> Fig. 1: Sampled-data system. Iwazumi [38] have developed a state-space-based design method of finite-dimensional controllers for output regulation of unstable continuous-time infinite-dimensional systems, by using residue mode filters proposed in [34]. On the other hand, we employ a frequency-domain technique based on coprime factorizations as in [9, 15, 16, 21, 22]. In particular, we extend a design of stabilizing controller in [18] to output regulation. Let A, B, C) and G be generating operators and a transfer function of a wellposed system Σ, respectively. The operator A is the generator of a strongly continuous semigroup T, which governs the dynamics of the system without control. In addition, B and C are control and observation operators, respectively. We consider only infinitedimensional systems that has finite-dimensional input and output spaces with the same dimension, namely, the transfer function G is a square-matrix-valued function. The well-posed system is connected with a discrete-time linear controller Σd through a zero-order hold Hτ and a generalized sampler Sτ, where τ > is a sampling period. Let u, y be the input and output of the well-posed system Σ and ud, yd be the input and output of the digital controller Σd, respectively. The generalized sampler Sτ is written as Z τ Sτ y)k) = wt)ykτ + t)dt k Z+, where w L2, τ ) is the scalar weighting function. Using the zero-order hold Hτ and the generalized sampler Sτ, we construct the sampled-data feedback in the following form: u = Hτ yd + v1r+ ud = yref 1Z+ Sτ y, where v1r+ and yref 1Z+ are constant external disturbance and reference signals, respectively. Fig. 1 illustrates the sampled-data system we study. Since the output y of well-posed systems belongs to L2loc, the output y is not guaranteed to converge to yref as t. In this paper, we therefore consider the convergence of the output in the following sense: There exist constants Γ > and α < such that! x) ky yref kl2α Γref + kvkcp + kyref kcp xd ) X X d for all t, all initial states x) X and xd ) Xd of Σ and Σd, and all v, yref Cp, where L2α is the L2 -space weighted by the exponential function e α. The above condition means that as t, the energy of the restricted tracking error y yref ) [t, ) exponentially converges to zero. If we add a precompensator between the plant and the zero-order hold as proposed in [18], then the output y converges to yref under a certain assumption on initial states.

3 SAMPLED-DATA OUTPUT REGULATION 3 Before studying sampled-data output regulation, we investigate discrete-time output regulation for infinite-dimensional systems. In the discrete-time setup, we propose a design method of finite-dimensional controllers that achieve output regulation. Although in the sampled-data setup, we consider constant signals for disturbance and reference signals, the proposed method in the discrete-time setup allows disturbance and reference signals that are finite superpositions of sinusoids. The controller design consists of two steps: First we design stabilizing controllers with a free parameter in H based on the results in [18]. Then we choose the free parameter so that the controller incorporates the internal model for output regulation. The design problem of the desired controller is reduced to a Nevanlinna-Pick interpolation problem with both interior and boundary conditions. Whereas the interior conditions are used for stabilization, output regulation leads to the boundary conditions. Our main result, Theorem 3.8, states that for a given unstable well-posed system, there exists a finite-dimensional digital controller that achieves output regulation for constant disturbance and reference signals if the following conditions are satisfied: i) The resolvent set of A contains. ii) det G). iii) The unstable part of the spectrum of A consists of finitely many eigenvalues with finite multiplicities. iv) The semigroup generated by the stable part of A is exponentially stable. v) The unstable part of A, B, C) is controllable and observable. vi) For every nonzero integer l, 2lπi/τ does not belong to the spectrum of A. vii) For every unstable eigenvalues λ of A, τ wt)eλt dt. viii) For every unstable eigenvalues λ, µ of A and nonzero integer l, τλ µ) 2lπi. ix) The multiplicities of all unstable eigenvalues of A are one. Assumptions iii) viii) are used for sampled-data stabilization in [18]. In fact, iii) vii) are sufficient for the existing of sampled-data stabilizing controller, and further, iii) viii) are necessary and sufficient in the single-input and single-output case. We place the remaining assumptions i), ii), and ix) for output regulation but can remove ix) in the single-input and single-output case. The paper is structured as follows. In section 2, we study an output regulation problem for discrete-time infinite-dimensional systems. In section 3, we obtain a sufficient condition for the existence of finite-dimensional sampled-data controllers achieving closed-loop stability, output tracking, and disturbance rejection for unstable well-posed systems with constant disturbance and reference signals. In section 4, we illustrate our results by applying to a system of retarded type with output delays. Notation and terminology. We denote by Z + and R + the set of nonnegative integers and the set of nonnegative real numbers, respectively. For α R, we define C α := {s C : Res > α}, and for η >, E η := {z C : z > η}. We also define D := {z C : z < 1} and T := {z C : z = 1}. For a set Ω C, its closure is denoted by clω). For an arbitrary set Ω, we define 1 Ω : Ω R by 1 Ω t) := 1 for every t Ω. For a matrix M, let us denote by M and M adj the conjugate transpose and the adjugate matrix of M, respectively. We also denote by M the matrix with complex conjugate entries. Let X and Y be Banach spaces. Let LX, Y ) denote the space of all bounded linear operators from X to Y. We set LX) := LX, X). An operator T LX) is said to be power stable if there exist Γ 1 and ρ, 1) such that T k LX) Γρ k for every k Z +. Let A be a linear operator from X to Y. The domain of A is denoted by doma). The spectrum and resolvent set of A are denoted by σa) and

4 4 MASASHI WAKAIKI AND HIDEKI SANO ϱa), respectively. For α R and a Banach space X, we define the weighted L 2 -space L 2 αr +, X) by L 2 αr +, X) := { f L 2 loc R +, X) : e α f ) L 2 R +, X) } and endow this space with the norm f L 2 α := e α f ) L 2. The space of all functions from Z + to X is denoted by F Z +, X), and f k) := fk + 1) for every k Z + and every f F Z +, X). Let Ω = C α, Ω = E η, or Ω = D. Let H Ω, C p m ) denote the space of all bounded holomorphic functions from Ω to C p m. The norm of H Ω, C p m ) is given by Φ := sup s Ω Φs). We write H Ω) for H Ω, C). 2. Discrete-time output regulation Assumptions and control objective. In this section, we consider the following discrete-time infinite-dimensional system: 2.1a) 2.1b) x k) = Axk) + Buk), yk) = Cxk) + Duk), x) = x X where the state space X is a separable complex Hilbert space, A LX), B LC p, X), C LX, C p ), and D LC p, C p ). We use a strictly causal controller 2.2a) 2.2b) x d k) = P x dk) + Qu d k), y d k) = Rx d k), x d ) = x d X d where the state space X d is a complex Hilbert space, P LX d ), Q LC p, X d ), and R LX d, C p ). The control objective is that the output y tracks a given reference signal y ref despite the external disturbance signal v, which are given by 2.3a) 2.3b) 2.3c) ξ k) = Sξk), vk) = Eξk) y ref k) = F ξk), ξ) = ξ C n where E LC p, C n ), F LC p, C n ), and S := diag e iθ1,..., e iθn) with θ 1,..., θ n [, 2π) distinct. The inputs u and u d of the plant and the controller are given by uk) = y d k) + vk), u d k) = y ref k) yk) =: ek). We can write the dynamics of the closed-loop system as 2.4a) x e k) = A e x e k) + B e ξk), x e ) = x e 2.4b) ek) = C e x e k) + D e ξk), [ ] [ ] xk) x where x e k) =, x x d k) e = x, and d 2.5) [ ] [ ] A BR BE A e :=, B QC P QDR e := QF DE) C e := [ C DR ], D e := F DE. Given a controller in 2.2) represented by the operators P, Q, R), we consider a set of perturbed plants and exosystems OP, Q, R) defined as follows.

5 SAMPLED-DATA OUTPUT REGULATION 5 Definition 2.1 Set of perturbed plants and exosystems). The operators Ã, B, C, D, Ẽ, F ) satisfy Ã, B, C, D, Ẽ, F ) OP, Q, R) if these operators satisfy the two conditions: 1. Ã LX), B LC p, X), C LX, C p ), D LC p, C p ), Ẽ LC p, C n ), and F LC p, C n ). 2. The perturbed operator Ãe defined by [ ] Ã BR Ã e := Q C P Q DR is power stable. If A e is power stable, the conditions above are satisfied for any bounded perturbations of sufficiently small norms. In this section, we consider a robust output regulation problem. Problem 2.2 Robust output regulation for discrete-time systems). Given the plant 2.1) and the exosystem 2.3), Find a controller 2.2) satisfying the following properties: Stability: The operator A e is power stable. Tracking: There exist M e > and ρ e > such that for every initial state x X, x d X d, and ξ C n, [ ) ek) C p M e ρ k x e xd] + ξ C n k Z +. X X d Robustness: If the operators A, B, C, D, E, F ) are changed to Ã, B, C, D, Ẽ, F ) OP, Q, R), then the above tracking condition still holds. Before proceeding to the design of finite-dimensional controllers for the robust output regulation, we recall the internal model principle. Definition 2.3 Definition 6.1 in [27]). A controller 2.2) is said to incorporate a p-copy internal model of the exosystem 2.3) if 2.6) dim kere iθ l I P ) p l {1,..., n}. Theorem 2.4 Theorem IV.5 in [26]). Suppose that A e is power stable. The controller 2.2) incorporates a p-copy internal model of the exosystem 2.3) if and only if it is a solution of Problem Design of finite-dimensional controllers. Throughout this section, we impose the following assumptions: a1 For every l {1,..., n}, e iθ l ϱa). a2 det Ge iθ l ) for every l {1,..., n}. a3 There exist subspaces X + and X with dim X + < such that X = X + X. a4 AX + X + and AX X. Let us denote the projection operator from X to X + by Π, and define A + := A X +, A := A X, B + := ΠB, B := I Π)B, C + := C X +, C := C X. We place the remaining assumptions.

6 6 MASASHI WAKAIKI AND HIDEKI SANO a5 σa) cle 1 ) consists of finitely many eigenvalues with finite algebraic multiplicities, σa + ) = σa) cle 1 ), and there exists η, 1) such that σa ) = σa) C \ cle η ) ). a6 A +, B +, C + ) is controllable and observable. a7 The zeros of detzi A + ) are simple. We place the conditions a1 and a2 for robust output regulation. The conditions a3 a6 are used for stabilization of discrete-time infinite dimensional systems; see, e.g., in [17]. We need a7 for a technical reason but can remove it in the case p = 1. By a4, we obtain where Gz) := CzI A) 1 B + D = G + z) + G z), G + z) := C + zi A + ) 1 B +, G z) := C zi A ) 1 B + D. By a5, G H E η, C p p ) with η, 1). By a6, the unstable part G + of the plant has no unstable pole-zero cancellations. There exist N +, D + with rational entries in H E 1, C p p ) such that G + = D 1 + N + and N +, D + are left coprime over the sets of rational functions in H E 1, C p p ). Choose such N +, D + arbitrarily, and let χ 1,..., χ Υ be the zeros of det D + in cle 1 ). Together with a6 and a7, Lemma A.7.39 of [5] shows that these zeros are equal to the eigenvalues of A + and that their multiplicities are one. The objective of this section is to prove the following theorem constructively: Theorem 2.5. Assume that a1 a7 hold. There exists a finite-dimensional controller 2.2) that is a solution of the robust output regulation problem, Problem Preliminary lemmas. We first show three preliminary results. The first lemma provides an upper bound of inverse matrices. Lemma 2.6. Let V, W C p p. If V is invertible and if then W is also invertible and V 1 C p p V W C p p < 1, 2.7) W 1 V 1 C p p C p p 1 V 1. C p p V W C p p Proof. Since I V 1 W C p p V 1 C p p V W C p p < 1, it follows that V 1 W and hence W are invertible. Using the identity we obtain V 1 W 1 = V 1 W V )W 1, V 1 W 1 C p p V 1 C p p V W C p p W 1 C p p.

7 SAMPLED-DATA OUTPUT REGULATION 7 This yields W 1 C p p V 1 C p p + V 1 W 1 C p p V 1 C p p + V 1 C p p V W C p p W 1 C p p. Thus, we obtain the desired inequality 2.7). The second preliminary result characterizes adjugate matrices. Lemma 2.7. For a region Ω C, consider a holomorphic function : Ω C p p. Suppose that the zero z of det is simple. Then dim ker z ) = 1. Furthermore, if a nonzero vector w C p satisfy ker z ) = {αw : α C}, then there exist α 1,..., α p C such that α l for some l {1,..., p} and adj z ) can be written as 2.8) adj z ) = α 1 w. α p w Proof. Suppose, to get a contradiction, that dim ker z ) 2. There exists nonzero vectors w 1, w 2 C p such that w 1, w 2 are linearly independent and z )w 1 =, z )w 2 =. Let e 1,..., e n be the standard basis of the p-dimensional Euclidean space. There exists an invertible matrix U C p p such that w 1 = Ue 1 and w 2 = Ue 2. Let us denote by l the lth column vector of the product U. Then. l z ) = z )Ue l = z )w l = l {1, 2}. Since each element of U is holomorphic, there exist vector-valued functions 1 and 2 with each entry holomorphic such that 1 z) = z z ) 1 z) and 2 z) = z z ) 2 z). Thus, det z) = det z)u) det U 1 = z z ) 2 det [ 1 z) 2 z) 3 z) p s) ] det U 1, which contradicts that z is a simple zero. To prove the second result, we employ Cramer s rule 2.9) adj = adj = det I. Immediately, we obtain adj z ) z ) = det z )I =. Since dim ker z ) = 1, it follows that all the row vector of adj z ) can be written as αw for some α C, where ker z ) = {αw C p : α C}. Thus 2.8) holds. Finally, let us show the existence of nonzero coefficients α l. By contradiction, assume that α l = in 2.8) for every l {1,..., n}, namely, that adj z ) =. Since adj and det are holomorphic, then there exist holomorphic functions F and f such that adj s) = s z )F and det s) = s z )f. Since z is a simple zero of det, it follows that fz ). Substituting these equations to Cramer s rule 2.9), we obtain F = fi. It follows that which contradicts fz ) and w. = w z )F z ) = fz )w,

8 8 MASASHI WAKAIKI AND HIDEKI SANO The third preliminary lemma provides a stabilizable and detectable realization of the product of two finite-dimensional systems. Lemma 2.8. For l 1, 2, consider the matrix pair P l, Q l, R l, S l ) and define the transfer function K l z) := R l zi P l ) 1 Q l + S l. Assume that σp 1 ) σp 2 ) cle 1 ) =. Assume also that K 1 λ) is full column rank for every λ σp 2 ) cle 1 ) and that K 2 λ) is full row rank for every λ σp 1 ) cle 1 ). If P l, Q l, R l ) is stabilizable and detectable for l {1, 2}, then the realization of K 1 K 2 given by ) [ ] R1 S 1 R 2, S1 S 2 2.1) [ ] P1 Q 1 R 2, P 2 is stabilizable and detectable. [ ] Q1 S 2, Q 2 Proof. It is well known that 2.1) is a realization of K 1 K 2 ; see, e.g., Section 3.6 of [43]. It suffices to show that the realization 2.1) is stabilizable and detectable. Assume, to reach a contradiction, that the realization 2.1) is not stabilizable. Then there exist a complex number λ σp 1 ) σp 2 ) with λ 1 and vectors w 1, w 2 such that [ ] w1, w 2 [ w 1 w2] [ ] λi P 1 Q 1 R 2 =, λi P 2 [ w 1 w2] [ ] Q 1 S 2 =. Q 2 For the case λ σp 2 ), since λ ϱp 1 ), we obtain w 1 = from w 1λI P 1 ) =. Therefore, w 2λI P 2 ) =, w 2Q 2 =. Using the stabilizability of P 2, Q 2 ), we find w 2 =. This is a contradiction. Suppose next that λ σp 1 ). Then λi P 2 is invertible. Therefore, 2.11) w 2 = w 1Q 1 R 2 λi P 2 ) 1. We obtain w 1Q 1 K 2 λ) = w 1Q 1 R 2 λi P 2 ) 1 Q 2 + S 2 ) =. Since K 2 λ) is full row rank, it follows that w 1Q 1 =. Together with w 1λI P 1 ) =, this implies that w 1 = by the stabilizability of P 1, Q 1 ). Hence w 2 = by 2.11). This is a contradiction. Thus, the realization 2.1) is stabilizable. The detectability of the realization 2.1) can be obtained in a similar way Proof of Theorem 2.5. Let us start to prove Theorem 2.5, by using Lemmas We obtain an upper bound on the norm of a certain matrix-valued function at z = e iθ l for every l {1,..., n}. Lemma 2.9. Assume that a1 a7 hold. Define { 2.12) δ := max D+ e iθ l )Ge iθ l ) ) 1 Cp p For every R H E 1, C p p ) satisfying 2.13) G 1 R H E 1) < 2δ, D + H E 1) we obtain 2.14) } : l {1,..., n}. N + e iθ l ) + D + e iθ l )Re iθ l ) ) 1 C p p < 2δ l {1,..., n}.

9 SAMPLED-DATA OUTPUT REGULATION 9 Proof. The assumption a1 leads to det De iθ l ) for every l {1,..., n}, which together with a2 implies that δ is well defined. Since G = G + + G = D 1 + N + + G, it follows that N + e iθ l ) + D + e iθ l )G e iθ l ) ) 1 C = D + e iθ l )Ge iθ l ) ) 1 C δ p p p p for every l {1,..., n}. Moreover, for every l {1,..., n}, D + e iθ l ) C p p D + H E 1), G e iθ l ) Re iθ l ) C p p G R H E 1). Thus we conclude from Lemma 2.6 that for every l {1,..., n}, the matrix N + e iθ l )+ D + e iθ l )Re iθ l ) is invertible and satisfies N + e iθ l ) + D + e iθ l )Re iθ l ) ) 1 C p p N + e iθ l ) + D + e iθ l )G e iθ l ) ) 1 C p p 1 ) N + e iθ l ) + D+ e iθ l )G e iθ 1 C l ) D + p p H E 1) G R H E 1) = 2δ. For the rational functions N +, D +, which are left coprime over the sets of rational functions in H E 1, C p p ), there exists a strictly proper rational function Y + H E 1, C p p ) and a rational function Z + H E 1, C p p ) such that the Bezout identity 2.15) N + Y + + D + Z + = I holds; see, e.g., Lemma of [39]. We provide interpolation conditions on Y + H E 1, C p p ), as in Theorem IV.3 of [4]. To that purpose, we see from Lemma 2.7 and a7 that, for every r {1,..., Υ}, there uniquely exists a nonzero vector w r C p such that w rd + χ r ) =. Lemma 2.1. Suppose that a1 a7 are satisfied. A rational function Y + H E 1, C p p ) is strictly proper and satisfies the Bezout identity 2.15) for some rational function Z + H E 1, C p p ) if and only if the interpolation conditions 2.16) Y + ) =, w rn + χ r )Y + χ r ) = w r r {1,..., Υ}. hold. Moreover, if the latter part of the interpolation conditions 2.16) holds, then a rational function 2.17) Z + := Dadj + D adj + N + Y + det D + satisfies Z + H E 1, C p p ) and the Bezout idendity 2.15).

10 1 MASASHI WAKAIKI AND HIDEKI SANO Proof. It is clear that the strictly properness of Y + is equivalent to Y + ) =. Suppose that a rational function Y + satisfies the Bezout identity 2.15) for some rational function Z + H E 1, C p p ). Using Cramer s rule for D +, we obtain 2.18) D adj + = D adj + N + Y + + D + Z + ) = D adj + N + Y + + det D + Z +. Since χ 1,..., χ r are the zeros of det D +, it follows that D adj + χ r ) = D adj + N + Y + )χ r ) r {1,..., Υ}. The second statement of Lemma 2.7 shows that wrn + χ r )Y + χ r ) = wr for every r {1,..., Υ}. Conversely, suppose that the interpolation condition 2.16) holds. To show that Y + H E 1, C p p ) satisfies the Bezout identity 2.15) for some rational function Z + H E 1, C p p ), it suffices to prove that Z + defined by 2.17) satisfies the Bezout identity 2.15) and Z + H E 1, C p p ). Using Cramer s rule for D +, we find that Z + satisfies the Bezout identity 2.15). By way of contradiction, assume that Z + H E 1, C p p ). Let the j, l)th entry Z j,l + of Z + satisfy Z j,l + H E 1 ). We obtain det D + Z + = D adj + I N + Y + ) H E 1, C p p ). Since a rational function det D + is not strictly proper by Theorem of [39], it follows that Z j,l + is proper. Therefore, there exists a pole of the rational function Z j,l + in cle 1 ) that is equal to a zero χ r of det D +. Since χ r is a simple zero, it follows that 2.19) det D+ Z j,l ) + χr ). However, by the latter part of the interpolation conditions 2.16) and Lemma 2.7, we obtain det D + Z + )χ r ) = D adj + χ r )I N + Y + )χ r ) ) =. This contradicts 2.19). Set M > as in 2.2) M > inf { Y H E 1) : Y + H E 1, C p p ) is rational and satisfies 2.16) }. Since there always exists a rational function Y + H E 1, C p p ) satisfying the interpolation conditions 2.16), the right side of 2.2) belongs to R +. Next, we prove the existence of rational functions satisfying a norm condition and boundary interpolation conditions in addition to the Bezout identity. Lemma Assume that a1 a7 hold, and define δ > by 2.12). For every R H E 1, C p p ) satisfying 2.14), there exist a strictly proper rational function Y + H E 1, C p p ) and a rational function Z + H E 1, C p p ) such that Y + satisfies the interpolation conditions 2.21a) 2.21b) Y + e iθ l ) = N + + D + R) 1 e iθ l ) l = {1,..., n} Y +e iθ l ) = N + + D + R) 1 e iθ l ) D + e iθ l ) + N + + D + R) e iθ l )N + + D + R) 1 e iθ l ) ) l = {1,..., n},

11 SAMPLED-DATA OUTPUT REGULATION 11 the norm condition 2.22) Y + H E 1) < max {2δ, M}, and the Bezout identity 2.15) hold. Proof. Lemma 2.1 shows that a rational function Y + H E 1, C p p ) is strictly proper and satisfies the Bezout identity 2.15) for some rational function Z + H E 1, C p p ) if and only if the interpolation conditions 2.16) hold. Hence the problem of finding the desired Y +, Z + H E 1, C p p ) is equivalent to that of finding a rational function Y + H E 1, C p p ) satisfying the interior interpolation conditions 2.16), the boundary interpolation conditions 2.21), and the norm condition 2.22), which is called the Nevanlinna-Pick interpolation problem with both interior and boundary conditions; see Appendix A for details. This interpolation problem is solvable if R satisfies 2.14). Once we obtain a solution Y + H E 1, C p p ) of the interpolation problem, Z + H E 1, C p p ) defined by 2.17) satisfies the Bezout identity 2.15). Lemma Assume that a1 a7 hold. Suppose that a rational function R H E 1, C p p ) satisfies 2.14) and that Y +, Z + H E 1, C p p ) satisfy the conditions in Lemma Then the following results holds: a) Y + and Z + RY + are right coprime over the set of rational functions in H E 1 ). b) The rational function defined by 2.23) K := Y + Z + RY + ) 1 is strictly proper and satisfies 2.24) K = Y + I N+ + D + R)Y + ) 1D+. c) There exists a rational function Ẑ+ H E 1, C p p ) such that 2.25a) 2.25b) n Z + RY + )z) = z e iθ l )Ẑ+z) l=1 det Ẑ+e iθ l ). Proof. a) By the Bezout identity 2.15), 2.26) N + + D + R)Y + + D + Z + RY + ) = I. Hence Y + and Z + RY + are right coprime over the sets of rational functions in H E 1 ). b) Since Y + ) =, it follows from the Bezout identity 2.15) that Z + ) is invertible. Therefore, K ) = and K is strictly proper. The Bezout identity 2.15) also leads to 2.27) Z + RY + = D 1 + I NY + ) RY + = D 1 + I N + D + R)Y + ). Therefore, we obtain 2.24). c) To show the existence of a rational function Ẑ+ H E 1, C p p ) satisfying 2.25), it suffices to prove 2.28) Z + RY + )e iθ l ) = l {1,..., n}

12 12 MASASHI WAKAIKI AND HIDEKI SANO and Z + RY + ) e iθ l ) is invertible for all l {1,..., n}. We immediately obtain 2.28) from 2.21a) and 2.27). We also have Z + RY + ) e iθ l ) = D 1 ) + e iθ l ) ) e iθ I N + D + R)Y l + ) l {1,..., n}. The interpolation condition 2.21b) yields I N+ + D + R)Y + ) e iθ l ) = D + e iθ l ) l {1,..., n}. Thus, Z + RY + ) e iθ l ) = I for all l {1,..., n}. This completes the proof. For δ in 2.12) and M in 2.2), define 2.29) M 1 := max {2δ, M}. The following lemma provides a sufficient condition for the robust output regulation to be solvable. Lemma Assume that a1 a7 hold. If we choose a rational function R H E 1, C p p ) so that 2.3) G 1 R H E 1) <, M 1 D + H E 1) then there exists a realization P, Q, R) of the rational function K defined by 2.23) such that the controller 2.2) with this P, Q, R) is a solution of Problem 2.2. Proof. Since 1 1 < M 1 D + H E 1) 2δ, D + H E 1) Lemma 2.9 shows that a rational function R satisfies 2.14). Due to Theorem 2.4, it is enough to prove that A e defined by 2.5) is power stable and that the controller 2.2) incorporates a p-copy internal model. To this end, we first decompose K. Fix a 1, 1). For Y +, Ẑ+ H E 1, C p p ) in Lemma 2.12, define 2.31) K 1 z) := n l=1 z a z e iθ I, K 2 z) := Y + z) z α) n Ẑ + z) ) 1. l Then K = K 1 K 2. For every l {1,..., n}, let c l C be the residue of n z a j=1 at z = e iθ l. z e iθ j Using the identity matrix I with dimension p, we define I 2.32) P 1 := diag e iθ1 I,..., e iθn I ), Q 1 :=, R 1 := [ c 1 I c n I ], S 1 := I. Then P 1, Q 1, R 1, S 1 ) is a minimal realization of K 1, and dim kere iθ l I P 1 ) p for every l {1,..., n}. Let P 2, Q 2, R 2, S 2 ) be a minimal realization of K 2. Since K 2 is strictly proper, it follows that S 2 =. In addition, the realizations P 1, Q 1, R 1, S 1 ) and P 2, Q 2, R 2, S 2 ) satisfy the conditions in Lemma 2.8. In fact, since 2.26) leads to N + + D + R)Y + + D + n l=1. I z e iθ l z a ) z a) n Ẑ + = I,

13 SAMPLED-DATA OUTPUT REGULATION 13 it follows that Y + and z a) n Ẑ + are right coprime. Lemma A.7.39 of [5] shows that every l {1,..., n} satisfies e iθl σp 2 ) by 2.25b). We obtain det K 1 λ) for every λ σp 2 ) cle 1 )). By 2.21a), det K 2 λ) for every λ σp 1 ) cle 1 ) = {e iθ1,..., e iθn }. Therefore, Lemma 2.8 shows that the realization P, Q, R) of K = K 1 K 2 in the form 2.1) is stabilizable and detectable. By construction, the controller 2.2) with P, Q, R) incorporates a p-copy internal model. We can see the power stability of A e from the same argument as in the proof of Theorem 9 in [18]. Since D + G R)Y + H E 1) D + H E 1) G R H E 1) Y + H E 1) 1 < D + H E 1) M 1 < 1. M 1 D + H E 1) Therefore U := D + G + N + )Y + + D + Z + RY + ) satisfies U I H E 1) = D + G R)Y + H E 1) < 1, which yields U, U 1 H E 1, C p p ). Since I + GK) 1 = Z + RY + )U 1 D +, G = D 1 + N + + D + G ), it follows that [ ] 1 [ ] I K I KI + GK) = 1 G KI + GK) 1 G I I + GK) 1 G I + GK) 1 [ ] I Y = + U 1 N + + D + G ) Y + U 1 D + Z + RY + )U 1 N + + D + G ) Z + RY + )U 1 D + H E 1, C 2p 2p ). Since [ I ] 1 K G I is a transfer function of the system [ ] [ ] [ ] ) 1 B R I A e,,,, Q C D I Theorem 2 of [17] shows that A e is power stable if 2.33) [ ]) B A e,, Q [ ] ) R, A C e is stabilizable and detectable, respectively, which is equivalent to the stabilizablity and detectablility of A, B, C) and P, Q, R). These properties of A, B, C) follow from a6, and we have already proved that P, Q, R) is stabilizable and detectable. This completes the proof. We are now in a position to prove Theorem 2.5. Proof of Theorem 2.5. Due to Lemma 2.13, it remains to show that the existence of a rational function R H E 1, C p p ) satisfying 2.3).

14 14 MASASHI WAKAIKI AND HIDEKI SANO Since G H E η, C p p ) with η, 1), the Taylor expansion of G at, G z) = G j z j, j= converges uniformly in E 1, i.e., lim sup N N z E 1 G z) G j z j Thus 2.3) holds with for a sufficiently large N N. R z) := j= N G j z j j= C p p We summarize the proposed method for the construction of finite-dimensional controllers. The problem of finding a rational function in the steps 2 and 5 of the procedure below is the Nevanlinna-Pick interpolation problem; see Appendix A for details. Design procedure of controllers 1. Obtain a left-coprime factorization D 1 + N + of a rational function G + over the set of rational functions in H E 1 ). 2. Find M > satisfying 2.2). 3. Set M 1 > as in 2.29). 4. Find a rational function R H E 1, C p p ) satisfying the norm condition 2.3). 5. Find a rational function Y + H E 1, C p p ) satisfying the interpolation conditions 2.16), 2.21) and the norm condition Y + H E 1) < M Define a rational function Z + H E 1, C p p ) by 2.17). 7. Find a rational function Ẑ+ H E 1, C p p ) satisfying 2.25). 8. Define the minimal realization P 1, Q 1, R 1, S 1 ) as in 2.32) and compute a minimal realization of K 2 defined by 2.31). 9. Calculate a realization 2.1), which is a realization of the desired controller. In the case p = 1, we can remove the assumption a7 and the redundant steps 6 8 in the above design procedure. To see this, let the multiplicity of the zeros χ 1,..., χ Υ in cle 1 ) of detsi A + ) be J r N for every r {1,..., Υ}. Then there exists a rational function Y + H E 1 ) satisfying the interpolation conditions = a) 2.16 b) 1 Y + ) =, Y + χ r ) = N + χ r ) Y j) + χ r ) = 1 j 1 ) j N j l) χ r )Y l) + χ r ) N + χ r ) l l= for all r {1,..., Υ}, j {1,..., J r } and 2.21 ) Y + e iθ l ) = 1 N + e iθ l ) + D+ e iθ l )Re iθ l)

15 SAMPLED-DATA OUTPUT REGULATION 15 for all l = {1,..., n} and the norm condition Y + H E 1) < M 1. Then Y + D + K := 1 N + + D + R)Y + is strictly proper, has a pole at z = e iθ l for all l {1,..., n}, and satisfies [ ] 1 1 K H E G 1 1, C 2 2 ). Lemma A.7.39 of [5] directly shows that a minimal realization P, Q, R) of K incorporates a p-copy internal model, without the decomposition K = K 1 K 2 as in Lemma Thus, P, Q, R) is the desired realization of K. Since a similar argument leads to this conclusion, we omit the details for the sake of brevity. 3. Sampled-data servo control for constant reference signals and disturbances Preliminaries on well-posed systems. We provide brief preliminaries on well-posed linear systems and refer the readers to the surveys [37,41] and the book [36] for more details. As a plant, we consider a well-posed system Σ with state space X, input space C p, and output space C p, generating operators A, B, C), transfer function G, and the input-output operator G. Here X is a separable complex Hilbert space with norm and A is the generator of a strongly continuous semigroup T = T t ) t on X. The spaces X 1 and X 1 are interpolation and extrapolation spaces associated with X, respectively. For λ ϱa), the space X 1 is defined as doma) endowed with the norm ζ 1 := λi A)ζ, and X 1 is the completion of X with respect to the norm ζ 1 := λi A) 1 ζ. Different choices of λ lead to equivalent norms on X 1 and X 1. The semigroup T restricts to a strongly continuous semigroup on X 1, and the generator of the restricted semigroup is a restriction of A. Similarly, T can be uniquely extended to a strongly continuous semigroup on X 1, and the generator of the extended semigroup is an extension of A. We denote the restrictions and extensions of T and A by the same symbols. The restriction and extension of T has the same exponential growth bound as the original semigroup T, and we denote the exponential growth bound by ωt), i.e., ln T t ωt) := lim. t t We place the following conditions for the operators A, B, C) and the transfer function G to be well posed: The control operator B satisfies B LC p, X 1 ) and is an admissible control operator for T, that is, for every t, there exists b t such that t T t s Bus)ds b t u L 2,t) u L 2 [, t], C p ). The observation operator C satisfies C LX 1, C p ) and is an admissible observation operator for T, that is, for every t, there exists c t such that t ) 1/2 CT s ζ 2 C pds c t ζ ζ X 1.

16 16 MASASHI WAKAIKI AND HIDEKI SANO The transfer function G : C ωt) C p p satisfies 1 ) Gs) Gλ) = CsI A) 1 λi A) 1 B s, λ C ωt) with s λ s λ and G H C α, C p p ) for every α > ωt). The transfer function G may have an analytic extension to a half plane C α with α < ωt). If it exists, we say that G is holomorphic meromorphic) on C α and use the same symbol G for an analytic extension to a larger right half plane. For every α > ωt), the input-output operator G : L 2 loc R +, C p ) L 2 loc R +, C p ) satisfies G L L 2 αr +, C p ), L 2 αr +, C p ) ) and LGu)s) ) = Gs) Lu) ) s) s Cα, u L 2 αr +, C p ), where L denotes the Laplace transform. The Λ-extension C Λ of C is defined by C Λ ζ := lim CssI s,s R A) 1 ζ with its domain domc Λ ) consisting of those ζ X for which the limit exists. For every ζ X, we obtain T t ζ domc Λ ) for a.e. t. By the admissibility of C, for ) every t, there exists c Λ t 1/2 t such that C ΛT s ζ 2 C pds c Λ t ζ for every ζ X. If we define the operator Ψ : X L 2 loc R +, C p ) by Ψζ)t) := C Λ T t ζ ζ X, a.e. t, then Ψ satisfies Ψ L X, L 2 αr +, C p ) ) for every α > ωt ). Fix λ C ωt) arbitrarily. Let x and y denote the state and output functions of Σ corresponding to the initial condition x) = x X and the input function u L 2 loc R +, C p ). Then we obtain 3.1) xt) = T t x + t T t s Bus)ds t, xt) λi A) 1 But) domc Λ ) for a.e. t, and 3.2a) 3.2b) ẋt) = Axt) + But), x) = x X a.e. t yt) = C Λ xt) λi A) 1 But) ) + Gλ)ut) a.e. t, where the differential equation 3.2a) is interpreted on X 1. We have from 3.1) and 3.2b) that for every u L 2 loc R +, C p ) and a.e. t, the input-output operator G satisfies t ) 3.3) Gu)t) = C Λ T t s Bus)ds λi A) 1 But) + Gλ)ut) Closed-loop system and control objective. Let τ > denote the sampling period. The zero-order hold operator H τ : F Z +, C p ) L 2 loc R +, C p ) is defined by H τ f)kτ + t) := fk) t [, τ), k Z +.

17 SAMPLED-DATA OUTPUT REGULATION 17 The generalized sampling operator S τ : L 2 loc R +, C p ) F Z +, C p ) is defined by S τ g)k) := τ wt)gkτ + t)dt k Z +, where w L 2, τ) is the scalar weighting function with τ wt)dt = 1. We connect the continuous-time system 3.2) and the discrete-time controller 2.2) via the following sampled-data feedback law: u = H τ y d + v1 R+, u d = y ref 1 Z+ S τ y, where v1 R+ and y ref 1 Z+ with v C p and y ref C p are the external disturbance and the reference signal whose values are constant but unknown. The dynamics of the sampled-data system is given by 3.4a) 3.4b) 3.4c) 3.4d) ẋ = Ax + BH τ y d + v1 R+ ), x) = x X y = C Λ x λi A) 1 BH τ y d + v1 R+ ) + Gλ)Hτ y d + v1 R+ ) x d = P x d + Qy ref 1 Z+ S τ y), x d ) = x d X d y d = Rx d. We define the exponential stability of this sampled-data system. Definition 3.1 Exponential stability). The sampled-data system 3.4) is exponentially stable if there exist Γ 1 and γ > such that [ ] [ ) xkτ + t) Γ e γkτ+t) x 3.5) x d k) x + v C p + y ref C p X X d d] X X d k Z +, t [, τ), x X, x d X d, v, y ref C p. We consider a set of perturbed plants O s P, Q, R) defined as follows. Definition 3.2 Set of perturbed plants). The operators Ã, B, C) and the transfer function G satisfy Ã, B, C, G) O s P, Q, R) if Ã, B, C, G) satisfy the two conditions: 1. The operators Ã, B, C) and the transfer function G generate a well-posed system. 2. The perturbed sampled-data system, in which A, B, C, G) are changed to Ã, B, C, G), is exponentially stable. In this section, we study the following sampled-data robust output regulation problem. Problem 3.3 Robust output regulation for sampled-data systems). Find a controller 2.2) such that the sampled-data system 3.4) satisfies the following three properties: Stability: The sampled-data system is exponentially stable. Tracking: There exist Γ ref > and α < such that [ ) x 3.6) y y ref L 2 α Γ ref xd] + v C p + y ref C p X X d x X, x d X d, v, y ref C p.

18 18 MASASHI WAKAIKI AND HIDEKI SANO Robustness: If the operators A, B, C) and the transfer function G are changed to Ã, B, C, G) O s P, Q, R), then the above tracking condition still holds Assumptions on well-posed systems. In what follows, we impose several assumptions on the well-posed system 3.2). b1 ϱa). b2 det G). b3 There exists ε > such that σa) clc ε ) consists of finitely many isolated eigenvalues of A with finite algebraic multiplicities. Under the assumption b3, we obtain the following spectral decomposition of X for A; see, e.g., Lemma of [5] and Proposition IV.1.16 of [6]. There exists a rectifiable, closed, simple curve Φ in C enclosing an open set that contains σa) clc ) in its interior and σa) C \ clc ) ) in its exterior. The operator Π := 1 si A) 1 ds 2πi Φ is a projection on X. Define X + := ΠX and X := I Π)X. Then X = X + X, dim X + <, and X + X 1. The subspaces X + and X are T t -invariant for all t. Define Then A + := A X +, A := A X1 X, T+ t := T t X +, T t := T t X. σa + ) = σa) clc ) σa + ) = σa) C \ clc ) ), and T + := T + t ) t and T := T t ) t are strongly continuous semigroups on X + and X with generators A + and A, respectively. The projection operator Π on X can be extended to a projection Π 1 on X 1, and Π 1 X 1 = ΠX = X +. We define B + := Π 1 B, B := I Π 1 )B, C + := C X +, C := C X1 X. We are now in a position to formulate the remaining assumptions. b4 The strongly continuous semigroup T = T t ) t b5 A +, B +, C + ) is controllable and observable. b6 2lπi/τ σa + ) for every l Z \ {}. b7 τ wt)eλt dt for every λ σa + ). b8 τλ µ) 2lπi for every λ, µ σa + ) and for every l Z \ {}. b9 The zeros of detsi A + ) are simple. As in the discrete-time case, we assume b1 and b2 for output regulation. For a technical reason, we place the assumption b9 but can remove it in the case p = 1, as commented in the last paragraph of Sec. 2. Proposition 5 and Lemma 9 of [18] show that for the existence of stabilizing controllers, the conditions b3 b7 are sufficient, and the conditions b3 b8 are necessary and sufficient in the case p = 1. Define the input-output operator G + of the system A +, B +, C + ) by G + u)t) := t C + e A+ t s) B + us)ds t, u L 2 locr +, C p ). and define G := G G +. We use the following result on the decomposition of the output:

19 SAMPLED-DATA OUTPUT REGULATION 19 Lemma 3.4 Lemma 4.2 in [2]). Assume that b3 holds. There exists a wellposed system Σ with generating operator A, B, C ) and input-output operator G. For every x X and every u L 2 loc R +, C p ), the output y of the well-posed system 3.2) can be written in the form 3.7) yt) = C + Πxt) + C ) Λ T t I Π)x + G u)t) a.e. t. The Λ-extension of C satisfies 3.8) C ) Λ ξ = C Λ z z dom C ) Λ ) = domcλ ) X Properties of discretized systems. We first recall the discrete-time dynamics of the plant combined with the zero-order hold and the sampler. Define A τ := T τ LX). By the admissibility of B, the operator B τ : L 2 [, τ], C p ) X defined by B τ g := τ T t Bgτ t)dt satisfies B τ LL 2 [, τ], C p ), X). Similarly, by the admissibility of C, the operator C τ : X C p defined by C τ ζ := τ wt)c Λ T t ζdt ζ X satisfies C τ LX, C p ). We define the operator D τ : L 2 [, τ], C p ) C p by D τ g := τ wt)gg)t)dt, which satisfies D τ L L 2 [, τ], C p ), C p). For ψ C p, we set B τ ψ := B τ ψ1,τ) ), Dτ ψ := D τ ψ1,τ) ). Lemma 3.5 Lemma 2 of [18]). Let u = H τ f + g, where f F Z +, C p ) and g L 2 loc R +, C p ), and let x X. Set xt) as in 3.1). Then x k + 1)τ ) = A τ xkτ) + B τ fk) + B τ L kτ g S τ y)k) = C τ xkτ) + D τ fk) + D τ L kτ g, where L kτ g L 2 [, τ], C p ) is defined by L kτ g)t) = gkτ + t) for all t [, τ]. Next we investigate relationships between the continuous-time dynamics and its discrete-time counterpart. Define A + τ := T + τ = A τ X +, A τ := T τ = A τ X. Lemma 3.6. If b1, b3, b4, and b6 hold, then 1 ϱa τ ). Proof. Since X + and X are A τ -invariant, it is enough to show that 1 ϱa + τ ) ϱa τ ). Since b1 holds, it follows that ϱa + ). Together with b6, this yields 2lπi/τ σa + ) for every l Z. By the spectral mapping theorem, 3.9) σ e τa+) = e τσa+). Therefore, 1 σ e τa+ ) = σa + τ ). On the other hand, b4 leads to the power stability of A τ, and hence 1 ϱa τ ). This completes the proof

20 2 MASASHI WAKAIKI AND HIDEKI SANO Define G τ z) := C τ zi A τ ) 1 B τ + D τ. For exponentially stable well-posed systems, G τ 1) = G) has been proved in Proposition 4.3 of [21] and Proposition 3.1 of [13]. Here we study the case where the plant has finitely many unstable poles. Lemma 3.7. If b1, b3, and b4 hold, then G τ 1) = G). Proof. Define G + s) := C + si A + ) 1 B +, G s) := Gs) G + s). Clearly, G + is the transfer function of a finite-dimensional system with generating matrices A +, B +, C + ) and input-output operator G +. Lemma 3.4 shows that G is the transfer function of an exponentially stable well-posed system with generating operators A, B, C ) and input-output operator G. Define D τ + : L 2 [, τ], C p ) C p by D + τ g := τ wt)g + g)t)dt. For ψ C p, we also set D + τ ψ := D + τ ψ1,τ) ) and regard D + τ as a matrix in C p p. We first show that 3.1) G + ) = C + τ A + ) 1 B + + D + τ, where A + is invertible by b1. Since if g ψ C p, then it follows from τ G + g)t) = C + e A+t I)A + ) 1 B + ψ, wt)dt = 1 that D + τ ψ = C + τ A + ) 1 B + ψ C + A + ) 1 B + ψ ψ C p. Thus, 3.1) holds. Since T is exponentially stable by b4, A is boundedly invertible. Next we shall prove that 3.11) G ) = C τ A ) 1 B + D τ D + τ. By definition, D τ g D + τ g = Similarly to 3.3), we obtain G g)t) = C ) Λ t τ wt)g g)t)dt g L 2 [, τ], C p ). ) T s B gt s)ds + A ) 1 B gt) + G )gt) g L 2 locr +, C p ), a.e. t, Using t T s B ψds = T t A ) 1 B ψ A ) 1 B ψ ψ C p,