# An Improved Internal Model Principle Based Multivariable Nonlinear Control Method with Multiclass Nonharmonic Disturbances and Its Application to Speed Control of a Motor Drive System

Authors: Yang Yu; Zengqiang Mi; Tong Zhao; Yikun Xu
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We study the global disturbance rejection problem for a class of general multivariable nonlinear systems with multiclass nonharmonic disturbances. The paper first introduces the importance and state of the art for disturbance rejection problem and describes the control problem in the form of mathematical expressions. It stresses the multiclass disturbances produced by the exosystem satisfying certain characteristic conditions. Then, the nonlinear internal models are designed in accordance with different characteristics of multiclass external disturbances. On the basis of introduction of the control law for disturbance-free system, a multivariable state feedback controller is devised in terms of the designed internal model equations and corresponding assumptions. A Lyapunov function is constructed to theoretically prove the global uniform boundness of all signals for the multivariable closed-loop system. Finally, the presented method is applied to implement the speed control and reject the multiclass nonharmonic disturbances for a two-input motor drive system. The simulation results testify correctness and effectiveness of the presented algorithm.


---

## Body

## 1. Introduction

Due to its importance, the problem of disturbance rejection generated by the exosystem has been paid great attention in the world over the past few decades [1–3]. Especially, this problem in nonlinear system has been concerned in a recent period of time [4–6]. In terms of internal model principle (IMP), one of the key technologies of rejection of exogenous disturbance is to construct reasonable internal model equations, which is able to duplicate the characteristics of the external disturbance.In addition, many industrial applications, such as servo system, wind power system, photovoltaic system, and other modern systems, which rely on power electronic equipments to work, [7–9], are frequently influenced by the external disturbances generated by exosystem [10, 11]. In the circumstance of the existence of external disturbances, the conventional PID controller with constant parameters is difficult to satisfy the requirement of high precision control for these systems [12, 13]. Therefore, rejection of external disturbances by means of new ways [14, 15] is an important work to improve the working performance of these industrial systems. From the state-of-the-art technologies in disturbance rejection, the researchers frequently assume the following:(1)
the external disturbances are frequently sinusoidal; that is to say, the exosystem generally satisfies the condition of neutral stability such as references [16, 17] discussing the suppression of sinusoidal interference with known and unknown frequency, respectively;
(2)
currently, the focus of the research on disturbance rejection is gradually converting from linear systems to nonlinear systems; nevertheless, the disturbance rejections with semiglobal stability [18] and single-input and single-output (SISO) [19] systems are concerned.Therefore, the disturbance rejection problem with nonharmonic interferences and multivariable nonlinear systems has not received adequate attention in current studies. However, the nonharmonic disturbances produced by nonlinear exosystem can induce noise and reduce accuracy to these industrial systems [20]. In addition, just a few articles on nonharmonic interference suppression aim at one class of disturbance [21]. Therefore, the research of multiclass disturbances rejection in multivariable nonlinear system can extend the application of present control theory to a generalized range [22]. Hence, in [23], a global multivariable control algorithm is proposed to reject multiclass external nonharmonic disturbances.In view of the issues described above, the purpose of the paper is to propose an improved internal model principle (IMP) based multivariable nonlinear control algorithm with multiclass nonharmonic disturbances and utilize the proposed method to control a motor drive system to run stably and suppress multiclass external nonharmonic disturbances. As written in [23], the critical points of the problem in disturbance rejection are to model the nonlinear exosystem and present a corresponding control method to suppress the disturbance. Therefore, the essential distinction with the work in [23], this paper addresses the problem of multivariable disturbances rejection with other multiclass exosystems, additionally, which are provided by matrix form. For the difference of external exosystems, the entirely different internal model and a matching control method should be proposed to model and reject the external disturbances in terms of the intrinsic characteristics of exosystem. In practical terms, the construction of internal model replicating with the exogenous interference in the paper is on the basis of self-contained concepts of expanded steady state generator and internal model proposed in [24]. Similarly, the problem of nonharmonic disturbances rejection in the paper is different from the work in [22], which aims to reject a constant disturbance.The innovative works in the paper are the following:(1)
the previous work in references [20, 21] is extended to a generalized field with multivariable inputs of nonlinear systems and multiclass nonlinear exosystems;
(2)
an improved IMP based multivariable nonlinear control algorithm with another multiclass nonharmonic disturbance is presented on the basis of certain characteristic conditions satisfied by exosystem and state information of closed-loop system;
(3)
a new type of nonlinear multivariable internal models for another multiclass nonlinear nonharmonic disturbance is constructed in terms of the expanded steady state generator and internal model.The structures of the paper are arranged as follows. It begins in Section1 with an analysis of the state of the art of disturbance rejection. The formulation of control problem in the paper is given in Section 2. The nonlinear internal models of multiclass external disturbances are designed in Section 3. The nonlinear multivariable state feedback controller is devised in Section 4, and the global convergence of the presented controller is theoretically proved. The implementation of the proposed algorithm to a two-input motor drive system is performed in Section 5 and the conclusions of the research are obtained in the last section.


