
2
membership approach is adopted to provide robust constraint
satisfaction with online parameter update [18]–[21]. However,
most of the existing works of robust adaptive MPC were
designed for the regulation problem and few works can be
found for the tracking problem of systems with constraints in
the presence of external disturbances and uncertain parameters.
Inspired by the timely needs of operating wheelchairs
safely in different application scenarios and for diverse users,
the main contribution of this work is to present a safety-
based constrained tracking control algorithm for wheelchair
systems with external disturbances and uncertain parameters
at the dynamics level. Under the assumption of unknown-but-
bounded disturbances, the set-membership approach is adopted
to provide an updated set of uncertain parameters which is
used in the subsequent MPC design. For any infeasible refer-
ence that violates the safety constraints, an optimization-based
method is utilized to compute the closest admissible reference
for tracking. The state and input constraints are robustly
satisfied in the proposed MPC framework with reference-
dependent and tube-based constraints in vertex representation
of states. The recursive feasibility and input-to-state stability
of the wheelchair system with the proposed MPC controller
are guaranteed. The effectiveness of the proposed safety-based
control algorithm is validated by two speed tracking tasks on
the high-fidelity model of a practical wheelchair.
The remainder of this paper is organized as follows: In
Section II, the dynamics of the wheelchair system at actuator
level is described and the system in the LPV form with
unknown parameters is discussed. The proposed robust adap-
tive tracking controller with details in parameter estimation,
feasible reference generation, robust constraint satisfaction and
tracking MPC formulation are presented in Section III. The
tracking results of the proposed safety-based control algorithm
are demonstrated in Section IV. Section V concludes this
work.
Notations: When defining the variable, we follow the rule
that capitalized letters are for matrices and small letters are
for vectors or scalars. Rand Zare the sets of real and integer
numbers. Given two integers a,b∈Z,Za+,{i∈Z|i≥a}
and Z[a,b],{i∈Z|a≤i≤b}. The i-th row of matrix
Xand the i-th element of vector xare represented by X[i]
and x[i], respectively. The i-th vertex of x∈ X is denoted
by x(i). A non-negative matrix is denoted by X≥0. The
positive definite and semi-definite matrices are represented as
X0and X0, respectively. For a matrix X∈Rn×n,
the smallest eigenvalue is denoted by λ(X). The identity
matrix of dimension nis denoted by Inand an m-dimension
vector with all elements as 1 is denoted by 1m. The m×n
matrix with all elements as zero is denoted by 0m,n. A
diagonal matrix with main diagonal elements a1, . . . , anis
denoted by diag(a1, . . . , an). The following sets are defined:
Sn={X∈Rn×n:X=X>},Sn
0={X∈Sn:X0}
and Sn
0={X∈Sn:X0}. A convex polyhedral set of x
is defined as Px(Fx, bx) = {x|Fxx≤bx}. With P∈Sn
0, an
ellipsoidal set of xis defined as E(P, 1) = {x|x>P x ≤1}.
For a vector x∈Rnwith a matrix Q,kxkdenotes the 2-norm
and kxkQstands for px>Qx. The vector xi|krepresents the
predicted value of xat a sampling time instant k+ibased on
measurement at k. The vector x(k)stands for the measured
value of xat a sampling time instant k. A continuous function
f: [0, a)→[0,∞)belongs to class K∞if a=∞and
f(r)→ ∞ as r→ ∞.
II. SYSTEM DESCRIPTION AND PROBLEM FORMULATION
A. Wheelchair Dynamics
With implicit variable change for rotary to translational
movement, the electrical and mechanical subsystems of elec-
tric powered wheelchair are given as follows [12]:
L˙
ic+Ric+Kev=u, (1a)
M˙v+Dv +wf=Ktic,(1b)
where v= [v1, v2]>is the vector of the linear velocities, u=
[u1, u2]>is the vector of the motor voltages, ic= [ic1, ic2]>
is the vector of the currents and wf= [w1, w2]>is the vector
of the lumped disturbance torques on the right and left wheels,
respectively. Mand Dare the equivalent mass and damping
coefficient matrices, which can be expressed as
M=m11 m12
m21 m22 , D = diag(d1, d2),
where m11,m12,m21 and m22 are scalars. The non-zero
matrix Mcouples the dynamics of left and right wheels,
and d1and d2are the damping coefficients for the right
and left wheels, respectively. Furthermore, L= diag(l1, l2)
and R= diag(r1, r2)are the inductance and resistance
matrices of system, respectively. Ke= diag(ke1, ke2)and
Kt= diag(kt1, kt2)are the back electromotive force constant
and torque constant matrices, respectively.
For system without current feedback ic, the dynamics in
(1a) and (1b) can be reformulated as
¨v+ (M−1D+L−1R) ˙v+ Γv=M−1L−1Ktu+w, (2)
where w=−M−1˙wf−M−1L−1Rwfand Γ =
M−1L−1RD +M−1L−1KtKeare the lumped terms.
Let x,[v>,˙v>]>be the state vector. By using the Euler
forward approximation, the discrete-time state-space model
with a sampling time interval Tscan be formulated as
x(k+ 1) = Ax(k) + Bu(k) + Ew(k)
=I2TsI2
−TsΓI2−Ts(M−1D+L−1R)x(k)
+02,2
TsM−1KtL−1u(k) + 02,2
TsI2w(k).
(3)
All the states x(k)and inputs u(k)are required to satisfy the
following constraints to guarantee the safety of the wheelchair:
Gx(k) + Hu(k)≤b, (4)
with given matrices G∈Rnc×4,H∈Rnc×2and b∈Rnc.
Remark 1. The safety constraints in (4) are defined in a
general form, where bis not necessarily a vector with all
elements of 1. This form offers more degrees of freedom