EN FR
EN FR


Section: Research Program

Background on non-standard analysis

Non-Standard analysis plays a central role in our research on hybrid systems modeling [3] , [6] , [15] , [14] . The following text provides a brief summary of this theory and gives some hints on its usefulness in the context of hybrid systems modeling. This presentation is based on our paper [3] , a chapter of Simon Bliudze's PhD thesis  [21] , and a recent presentation of non-standard analysis, not axiomatic in style, due to the mathematician Lindström  [41] .

Non-standard numbers allowed us to reconsider the semantics of hybrid systems and propose a radical alternative to the super-dense time semantics developed by Edward Lee and his team as part of the Ptolemy II project, where cascades of successive instants can occur in zero time by using +× as a time index. In the non-standard semantics, the time index is defined as a set 𝕋={nn*}, where is an infinitesimal and * is the set of non-standard integers. Remark that 1/ 𝕋 is dense in +, making it “continuous”, and 2/ every t𝕋 has a predecessor in 𝕋 and a successor in 𝕋, making it “discrete”. Although it is not effective from a computability point of view, the non-standard semantics provides a framework that is familiar to the computer scientist and at the same time efficient as a symbolic abstraction. This makes it an excellent candidate for the development of provably correct compilation schemes and type systems for hybrid systems modeling languages.

Non-standard analysis was proposed by Abraham Robinson in the 1960s to allow the explicit manipulation of “infinitesimals” in analysis  [48] , [35] , [10] . Robinson's approach is axiomatic; he proposes adding three new axioms to the basic Zermelo-Fraenkel (ZFC) framework. There has been much debate in the mathematical community as to whether it is worth considering non-standard analysis instead of staying with the traditional one. We do not enter this debate. The important thing for us is that non-standard analysis allows the use of the non-standard discretization of continuous dynamics “as if” it was operational.

Not surprisingly, such an idea is quite ancient. Iwasaki et al.  [37] first proposed using non-standard analysis to discuss the nature of time in hybrid systems. Bliudze and Krob  [22] , [21] have also used non-standard analysis as a mathematical support for defining a system theory for hybrid systems. They discuss in detail the notion of “system” and investigate computability issues. The formalization they propose closely follows that of Turing machines, with a memory tape and a control mechanism.

The introduction to non-standard analysis in  [21] is very pleasant and we take the liberty to borrow it. This presentation was originally due to Lindstrøm, see  [41] . Its interest is that it does not require any fancy axiomatic material but only makes use of the axiom of choice — actually a weaker form of it. The proposed construction bears some resemblance to the construction of as the set of equivalence classes of Cauchy sequences in modulo the equivalence relation (un)(vn) iff limn(un-vn)=0.

Motivation and intuitive introduction

We begin with an intuitive introduction to the construction of the non-standard reals. The goal is to augment {±} by adding, to each x in the set, a set of elements that are “infinitesimally close” to it. We will call the resulting set *. Another requirement is that all operations and relations defined on should extend to *.

A first idea is to represent such additional numbers as convergent sequences of reals. For example, elements infinitesimally close to the real number zero are the sequences un=1/n, vn=1/n and wn=1/n2. Observe that the above three sequences can be ordered: vn>un>wn>0 where 0 denotes the constant zero sequence. Of course, infinitely large elements (close to +) can also be considered, e.g., sequences xu=n, yn=n, and zn=n2.

Unfortunately, this way of defining * does not yield a total order since two sequences converging to zero cannot always be compared: if un and un' are two such sequences, the three sets {nun>un'}, {nun=un'}, and {nun<un'} may even all be infinitely large. The beautiful idea of Lindstrøm is to enforce that exactly one of the above sets is important and the other two can be neglected. This is achieved by fixing once and for all a finitely additive positive measure μ over the set of integers with the following properties: (The existence of such a measure is non trivial and is explained later.)

  1. μ:2{0,1};

  2. μ(X)=0 whenever X is finite;

  3. μ()=1.