## 2. Formulation of Control Problem

The multivariable nonlinear system with multiclass external nonharmonic disturbances studied in the paper is described in an affine form shown as below [23]:
(1)
x
˙
=
f
(
x
)
+
∑
i
=
1
m
g
i
(
x
)
(
u
i
+
v
i
(
w
)
)
,
1
≤
i
≤
m
,
where x
∈
R
n indicate the state vectors, u
i
∈
R describe the control inputs, f
(
x
) and g
i
(
x
) denote the known smooth vector fields, v
i
(
w
) represent the external nonlinear disturbances, w
∈
R
q are the exogenous vectors generated by the nonlinear exosystem and shown as follows [24]:
(2)
w
˙
=
s
i
(
w
)
=
A
i
1
w
+
∑
k
=
2
n
A
i
k
w
s
i
k
(
w
)
,
1
≤
i
≤
m
,
where matrices A
i
1 and A
i
k
∈
R
n
q
×
n
q
,
k
=
2
,
…
,
n, and s
i
k
(
w
)
:
R
n
q
→
R are sufficient smooth functions, and s
i
k
(
0
)
=
0.For the stability problems of multivariable input system, the key of the study is to convert multivariable nonlinear system into multiple single-input systems [25].Assumption 1.
For the multivariable nonlinear disturbance-free system shown below:(3)
x
˙
=
f
(
x
)
+
∑
i
=
1
m
g
i
(
x
)
u
i
,
1
≤
i
≤
m
,
there exists control law of sate feedback α
i
(
x
), which can make the system
(4)
x
˙
=
f
(
x
)
+
∑
i
=
1
m
g
i
(
x
)
α
i
(
x
)
asymptotically converge to the origin. In consequence, a Lyapunov function V
(
x
) exists and makes the following inequalities hold [20, 23]:
(5)
d
_
(
∥
x
∥
)
≤
V
(
x
)
≤
d
-
(
∥
x
∥
)
,
∂
V
(
x
)
∂
x
(
f
(
x
)
+
∑
i
=
1
m
g
i
(
x
)
α
i
(
x
)
)
≤
-
d
0
(
∥
x
∥
)
,
|
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
|
2
≤
d
0
(
∥
x
∥
)
,
where d
_, d
-, and d
0 belong to class K
∞ functions.Assumption 2.
The flows of the vector field for nonlinear exosystem (2) are bounded.Remark 3.
Generally speaking, there exist numerous nonlinear systems satisfying Assumption2, such as harmonic functions and limit cycles of nonlinear dynamic system. The well-known Van der Pol oscillator [20] is selected as the exosystem in the paper and can be expressed as below:
(6)
w
˙
1
=
w
2
,
w
˙
2
=
-
a
w
1
+
b
(
1
-
w
1
2
)
w
2
,
where a and b are constants greater than zero. Under such circumstances, the limit cycle generated by Van der Pol oscillator is stable. Consequently, (6) is written in the form of (2), and we can obtain
(7)
w
˙
=
A
11
w
+
A
12
w
s
12
(
w
)
,
where A
11
=
[
0
1
-
a
b
], A
12
=
[
0
0
0
-
b
], s
12
(
w
)
=
w
1
2.Assumption 4.
There exist a constantr
i and a set of real numbers a
i
0
,
a
i
1
,
…
,
a
i
(
r
i
-
1
) satisfying the following equation:
(8)
L
A
i
1
w
r
i
v
i
(
w
)
=
a
i
0
v
i
(
w
)
+
a
i
1
L
A
i
1
w
v
i
(
w
)
+
⋯
+
a
i
(
r
i
-
1
)
L
A
i
1
w
r
i
-
1
v
i
(
w
)
,
where L
A
i
1
w
v
i
(
w
)
=
(
∂
v
i
(
w
)
/
∂
w
)
A
i
1
w, L
A
i
1
w
m
v
i
(
w
)
=
(
∂
L
A
i
1
w
m
-
1
v
i
(
w
)
/
∂
w
)
A
i
1
w, m
=
2,3
,
…
,
r
i, and L is an operator of Lee derivative.Apparently, ifv
i
(
w
) denotes a linear expression of w, Assumption 4 will hold automatically.Assumption 5.
There exists a smooth functionh
i
(
x
)
:
R
n
→
R
n
q making
(9)
∂
h
i
(
x
)
∂
x
g
i
(
x
)
=
G
i
,
1
≤
i
≤
m
,
where G
i is a nonzero constant vector defined in R
n
q.The questions resolved in the paper can be depicted as below.Definition 6.
For an arbitrarily given compact subsetD
w
∈
R
n
q, a state feedback controller u
i can always be found to guarantee the solution of original multivariable closed-loop system (1) existence and boundness under arbitrary initial conditions for all w
(
0
)
∈
D
w, when t
≥
0, and even lim
⁡
t
→
∞
x
(
t
)
=
C, where C is a constant vector representing the reference value.


