# Robust Nonlinear Control with Compensation Operator for a Peltier System

Authors: Sheng-Jun Wen; Dong-Yun Wang; Lei Zhang; Ming-Cong Deng
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Robust nonlinear control with compensation operator is presented for a Peltier actuated system, where the compensation operator is designed by using a predictive model on heat radiation. For the Peltier system, the heat radiation is related to the fourth power of temperature. So, the heat radiation is affected evidently by the temperature when it is high and temperature difference between the system and environment is large. A new nonlinear model with the heat radiation is set up for the system according to some thermal conduction laws. To ensure robust stability of the nonlinear system, operator based robust right coprime factorization design is considered. Also, a compensation operator based on a predictive model is proposed to cancel effect of the heat radiation, where the predictive model is set up by using radial basis kernel function based SVM (support vector machine) method. Finally, simulation results are given to show the effectiveness of the proposed scheme.


---

## Body

## 1. Introduction

With the rapid development of semiconductor materials, semiconductor refrigeration has been paid more and more attention in refrigeration research and industrial fields. And because of some advantages, such as no refrigerant, no pollution, no vibration, and no noise, semiconductor refrigerator is gradually going from lab to practical engineering, such as national defense industry, agriculture, and medical treatment [1]. Semiconductor refrigeration is also called electronic refrigeration or thermoelectric refrigeration, which uses special semiconductor material with composition ofN-type semiconductor andP-type semiconductor. That is, it is a new type of refrigeration method through the input direct current flow, and it has been one of the three typical refrigeration ways in the world. In the semiconductor refrigeration, Peltier device is usually used as one of the refrigeration components and aluminum plate is considered as the thermal sink. Due to the Peltier effect, heat is absorbed at one end and pumped out through the other end when an external voltage is applied to aN-type/P-type thermal couple. Therefore, a Peltier actuated aluminum plate thermal system is often utilized to analyze the refrigeration performance of the Peltier device.In the Peltier actuated system, heat exchange mainly includes three forms: heat convection, heat conduction, and heat radiation. By modeling using some thermal exchange laws, it shows that the Peltier actuated system is a typical nonlinear affine system, where input and output are current and temperature, respectively. Not only the system temperature depends on input current and its square, but also there exists the fourth power of the temperature in the model. For the nonlinear system ignoring heat radiation, operator theory based nonlinear control method is confirmed to be effective to ensure robust stability of the controlled system [2–9]. Operator theory is a control theory based on an idea that a signal in the input space is mapped to the output space, which can be applied in linear system and nonlinear system [3]. So far, a lot of researchers have paid attention to the nonlinear control in regard to the theory, and some interesting results have been concluded [2–9]. An outstanding example is that robust nonlinear control based on robust right coprime factorization is presented for the Peltier actuated system in [4], where the plant and each controller are regarded as an operator, and signals are transferred from one space into another. Based on the design scheme, the realization of right coprime factorization controllers and design of robust tracking controller are discussed in [5]. Further, a control scheme is considered when the thermal system exist input constraints, where the designed controller is required to satisfy some robust right coprime factorization conditions [6]. Also, control and fault detection by using operator based nonlinear design scheme are investigated in [6–8]. Moreover, in [9, 10], operator based robust nonlinear control is designed for a multi-input multioutput thermal process with perturbation. However, the heat radiation will affect evidently the temperature of the Peltier actuated system when the temperature is high and temperature difference between the system and environment is large. A new nonlinear model with heat radiation needs to be built to design operator based robust right coprime factorization controllers. Also, in order to cancel the effect of the heat radiation, a compensation operator based on predictive model is investigated, where the model is considered by using radial basis kernel function based SVM method. SVM is a type of machine learning algorithm which is based on statistical learning theory and arises in recent years [11]. Comparing with SVM method and other machine learning methods, SVM can effectively improve generalization abilities of the algorithm, while not having a local minimum solution but a unique global optimal solution, and there are some characteristics, such as small sample learning ability and automatic identification of the model structure and the parameters [11–15]. It makes SVM very suitable for the modeling of the system with unknown structure and parameters according to the input and output data. Therefore, operator based robust nonlinear control with a compensation operator is put forward for the Peltier actuated system considering heat radiation in this paper.The outline of this paper is given as follows. Section2 begins with the introduction of the Peltier system, and modeling on the Peltier system and SVM based modeling on the heat radiation are also shown in this part. Design of operator based robust nonlinear controllers and compensation operator are proposed for the Peltier system in Section 3. Simulation results are given to illustrate the effectiveness of the modeling and the proposed scheme in Section 4. The conclusion is drawn in Section 5.


## 2. Modeling

### 2.1. Peltier Refrigeration System

The diagram of Peltier refrigeration system is shown in Figure1, which mainly includes Advantech PCL812PG board, PCI1760U board, aluminium plate, Peltier actuator, temperature sensor, and PC. LM35 is used to measure the temperature of the aluminium plate, and its output voltage range is 0~1 V. Peltier actuator is a kind of refrigeration device, which can produce cooling effect through changing duty cycle of PWM wave signal. PCL812PG board and PCI1760U board are inserted in PC, which are connected to the signal processing board by terminal boards. PCL812PG board with 16 A/D channels and 4 D/A channels is selected to get temperature from sensor, whose available voltage range is 0~5 V with 12-bit precision. PCI1760U board has 10 output channels with PWM wave, whose duty cycle is limited between 0% and 100%. Signal processing board is designed to amplify voltage signal of the temperature sensor from 0~1 V to 0~5 V and convert the control signal (5~12 V) into PWM wave with different duty cycle. Control algorithm is realized and control signal is produced by PC.Figure 1
The diagram of Peltier refrigeration system.In the refrigeration system, Peltier actuator is a semiconductor refrigeration device with composition ofN-type semiconductor andP-type semiconductor. The heat flux and direction are controlled by the input current flow. Four main thermoelectric effects govern the thermodynamics of a Peltier actuator: the Seebeck effect, the Peltier effect, the Thomson effect, and the Joule effect. Due to the Peltier effect, heat will be absorbed at one end and pumped out through the other end when an external voltage is applied to aN-type/P-type thermal couple. Thus, effectiveness of refrigeration can be realized by the Peltier actuator.


### 2.2. Modeling on the Peltier System

