• No results found

Jan C. Willems Department of Electrical Engineering - SCD (SISTA) University of Leuven Kasteelpark Arenberg 10 B-3001 Leuven-Heverlee Belgium

N/A
N/A
Protected

Academic year: 2021

Share "Jan C. Willems Department of Electrical Engineering - SCD (SISTA) University of Leuven Kasteelpark Arenberg 10 B-3001 Leuven-Heverlee Belgium"

Copied!
5
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

STATE AND FIRST ORDER REPRESENTATIONS

Jan C. Willems

Department of Electrical Engineering - SCD (SISTA) University of Leuven Kasteelpark Arenberg 10 B-3001 Leuven-Heverlee Belgium Jan.Willems@esat.kuleuven.ac.be

0.1

Abstract

We conjecture that the solution set of a system of linear constant coefficient PDE’s is Markovian if and only if it is the solution set of a system of first order PDE’s. An analogous conjecture regarding state systems is also made. Keywords: Linear differential systems, Markovian systems, state systems, kernel representations.

0.2

Description of the problem

0.2.1 Notation

First, we introduce our notation for the solution sets of linear PDE’s in the n real independent variables x = (x1, . . . , xn). Let D0n denote, as usual, the set of real distributions on Rn, and Lw

n the linear subspaces of (D 0

n)w consisting of the solutions of a system of linear constant coefficient PDE’s in the w real-valued dependent variables w = col(w1, . . . , ww). More precisely, each element B ∈ Lw

n is defined by a polynomial matrix R ∈ R•×w[ξ1, ξ2, . . . , ξn], with w columns, but any number of rows, such that

B= {w ∈ (D0n)w | R( ∂ ∂x1 , ∂ ∂x2 , . . . , ∂ ∂xn )w = 0}. We refer to elements of Lw

n as linear differential n-D systems. The above PDE is called a kernel representation of B ∈ Lw

n. Note that each B ∈ Lwn has many kernel representations. For an in depth study of Lw

n, see [1] and [2].

(2)

2 CHAPTER 0. Next, we introduce a class of special three-way partitions of Rn. Denote by P the following set of partitions of Rn:

[(S−, S0, S+) ∈ P] :⇔ [(S−, S0, S+ are disjoint subsets of Rn) ∧ (S−∪ S0∪ S+= Rn) ∧ (S− and S+ are open, and S0 is closed)]. Finally, we define concatenation of maps on Rn. Let f

−, f+ : Rn → F, and let π = (S−, S0, S+) ∈ P. Define the map f−∧πf+: Rn → F, called the concatenation of (f−, f+) along π, by (f−∧πf+)(x) :=  f−(x) for x ∈ S− f+(x) for x ∈ S0∪ S+ 0.2.2 Markovian systems Define B ∈ Lw n to be Markovian :⇔ [(w−, w+∈ B ∩ C∞(Rn, Rw)) ∧ (π = (S−, S0, S+) ∈ P) ∧ (w−|S0 = w+|S0)] ⇒ [(w−∧π w+ ∈ B].

Think of S− as the ‘past’, S0 as the ‘present’, and S+ as the ‘future’. Markovian means that if two solutions of the PDE agree on the present, then their pasts and futures are compatible, in the sense that the past (and present) of one, concatenated with the (present and) future of the other, is also a solution. In the language of probability: the past and the future are independent given the present.

We come to our first conjecture: B∈ Lw

n is Markovian if and only if

it has a kernel representation that is first order.

I.e., it is conjectured that a Markovian system admits a kernel representation of the form R0w + R1 ∂ ∂x1 w + R2 ∂ ∂x2 w + · · · Rn ∂ ∂xn w = 0.

Oberst [2] has proven that there is a one-to-one relation between Lw n and the submodules of Rw

1, ξ2, . . . , ξn], each B ∈ Lwn being identifiable with its set of annihilators NB:= {n ∈ Rw[ξ1, ξ2, . . . , ξn] | n>( ∂ ∂x1 , ∂ ∂x2 , . . . , ∂ ∂xn )B = 0}.

Markovianity is hence conjectured to correspond exactly to those B ∈ Lw n for which the submodule NB has a set of first order generators.

(3)

0.2.3 State systems

In this section we consider systems with two kind of variables: w real-valued manifest variables, w = col(w1, . . . , ww), and z real-valued state variables, z = col(z1, . . . , zz). Their joint behavior is again assumed to be the solu-tion set of a system of linear constant coefficient PDE’s. Thus we consider behaviors in Lw+z