## 3. Design of Multiclass Nonlinear Internal Models

Let(10)
τ
i
(
w
)
=
col
(
v
i
(
w
)
,
L
A
i
1
w
v
i
(
w
)
,
…
,
L
A
i
1
w
r
i
-
1
v
i
(
w
)
)
;
then there exist matrices Φ
i
=
[
0
(
r
i
-
1
)
×
1
I
r
i
-
1
a
i
0
[
a
i
1
,
…
,
a
r
i
-
1
]
] and ψ
i
=
[
1,0
,
…
,
0
] making
(11)
∂
τ
i
(
w
)
∂
w
A
i
1
w
=
Φ
i
τ
i
(
w
)
,
v
i
(
w
)
=
ψ
i
τ
i
(
w
)
,
where matrix pair (
Φ
i
,
ψ
i
) is observable and I
r
i
-
1 is (
r
i
-
1
)
×
(
r
i
-
1
) unit matrix.For the sake of establishment of nonlinear internal model equation, the following assumption is brought up.Assumption 7.
There exists a matrixΦ
i
k,  k
=
2
,
…
,
n, making
(12)
∂
τ
i
(
w
)
∂
w
A
i
k
w
=
Φ
i
k
τ
i
(
w
)
.
Assume τ
^
i
(
w
)
=
T
i
τ
i
(
w
) and T
i
∈
R
r
i
×
r
i to be nonsingular matrices. With derivative of τ
^
i
(
w
) along with (2), we can obtain
(13)
∂
τ
^
i
(
w
)
∂
w
s
i
(
w
)
=
∂
τ
^
i
(
w
)
∂
w
(
A
i
1
w
+
∑
k
=
2
n
A
i
k
w
s
i
k
(
w
)
)
∂
τ
^
i
(
w
)
∂
w
s
i
(
w
)
=
T
i
Φ
i
T
i
-
1
τ
^
i
(
w
)
+
∑
k
=
2
n
T
i
Φ
i
k
s
i
k
(
w
)
T
i
-
1
τ
^
i
(
w
)
∂
τ
^
i
(
w
)
∂
w
s
i
(
w
)
=
T
i
ϕ
i
(
w
)
T
i
-
1
τ
^
i
(
w
)
,
w
w
w
w
v
i
(
w
)
=
ψ
i
T
i
-
1
τ
^
i
(
w
)
,
where ϕ
i
(
w
)
=
Φ
i
+
φ
i
(
w
),  φ
i
(
w
)
=
∑
k
=
2
n
Φ
i
k
s
i
k
(
w
), and φ
i
(
0
)
=
0,  ϕ
i
(
0
)
=
Φ
i.
In terms of linear observer theory, a Hurwitz matrixF
i is chosen to make matrix pair (
F
i
,
G
i
) be controllable for nonzero constant vector G
i defined in (9). Due to the observability of matrix pair (
Φ
i
,
ψ
i
) and the fact that F
i and Φ
i have nonintersecting frequency spectrum, hence Sylvester equation T
i
Φ
i
-
F
i
T
i
=
G
i
ψ
i has a unique nonsingular solution T
i.
Letq
i
=
ψ
i
T
i
-
1, the nonlinear exosystem can immerse into the system shown as below:
(14)
η
˙
i
=
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
)
η
i
+
G
i
q
i
η
i
;
v
i
(
w
)
=
q
i
η
i
.
In consequence, multiclass nonlinear internal models can be designed as follows:(15)
ξ
˙
i
=
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
)
(
ξ
i
+
h
i
(
x
)
)
-
∂
h
i
(
x
)
∂
x
(
f
i
(
x
)
+
g
i
(
x
)
u
i
)
.
Define an auxiliary errore
i as
(16)
e
i
=
η
i
-
ξ
i
-
h
i
(
x
)
,and with derivative of (16) along with (1), (14), and (15), we can get
(17)
e
˙
i
=
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
)
η
i
+
G
i
q
i
η
i
-
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
)
(
ξ
i
+
h
i
(
x
)
)
+
∂
h
i
(
x
)
∂
x
(
f
i
(
x
)
+
g
i
(
x
)
u
i
)
-
∂
h
i
(
x
)
∂
x
(
f
i
(
x
)
+
g
i
(
x
)
(
u
i
+
v
i
(
w
)
)
)
=
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
)
(
η
i
-
ξ
i
-
h
i
(
x
)
)
=
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
)
e
i
.