In order to set up the model of the Peltier system, some laws related to thermal exchange are utilized as follows [1].(1) Fourier’s law concerning thermal conduction is(1)
q
[
W
/
m
2
]
=
-
λ
[
W
/
mK
]
(
d
θ
d
n
)
[
K
/
m
]
,
where q is heat flux, λ is thermal conductivity, and d
θ
/
d
n is temperature gradient.(2) Newton’s cooling law is(2)
q
[
W
/
m
2
]
=
α
[
W
/
mK
]
(
Δ
T
)
[
K
]
,
where α is thermal conductivity (flowing air: 10~600 [
W
/
mK
]) and Δ
T is temperature difference.(3) Stefan-Boltzmann law is(3)
Φ
[
W
]
=
ε
σ
[
W
/
m
2
K
4
]
A
S
[
m
2
]
T
4
[
K
4
]
,
where Φ is radiation energy, ε is emissivity (ε
<
1), σ is radiation coefficient, A
S is area, and T is temperature.(4) Thermal conduction law is(4)
d
Q
[
J
]
=
c
[
J
/
kgK
]
m
[
kg
]
d
θ
[
K
]
,
where Q is Calorie heat, c is specific heat, and m is mass of the aluminum plate.(5) Electrothermal amount on Peltier effect is(5)
Q
[
J
]
=
S
C
[
J
/
KA
]
T
[
K
]
I
[
A
]
,
where S
C is Seebeck coefficient and I is current.(6) Joule exothermic heat on current is(6)
Q
[
J
]
=
R
P
[
Ω
]
I
2
[
A
2
]
,
where R
P is Peltier’s resistance.The aluminum plate thermal process with a Peltier device is shown in Figure2.Figure 2
The aluminum plate thermal process with a Peltier device.A differential equation based on the above laws is obtained for the Peltier system:(7)
c
m
d
(
T
0
-
T
x
)
d
t
=
u
d
-
α
(
T
0
-
T
x
)
(
2
S
4
+
2
S
5
-
S
6
)
-
2
λ
(
T
0
-
T
x
)
S
2
d
1
-
ε
σ
(
4
S
1
+
2
S
2
+
4
S
3
+
2
S
4
+
2
S
5
-
S
6
)
(
T
0
4
-
T
x
4
)
,
where T
0 is the environmental temperature, T
x is the temperature, and u
d denotes the heat absorbed by the Peltier device. Some parameters of the Peltier system are given in Table 1.Table 1
Parameters of the Peltier system.
I
Current: 0~2.2[
A
]
u
d
Endothermic amount[
W
]
c
900[
J/kg/K
]
α
15[
ω
/
m
2
/
K
]
λ
238[
ω
/
m
/
K
]
σ
5.6697
×
1
0
-
8  [
W
/
m
2
/
K
4
]
K
0.63[
W
/
K
]
m
d
3
×
d
4
×
d
5
×
d  [
kg
]
d
1
0.086[
m
]
d
2
0.068[
m
]
d
3
0.1[
m
]
d
4
0.24[
m
]
d
5
0.005[
m
]
d
2700[
kg
/
m
3
]
S
C
0.053[
V
/
K
]
S
1
d
1
×
d
3  [
m
2
]
S
2
d
3
×
d
5  [
m
2
]
S
3
d
1
×
d
5  [
m
2
]
S
4
d
2
×
d
3  [
m
2
]
S
5
d
2
×
d
5  [
m
2
]
S
6
9
×
1
0
-
4  [
m
2
]
R
P
5.5[
Ω
]In the model,α
(
T
0
-
T
x
)
(
2
S
4
+
S
5
-
S
6
) denotes the energy produced by thermal convection between the middle aluminum plate and the air. 2
λ
(
T
x
-
T
0
)
S
2
/
d
1 is the conducted energy from the other two parts to the middle part. And ε
σ
(
4
S
1
+
2
S
2
+
4
S
3
+
2
S
4
+
2
S
5
-
S
6
)
(
T
0
4
-
T
x
4
) is the energy produced by the thermal radiation. The endothermic amount of the Peltier device u
d is shown in
(8)
u
d
=
S
C
T
1
I
-
K
(
T
h
-
T
1
)
-
R
P
I
2
2
,
where T
1 is the temperature of endothermic side, T
h is the temperature of radiation side, and I is the current flow. S
C
T
1
I denotes the amount of heat that moves from the endothermic side to the radiation side by Peltier effect, K
(
T
h
-
T
1
) shows the movement heat caused by temperature difference between the two sides of the device, and R
P
I
2
/
2 is Joule heat which is generated by the input current through the device.


### 2.3. SVM Based Modeling on the Heat Radiation Part