n , whence each element B ∈ Lw+zn is described in terms of two polynomial matrices (R, M ) ∈ R•×(w+z)[ξ1, ξ2, . . . , ξn] by

B= {(w, z) ∈ (D0n)w× (D0n)z | R( ∂ ∂x1 , ∂ ∂x2 , . . . , ∂ ∂xn )w + M ( ∂ ∂x1 , ∂ ∂x2 , . . . , ∂ ∂xn )z = 0}. Define B ∈ Lw+z

n to be a state system with state z :⇔

[((w−, z−), (w+, z+) ∈ B ∩ C∞(Rn, Rw+z)) ∧ (π = (S−, S0, S+) ∈ P) ∧ (z−|S0 = z+|S0)] ⇒ [(w−, z−) ∧π (w+, z+) ∈ B].

Think again of S− as the ‘past’, S0 as the ‘present’, S−+ as the ‘fu-ture’. State means that if the state components of two solutions agree on the present, then their pasts and futures are compatible, in the sense that the past of one solution (involving both the manifest and the state vari-ables), concatenated with the present and future of the other solution, is also a solution. In the language of probability: the present state ‘splits’ the past and the present plus future of the manifest and the state trajectory combined.

We come to our second conjecture: B∈ Lw+z

n is a state system if and only if

it has a kernel representation that is first order in the state variables z and zero-th order in the manifest variables w.

I.e., it is conjectured that a state system admits a kernel representation of the form R0w + M0z + M1 ∂ ∂x1 z + M2 ∂ ∂x2 z + · · · Mn ∂ ∂xn z = 0.

(4)

4 CHAPTER 0.

0.3

Motivation and history of the problem

These open problems aim at understanding state and state construction for n-D systems.

Maxwell’s equations constitute an example of a Markovian system. The diffusion equation and the wave equation are non-examples.

0.4

Available results

It is straightforward to prove the ‘if’-part of both conjectures. The conjec-tures are true for n = 1, i.e. in the ODE case, see [3].

Acknowledgement This research is supported by grants from several fund-ing agencies and sources: Research Council KUL: Concerted Research Ac-tion GOA-Mefisto 666 (Mathematical Engineering); Flemish Government: Fund for Scientific Research Flanders, project G.0256.97 (subspace); Re-search communities ICCoS, ANMMM, IWT (Soft4s, softsensors), Eureka-Impact (MPC-control), Eureka-FLiTE (flutter modeling); Belgian Federal Government: DWTC IUAP V-22 (2002-2006): Dynamical Systems and Control: Computation, Identification & Modelling).

(5)

[1] H.K. Pillai and S. Shankar, A behavioral approach to control of dis-tributed systems, SIAM Journal on Control and Optimization, vol-ume 37, pages 388-408, 1999.

[2] U. Oberst, Multidimensional constant linear systems, Acta Applicandae Mathematicae, volume 20, pages 1-175, 1990.

[3] P. Rapisarda and J.C. Willems, State maps for linear systems, SIAM Journal on Control and Optimization, volume 35, pages 1053-1091, 1997.

Referenties

GERELATEERDE DOCUMENTEN

Figure 4.3: 4x4 Torus Topology Throughput VS Injection rate (flits/cycle) In Figure 4.4 a comprehensive comparison of average packet latency vs injection rate has been obtained

Figure 1: Latent traversals along the first principal component from a trained St-RKM model on 3Dshapes (left) and Dsprites (right) dataset.. It has been argued that

A gossip algorithm is a decentralized algorithm which com- putes a global quantity (or an estimate) by repeated ap- plication of a local computation, following the topology of a

Keywords: Dominant eigenspace, eigenvalues, updating, tracking, kernel Gram matrix, Prinicipal Components, large scale

Department of Electrical Engineering, ESAT- SCD(SISTA), Katholieke Universiteit Leuven, Leuven, Belgium; and *Department of Neurology, University Hospital Gasthuisberg, 3000

As the VDSL reach performance is upstream limited, the results of optimal power allocation can be used to improve reach performance by improving the upstream bit

-DATA4s delivers end-to-end solutions to financial institutions and telecom operators for improved risk analysis and management. of their customer and

CANONICAL CONTROLLER to−be−controlled to−be−controlled variables control DESIRED CONTROLLED BEHAVIOR variables PLANT variables