## 4. Design of Multivariable Nonlinear State Feedback Controller

On the basis of the established multiclass nonlinear internal models (15) and Assumption 1, the multivariable nonlinear state feedback controller can be designed as follows:
(18)
u
i
=
α
i
(
x
)
-
q
i
(
ξ
i
+
h
i
(
x
)
)
.Construct a Lyapunov functionW as below:
(19)
W
=
V
(
x
)
+
∑
i
=
1
m
γ
e
i
T
e
i
,
where γ is a positive real constant. With derivative of function W along with originally nonlinear system (1) and auxiliary error (17), we can obtain
(20)
W
˙
=
∂
V
(
x
)
∂
x
(
f
(
x
)
+
∑
i
=
1
m
g
i
(
x
)
(
u
i
+
v
i
(
w
)
)
)
+
∑
i
=
1
m
(
γ
e
i
T
(
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
+
F
i
T
+
T
i
-
T
φ
i
T
(
w
)
T
i
T
)
e
i
)
=
∂
V
(
x
)
∂
x
(
∑
i
=
1
m
f
i
(
x
)
w
w
w
w
i
+
∑
i
=
1
m
g
i
(
x
)
(
α
i
(
x
)
-
q
i
(
ξ
i
+
h
i
(
x
)
)
+
q
i
η
i
)
∑
i
=
1
m
)
+
∑
i
=
1
m
γ
e
i
T
(
S
(
w
)
)
e
i
=
∂
V
(
x
)
∂
x
(
f
i
(
x
)
+
∑
i
=
1
m
g
i
(
x
)
α
i
(
x
)
)
+
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
q
i
e
i
+
∑
i
=
1
m
γ
e
i
T
(
S
(
w
)
)
e
i
≤
-
d
0
(
∥
x
∥
)
+
|
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
|
∥
q
i
e
i
∥
+
∑
i
=
1
m
γ
e
i
T
(
S
(
w
)
)
e
i
,
where S
(
w
)
=
F
i
+
T
i
φ
i
(
w
)
T
i
-
1
+
F
i
T
+
T
i
-
T
φ
i
T
(
w
)
T
i
T. By application of inequality 2
a
b
≤
c
a
2
+
c
-
1
b
2 (choose c
=
1) into the second item of (20), we can get
(21)
|
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
|
∥
q
i
e
i
∥
≤
1
2
|
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
|
2
+
∑
i
=
1
m
1
2
∥
q
i
∥
2
∥
e
i
∥
2
≤
1
2
d
0
(
∥
x
∥
)
+
∑
i
=
1
m
1
2
∥
q
i
∥
2
∥
e
i
∥
2
.Assume that there exists a compact setS
-
(
w
) making S
(
w
) be a symmetric matrix of negative definiteness for all w
∈
S
-
(
w
). Hence, there exists a positive real number f satisfying the following inequality for all w
∈
S
-
(
w
):
(22)
e
i
T
(
S
(
w
)
)
e
i
≤
-
f
∥
e
i
∥
2
.Substitute (21) and (22) into (20); we can get
(23)
W
˙
≤
-
d
0
(
∥
x
∥
)
+
|
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
|
∥
q
i
e
i
∥
+
∑
i
=
1
m
γ
e
i
T
(
S
(
w
)
)
e
i
≤
-
d
0
(
∥
x
∥
)
+
1
2
d
0
(
∥
x
∥
)
+
∑
i
=
1
m
1
2
∥
q
i
∥
2
∥
e
i
∥
2
-
∑
i
=
1
m
γ
f
∥
e
i
∥
2
≤
-
1
2
d
0
(
∥
x
∥
)
+
∑
i
=
1
m
(
1
2
∥
q
i
∥
2
-
γ
f
)
∥
e
i
∥
2
.
Choose appropriate γ and f to satisfy
(24)
θ
=
1
2
∥
q
i
∥
2
-
γ
f
<
0
;
then, we can obtain
(25)
W
˙
≤
-
1
2
d
0
(
∥
x
∥
)
+
∑
i
=
1
m
θ
∥
e
i
∥
2
<
0
.
Hence, all variables are bounded. Combined with the utilization of invariant set theorem, we can get lim
⁡
t
→
∞
x
(
t
)
=
C and lim
⁡
t
→
∞
e
i
=
0. The following theorem can be obtained.Theorem 8.
For the multivariable nonlinear system (1) and nonlinear exosystem (2) satisfying Assumption 1 to Assumption 7, multiclass nonlinear internal models (15) and multivariable nonlinear state feedback controller (18) are able to make the closed-loop system globally bounded, and lim
⁡
t
→
∞
x
(
t
)
=
C.


