Convergence properties of indefinite linear quadratic problems
with receding horizon
Citation for published version (APA):
Trentelman, H. L., & Soethoudt, J. M. (1988). Convergence properties of indefinite linear quadratic problems with receding horizon. (Memorandum COSOR; Vol. 8825). Technische Universiteit Eindhoven.
Document status and date: Published: 01/01/1988
Document Version:
Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers)
Please check the document version of this publication:
• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.
• The final author version and the galley proof are versions of the publication after peer review.
• The final published version features the final layout of the paper including the volume, issue and page numbers.
Link to publication
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal.
If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:
www.tue.nl/taverne Take down policy
If you believe that this document breaches copyright please contact us at: openaccess@tue.nl
EINDHOVEN UNIVERSITY OF TECHNOLOGY
Department of Mathematics and Computing Science
Memorandum COSOR 88-25
CONVERGENCEPROPERTffiSOF
INDEFINITE LINEAR QUADRATIC
PROBLEMS WITH RECEDING HORIZON
HL Trentelman and
J.M. Soethoudt
Eindhoven University of Technology
Department of Mathematics and Computing Science
P.O. Box 513
5600 MB Eindhoven
The Netherlands
Eindhoven, October 1988
The Netherlands
CONVERGENCE PROPERTIES OF INDEFINITE LINEAR QUADRATIC PROBLEMS WITH RECEDING HORIZON
H.L. Trentelman & J.M. Soethoudt
Faculty of Mathematics & Computer Science Eindhoven University of Technology,
P.O. Box 513, 5600 MB Eindhoven,
The Netherlands.
Abstract: In this paper we study the following question: given a finite dimensional linear system together with a finite horizon (possibly indefinite) quadratic cost functional, when does the corresponding optimal cost converge to the optimal cost of the corresponding infinite horizon problem, as the length of the horizon tends to infinity? For the case that the linear quadratic problems are regular we establish necessary and sufficient conditions for this convergence to hold.
Keywords: Linear quadratic optimal control, indefinite cost, large horizon, Riccati equation.
1.
INTRODUCTIONThe finite horizon linear quadratic control problem for the linear time-invariant system
(1.1) x(t)
.
Ax(t) + Bu(t), x(O)is concerned with choosing a control function u such that the cost functional
(1. 2)
T
JT(xo'u) :=
J
w(x(t),u(t»dt ois minimized. Here,
w
is a real quadratic form on Rnx Rmgiven by(1.3) w(x,u) xTQ x
+
2u S xT+
u R uT .h b . IOnxn IOnxm Q
e
..,nxn. Se
..,mxn and Re
..,mxm. In t e a ove express~ons, Ae l l ' ,
Be
1 1 \ , lI' lI' lI'It is assumed that Q is symmetric and that R is positive definite. Apart from this definiteness assumption on R, we allow
w
to be indefinite.In order for the integral in (1.2) to be well-defined, we restrict the control functions to be elements of the class
functions that are square integrable over [0, T] . above problem is then defined as
L [O,T] of all Rm-valued 2
The optimal cost for the
(1.4)
I
ue
L [0, T] }.2
In addition to the above, the infinite horizon linear quadratic problem deals with minimizing the indefinite integral
(1.5) J(x ,u) := limJT(xO'u).
o T~oo
on
~+.
The functional J(x ,u) is defined for all owhich the limit in (1.5) exists in the sense that Let L (R+):=
n
L2[O,T] denote the space of all Rm-valued 2.1oc
T~O square integrable functions
• ..,+
funct~ons u
e
L (lI') for2.1oc
i t is finite or infinite. This class of functions is denoted by
locally
(1. 6) U(x) := { u e L (~+)
o 2.1oc lim JT~oo T(x ,u)0 exists in
The optimal cost for the above infinite horizon linear quadratic problem is
(1.7) V (x )+
.=
inf { J(x ,u)I
ue
U(x ) }.f 0 0 0
for all x (provided that (A,B) is stabilizable). I t o
if
w
is indefinite then this convergence no longer In this paper we are interested in the question whether the optimal cost for the finite horizon problem (1.4) converges to the optimal cost for the infinite horizon problem (1.7) as T tends to infinity. Of course, if the quadratic formw
is positive semi-definite then it is well-known that indeed+
~ V ( x ) (T-+co) f 0
out however that
holds in general. In this paper we establish necessary and sufficient conditions for convergence to hold.
Actually, we shall treat the above question of convergence in the
IRnxn •
following, more general context. Let N
e
be a symmetr~c matrix, let T>
0 and consider the finite horizon problem with cost functional(1. 8)
The second term in the above represents a penalty on the terminal state. The matrix N is allowed to be indefinite. The optimal cost associated with the
latter problem is given by
(1.9)
In addition to this finite horizon problem we consider the infinite horizon problem of infimizing (1.5) under the constraint that Nx(t) converges to zero as t tends to infinity. More specifically, let
{ u e
U(x )o lim Nx(t)T-+CO
o }
and consider the problem of inf imizing (1.5) over the class UN (x0)' The optimal cost for this problem is given by
(1.10) V (x)
+
:= inf { J(x ,u)I
u e UN(x o) }.N 0 0
The latter optimization problem was studied in detail in [8]. Of course, the problem (1.7) can be reobtained from this formulation as a special case by taking N = 0 (see also [9]).
I , the identity Another special case of (1.10) is obtained by taking N
matrix. This special case was treated in [10].
Now, in this paper we shall ask ourselves the question: when does the optimal cost for the finite horizon problem (1.9) converge to the optimal cost for the infinite horizon problem (1.10) as T ~ m ?
We conclude this introduction by noting that the questions to be studied here have been studied before in [2] and [11]. However, in these references only the case that both
w
~ 0 and N ~ 0 was considered, while we intend to treat the most general case that wand N are allowed to be indefinite.2
THE
ALGEBRAIC RICCATI EQUATIONThe characterization of the optimal costs for the infinite horizon problems (1.7) and (1.10) centers around the algebraic Riccati equation (ARE):
(2.1 )
o.
Let
r
be the set of all real sYmmetric solutions of (2.1). According to [10] , if (A,B) is controllable andr
¢ 0 then there is exactly one K Er
such that. -1 T -
°
the matr~x A
K := A - BR (B K
+
S) has all its eigenvalues in C v C and exactly one K Er
such that A has all its eigenvalues in C+v Co. Here, weK
- °
+define C (C ,C) := { S E C
I
Re s < 0 (Re s = 0, Re s>
0 ). These +elements of
r
are denoted by K and K, respectively. It can be shown that if+ +
K E
r
then K ~ K ~ K . We denote AK- by A and AK+ by A .
If M E Rnxn then
X-(M)
( XO(M), X+(M) )
denotes the span of all generalized-
°
+eigenvectors of M corresponding to its eigenvalues in C ( C , C ).
- + - +
Let Q denote the set of all A -invariant subspaces of
X
(A ). Let ~ := K K (the 'gap' of the ARE ). The following result states that there exists a bijection between Q andr :
Theorem 2.1 ( [10], [3], [7] ). Let (A,B) be controllable and assume that
r
¢ 0. IfV
E Q thenR
n =V
@~-lV~.
There exists a bijectionr :
Q~
r
defined by-1 ~ n ~
where P
v
is the projector along ~ V := { x E RI
~x E V }. If KX+
(AK)
=
V , XO
(AK)=
ker ~ andX-
(AK)=
X-
(A+) f"I ~-lV~. 0r
(V) then3
THE INFINITE HORIZON PROBLEM WITH ASYMPTOTIC CONSTRAINTSIn this section we briefly recall the results from [8] on the infimization problem (1.10). Let
1£
be a subspace of IRn• A symmetric matrix K e IRnxn is called negative semi-definite on1£
if the following two conditions hold: (i) V xe
1£ :
xTK x S 0 (ii)V
xe
1£ :
xTK x = 0*
Kx =O.
If
V
S;; IRnand Me
IRnxn then<V
1M> will denote the largest M-invariant + subspace ofV.
A key role in the characterization of the optimal cost VN(xO) is played by the subspace
(3.1 )
Observe that
V
N e Q. Thus, with
V
N there corresponds exactly one solution ofthe ARE. This solution
~(VN)
is denoted byK~.
The following theorem is the main result of [8]:Theorem 3.1 Let (A, B) be controllable. Assume that
r
*
l2l and that K is negative semi-definite on ker N. Then we have. + T + n
(~) VN(x ) = x K x for all x e IR .
o 0 N 0 0
(ii) For all x e IRn there exists an optimal o
ker ~ S;; ker N n ker K
u
*
i f and only if(iii) If ker ~ S;; ker N n ker K then there exists exactly one optimal input and, moreover, this input is given by the feedback control law u
-1 T + -R (B K
N + S)x. 0
4
THE FINITE HORIZON PROBLEM WITH ENDPOINT PENALTYIn this section we consider the finite horizon problem of infimizing the cost functional (1.8). We note that for the case that both W ~ 0 and N ~ 0 this problem is quite standard and is treated, for example, in [6]. The general case however is slightly more complicated. The following result can be found in [1,p.13l]:
Lemma 4.1 Suppose there exists on the interval [0,00) a solution K(t) of the Riccati differential equation (ROE):
(4.1) K(t)
T T -1 T
A K(t) + K(t)A + Q - (K(t)B + S)R (B K(t) + S), K(O) N.
Then for all T
>
0 and for all x e Rn we have oT +
x KN(T)X .
o 0
* *
For all T
>
0 and x e Rn there is exactly one u e L[O,T] such that JT N(xO'u )
+ 0 -1' T
VT,N(x
O). This input is given by the feedback control law u = -R (B K(T-t) + S)x; te[O,T]. [J
It is well-known that if
w
~ 0 and N ~ 0 then (4.1) indeed has a unique solution on [0,00). In the general case, the ROE need not have a solution on a given interval. We do have the following:Lemma 4.2 Suppose that (A,B) is controllable and that
r
~ 0. If N - K ~ 0 then (4.1) has a unique solution KN(t) on [0,00). In fact, KN(t) = K + D(t), where D(t) is the unique solution on [0,00) of(4.2)
.
D(t) A-TD (t ) + D (t ) A- - D (t ) BR-1B D (t) ,T D(0) N - K
(with A = A -, see section 2) . K
Proof The fact that (4.2) has a unique solution D(t) on [0,00) is standard (see for example [4 ,cor. 2.4.4]). It is then a matter of straightforward calculation to show that K-+ D(t) satisfies (4.1). Uniqueness of K + D(t) then follows from the uniqueness of D(t). [J
5
CONVERGENCE OF THE OPTIMAL COSTIn this section we shall give a formulation of our main result. Before doing t h i s , we state the following lemma:
n lRn= . , Lemma 5.1 Let
:e
be a subspace of lR and let K E be a symmetr~c matr~x.Then K is negative semi-definite on
:e
if and only if there exists a symmetric matrix N such that ker N =:e
and N - K ~ O.Proof A proof of this can be given similar to the proof of [8,lemma 3.3]. c
Consider the problems (1.8) and (1.10). In the remainder of this section we assume that (A,B) is controllable and that
r
~ 0. According to the previous lemma, if N - K- ~ 0 then K is negative semi-definite on ker N. Conversely, if K- is negative semi-definite on ker N then. one can always find a symmetric matrix N such that ker N=
ker Nand N - K ~ O.1 1 1
We now formulate our main result:
Theorem 5.2 Assume that (A,B) is controllable,
r
~ 0, N-
K ~ 0 and ker 6 S;;+
+
all lRni f and only i f
ker N ('\ ker K
.
Then VT,N(xo ) -+ VN(xo) (T-+oo) for x0E ker N ('\ ker K is A--invariant.
A proof of this theorem is given in section 6. In the remainder of this section we shall discuss some important special cases of this theorem. First
+ we consider the case that N = 0, that is , the question whether VT(x
o) given
+
by (1.4) converges to the optimal cost Vf(x
o) given by (1.7). (T-+oo ) Corollary 5.3 ker K-. Then A- -invariant. that (A,B)
+
-+ Vf(x o)is controllable,
r
~ 0, K s 0 and ker6
S;; for all x if and only if ker K- iso
It is well-known that if, in addition,
w
~ 0 then if (A,B) is controllable+
+
we always have convergence of V (x ) to Vf(x ). This result can be reobtained T o o
using the previous corollary in combination with the following lemma:
Lemma 5.4 Assume that
w
~ O. If (A, B) is controllable thenr
~ 0, K SO, ker 6 S;; ker K-and ker K- is A -invariant.Proof A proof of these claims can be found in [9]. c
Next, we consider the special case that N ~ 0 and
w
~ O. It turns out that ifw
~ 0 then the subspace ker K- can be calculated in terms of the matrices A, B, Q, Sand R. In fact, i t can be shown that if (A,B) is controllable thenker K <ker (Q - S TR-1S)
I
A - BR-1S>+ -1 0 -1
() { X
(A - BR S) (t)X
(A - BR S) },exactly the undetectable subspace of the respect to the stability set
C-).
Denote-1 A - BR S).
[9]). The latter subspace is
(Q - S TR-1 S , A - BR-1S) (with T -1 subspace by HV«Q - S R S, (see pair this
Corollary 5.5 Assume
w
~0
andN
~O.
If (A,B) is controllable and kerA
~T -1 -1 + +
ker N " HV( (Q - S R S, A - BR S). Then V (x) -+ Vf(x) (T-+oo) for all x
T,N 0 0 0
T -1 -1 -1
if and only if ker N " HV«Q - S R S, A - BR S) is (A - BR S)-invariant.
c
A similar result was obtained in [2], with' (A,B) controllable' replaced by
T -1 -1
, (A, B) stabilizable' and 'ker A S; ker N () HV(Q - S R S, A - BR S)' replaced by 'A > 0'. The methods used in [11] ( and [2]) to prove this result will however turn out to be completely different from ours.
Note that the results presented up to now all require N - K ~
O.
We now ask ourselves the question: what can we do if this is not the case'? Again consider the problems (1.9) and (1.10) and assume that K is negative semi-definite on ker N. (see th.3.1). Now, the latter does not imply that N -K ~ 0 but according to lemma 5.1 we can find a new weighting matrix N such1
and x
o with N and N - K ~ O. Clearly, U (x ) =' U (x) for all
1 N o N 0
1
for all x . We can then apply theorem 5.2 o = ker N 1
+
=' V (x ) N 0 + hence V (x) N 0 1+
+
replaced by N1 and study the question whether VT,N (xo) converges to VN (xo)
+ 1 1
VN(x
o)' An example of this is following: suppose we consider the problem with N = I. There always exist a real number k
*
0 such that kI - K ~O.
Of course,+ +
VI (x
o) = VkI(xo) for all xo' We can then apply theorem 5.2 to obtain: that ker N
+ exists k
*
0 such that VT,kI(xO) convergence holds for all k
*
0 suchCorollary 5.6 Assume that (A,B) is controllable,
r
*
+
-+ Vr(x
o) (T-+oo) that kI ~ K . c
o and
A
>
O.
Then there for all x . In fact,o
6
PROOF OF THE MAIN RESULTOur proof of theorem 5.2 runs along a series of lemmas that we consider to be interesting in their own right. Due to lack of space the proofs of these lemmas are deferred to a future paper. Our first lemma deals with an arbitrary system (A, B) and an arbitrary matrix R
>
0, independent of the previous context. Consider the standard Riccati differential equation(6.1) P (t) P (0) T -1 T A P (t)
+
P (t) A - P (t) BR B P (t) , P , otogether with the standard algebraic Riccati equation
(6.2)
Recall that if P ?; 0 and (A, B) is controllable, then (6.1) has a unique
o
solution P (t) ?; 0 on [0,(0) (see [4]). Also, (6.2) has at least one solution
+
([5]). Let P be the largest real symmetric solution of (6.2).
Lemma 6.1 Assume that (A,B) is controllable and ~(A) ~ C+. Then p+
>
O.
For+
any P
>
0 we have lim P(t) = P . c o t-+CX1Our following result again deals with the Riccati differential equation (6.1). The result
is,
in a sense, the converse of the previous lemma:Lemma 6.2 Assume P ?; 0 and assume that lim P(t) =: P E R
nxn
exists. If P
>
0o t-+CX1
then P
> O.
c oWe now return to the original context of this paper. Consider the system (1.1), together with the quadratic form
w
given by (1. 3). Recall that we denote AK- by A (see section 2). An important role
is
played by the following algebraic Riccati equation in the unknown D:(6.3)
o
We make the following observation:
Dnxn .
Lemma 6.3 Assume (A,B)
is
controllable andr
*
0. Let K E Il' be symmetr~c. Then Kis
a solution of (2.1) if and only if D = K - Kis
a solution of+
(6.3). In particular D = 0 and D = f!. (= K - K) are solutions of (6.3). In fact, 0 and f!. are the extremal solutions of (6.3) in the sense that any solution of (6.3) satisfies 0 ~ D ~ f!.. c
Finally, we shall need the following result:
Lemma 6.4 Assume that (A,B) is controllable and
r
*
0. Assume that N - K?;O.
We have now collected the most important ingredients that will be used in our proof of theorem 5.2. In order to give this proof we shall make a suitable decomposition of the state space. Let VN be the subspace defined in (3.1). Define X := V 1 N X 2:= ker t:. , X
.=
X-(A+) " t:.-1V.L 3· N Denote A N:= + 0AK+.According to theorem 2.1 we have
X
=X
(A ),X
=X
(AN) and1 N 2
N
Hence IRn = Xl $
X
2 $
X
3 • With respect to this decomposition we(6. 4)
have
This follows from the fact that both
X
1
s;; C+. According to theorem 2.1, for given matrices A
ij
A- -invariant. Note that IT(A )
33
t:.-1X.L. Since also
X
= ker t:., we1 2 (6.5)
o
o
o
with t:.ll
>
0 andt:.
33>
O. Finally, we partitionB [ : : ] •
Proof of theorem 5.2 From the assumption ker t:. S;; ker N " ker K i t follows that Xl $
X
2 S;; ker N " ker K- Hence Nand K have the form
~ ~-],
33with N - K >
o.
Since K+= K + t:. , we have[
0 0 0J.
+ 0 0 0 K N 0 0 K+
33(Recall + + is the projector
X=
VN along that K
N KP
+
K (I - P), where P onto 1l1-1V+
X
e
X
3) • By combining the above we see that
N 2
Let D(t) be the solution of the Riccati differential equation (4.2). By lemma 6.4 and the fact that
X
= ker £1 is A--invariant, we have2
o
o
o
for some D (t). By writing out (4.2) in the decomposition employed, we see 33
that D (t) is the unique solution of 33 (6.6) D (t) 33 D (0) 33 AT D (t) 33 33 + N - K 33 33 + also, K
N - K
is
a solution of the algebraic Riccati equation (6.3) (see lemma (6.3). This implies that ~33 is a solution to(6.7)
In fact, since ~ is the largest solution·of (6.3), ~33 is the largest solution of (6.7). that , x =
>
o.
+ -+ K N+
+
('~ ') Assume that VT,N(xO) ~ VN(xO) for all Xo or, equivalently, KN(t)
(t~oo). Here, KN(t) is the unique solution of (4.1). Then we have KN(t) - K
+
K
N - K (t~oo). Now, the point is that KN(t) - K = D (t). Hence we find D (t) ~ £1 (t~oo). Since £1
>
0 i t follows from lemma 6.2 that N - K33 33 33 33 33
We now prove the A--invariance of ker N n ker K . Let x e ker N n ker K (x ,x ,x ). Then (N - K )x = 0
1 2 3
(A x , A x , 0). Since
X
eX
11 1 22 2 1 2
('~ ') Assume ker N n ker K ker N n ker K-. Indeed,
whence (N - K- ) x
=
0 so x=
o.
Thus A x =33 33 3 3
~ ker N n ker K-, the claim follows.
invariant under A-. We then claim that Xl e
X
- - +
-t:..
X @ X <ker N f'I ker K
IA
>
f'I X (A ) } @ ker1 2
+
-ker
t:..
(ker N f'I ker K f'I X (A ) } @
+
-ker N f'I ker K f'I ( ker
t:..
@ X (A ) ),
where the last equality again uses the assumption ker
t:..
~ ker K f'I ker N. Now, kert:..
@ X+(A-)=
XO(A-) $ X+(A-)=
lRn• This proves the claim. It follows from this that N - K>
O. Hence, since cr(A ) ~c+,
the solution D (t) of (6.6)33 33 33 33
converges to t:..
33, the largest solution of (6.7) (see lemma 6.1). In turn this
+ + +
+
implies that D
33(t) ~ KN - K or, equivalently, that KN(t) ~ KN· Thus VT,N(xO)
+
~ V (x) (T~~) for all x . This completes our proof of theorem 5.2. c
N 0 0
REFERENCES
[1] R.W. Brockett, Finite Dimensional Linear Systems, New York, Wiley, 1970. [ 2 ] F.M. Callier & J. L. Willems, , Criterion for the convergence of the
solution of the Riccati differential equation' IEEE Trans. Aut. Control,
AC-26, p. 1232-1242, 1981.
[3] W.A. Coppel, 'Matrix quadratic equation', Bull. Austral. Math Soc., Vol. 10, p. 377-401, 1974.
[4] O. Hijab, Stabilization of Control Systems, New York, Springer Verlag, 1987.
[5] R.E. Kalman, 'Contributions to the theory of optimal control', Bol. Soc. Mat. Mex. Vol.5, p. 102-199, 1960.
[6] H. Kwakernaak & R. Sivan, Linear Optimal Control Systems, Wiley, New York, 1972.
[7] M Shayman, 'Geometry of the algebraic Riccati equation - part 1', SIAM Journal of. Contr. &Opt., Vol. 21, no 3, p. 375-393, 1983.
[8] J.M. Soethoudt & H.L. Trentelman, 'The regular indefinite linear quadratic problem with linear endpoint constraints', Memorandum Eindhoven University of Technology, 1988, to appear in Systems & Contr. Letters.
[9] H.L. Trentelman 'The regular free-endpoint linear quadratic problem with indefinite cost', Memorandum Eindhoven University of Technology, 1987 , to appear in SIAM Journal of Cont.
&
Opt ..[10] J.C. Willems, 'Least squares stationary optimal control and the algebraic Riccati equation', IEEE Trans. Aut. Control, AC 16, no.6, p. 621-634, 1971.
[11] J.L. Willems & F.M. Callier, 'Large finite horizon and infinite horizon LQ-optimal control problems', Optimal Contr. Appl. & Meth., Vol. 4, P 31-45, 1983.
EINDHOVEN UNIVERSITY OF TECHNOLOGY DepaI1Inent of Mathematics and Computing Science
PROBABILITY
THEORY,
STATISTICS,
OPERATIONS
RESEARCH
AND
SYSTEMS
THEORY
P.O.
Box 5135600 MB Eindhoven - The Netherlands Secretariate: Domme1building 0.02 Telephone: 040 - 473130
List of COSOR-memoranda - 1988
Number Month Author Title
M 88-01 January F.W. Steutel, Haight's distribution and busy periods. B.G. Hansen
M 88-02 January 1. ten Vregelaar On estimating the parameters of a dynamics model from noisy input and output measurement.
M 88-03 January B.G. Hansen, The generalized logarithmic series distribution. E. Willekens
M 88-04 January J. van Geldrop, A general equilibrium model of international trade with C. Withagen exhaustible natural resourse commodities.
M 88-05 February A.H.W. Geerts A note on "Families oflinear-quadratic problems": continuity properties.
M 88-06 February Siquan, Zhu A continuity property of a parametric projection and an iterative process for solving linear variational inequalities.
M 88-07 February
1.
Beirlant, Rapid variation with remainder and rates of convergence. WillekensM 88-08 April Jan v. Doremalen, A recursive aggregation-disaggregation methodto approxi-J. Wessels mate large-scale closed queuing networks with multiple job
Number Month Author
-
2-Title
M 88-09 April J. Hoogendoom, The Vax/YMS Analysis and measurement packet (VAMP): R. C. Marcelis, a case study.
AP. de Grient Dreux, J. v.d. Wal,
R.I.
WijbrandsM 88-10 April E.Omey, Abelian and Tauberian theorems for the Laplace transform E. Willekens of functions in several variables.
M 88-11 April E. Willekens, Quantifying closeness of distributions of sums and maxima S.I. Resnick when tails are fat.
M 88-12 May E.E.M. v. Berkum Exact paired comparison designs for quadratic models.
M 88-13 May J. ten Vregelaar Parameter estimation from noisy obseIVations of inputs and outputs.
M 88-14 May L. Frijters, Lot-sizing and flow production in an MRP-environment. T. de Kok,
J. Wessels
M 88-15 June J.M. Soethoudt, The regular indefinite linear quadratic problem with linear H.L. Trentelman endpoint constraints.
M 88-16 July lC. Engwerda Stabilizability and detectability of discrete-time time-varying systems.
M 88-17 August AH. W. Geerts Continuity properties of one-parameter families of linear-quadratic problems without stability.
M 88-18 September W.E.J.M. Bens Design and implementation of a push-pull algorithm for manpower planning.
M 88-19 September AJ.M. Driessens Ontwikkeling van een informatie systeem voor het werken met Markov-modellen.