In the above model, the heat radiation part contains the fourth power of the temperature, and it is not easy to be solved directly. Also, the heat radiation affects evidently the temperature of the Peltier actuated system when the temperature is high and temperature difference is large between the system and environment. In order to cancel the effect of the heat radiation, a compensation operator on the heat radiation is designed, which is based on a predictive model of the heat radiation. In this paper, the predictive model is set up by using SVM method and used to estimate the heat radiation. That is, SVM based predictive model on the heat radiation needs to be built.SVM is firstly put forward by Vapnik, like multilayer perception network and radial basis function network, which can be used for pattern classification and nonlinear regression [11]. The basic idea of SVM is to create a classification hyperplane as decision surface which makes the isolation edge between positive examples and counter examples maximized. Its regression can construct an identification model of the system according to input and output data under the condition of some unknown parameters. SVM theory is based on Vapnik-Chervonenkis’ statistical learning theory, where a structural risk minimization principle is utilized to get a unique global optimal solution. It can avoid local minimum problem and there are some good promotion characteristics, such as learning ability of small sample and automatic identification of structure and parameters of the model [12–15]. It is suitable for building the nonlinear predictive model of the heat radiation.The detailed modeling on the heat radiation is shown as follows. Firstly, we define(9)
A
=
[
α
(
2
S
4
+
2
S
5
-
S
6
)
+
2
λ
S
2
/
d
1
]
(
c
m
)
,
M
svm
=
ε
σ
(
4
S
1
+
2
S
2
+
4
S
3
+
2
S
4
+
2
S
5
-
S
6
)
(
T
0
4
-
T
x
4
)
(
c
m
)
.
Then (7) is transformed as the following form:
(10)
d
(
T
0
-
T
x
)
d
t
=
u
d
c
m
-
A
(
T
0
-
T
x
)
-
M
svm
.
It shows that the term M
svm denotes the heat radiation part, which contains the fourth power of the temperature. Therefore, SVM method is considered to set up the model of M
svm in this paper, and it is used to estimate the heat radiation real time. SVM concerning the estimation of a real-valued function can be described as follows:
(11)
f
(
x
~
)
=
w
T
·
ϕ
(
x
~
)
+
b
,
w
∈
R
d
,
b
∈
R
.The model is based on a finite number set of independent and distributed data(
x
~
i
,
y
~
i
), where x
i and y
i denote input and output, respectively, w is weight vector, and b is offset. In this paper, we make u
d and T
x as the input and M
svm as the output when we build the SVM modeling, and w and b are identification model parameters after training. In Vapnik’s ε-insensitive support vector regression, the aim is to find a function f
(
x
~
) which allows error of y
i to be no more than ε and makes y
i flat for all the training data. Considering more interferential errors, nonnegative slack variables ξ and ξ
* are introduced. Then, the optimization problem on the modeling can be described as the following form:
(12)
min
⁡
w
,
b
⁡
1
2
w
T
·
w
+
C
∑
i
=
1
l
(
ξ
i
+
ξ
i
*
)
s
.
t
.
{
y
~
i
-
(
w
T
·
ϕ
(
x
~
i
)
+
b
)
≤
ε
+
ξ
i
(
w
T
·
ϕ
(
x
~
i
)
+
b
)
-
y
~
i
≤
ε
+
ξ
i
*
ξ
i
,
ξ
i
*
≥
0
,
i
=
1
,
…
,
l
,
where C is a positive constant and used to control the punishment to the samples beyond the error ε. Generalizing to kernel function based regression estimation by introducing Lagrange multipliers, we can arrive at the following optimization problem [15]:
(13)
max
⁡
⁡
α
,
α
*
W
(
α
,
α
*
)
=
-
1
2
∑
i
=
1
l
∑
j
=
1
l
Q
i
j
(
α
i
-
α
i
*
)
(
α
i
-
α
j
*
)
max
⁡
⁡
α
,
α
*
W
(
α
,
α
*
)
=
+
∑
i
=
1
l
y
~
i
(
α
i
-
α
i
*
)
-
∑
i
=
1
l
ε
(
α
i
+
α
i
*
)
s
.
t
.
{
α
i
,
α
i
*
∈
[
0
,
C
]
,
i
=
1
,
…
,
l
∑
i
=
1
l
(
α
i
-
α
i
*
)
=
0
K
(
x
~
i
,
x
~
j
)
=
ϕ
(
x
~
i
)
·
ϕ
(
x
~
j
)
=
Q
i
j
.Ultimately, the regression estimation takes the form(14)
f
(
x
~
)
=
∑
i
=
1
l
(
α
i
-
α
i
*
)
(
x
~
i
·
x
~
)
+
b
=
∑
i
=
1
l
(
α
i
-
α
i
*
)
K
(
x
~
i
·
x
~
)
+
b
,
where K
(
x
~
i
,
x
~
j
)
=
ϕ
(
x
~
i
)
·
ϕ
(
x
~
j
) is kernel function that returns the dot product in feature space by using a nonlinear mapping ϕ when two vectors are given in input space. That is, linear regression in a high dimensional space corresponding to nonlinear regression in the low dimensional input space can be realized by mapping the input vectors into a feature space. That is, nonlinear regression problem is resolved by linear regression form.There are many kernel functions used in SVM, such as linear kernel function, polynomial function, radial basis function, and sigmoid function. Linear kernel can be looked as a special kind of form of RBF kernel function. The characteristic of sigmoid function is similar to RBF. In the radial basis function case, the support vector algorithm automatically determines centers, weights, and threshold such that radial basis kernel function is used in this paper, which is shown as follows:(15)
K
(
x
~
i
,
x
~
)
=
exp
⁡
⁡
(
-
∥
x
~
i
-
x
~
∥
2
2
σ
2
)
.Therefore, under the above case, some parameters of support vector algorithm can be automatically determined, such as weight vectorw, offset b, slack variables ξ
i,  ξ
i
*, and Lagrange multipliers α
i and α
i
*. Then, SVM model based estimation accuracy just depends on a setting of parameters C, ε, and σ. C is penalizing constant, which determines the tradeoff between the model complexity and fatness. ε is error-accuracy parameter, and it controls the width of the ε-insensitive zone for fitting the training data. σ is kernel width parameter, which is appropriately selected to reflect the input range of the training data. By selecting the appropriate parameters C, ε, and σ, the predictive model of the heat radiation part M
svm can be obtained and the heat radiation can be estimated in real time by the built predictive model. Then, operator based robust nonlinear control with compensation operator is presented for the Peltier system by the following section, where the compensation operator is proposed by using the SVM based predictive model.


## 2.1. Peltier Refrigeration System

The diagram of Peltier refrigeration system is shown in Figure1, which mainly includes Advantech PCL812PG board, PCI1760U board, aluminium plate, Peltier actuator, temperature sensor, and PC. LM35 is used to measure the temperature of the aluminium plate, and its output voltage range is 0~1 V. Peltier actuator is a kind of refrigeration device, which can produce cooling effect through changing duty cycle of PWM wave signal. PCL812PG board and PCI1760U board are inserted in PC, which are connected to the signal processing board by terminal boards. PCL812PG board with 16 A/D channels and 4 D/A channels is selected to get temperature from sensor, whose available voltage range is 0~5 V with 12-bit precision. PCI1760U board has 10 output channels with PWM wave, whose duty cycle is limited between 0% and 100%. Signal processing board is designed to amplify voltage signal of the temperature sensor from 0~1 V to 0~5 V and convert the control signal (5~12 V) into PWM wave with different duty cycle. Control algorithm is realized and control signal is produced by PC.Figure 1
The diagram of Peltier refrigeration system.In the refrigeration system, Peltier actuator is a semiconductor refrigeration device with composition ofN-type semiconductor andP-type semiconductor. The heat flux and direction are controlled by the input current flow. Four main thermoelectric effects govern the thermodynamics of a Peltier actuator: the Seebeck effect, the Peltier effect, the Thomson effect, and the Joule effect. Due to the Peltier effect, heat will be absorbed at one end and pumped out through the other end when an external voltage is applied to aN-type/P-type thermal couple. Thus, effectiveness of refrigeration can be realized by the Peltier actuator.


## 2.2. Modeling on the Peltier System

In order to set up the model of the Peltier system, some laws related to thermal exchange are utilized as follows [1].(1) Fourier’s law concerning thermal conduction is(1)
q
[
W
/
m
2
]
=
-
λ
[
W
/
mK
]
(
d
θ
d
n
)
[
K
/
m
]
,
where q is heat flux, λ is thermal conductivity, and d
θ
/
d
n is temperature gradient.(2) Newton’s cooling law is(2)
q
[
W
/
m
2
]
=
α
[
W
/
mK
]
(
Δ
T
)
[
K
]
,
where α is thermal conductivity (flowing air: 10~600 [
W
/
mK
]) and Δ
T is temperature difference.(3) Stefan-Boltzmann law is(3)
Φ
[
W
]
=
ε
σ
[
W
/
m
2
K
4
]
A
S
[
m
2
]
T
4
[
K
4
]
,
where Φ is radiation energy, ε is emissivity (ε
<
1), σ is radiation coefficient, A
S is area, and T is temperature.(4) Thermal conduction law is(4)
d
Q
[
J
]
=
c
[
J
/
kgK
]
m
[
kg
]
d
θ
[
K
]
,
where Q is Calorie heat, c is specific heat, and m is mass of the aluminum plate.(5) Electrothermal amount on Peltier effect is(5)
Q
[
J
]
=
S
C
[
J
/
KA
]
T
[
K
]
I
[
A
]
,
where S
C is Seebeck coefficient and I is current.(6) Joule exothermic heat on current is(6)
Q
[
J
]
=
R
P
[
Ω
]
I
2
[
A
2
]
,
where R
P is Peltier’s resistance.The aluminum plate thermal process with a Peltier device is shown in Figure2.Figure 2
The aluminum plate thermal process with a Peltier device.A differential equation based on the above laws is obtained for the Peltier system:(7)
c
m
d
(
T
0
-
T
x
)
d
t
=
u
d
-
α
(
T
0
-
T
x
)
(
2
S
4
+
2
S
5
-
S
6
)
-
2
λ
(
T
0
-
T
x
)
S
2
d
1
-
ε
σ
(
4
S
1
+
2
S
2
+
4
S
3
+
2
S
4
+
2
S
5
-
S
6
)
(
T
0
4
-
T
x
4
)
,
where T
0 is the environmental temperature, T
x is the temperature, and u
d denotes the heat absorbed by the Peltier device. Some parameters of the Peltier system are given in Table 1.Table 1
Parameters of the Peltier system.
I
Current: 0~2.2[
A
]
u
d
Endothermic amount[
W
]
c
900[
J/kg/K
]
α
15[
ω
/
m
2
/
K
]
λ
238[
ω
/
m
/
K
]
σ
5.6697
×
1
0
-
8  [
W
/
m
2
/
K
4
]
K
0.63[
W
/
K
]
m
d
3
×
d
4
×
d
5
×
d  [
kg
]
d
1
0.086[
m
]
d
2
0.068[
m
]
d
3
0.1[
m
]
d
4
0.24[
m
]
d
5
0.005[
m
]
d
2700[
kg
/
m
3
]
S
C
0.053[
V
/
K
]
S
1
d
1
×
d
3  [
m
2
]
S
2
d
3
×
d
5  [
m
2
]
S
3
d
1
×
d
5  [
m
2
]
S
4
d
2
×
d
3  [
m
2
]
S
5
d
2
×
d
5  [
m
2
]
S
6
9
×
1
0
-
4  [
m
2
]
R
P
5.5[
Ω
]In the model,α
(
T
0
-
T
x
)
(
2
S
4
+
S
5
-
S
6
) denotes the energy produced by thermal convection between the middle aluminum plate and the air. 2
λ
(
T
x
-
T
0
)
S
2
/
d
1 is the conducted energy from the other two parts to the middle part. And ε
σ
(
4
S
1
+
2
S
2
+
4
S
3
+
2
S
4
+
2
S
5
-
S
6
)
(
T
0
4
-
T
x
4
) is the energy produced by the thermal radiation. The endothermic amount of the Peltier device u
d is shown in
(8)
u
d
=
S
C
T
1
I
-
K
(
T
h
-
T
1
)
-
R
P
I
2
2
,
where T
1 is the temperature of endothermic side, T
h is the temperature of radiation side, and I is the current flow. S
C
T
1
I denotes the amount of heat that moves from the endothermic side to the radiation side by Peltier effect, K
(
T
h
-
T
1
) shows the movement heat caused by temperature difference between the two sides of the device, and R
P
I
2
/
2 is Joule heat which is generated by the input current through the device.


## 2.3. SVM Based Modeling on the Heat Radiation Part

In the above model, the heat radiation part contains the fourth power of the temperature, and it is not easy to be solved directly. Also, the heat radiation affects evidently the temperature of the Peltier actuated system when the temperature is high and temperature difference is large between the system and environment. In order to cancel the effect of the heat radiation, a compensation operator on the heat radiation is designed, which is based on a predictive model of the heat radiation. In this paper, the predictive model is set up by using SVM method and used to estimate the heat radiation. That is, SVM based predictive model on the heat radiation needs to be built.SVM is firstly put forward by Vapnik, like multilayer perception network and radial basis function network, which can be used for pattern classification and nonlinear regression [11]. The basic idea of SVM is to create a classification hyperplane as decision surface which makes the isolation edge between positive examples and counter examples maximized. Its regression can construct an identification model of the system according to input and output data under the condition of some unknown parameters. SVM theory is based on Vapnik-Chervonenkis’ statistical learning theory, where a structural risk minimization principle is utilized to get a unique global optimal solution. It can avoid local minimum problem and there are some good promotion characteristics, such as learning ability of small sample and automatic identification of structure and parameters of the model [12–15]. It is suitable for building the nonlinear predictive model of the heat radiation.The detailed modeling on the heat radiation is shown as follows. Firstly, we define(9)
A
=
[
α
(
2
S
4
+
2
S
5
-
S
6
)
+
2
λ
S
2
/
d
1
]
(
c
m
)
,
M
svm
=
ε
σ
(
4
S
1
+
2
S
2
+
4
S
3
+
2
S
4
+
2
S
5
-
S
6
)
(
T
0
4
-
T
x
4
)
(
c
m
)
.
Then (7) is transformed as the following form:
(10)
d
(
T
0
-
T
x
)
d
t
=
u
d
c
m
-
A
(
T
0
-
T
x
)
-
M
svm
.
It shows that the term M
svm denotes the heat radiation part, which contains the fourth power of the temperature. Therefore, SVM method is considered to set up the model of M
svm in this paper, and it is used to estimate the heat radiation real time. SVM concerning the estimation of a real-valued function can be described as follows:
(11)
f
(
x
~
)
=
w
T
·
ϕ
(
x
~
)
+
b
,
w
∈
R
d
,
b
∈
R
.The model is based on a finite number set of independent and distributed data(
x
~
i
,
y
~
i
), where x
i and y
i denote input and output, respectively, w is weight vector, and b is offset. In this paper, we make u
d and T
x as the input and M
svm as the output when we build the SVM modeling, and w and b are identification model parameters after training. In Vapnik’s ε-insensitive support vector regression, the aim is to find a function f
(
x
~
) which allows error of y
i to be no more than ε and makes y
i flat for all the training data. Considering more interferential errors, nonnegative slack variables ξ and ξ
* are introduced. Then, the optimization problem on the modeling can be described as the following form:
(12)
min
⁡
w
,
b
⁡
1
2
w
T
·
w
+
C
∑
i
=
1
l
(
ξ
i
+
ξ
i
*
)
s
.
t
.
{
y
~
i
-
(
w
T
·
ϕ
(
x
~
i
)
+
b
)
≤
ε
+
ξ
i
(
w
T
·
ϕ
(
x
~
i
)
+
b
)
-
y
~
i
≤
ε
+
ξ
i
*
ξ
i
,
ξ
i
*
≥
0
,
i
=
1
,
…
,
l
,
where C is a positive constant and used to control the punishment to the samples beyond the error ε. Generalizing to kernel function based regression estimation by introducing Lagrange multipliers, we can arrive at the following optimization problem [15]:
(13)
max
⁡
⁡
α
,
α
*
W
(
α
,
α
*
)
=
-
1
2
∑
i
=
1
l
∑
j
=
1
l
Q
i
j
(
α
i
-
α
i
*
)
(
α
i
-
α
j
*
)
max
⁡
⁡
α
,
α
*
W
(
α
,
α
*
)
=
+
∑
i
=
1
l
y
~
i
(
α
i
-
α
i
*
)
-
∑
i
=
1
l
ε
(
α
i
+
α
i
*
)
s
.
t
.
{
α
i
,
α
i
*
∈
[
0
,
C
]
,
i
=
1
,
…
,
l
∑
i
=
1
l
(
α
i
-
α
i
*
)
=
0
K
(
x
~
i
,
x
~
j
)
=
ϕ
(
x
~
i
)
·
ϕ
(
x
~
j
)
=
Q
i
j
.Ultimately, the regression estimation takes the form(14)
f
(
x
~
)
=
∑
i
=
1
l
(
α
i
-
α
i
*
)
(
x
~
i
·
x
~
)
+
b
=
∑
i
=
1
l
(
α
i
-
α
i
*
)
K
(
x
~
i
·
x
~
)
+
b
,
where K
(
x
~
i
,
x
~
j
)
=
ϕ
(
x
~
i
)
·
ϕ
(
x
~
j
) is kernel function that returns the dot product in feature space by using a nonlinear mapping ϕ when two vectors are given in input space. That is, linear regression in a high dimensional space corresponding to nonlinear regression in the low dimensional input space can be realized by mapping the input vectors into a feature space. That is, nonlinear regression problem is resolved by linear regression form.There are many kernel functions used in SVM, such as linear kernel function, polynomial function, radial basis function, and sigmoid function. Linear kernel can be looked as a special kind of form of RBF kernel function. The characteristic of sigmoid function is similar to RBF. In the radial basis function case, the support vector algorithm automatically determines centers, weights, and threshold such that radial basis kernel function is used in this paper, which is shown as follows:(15)
K
(
x
~
i
,
x
~
)
=
exp
⁡
⁡
(
-
∥
x
~
i
-
x
~
∥
2
2
σ
2
)
.Therefore, under the above case, some parameters of support vector algorithm can be automatically determined, such as weight vectorw, offset b, slack variables ξ
i,  ξ
i
*, and Lagrange multipliers α
i and α
i
*. Then, SVM model based estimation accuracy just depends on a setting of parameters C, ε, and σ. C is penalizing constant, which determines the tradeoff between the model complexity and fatness. ε is error-accuracy parameter, and it controls the width of the ε-insensitive zone for fitting the training data. σ is kernel width parameter, which is appropriately selected to reflect the input range of the training data. By selecting the appropriate parameters C, ε, and σ, the predictive model of the heat radiation part M
svm can be obtained and the heat radiation can be estimated in real time by the built predictive model. Then, operator based robust nonlinear control with compensation operator is presented for the Peltier system by the following section, where the compensation operator is proposed by using the SVM based predictive model.


## 3. Robust Nonlinear Control with Compensation Operator

### 3.1. Operator Based Robust Right Coprime Factorization

Robust right coprime factorization of a nonlinear system is to ensure that the system can be right coprime factorized and the system remains its right coprime factorization while it is affected by external perturbation. More specifically, we suppose that the nonlinear systemP has right factorization P
=
N
D
-
1, where N and D
-
1 are right factorization operators and the operation between N and D
-
1 means that the output of D
-
1 is the input of N. Then, if there are two operators S and R satisfying Bezout identity S
N
+
R
D
=
I, we can say that right coprime factorization factors of the system exist and the factorization is right coprime. Assume that the system P is disturbed by a perturbation Δ
P. If P
+
Δ
P has a right coprime factorization and satisfies some conditions, it shows that the system P has robust right coprime factorization shown in Figure 3.Figure 3
Robust right coprime factorization.The systemP without disturbance has the decomposition P
=
N
D
-
1. We define a disturbance operator as follows while P is disturbed by Δ
P:
(16)
Δ
N
=
Δ
P
D
.So, for the stable operatorN
+
Δ
N,  P
+
Δ
P can be factorized into
(17)
P
+
Δ
P
=
(
N
+
Δ
N
)
D
-
1
.Therefore, the systemP has the robust right coprime factorizations when the following conditions of robust stability are satisfied [3]:
(18)
S
N
+
R
D
=
L
,
where
L
is
unimodular
operator
,
∥
(
S
(
N
+
Δ
N
)
-
S
N
)
L
-
1
∥
<
1
.Based on the above robust right coprime factorization technology, a robust nonlinear control method with compensation operator is proposed for the Peltier actuated system.


### 3.2. Robust Right Coprime Factorization with Compensation Operator

For the Peltier system shown in (10), we define y
(
t
)
=
T
0
-
T
x. The model can be reexpressed as
(19)
d
y
d
t
=
u
d
c
m
-
A
y
-
M
svm
,
where the term M
svm denotes the heat radiation part which can be estimated in real time by using the SVM based predictive model. In a word, the M
svm can be regarded as a predictive operator of the heat radiation part.Considering the heat radiation as the perturbation, the Peltier actuated system can be expressed by (20):
(20)
y
=
1
c
m
(
e
-
A
t
+
Δ
)
∫
u
d
(
t
)
e
A
t
d
t
.In order to cancel the effect of the heat radiation, a compensation operatorP
svm is designed according to the predictive operator M
svm and the transformer operator T. Robust right coprime factorization with the compensation operator is shown in Figure 4 for the Peltier system.Figure 4
Robust nonlinear control scheme with compensation operator.In Figure4, U,  Y, and W are the input space, output space, and quasistate space, respectively. The operator M
svm is designed by using the SVM based predictive model and T is designed to transform the signal y
svm from output space Y to input space U.Theorem 1.
Assume thatΔ
N
~
=
S
-
1
T
M
(
D
,
N
) is stable, where M
(
D
,
N
) denotes M
svm
(
D
(
ω
)
,
N
(
ω
)
). Let S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N be stable and invertible operator (
W
→
U
). The robust stability of the Peltier system considering heat radiation is ensured if the following equations and inequality are satisfied:
(21)
S
N
+
R
D
=
I
,
(22)
S
(
N
+
Δ
N
-
Δ
N
~
)
+
R
D
=
L
~
,
(23)
∥
(
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
)
∥
<
1
.Proof.
Under the assumptionΔ
N
~
=
S
-
1
T
M
(
D
,
N
), the equivalent diagram of Figure 4 is obtained as shown in Figure 5. From (21) and (22), we have
(24)
L
~
=
I
+
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
.
SinceS
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N is stable operator and ∥
(
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
)
∥
<
1, then L
~ is invertible based on [3]. That is, L
~
-
1
=
(
I
+
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
)
-
1. Since S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N is stable and invertible operator, and I is identity operator, we have that L
~ is also a unimodular operator. Therefore, the stability of operators S, N, R, D, Δ
N, and Δ
N
~ can ensure that the outputs are bounded when the input signals are bounded in Figure 4.
So, the overall stability of the system is proved.
The nominal Peltier system is right divided into two partsN and D satisfying P
=
N
D
-
1:
(25)
D
-
1
(
u
d
)
=
u
d
(
t
)
c
m
,
(26)
N
(
ω
)
(
t
)
=
e
-
A
t
∫
e
A
t
ω
(
τ
)
d
τ
.
Considering the perturbationΔ
P caused by the heat radiation, the corresponding perturbed operator N
+
Δ
N can be obtained as
(27)
(
N
+
Δ
N
)
(
ω
)
(
t
)
=
(
e
-
A
t
+
Δ
)
∫
e
A
t
ω
(
τ
)
d
τ
which is bounded and stable. The uncertain term Δ is related to the heat radiation.
For the nonlinear system, nonlinear operator controllersS and R are designed to satisfy the Bezout identity (21) and inequality (23). It is worth to mention that the initial state should also be considered; that is, S
N
(
ω
0
)
(
t
0
)
+
R
D
(
ω
0
)
(
t
0
)
=
I
(
ω
0
)
(
t
0
) should be satisfied. In this paper, we select the initial time t
0
=
0 and ω
0
=
ω
0
(
t
0
). By considering the structure of D, the controller R is designed as
(28)
R
(
u
d
)
=
B
u
d
(
t
)
c
m
,
where B is a constant. Then we can get
(29)
R
D
(
ω
)
=
R
(
D
)
=
R
(
c
m
ω
(
t
)
)
=
B
ω
(
t
)
which shows that the expression of the R
D is simple. After that, it will be easy to get the expression of controller S. From Figure 4, we know that
(30)
S
N
=
S
(
N
)
=
(
1
-
B
)
ω
(
t
)
.
Based on N
(
ω
)
(
t
)
=
e
-
A
t
∫
e
A
t
ω
(
τ
)
d
τ, we can get the following operator:
(31)
S
(
y
)
=
(
1
-
B
)
(
A
y
(
t
)
+
d
y
(
t
)
d
t
)
.
The predictive operatorM
svm is realized by using the SVM based predictive model introduced in Section 2. The transformer operator T is
(32)
T
(
y
svm
)
=
(
1
-
B
)
c
m
y
svm
(
t
)
.
Furthermore, the tracking filterM is designed to make the system output y track the reference input r shown in Figure 6. According to the Bezout identity (21), Figure 6 is equivalent to Figure 7. Using the result in [5], the tracking filter M is designed as the following operator:
(33)
M
(
r
)
(
t
)
=
(
N
+
Δ
N
)
-
1
(
r
)
(
t
)
.
In the following section, the proposed control scheme is confirmed by simulation.Figure 5
The equivalent diagram of Figure4.Figure 6
Robust tracking control scheme.Figure 7
The equivalent diagram of Figure6.


## 3.1. Operator Based Robust Right Coprime Factorization

Robust right coprime factorization of a nonlinear system is to ensure that the system can be right coprime factorized and the system remains its right coprime factorization while it is affected by external perturbation. More specifically, we suppose that the nonlinear systemP has right factorization P
=
N
D
-
1, where N and D
-
1 are right factorization operators and the operation between N and D
-
1 means that the output of D
-
1 is the input of N. Then, if there are two operators S and R satisfying Bezout identity S
N
+
R
D
=
I, we can say that right coprime factorization factors of the system exist and the factorization is right coprime. Assume that the system P is disturbed by a perturbation Δ
P. If P
+
Δ
P has a right coprime factorization and satisfies some conditions, it shows that the system P has robust right coprime factorization shown in Figure 3.Figure 3
Robust right coprime factorization.The systemP without disturbance has the decomposition P
=
N
D
-
1. We define a disturbance operator as follows while P is disturbed by Δ
P:
(16)
Δ
N
=
Δ
P
D
.So, for the stable operatorN
+
Δ
N,  P
+
Δ
P can be factorized into
(17)
P
+
Δ
P
=
(
N
+
Δ
N
)
D
-
1
.Therefore, the systemP has the robust right coprime factorizations when the following conditions of robust stability are satisfied [3]:
(18)
S
N
+
R
D
=
L
,
where
L
is
unimodular
operator
,
∥
(
S
(
N
+
Δ
N
)
-
S
N
)
L
-
1
∥
<
1
.Based on the above robust right coprime factorization technology, a robust nonlinear control method with compensation operator is proposed for the Peltier actuated system.


## 3.2. Robust Right Coprime Factorization with Compensation Operator

For the Peltier system shown in (10), we define y
(
t
)
=
T
0
-
T
x. The model can be reexpressed as
(19)
d
y
d
t
=
u
d
c
m
-
A
y
-
M
svm
,
where the term M
svm denotes the heat radiation part which can be estimated in real time by using the SVM based predictive model. In a word, the M
svm can be regarded as a predictive operator of the heat radiation part.Considering the heat radiation as the perturbation, the Peltier actuated system can be expressed by (20):
(20)
y
=
1
c
m
(
e
-
A
t
+
Δ
)
∫
u
d
(
t
)
e
A
t
d
t
.In order to cancel the effect of the heat radiation, a compensation operatorP
svm is designed according to the predictive operator M
svm and the transformer operator T. Robust right coprime factorization with the compensation operator is shown in Figure 4 for the Peltier system.Figure 4
Robust nonlinear control scheme with compensation operator.In Figure4, U,  Y, and W are the input space, output space, and quasistate space, respectively. The operator M
svm is designed by using the SVM based predictive model and T is designed to transform the signal y
svm from output space Y to input space U.Theorem 1.
Assume thatΔ
N
~
=
S
-
1
T
M
(
D
,
N
) is stable, where M
(
D
,
N
) denotes M
svm
(
D
(
ω
)
,
N
(
ω
)
). Let S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N be stable and invertible operator (
W
→
U
). The robust stability of the Peltier system considering heat radiation is ensured if the following equations and inequality are satisfied:
(21)
S
N
+
R
D
=
I
,
(22)
S
(
N
+
Δ
N
-
Δ
N
~
)
+
R
D
=
L
~
,
(23)
∥
(
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
)
∥
<
1
.Proof.
Under the assumptionΔ
N
~
=
S
-
1
T
M
(
D
,
N
), the equivalent diagram of Figure 4 is obtained as shown in Figure 5. From (21) and (22), we have
(24)
L
~
=
I
+
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
.
SinceS
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N is stable operator and ∥
(
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
)
∥
<
1, then L
~ is invertible based on [3]. That is, L
~
-
1
=
(
I
+
S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N
)
-
1. Since S
(
N
+
Δ
N
-
Δ
N
~
)
-
S
N is stable and invertible operator, and I is identity operator, we have that L
~ is also a unimodular operator. Therefore, the stability of operators S, N, R, D, Δ
N, and Δ
N
~ can ensure that the outputs are bounded when the input signals are bounded in Figure 4.
So, the overall stability of the system is proved.
The nominal Peltier system is right divided into two partsN and D satisfying P
=
N
D
-
1:
(25)
D
-
1
(
u
d
)
=
u
d
(
t
)
c
m
,
(26)
N
(
ω
)
(
t
)
=
e
-
A
t
∫
e
A
t
ω
(
τ
)
d
τ
.
Considering the perturbationΔ
P caused by the heat radiation, the corresponding perturbed operator N
+
Δ
N can be obtained as
(27)
(
N
+
Δ
N
)
(
ω
)
(
t
)
=
(
e
-
A
t
+
Δ
)
∫
e
A
t
ω
(
τ
)
d
τ
which is bounded and stable. The uncertain term Δ is related to the heat radiation.
For the nonlinear system, nonlinear operator controllersS and R are designed to satisfy the Bezout identity (21) and inequality (23). It is worth to mention that the initial state should also be considered; that is, S
N
(
ω
0
)
(
t
0
)
+
R
D
(
ω
0
)
(
t
0
)
=
I
(
ω
0
)
(
t
0
) should be satisfied. In this paper, we select the initial time t
0
=
0 and ω
0
=
ω
0
(
t
0
). By considering the structure of D, the controller R is designed as
(28)
R
(
u
d
)
=
B
u
d
(
t
)
c
m
,
where B is a constant. Then we can get
(29)
R
D
(
ω
)
=
R
(
D
)
=
R
(
c
m
ω
(
t
)
)
=
B
ω
(
t
)
which shows that the expression of the R
D is simple. After that, it will be easy to get the expression of controller S. From Figure 4, we know that
(30)
S
N
=
S
(
N
)
=
(
1
-
B
)
ω
(
t
)
.
Based on N
(
ω
)
(
t
)
=
e
-
A
t
∫
e
A
t
ω
(
τ
)
d
τ, we can get the following operator:
(31)
S
(
y
)
=
(
1
-
B
)
(
A
y
(
t
)
+
d
y
(
t
)
d
t
)
.
The predictive operatorM
svm is realized by using the SVM based predictive model introduced in Section 2. The transformer operator T is
(32)
T
(
y
svm
)
=
(
1
-
B
)
c
m
y
svm
(
t
)
.
Furthermore, the tracking filterM is designed to make the system output y track the reference input r shown in Figure 6. According to the Bezout identity (21), Figure 6 is equivalent to Figure 7. Using the result in [5], the tracking filter M is designed as the following operator:
(33)
M
(
r
)
(
t
)
=
(
N
+
Δ
N
)
-
1
(
r
)
(
t
)
.
In the following section, the proposed control scheme is confirmed by simulation.Figure 5
The equivalent diagram of Figure4.Figure 6
Robust tracking control scheme.Figure 7
The equivalent diagram of Figure6.


## 4. Simulation Results

The system with Peltier device is a typical nonlinear control system considering the heat radiation. Because the heat radiation is related to the fourth power of temperature and is affected evidently when the temperature is high and temperature difference between the system and environment is large, a high initial temperature (50°C) is settled to show the effectiveness of the control scheme. The input current is limited between 0.0[A] and 2.2[A]. The parameters used in the simulation are shown in Table2. In order to show the simulation results clearly, we change the unit of temperature “K” into “°C” in the simulation.Table 2
Simulation parameters.
Reference input (desired temperature)
r
=
6  [°C]
Parameter
B
=
0.99
Current
I = 0.0~2.2 [A]
Initial temperature
50[°C]
Simulation time
600[sec]
Sampling time
100[msec]


### 4.1. Simulation on the Prediction Model

In this paper, the radial basis function is used to establish the predictive modelM
svm according to the SVM method, where some parameters of SVM are automatically determined, such as weight vector w, offset b, slack variables ξ
i, ξ
i
*, and Lagrange multipliers α
i, α
i
*. For the important design parameters C, ε, and σ, we set 2
σ
2 as 0.1, and C and ε are 100 and 0.1, respectively.According to the real data of the physical system, the predictive model of the heat radiation is obtained by using the SVM method. Figure8 shows the prediction results of the heat radiation using the built predictive model. The black line is the real value of the heat radiation and the red one is the predicting result of the heat radiation using the SVM based predictive model. Its standard deviation and maximum deviation are 4.3815
×
1
0
-
11 and 1.005
×
1
0
-
7, respectively.Figure 8
The prediction result of the heat radiation using SVM method.It is clear to see from Figure8 that the two curves are almost coincident. Moreover, the heat radiation part is estimated by using the least square method based predictive model shown in Figure 9. Its standard deviation and maximum deviation are 1.2905
×
1
0
-
8 and 3.9561
×
1
0
-
6, respectively. So it can be concluded that the SVM based predictive model of the heat radiation is more efficient and exact. The obtained model can be used to design the compensation operator.Figure 9
The prediction result of the heat radiation using least square method.


### 4.2. Simulation on Robust Nonlinear Control with Compensation Operator

Firstly, in order to show effect of the heat radiation in the refrigeration system, the simulations are carried out for the system considering the heat radiation or not, which is simulated by the model. The parameters used in the two simulations are the same and the controllers without compensation operator are designed. Figure10 shows the simulation result where the heat radiation has not been considered. And Figure 11 shows the simulation result where the heat radiation has been included.Simulation result without heat radiation.
(a)
(b)Simulation result with heat radiation.
(a)
(b)From Figure10, we can see that the output temperature of the system declined from 50°C to 44°C and finally falls in 44°C. But in Figure 11, the output temperature of the system cannot decline to the 44°C, and it just falls in 44.3°C at last. So it can be concluded that the heat radiation of the aluminum plate influences the output temperature of the system when the temperature is much higher than the environment temperature.In order to cancel the effect of the heat radiation, robust right coprime factorization control with compensation operator is designed to make the system reduce to 44°C. The simulation results are shown in Figure12, which is similar to Figure 10. It shows that the temperature of the aluminum plate considering the heat radiation can be fallen in 44°C through using the proposed controllers. Comparing with Figures 10 and 12, it also can be concluded that the designed compensation operator based on the predictive model of the heat radiation and the proposed scheme are efficient.Simulation results using the proposed scheme.
(a)
(b)Figure13 shows the result of robust stability analysis considering the heat radiation. In Figure 13, the ordinate axis denotes robust stability magnitude based on a time sequence, which is calculated according to the robust stability condition (23). In Figure 13, the result shows that the robust stability magnitude is less than 1. That is, the robust stability condition (23) is satisfied and the control system is robust stable.Figure 13
Result of robust stability analysis.


## 4.1. Simulation on the Prediction Model

In this paper, the radial basis function is used to establish the predictive modelM
svm according to the SVM method, where some parameters of SVM are automatically determined, such as weight vector w, offset b, slack variables ξ
i, ξ
i
*, and Lagrange multipliers α
i, α
i
*. For the important design parameters C, ε, and σ, we set 2
σ
2 as 0.1, and C and ε are 100 and 0.1, respectively.According to the real data of the physical system, the predictive model of the heat radiation is obtained by using the SVM method. Figure8 shows the prediction results of the heat radiation using the built predictive model. The black line is the real value of the heat radiation and the red one is the predicting result of the heat radiation using the SVM based predictive model. Its standard deviation and maximum deviation are 4.3815
×
1
0
-
11 and 1.005
×
1
0
-
7, respectively.Figure 8
The prediction result of the heat radiation using SVM method.It is clear to see from Figure8 that the two curves are almost coincident. Moreover, the heat radiation part is estimated by using the least square method based predictive model shown in Figure 9. Its standard deviation and maximum deviation are 1.2905
×
1
0
-
8 and 3.9561
×
1
0
-
6, respectively. So it can be concluded that the SVM based predictive model of the heat radiation is more efficient and exact. The obtained model can be used to design the compensation operator.Figure 9
The prediction result of the heat radiation using least square method.


## 4.2. Simulation on Robust Nonlinear Control with Compensation Operator

Firstly, in order to show effect of the heat radiation in the refrigeration system, the simulations are carried out for the system considering the heat radiation or not, which is simulated by the model. The parameters used in the two simulations are the same and the controllers without compensation operator are designed. Figure10 shows the simulation result where the heat radiation has not been considered. And Figure 11 shows the simulation result where the heat radiation has been included.Simulation result without heat radiation.
(a)
(b)Simulation result with heat radiation.
(a)
(b)From Figure10, we can see that the output temperature of the system declined from 50°C to 44°C and finally falls in 44°C. But in Figure 11, the output temperature of the system cannot decline to the 44°C, and it just falls in 44.3°C at last. So it can be concluded that the heat radiation of the aluminum plate influences the output temperature of the system when the temperature is much higher than the environment temperature.In order to cancel the effect of the heat radiation, robust right coprime factorization control with compensation operator is designed to make the system reduce to 44°C. The simulation results are shown in Figure12, which is similar to Figure 10. It shows that the temperature of the aluminum plate considering the heat radiation can be fallen in 44°C through using the proposed controllers. Comparing with Figures 10 and 12, it also can be concluded that the designed compensation operator based on the predictive model of the heat radiation and the proposed scheme are efficient.Simulation results using the proposed scheme.
(a)
(b)Figure13 shows the result of robust stability analysis considering the heat radiation. In Figure 13, the ordinate axis denotes robust stability magnitude based on a time sequence, which is calculated according to the robust stability condition (23). In Figure 13, the result shows that the robust stability magnitude is less than 1. That is, the robust stability condition (23) is satisfied and the control system is robust stable.Figure 13
Result of robust stability analysis.


## 5. Conclusions

For the Peltier actuated aluminum plate thermal system, it is certified that heat radiation influences evidently the temperature when the temperature is high and temperature difference is large between the system and environment. The model of the Peltier actuated system with heat radiation is built, and a predictive model on the heat radiation is obtained by using SVM method. Based on the predictive model, a compensation operator is designed to cancel the effect of the heat radiation. Furthermore, operator based robust nonlinear control with the compensation operator is presented to guarantee the robust stability and tracking performance for the Peltier actuated system. Finally, simulation results show the effectiveness of the proposed scheme.



---

*Source: 597236-2014-07-09.xml*