Now, once μ is fixed, one can compare any two sequences: for the above case, exactly one of the three sets must have μ-measure 1 and the others must have μ-measure 0. Thus, say that u>u',u=u', or u<u', if μ({nun>un'}=1), μ({nun=un'})=1, or μ({nun<un'})=1, respectively. Indeed, the same trick works for many other relations and operations on non-standard real numbers, as we shall see. We now proceed with a more formal presentation.

Construction of non-standard domains

For I an arbitrary set, a filter over I is a family of subsets of I such that:

  1. the empty set does not belong to ,

  2. P,Q implies PQ, and

  3. P and PQI implies Q.

Consequently, cannot contain both a set P and its complement Pc. A filter that contains one of the two for any subset PI is called an ultra-filter. At this point we recall Zorn's lemma, known to be equivalent to the axiom of choice:

Lemma 1 (Zorn's lemma) Any partially ordered set (X,) such that any chain in X possesses an upper bound has a maximal element.

A filter over I is an ultra-filter if and only if it is maximal with respect to set inclusion. By Zorn's lemma, any filter over I can be extended to an ultra-filter over I. Now, if I is infinite, the family of sets = {PIPcisfinite} is a free filter, meaning it contains no finite set. It can thus be extended to a free ultra-filter over I:

Lemma 2 Any infinite set has a free ultra-filter.

Every free ultra-filter over I uniquely defines, by setting μ(P)=1 if P and otherwise 0, a finitely additive measure (Observe that, as a consequence, μ cannot be sigma-additive (in contrast to probability measures or Radon measures) in that it is not true that μ(nAn)=nμ(An) holds for an infinite denumerable sequence An of pairwise disjoint subsets of .) μ:2I{0,1}, which satisfies

μ ( I ) = 1 and, if P is finite, then μ ( P ) = 0 .

Now, fix an infinite set I and a finitely additive measure μ over I as above. Let 𝕏 be a set and consider the Cartesian product 𝕏I=(xi)iI. Define (xi)(xi') iff μ{iIxixi'}=0. Relation is an equivalence relation whose equivalence classes are denoted by [xi] and we define:

* 𝕏 = 𝕏 I / (1)

𝕏 is naturally embedded into *𝕏 by mapping every x𝕏 to the constant tuple such that xi=x for every iI. Any algebraic structure over 𝕏 (group, ring, field) carries over to *𝕏 by almost point-wise extension. In particular, if [xi]0, meaning that μ{ixi=0}=0 we can define its inverse [xi]-1 by taking yi=xi-1 if xi0 and yi=0 otherwise. This construction yields μ{iyixi=1}=1, whence [yi][xi]=1 in *𝕏. The existence of an inverse for any non-zero element of a ring is indeed stated by the formula: x(x=0y(xy=1)). More generally:

Lemma 3 (Transfer Principle) Every first order formula is true over *𝕏 iff it is true over 𝕏.

The above general construction can simply be applied to 𝕏= and I=. The result is denoted *; it is a field according to the transfer principle. By the same principle, * is totally ordered by [un][vn] iff μ{nun>vn}=0. We claim that, for any finite [xn]*, there exists a unique 𝑠𝑡([xn]), call it the standard part of [xn], such that

𝑠𝑡 ( [ x n ] ) and 𝑠𝑡 ( [ x n ] ) [ x n ] . (2)

To prove this, let x=sup{u[u][xn]}, where [u] denotes the constant sequence equal to u. Since [xn] is finite, x exists and we only need to show that [xn]-x is infinitesimal. If not, then there exists y,y>0 such that y<|x-[xn]|, that is, either x<[xn]-[y] or x>[xn]+[y], which both contradict the definition of x. The uniqueness of x is clear, thus we can define 𝑠𝑡([xn])=x. Infinite non-standard reals have no standard part in .

It is also of interest to apply the general construction (1 ) to 𝕏=I=, which results in the set * of non-standard natural numbers. The non-standard set * differs from by the addition of infinite natural numbers, which are equivalence classes of sequences of integers whose essential limit is +.