I did not say this in the question, but I have a very specific system of this form which I want to analyze, and whose behavior I do understand when $x,A,B$ are in $R^3$. I'm looking for any related work out there. Its true that a complete classification of the behavior of such systems for all dimensions seems like too much to hope for.
–
alexOct 29 '10 at 0:04

What are the restrictions on $A, B,$ and $x(0)$ in the type of system which you do already understand? What physical system are you modeling? (My question is merely out of curiousity; sometimes, however, knowing the underlying system being modeled helps to constrain some of the equations and assumptions about the dynamical equations...)
–
sleepless in beantownOct 29 '10 at 0:43

One keyphrase that might be remotely related is iterative rounding, used in multiobjective optimization problems to obtain approximations.
–
Joseph O'RourkeOct 29 '10 at 0:57

3 Answers
3

The dynamical system you have presented is in the class of discrete-time piecewise-affine (PWA) dynamical systems, and they are notoriously hard to analyze. The general form of a discrete-time PWA system is:

where $x \in \mathbb{R}^n$ is the state vector, $u \in \mathbb{R}^m$ is the input vector, $c_i$ is a constant vector, and $\Omega_i$ is a region of the state-space. Note that $x^+$ is the updated state vector. Essentially, we have different affine dynamical systems on different "patches" of the state-space $\mathbb{R}^n$. This has been studied quite a bit in the past couple of decades by the control theory community. A good paper for you to start would be:

Just to add to this: rounding is a bit of a tricky function to handle though. Although I suppose one could obtain rounding behavior (en.wikipedia.org/wiki/Floor_and_ceiling_functions#Rounding) through the use of integer variables (switch between floor/ceiling functions depending on sign), which leads to MLD (Mixed Logical Dynamical) systems, which has been shown to be equivalent to PWAs in W. Heemels, B. D. Schutter, and A. Bemporad, "Equivalence of hybrid dynamical models", Automatica, vol. 37, no. 7, pp. 1085-1091, July 2001
–
GileadNov 10 '10 at 6:14

I would suggest searching for round-off error, floating-point representation and floating-point precision and numerical simulation of dynamical systems.

This is equivalent to doing numerical simulation computations using integers for all of the variables. Integers, short and long, signed and unsigned, are represented using a fixed number of bits in computational systems. Floating point numbers are also represented with a fixed number of bits, with a certain number allocation to the mantissa, and the remainder of bits allocated to the exponent in base $2$, representing what the mantissa is multiplied by

This is also equivalent to what must happen in fixed-digit or limited precision representation of floating point digits in any computerized simulation of any dynamical systems. The numerical simulation is only correct up to a certain number of digits, and the imprecision can build up rather quickly depending upon the variance of the elements of the matrix. Rounding occurs in every floating-point representation using a fixed number of bits; it just occurs at smaller magnitudes with larger number of bits used for the mantissa.

There are two points to consider:

fixed point errors in the numerical simulation (or round-off error), due to limited precision representation of floating point numbers or fixed-number-of-bits representation of integers

the dynamics of the system itself. The dynamics may be such that the sequence is diverging, and that the series will quicly go beyond the limit of what the numerical representation system is capable of dealing with: effectively hitting infinity (for the representation system) and not being able to go further. Even if the dynamics are such that the system is oscillatory or chaotic around an attractor, the round-off error in the numerical simulation may begin to dominate after a set number of iterations.

You will find a lot of literature about this from the 1950's and 1960's about matrix representations and matrix manipulations with computers, and a lot about this if you research floating-point representation standards.

A simple example which leads to open questions (afaik) is $(x,y)->[(\frac{x}{2}-y,x)]$. Without the rounding, the orbit of a given point belongs to an ellipse $x^2-\frac{xy}2+y^2=k$. This rounded version is a bijection on $\mathbb{Z} \times \mathbb{Z}$. It seems likely that all orbits are finite, but afaik this is not known to be true. A relevant reference seems to be The arithmetic of discretized rotations .

Notes:

This is not exactly in the form you request but it is equivalent since $A$ is the $0$ matrix.

I recall a problem like this floating around Ohio State University in the late 80's although it may not have been exactly this.