## 5. Application to Speed Control of a Two-Input Nonlinear System

Speed control of a permanent magnet synchronous motor drive system is widely utilized in various industrial fields. Hence, consider a two-input variable motor drive system shown as below [25], which is not able to tackle with a single-input method of previous disturbance rejection:
(26)
d
i
d
d
t
=
-
R
s
L
d
i
d
+
p
i
q
ω
m
+
1
L
d
u
d
,
d
ω
m
d
t
=
p
φ
f
J
m
i
q
-
B
m
J
m
ω
m
-
1
J
m
T
L
,
d
i
q
d
t
=
-
R
s
L
q
i
q
-
p
i
d
ω
m
-
p
φ
f
L
q
ω
m
+
1
L
q
u
q
,
where i
d and i
q represent the current in d-axis and q-axis, respectively, u
d and u
q indicate the voltage in d-axis and q-axis, respectively, R
s and φ
f denote the resistance of stator and flux linkage of permanent magnet, ω
m is the speed, p and J
m denote number of pole pairs and moment of inertia of rotor, respectively, and T
L is an increasing load.With conversion of (26) into an affine form as (1), we can get
(27)
x
˙
=
f
(
x
)
+
∑
i
=
1
2
g
i
(
x
)
(
u
i
+
v
i
(
w
)
)
,
where x
=
[
x
1
x
2
x
3
]
T
=
[
i
d
ω
m
i
q
]
T,
(28)
f
(
x
)
=
[
-
R
s
L
d
x
1
+
p
x
2
x
3
p
φ
f
J
m
x
3
-
B
m
J
m
x
2
-
1
J
m
T
L
-
R
s
L
q
x
3
-
p
x
2
x
1
-
p
φ
f
L
q
x
2
]
,
g
1
(
x
)
=
[
1
L
d
0
]
T
,
g
2
(
x
)
=
[
1
L
q
]
,
and control inputs u
=
[
u
1
u
2
]
T
=
[
u
d
u
q
]
T.For the sake of understandability, the external disturbances inputs,v
1 and v
2, which denote different nonlinear nonharmonic disturbance input signals, are both produced by Van der Pol oscillator shown in (6), and let a
=
b
=
1. Hence, A
11
=
[
0
1
-
1
1
], A
12
=
[
0
0
0
-
1
]. Evidently, A
21
=
A
11 and A
22
=
A
12. It is important to note that A
11
,
A
21
,
A
12, and A
22 are not elements in matrix A
1 or A
2, which are all independent matrices defined in (2). For the sake of convenience, let A
11
=
A
21
=
A
1, A
12
=
A
22
=
A
2, and s
12
(
w
)
=
s
22
(
w
)
=
w
1
2. The bounded limit cycle will be generated by the oscillator; hence, Assumption 2 holds.Letc
1, c
2, and c
3 all be certain positive constants; the control law of disturbance-free system for (27) can be obtained by means of backstepping control and can be shown as follows:
(29)
α
(
x
)
=
[
α
1
(
x
)
α
2
(
x
)
]
=
[
-
L
d
(
c
1
x
1
+
L
q
L
d
p
x
2
x
3
)
L
q
(
-
p
φ
f
J
m
c
2
(
x
2
-
ω
ref
)
+
R
s
L
q
x
3
+
L
d
L
q
p
x
2
x
1
+
p
φ
f
L
q
x
2
-
c
3
(
x
3
-
J
m
p
φ
f
(
φ
f
J
m
ω
ref
+
1
J
m
T
L
)
)
)
]
,
where ω
ref is the reference speed. It can be verified that disturbance-free system (3) can be stabilized by α
(
x
). Owing to its unimportance for the study of the paper, the derivation process of the controller for the disturbance-free system is omitted.Let(30)
V
(
x
)
=
1
2
x
1
2
+
1
2
c
2
(
x
2
-
ω
ref
)
2
+
1
2
(
x
3
-
J
m
p
φ
f
(
φ
f
J
m
ω
ref
+
1
J
m
T
L
)
)
2
.By means of some calculations and simplifications, we can get(31)
∂
V
(
x
)
∂
x
(
f
(
x
)
+
∑
i
=
1
2
g
i
(
x
)
α
i
)
=
-
(
c
1
+
R
s
L
d
)
x
1
2
-
φ
f
J
m
c
2
(
x
2
-
ω
ref
)
2
-
c
3
(
x
3
-
J
m
p
φ
f
(
φ
f
J
m
ω
ref
+
1
J
m
T
L
)
)
2
,
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
=
1
L
d
x
1
+
1
L
q
(
x
3
-
J
m
p
φ
f
(
φ
f
J
m
ω
ref
+
1
J
m
T
L
)
)
.To be more specific, the relevant parameters of the motor drive system are depicted as follows: the rated torque of motorT
e
=
5.0 N·m, number of rotor pole pairs p
=
4, flux linkage of permanent magnet φ
f
=
0.18 Wb, resistance of stator R
s
=
1.95 Ω, inductances of stator in d-axis and q-axis L
d
=
L
q
=
0.0115 H, moment of inertia of rotor J
m
=
0.008 kg·m2, and damping coefficient of motor B
m
=
0.01
N/rad/s.In addition, we assumex
′
=
[
x
1
,
(
x
2
-
ω
ref
)
,
(
x
3
-
(
J
m
/
p
φ
f
)
(
(
φ
f
/
J
m
)
ω
ref
+
(
1
/
J
m
)
T
L
)
)
]
T; by application of (30) and (31) and selection of c
1
=
8000
,
c
2
=
40, and c
3
=
8000, we can get
(32)
1
2
∥
x
′
∥
2
≤
V
(
x
)
≤
20
∥
x
′
∥
2
,
∂
V
(
x
)
∂
x
(
f
(
x
)
+
∑
i
=
1
2
g
i
(
x
)
α
i
)
≤
-
7562
∥
x
′
∥
2
,
|
∂
V
(
x
)
∂
x
∑
i
=
1
m
g
i
(
x
)
|
2
≤
7562
∥
x
′
∥
2
.
Hence, Assumption 1 is satisfied.Let(33)
h
1
(
x
)
=
[
0
3
L
d
x
1
]
T
,
h
2
(
x
)
=
[
0
3
L
q
x
3
]
T
.
Hence,
(34)
G
1
=
∂
h
1
(
x
)
∂
x
g
1
(
x
)
=
[
0
3
]
T
,
G
2
=
∂
h
2
(
x
)
∂
x
g
2
(
x
)
=
[
0
3
]
T
.
Assumption 5 holds.Assume thatv
1 and v
2 denote different nonlinear disturbance input signals and act on the currents of d-axis and q
-axis, respectively. Choose v
1
=
w
1 and v
2
=
w
2. In other words, v
1 and v
2  represent that the system are immersed into two differently external disturbances.Forv
1
=
w
1, by means of some mathematical calculations, it can be obtained that
(35)
L
A
1
w
v
1
(
w
)
=
∂
v
1
(
w
)
∂
w
A
1
w
=
[
1
0
]
[
0
1
-
1
1
]
[
w
1
w
2
]
=
w
2
,
L
A
1
w
2
v
1
(
w
)
=
∂
L
A
1
w
v
1
(
w
)
∂
w
A
1
w
=
[
0
1
]
[
0
1
-
1
1
]
[
w
1
w
2
]
=
-
w
1
+
w
2
=
-
a
10
v
1
(
w
)
+
a
11
L
A
1
w
v
1
(
w
)
.
Hence, there exist constants r
1
=
2
,
a
10
=
-
1, and a
11
=
1 satisfying Assumption 4.Forv
2
=
w
2, by means of similar mathematical calculations, we can get
(36)
L
A
1
w
v
2
(
w
)
=
∂
v
2
(
w
)
∂
w
A
1
w
=
[
0
1
]
[
0
1
-
1
1
]
[
w
1
w
2
]
=
-
w
1
+
w
2
,
L
A
1
w
2
v
2
(
w
)
=
∂
L
A
1
w
v
2
(
w
)
∂
w
A
1
w
=
[
-
1
1
]
[
0
1
-
1
1
]
[
w
1
w
2
]
=
-
w
1
=
-
a
20
v
1
(
w
)
+
a
21
L
A
1
w
v
1
(
w
)
.
Therefore, there exist constants r
1
=
2
,
a
20
=
1, and a
21
=
1 making Assumption 4 be satisfied.In terms of (35), let
(37)
τ
1
(
w
)
=
col
(
v
1
(
w
)
,
L
A
1
w
v
1
(
w
)
)
=
col
(
w
1
,
w
2
)
.Then, there exist matricesΦ
1
=
[
0
(
r
1
-
1
)
×
1
I
r
1
-
1
a
10
[
a
11
,
…
,
a
r
1
-
1
]
]
=
[
0
1
-
1
1
] and ψ
1
=
[
1,0
] making (11) hold and rank
⁡
[
ψ
1
ψ
1
Φ
1
]
T
=
2.According to (36), let
(38)
τ
2
(
w
)
=
col
(
v
2
(
w
)
,
L
A
1
w
v
2
(
w
)
)
=
col
(
w
2
,
-
w
1
+
w
2
)
.In that way, there exist matricesΦ
2
=
[
0
(
r
2
-
1
)
×
1
I
r
2
-
1
a
20
[
a
21
,
…
,
a
r
2
-
1
]
]
=
[
0
1
1
1
] and ψ
2
=
[
1,0
] making (11) satisfy and rank
⁡
[
ψ
2
ψ
2
Φ
2
]
T
=
2.ChooseΦ
12
=
A
2; after some calculations it can be gotten that
(39)
∂
τ
1
(
w
)
∂
w
A
2
w
=
Φ
12
τ
1
(
w
)
.SelectΦ
22
=
[
-
1
0
-
1
0
]; after some calculations we can get
(40)
∂
τ
2
(
w
)
∂
w
A
2
w
=
Φ
22
τ
2
(
w
)
.Consequently, Assumption7 is satisfied.Assumeτ
^
1
(
w
)
=
T
1
τ
1
(
w
), where T
1
∈
R
r
1
×
r
1 is a nonsingular matrix, and by derivation of τ
^
1
(
w
) along with (2), we can get
(41)
∂
τ
^
1
(
w
)
∂
w
s
1
(
w
)
=
T
1
ϕ
1
(
w
)
T
1
-
1
τ
^
1
(
w
)
,
where ϕ
1
(
w
)
=
Φ
1
+
φ
1
(
w
)
,
φ
1
(
w
)
=
Φ
12
s
12
(
w
).Supposeτ
^
2
(
w
)
=
T
2
τ
2
(
w
), where T
2
∈
R
r
2
×
r
2 is a nonsingular matrix and by derivation of τ
^
2
(
w
) along with (2), we can get
(42)
∂
τ
^
2
(
w
)
∂
w
s
2
(
w
)
=
T
2
ϕ
2
(
w
)
T
2
-
1
τ
^
2
(
w
)
,
where ϕ
2
(
w
)
=
Φ
2
+
φ
2
(
w
)
,
φ
2
(
w
)
=
Φ
22
s
22
(
w
).ChooseF
1
=
[
-
2
1
0
-
6
]; we can verify rank
⁡
[
G
1
G
1
F
1
]
T
=
2. Solving Sylvester equation T
1
Φ
1
-
F
1
T
1
=
G
1
ψ
1, it can be gotten that
(43)
T
1
=
[
0.1993
-
0.0897
0.4884
-
0.0698
]
.SelectF
2
=
[
-
2
1
0
-
6
]. Similarly, we are able to obtain rank
⁡
[
G
2
G
2
F
2
]
T
=
2. Solving Sylvester equation T
2
Φ
2
-
F
2
T
2
=
G
2
ψ
2, we can get
(44)
T
2
=
[
0.3220
-
0.1317
0.5122
-
0.0732
]
.It has already been verified that all assumed conditions required for systems (27) and (6) are satisfied. Following the design procedure given in Sections 3 and 4, the corresponding nonlinear internal models and state feedback controller are able to be designed as below:
(45)
ξ
˙
1
=
(
-
2
-
1.4651
w
1
2
)
ξ
1
+
(
1
+
0.598
w
1
2
)
(
ξ
2
+
0.0345
x
1
)
,
ξ
˙
2
=
-
1.1395
w
1
2
ξ
1
+
(
-
6
+
0.4651
w
1
2
)
(
ξ
2
+
0.0345
x
1
)
+
5.85
x
1
-
0.138
x
2
x
3
-
3
u
1
,
ξ
˙
3
=
(
-
2
+
0.3171
w
1
2
)
ξ
3
+
(
1
-
0.5707
w
1
2
)
(
ξ
4
+
0.0345
x
1
)
,
ξ
˙
4
=
0.7317
w
1
2
ξ
3
+
(
-
6
-
1.3171
w
1
2
)
(
ξ
4
+
0.0345
x
1
)
-
5.85
x
3
-
0.138
x
2
x
1
-
2.16
x
2
+
3
u
2
,
u
1
=
-
92.046
x
1
-
0.046
x
2
x
3
+
2.3333
ξ
1
-
3
ξ
2
,
u
2
=
-
413.28
x
2
-
90.096
x
3
+
249550
+
0.046
x
2
x
1
+
127.7778
T
L
+
1.6667
ξ
3
-
3
ξ
4
.In order to evaluate the validity of the proposed method, the numerical simulations are performed in Matlab software with the sampling intervalΔ
t
=
0.001 s and initial conditions x
(
0
)
=
[
0.1
0
1.0
]
,
η
^
(
0
)
=
[
0
0
0
0
], and w
(
0
)
=
[
1
-
1
]. In addition, the reference values of d-axis current i
d and speed ω
m are set to zero and 600 r/min, respectively, and q-axis current i
q tracks the linear variation value of load T
L. Ultimately, the simulation results are depicted in Figures 1, 2, and 3. Figure 1 describes the multiclass nonlinear nonharmonic disturbances v
1 and v
2, in which we can see that the exogenous disturbances immersing into the system are bounded and different, and they can originate from different exosystems so long as the exosystems can be written in the form of (2). Figure 2 demonstrates the control inputs u
d and u
q, which matches the external disturbances and reference values; and their combined actions make the different disturbances be rejected and ensure the system outputs to track the reference values. Figure 3 indicates the state outputs of the system, in which we can see that the d-axis current i
d and speed ω
m are able to track their own reference values in high precision. In other words, the simulation results signify that the proposed algorithm is able to reject the multiclass external disturbance effectively and ensures that the closed-loop system rapidly converge to the reference values. In consequence, the designed internal models can replicate with the different exogenous disturbances and the proposed controller has an excellent control performance.Figure 1
Multiclass nonlinear disturbance inputs.Figure 2
System control inputs.Figure 3
System state outputs.


## 6. Conclusions

In the paper, an IMP based multivariable nonlinear control algorithm with multiclass nonharmonic disturbances is proposed to suppress the differently external nonharmonic disturbances and control the closed-loop system to track the reference values. The major conclusions of the research are summarized as follows:(1)
in light of the multiclass nonharmonic disturbances, a class of new multiclass nonlinear internal models can be constructed on the basis of the definitions of expanded steady state generator and IMP, and the characteristics of which depend on the structural information of exosystem and state information of closed-loop system;
(2)
a multivariable nonlinear state feedback controller associated with the designed internal model equations is designed; the application of Lyapunov theory can justify the convergence of the proposed control algorithm;
(3)
the application of the algorithm into a two-input variable motor drive system demonstrates that the presented algorithm is able to suppress the multiclass disturbances and guarantee the multivariable closed-loop system global uniform convergence.



---

*Source: 570542-2014-06-19.xml*