# Stabilization Using a Discrete Fuzzy PDC Control with  PID Controllers and Pole Placement: Application to an Experimental Greenhouse

Authors: Amine Chouchaine; Elyes Feki; Abdelkader Mami
Journal: Journal of Control Science and Engineering
 (2011)

---

## Abstract

This paper proposes a control strategy for complex and nonlinear systems, based on a parallel distributed compensation (PDC) controller. A solution is presented to solve a stability problem that arises when dealing with a Takagi-Sugeno discrete system with great numbers of rules. The PDC controller will use a classical controller like a PI, PID, or RST in each rule with a pole placement strategy to avoid causing instability. The fuzzy controller presented combines the multicontrol approach and the performance of the classical controllers to obtain a robust nonlinear control action that can also deal with time-variant systems. The presented method was applied to a small greenhouse to control its inside temperature by variation in ventilation rate inside the process. The results obtained will show the efficiency of the adopted method to control the nonlinear and complex systems.


---

## Body

## 1. Introduction

In the last few decades the Takagi-Sugeno (TS) fuzzy systems have become an important means for both modeling and control, and their performance in these two domains is proved in the research especially in the application of complex and nonlinear systems. TS fuzzy systems use a multimodel approach by fitting enough local linear models each describing an operating functional zone of the process [1]. The final output of the fuzzy system is the fuzzy weighted contribution of all the local output models guaranteeing precision and stability.Such a tool can be used for control purposes; in fact when dealing with complex systems classical methods become inefficient and so emerges the need for other control techniques. The PDC controllers [2] are a convenient solution, having the same structure as a TS fuzzy system, where every rule is associated with a local linear controller synthesized from the local model with the same rule [3]. The result is a nonlinear control action, which is a fuzzy blending of each individual linear controller. Many applications using TS model and PDC controller [3, 4] were presented, and several propositions of new techniques of PDC control appear in the literature [5]. Where the most proposed approaches present a state feedback for each rule, the stability is proved by finding a common Lyapunov function which can satisfy all the fuzzy subsystems [6]. However, in case of a large number of rules describing the system it is very difficult to apply this approach [7].In the continuous domain, a solution was presented using a PD or a PID controller for each rule. The denominator of the closed loop system is treated as an uncertain polynomial with affine linear uncertainty structure where the stability can be verified using a frequency domain criterion [7]. In the discrete domain, there is not an efficient solution to deal with the stability problem which is one of the objectives of this paper that will be discussed in Section 3. Another goal that will be described in the same section is to find the appropriate parameters of the controller for each rule using pole placement approach. But before that, in Section 2, an introduction to the TS fuzzy model and the PDC controller will be presented. In Section 4, a practical application will be presented where the proposed control strategy is used in an adaptive fuzzy control scheme. The process is an experimental greenhouse that was fuzzy identified and controlled in real time. Finally, Section 5 will conclude this paper.


## 2. Takagi-Sugeno Fuzzy Systems

There are different classes of fuzzy systems; the most often used are Mamdani fuzzy systems [8] and Takagi-Sugeno fuzzy systems [1]. The latter differ from the former in the rules consequents: it is not a fuzzy set but it is a local model of the system (submodel related to the rule that describes an operative zone of the process) to be approximated. For the jth rule a TS fuzzy system has the following form:(1)R(j): ifz1(k)isΩj1and⋯andzn(k)isΩjn,then{xj(k+1)=Ajx(k)+Bjuc(k)+Ljv(k),yj(k)=Cjxj(k),j=1,…,N.The fuzzy proposition “z is Ω” is the antecedent of the rule, the system of equation in the second part of (1) is the consequent, “z=(z1,…,zn)” are the inputs of the TS fuzzy system, they can be the states “x=(x1,…,xna),” the input “uc=(u1,…,una),” or the disturbances inputs “v=(v1,…,vnv),” “Ωji” is the membership function representing the fuzzy subset with a corresponding membership value “Ωji(zi),” and N is the total number of rules.Assume that all the subsystems considered are completely controllable and completely observable. Also, denote the following states and inputs variables:(2)x1(k)=y(k),x2(k)=y(k-1),…,xna(k)=y(k-na+1),u1(k)=u(k),u2(k)=u(k-1),…,una(k)=u(k-na+1).The matrix (Aj∈ℝna×na, Bj∈ℝnb×1, and Dj∈ℝnd×nd), represents the parameters of the TS fuzzy system (1) with the following Frobenius canonical structure:(3)Aj=(010⋯0⋮01⋱0⋮⋱⋱⋱⋮⋮⋯⋅⋯⋅aj0aj1⋯⋯ajna),Bj=(00⋯000⋯0⋮⋱⋱⋮00⋯0bj1bj1⋯bjna),Lj=(00⋯000⋯0⋮⋱⋱⋮00⋯0lj0lj1⋯ljnv),Cj=(10⋯0).The output levelxj of each local model is weighted by the firing strength μj(z(k))=∏i=1n(Ωji(zi(k))). The final output of the system is the weighted average of all the rule's outputs, the computed as (4)x(k+1)=∑j=1Nμj(z(k))yj(k+1)∑j=1Nμj(z(k)).Assuming that(5)μj(z(k))∑j=1Nμj(z(k))=βj(z(k)),
the output of the TS model becomes(6)x(k+1)=∑j=1Nβj(z(k))xj(k+1)
with 0<βj(z(k))<1 and ∑j=1Nβj(z(k))=1.On the other hand, the controller is a TS fuzzy system, having the same antecedent of the fuzzy model and differing from it in its consequent. A linear controller is developed for each rule, and the global control action is synthesized in the same way as the TS fuzzy model output [2] as follows:(7)u(k)=∑j=1Nβj(z(k))uj(k)
with uj representing the local output control.


## 3. Stability of a Closed-Loop TS System

The stability problem of (1) can be solved in general using the Lyapunov approach [9–11]. In this case, the proposed control law of the PDC controller is a state feedback having the following expression:(8)u=-∑j=1Nβj(z)Kjx.
Replacing (8) into (1) without considering the disturbance inputs leads to the following expression:(9)x(k+1)=∑j=1Nβj(z(k))[Ajx(k)-Bj∑i=1Nβi(z(k))Kix(k)]=∑j=1Nβj(z(k))[Aj∑i=1Nβi(z(k))x(k)-Bj∑i=1Nβi(z(k))Kix(k)]=∑j=1N∑i=1Nβi(z(k))βj(z(k))[Aj-BjKi]x(k),(10)x(k+1)=∑j=1Nβj(z(k))2Gjjx(k)+2∑j<iNβi(z(k))βj(z(k))Gji+Gij2x(k)
with Gji=Aj-BjKi.The following theorem gives a sufficient stability condition for a system of the form of (11) [12].Theorem 1.
The closed-loop TS fuzzy system of the form of (10) is quadratically stable for some state feedback Kj (via PDC scheme) if there exists a common positive definite matrix P such that
(11)GjjTPGjj-P<0forj=1,…,N,(12)(Gji+Gij2)TP(Gji+Gij2)-P<0forj<i.The problem arises when the TS fuzzy system has a great number of rules, for it is difficult to find matrixP that verifies the stability condition (see (11) and (12)) [7].The solution presented in this paper deals with a transfer function form of the TS fuzzy system. Combining (1) and (6) without considering the disturbance inputs, we obtain the following NARX form:(13)y(k)=∑j=1Nβj(z)(∑m=1naajmy(k-m)+∑l=1nabjlu(k-l)).
Since ∑j=1Nβj(z)=1, the following expression can be obtained:(14)∑j=1Nβj(z)[y(k)-∑m=1naajmy(k-m)]=∑j=1Nβj(z)[∑l=1nabjlu(k-l)].
The Z transformation leads to the new form:(15)∑j=1Nβj(z)[1-∑m=1naajmq-m]Y(q-1)=∑j=1Nβj(z)[∑l=1nabjlq-lU(q-1)].Let(16)Aj(q-1)=1-∑m=1naajmq-m,Bj(q-1)=∑m=1nabjmq-m,
and so the transfer function of the TS fuzzy system is as follows:(17)H(q-1)=Y(q-1)U(q-1)=∑j=1Nβj(z)Bj(q-1)∑j=1Nβj(z)Aj(q-1).
Let us also consider the transfer function of the PDC controller to be(18)P(q-1)=∑j=1Nβj(z)Dj(q-1)∑j=1Nβj(z)Cj(q-1),
where Dj(q-1)=∑m=0nadjmq-m and Cj(q-1)=∑m=0nacjmq-m are the numerator and the denominator of local controllers correspondent to each rule j of the PDC controller.The transfer function of the closed-loop fuzzy system will be(19)G(q-1)=(∑j=1NβjBj)(∑i=1NβiDi)(∑j=1NβjAj)(∑i=1NβiCi)+(∑j=1NβjBj)(∑i=1NβiDi).
This leads to the following expression:(20)G(q-1)=(∑j=1NβjBj)(∑i=1NβiDi)∑j=iNβj2(AjCj+BjDj)+∑j≠iNβjβi(AjCi+BjDi).From (20) we can see that the controller related to a rule j, initially designed to stabilize the local model related to the same rule, has an influence on the stability of the other submodel related to different rules as indicated by the expression ∑j≠iNβjβi(AjCi+BjDi) located in the denominator of the closed-loop transfer function. There is a high probability that the local controller for a particular rule is not suitable for a local model associated with a different rule. Thus, the stability of the overall closed-loop system with the current structure cannot be guaranteed especially when the TS fuzzy system has a great number of rules.In fact, the cross-coupling and possible instability are created by the two multiplications(∑j=1NβjAj)(∑i=1NβiCi) and (∑j=1NβjBj)(∑i=1NβiDi); if we manage to eliminate the multiplication, then the problem can be solved.In order to fulfill this objective the chosen local controllers must have the same order as the local models. For instance, if the system is a first-order one, then the local controllers must be a PI, and if the system is a second-order one, then the correspondent local controllers should be a PID. For higher-order systems, one can choose an RST controller to each local model with an equivalent order having the polynomialR(q-1) equal to the polynomial T(q-1).Also, assuming all stable zeros of the transfer function (17), the denominator Cj(q-1) of each local controller must include the numerator of the local model having the same rule j. So let us consider the denominator of a local controller related to a rule j:(21)Cj(q-1)=e(q-1)B̃j(q-1)withe(q-1)=1-q-1andB̃j(q-1)=∑m=1nabjmq-(m-1).Here all the local controllers share the same polynomiale(q-1)=1-q-1 representing the numerical expression of an integrator that will allow the rejection of disturbances. This expression already exists in a PI, a PID, or an RST controller having an integrator in the open loop [13].We have also(22)Bj(q-1)=q-1B̃j(q-1).
Replacing (21) and (22) in (19) leads to the following transfer function of the closed-loop fuzzy system: (23)G(q-1)=q-1(∑j=1Nβj∑m=0nadjmq-m)(1-q-1)(∑j=1Nβj(1-∑m=1naajmq-m))+(q-1∑j=1Nβj∑m=0nadjmq-m)=q-1∑j=1Nβj∑m=0nadjmq-m∑j=1Nβj(1+(dj0-1-aj1)q-1+∑m=1na(ajm-aj(m+1)+djm)q-(m+1)).In the transfer function (23) of the closed-loop system, the multiplication was eliminated and so the cause of the interaction or cross-coupling is removed. However, the stability is not yet guaranteed. We still need to find the appropriate numerator's parameters Dj(q-1) of each local controller. A pole placement strategy is adopted with an appropriate polynomial F(q-1)=1+∑m=1na+1lmq-m, having its roots inside the unit circle, to be identified with the denominator of every closed-loop subsystem as follows:(24)1+(dj0-1-aj1)q-1+∑m=1na(ajm-aj(m+1)+djm)q-(m+1)=1+∑m=1na+1lmq-m,∀j=1,…,N.
The result is the following system of equations:(25)dj0=l1+1+aj1,dj1=l2+aj2-aj1,⋮djna=lna+1+ana+1-ajna,∀j=1,…,N.
The denominator of the transfer function (23) becomes invariant whatever the rule j, and since we have ∑j=1Nβj=1, the new form of the transfer function (23) becomes(26)G(q-1)=∑j=1Nβj∑m=0nadjmq-(m+1)F(q-1),
with the poles located inside the unit circle which guarantees the stability of the system. Also, with the described method, whatever the number of rules is, it is possible to take advantage of both: the performance of a classical controller like a PI, a PID, or an RST and in the same time the multicontrol approach given by the PDC controller. The result is a robust nonlinear controller that can even deal with complex system with time-variant parameters when using a recursive fuzzy identification of the TS fuzzy model. A simulation of the proposed method will be presented in the next section to control such a system.


## 4. Application to an Experimental Greenhouse

The presented method was applied to control the inside temperature of a small greenhouse (Figure1). The inside climate is subjected to various phenomena like the solar radiation transfer, heat and mass transfer between the soil, the cover, the canopy and the outside climate. There is also the ventilation conducted by two three-phase motors that interfere to cool the inside temperature.Figure 1
Greenhouse process.The effect of the solar radiation is preponderant during daylight; in fact, the components of the greenhouse absorb the radiation energy and convert it to heat energy released in the air by heat transfer [14], which enhances the inside temperature that will surpass the outside temperature (Figures 2 and 3).Figure 2
Difference between inside and outside temperatures during  daylight.Figure 3
Solar radiation intercepted during daylight.All these factors result in a nonlinear evolution with time-varying parameters of the climate inside the greenhouse [15]. The use of ventilation to cool down the climate during  daylight is very common. In fact, with an adequate and skilful management it can be possible to maintain suitable climate-state variables (temperature, humidity, and CO2 concentration) [16]. Also, the ventilation represents an alternative and cheap solution to the air conditioner that consumes an important electrical power.But the complexity of the system makes it hard to achieve successful results, and so most of the farmers use an on-off control action between two boundary temperature values. But, this kind of control is unhealthy to the canopy. On the other hand, the time where the ventilation is effective is related to the weather. In fact, the possible constraints that oppose the use of this method all along the year are(i)
the presence of solar radiation that will enhance the inside temperature and make it surpass the outside temperature,(ii)
the desired temperature of the climate inside the greenhouse that must be greater with several degrees than the outside temperature. The choice of the desired temperature depends on the requirement of the cultivation inside the greenhouse.These conditions can be met in 3 to 4 months of the year. When the outside temperature is near the desired inside temperature, the farmers use the pad cooling system; it is a combination between fans and a wet pad that decreases the introduced air with several degrees [17–19]. Thus, the ventilation use is effective for another one or two months depending on the quality of the pad.In this paper, we focus on the evolution of the inside temperature by manipulating the ventilation rate to obtain a desired behavior of the process. As described in Section2, the TS fuzzy system is able to describe this kind of system and can be used to control the inside temperature using the fans during daylight.


### 4.1. Materials and Measurements

The process is a greenhouse with the following geometrical characteristics: 1.5 m height, total width of 1 m, and total length of 1.5 m.The inputs of the system are obtained from several sensors located inside and outside the greenhouse as follows.(i)
The inside and outside temperatures were measured by two LM35 transistors; both have an AD 620 amplifier to amplify the delivered signal.(ii)
The inside humidity was measured using a humidity sensor (SY-HS-230BT).(iii)
The solar radiation was obtained by a pyranometer type LP-PYRA 03 module.In order to control the ventilation rate, the greenhouse was equipped with two fans driven by two three-phase motors. These engines have a power supply delivered by a frequency converter (microdrive FC 51 Danfoss). The rotational speed of the fans is proportional to the frequency of the three-phase power supply and so the ventilation rate inside the greenhouse depends on the frequency of the power supply. Thus, the control input will be the frequency that will be computed by the fuzzy controller.The communication between MATLAB, the sensors, and the frequency converter is carried by the data acquisition module (KUSB-3100).


### 4.2. Fuzzy Identification

The structure of the discrete TS fuzzy system that will be approximated will have a local model for every rulej having the following form [20]:(27)R(j):Tj(k+1)=αj1(k)T(k)+αj2(k)u(k)+αj3(k)To(k)+αj4(k)Rs(k)+αj5(k)Rhi(k)+αj6(k)forj=1,…,N,
where u is the frequency delivered by the frequency converter to the motors equipped with fans to control the ventilation rate inside the greenhouse.The fuzzy identification is performed in two steps. The first one is the determination of the appropriate membership functions of each input to set the antecedent part of the TS fuzzy system, which is done using fuzzy clustering based on C-means algorithm. The membership functions were fixed andwould not change during the control application. For that purpose, on the day of 25-02-2011, we gave a random action control to collect the data needed under the sample time of 15 seconds; after that we set the membership functions in Figures 4, 5, 6, 7, and 8.Figure 4
Membership function of the inside temperature.Figure 5
Membership function of the delivered frequency.Figure 6
Membership function of the outside temperature.Figure 7
Membership function of the inside humidity.Figure 8
Membership function of the intercepted solar radiation.The second step is the estimation of (27) that represents the consequent part of the TS fuzzy system. For a time-variant system like a greenhouse, the more convenient way is an online fuzzy identification using the following recursive equations based on the recursive least square algorithm [20]:(28)Kj(k)=Pj(k-1)ze(k-1)λ/βj(z(k-1))+zeT(k-1)Pj(k-1)ze(k-1),θ̂j(k)=θ̂j(k-1)+Kj(k)(y(k)-zeT(k-1)θ̂j(k-1)),Pj(k)=1λ(Pj(k-1)-Kj(k)zeT(k-1)Pj(k-1)),
where λ is a forgetting factor, zeT=(TuToRsRhi1) is the regression vector of inputs, θ̂jT=(α̂j1α̂j2α̂j3α̂j4α̂j5α̂j6) is the estimated vector of parameters, and y is the output representing the inside temperature T.Pj(0) and θ̂j(0) are obtained from an offline fuzzy identification based on an ordinary weighted least square method [20], where the data used are those collected during the day before the recursive identification and also served to build the membership functions.The estimated vectorθ̂ will be used in an adaptive fuzzy control scheme (see Figure 9) to bring the inside temperature T to the desired temperature r using the fuzzy controller output u.Figure 9
Adaptive fuzzy control scheme.


### 4.3. Fuzzy Control

The system is a first-order one so the local controller for each rule will be a PI controller. From the previous section, the transfer function of the PDC controller can be deduced as follows:(29)P(q-1)=∑j=1Nβj(z)(dj0+dj1q-1)∑j=1Nβj(z)(1-q-1)bj1,withbj1=αj2.
Consider the following polynomial having its roots inside the unit circle of the complex plane:(30)F(q-1)=1-1.0164q-1+0.5102q-2,
where its poles, having a complex form, are the following:(31)z1=0.5082+0.5020i,z2=0.5082-0.5020i.
The chosen polynomial is identified with the denominator of every closed-loop subsystem which leads to the following equality:(32)1+(dj0-1-αj1)q-1+(αj1+dj1)q-2=1-1.0164q-1+0.5102q-2,∀j=1,…,N.
The solution of (32) gives us the numerator's parameters of the PDC controller as follows: (33)dj0=-0.0164+αj1,dj1=0.5102-αj1,∀j=1,…,N.
The transfer function of the PDC controller (30) leads to the following equalities:(34)∑j=1Nβj(1-q-1)αj2u(q-1)=∑j=1Nβj(dj0+dj1q-1)E(q-1),(35)∑j=1Nβj(αj2u(q-1)-αj2q-1u(q-1))=∑j=1Nβj(dj0E(q-1)+dj1q-1E(q-1)),(36)∑j=1Nβjαj2u(q-1)=∑j=1Nβj(αj2q-1u(q-1)+dj0E(q-1)+dj1q-1E(q-1)).
From (36) a recurrent form of the control law can be deduced, which will be(37)u(k)=1∑j=1Nβj(k)αj2(k)×∑j=1Nβj(k)(αj2(k)u(k-1)+dj0(k)e(k)+dj1(k)e(k-1)),
where e represents the error between the reference output and the system output,
(38)u(k)=∑j=1Nβj(k)(α̃j2(k)u(k-1)+d̃j0(k)e(k)+d̃j1(k)e(k-1)),
whereα̃j2=αj2(k)/∑j=1Nβj(k)αj2(k),d̃j0=dj0(k)/∑j=1Nβj(k)αj2(k),d̃j1=dj1(k)/∑j=1Nβj(k)αj2(k).With this expression it is also possible to take into account the saturation of the control input since we deal with local classical controllers with recurrent form [21], and so the new expression of the control law will be(39)us(k)=∑j=1Nβj(k)(α̃j2(k)us(k-1)+d̃j0(k)e(k)+d̃j1(k)e(k-1))
With
(40)us={50Hzifu≥50Hz,uif0<u<50Hz,0ifu≤0
represent the saturated control input.


### 4.4. Practical Results

The desired outputr is 27°C. The application was conducted with a sample time equal to 15 seconds on the day of 24 March 2011. It started at 8:39 AM and ended at 6:23 PM. The results of the proposed method are shown in Figures 10 and 11. Figures 12 and 13 represent the measured disturbance inputs of the system during the application. In the beginning and in the end of the application, where the values of the intercepted solar radiation and outside temperature are very low (corresponding to the morning and later in the evening before the sun set), the controller does not take action. This is because the value of the inside temperature is below the desired temperature.Figure 10
Practical results: outside temperature (blue line) and inside temperature (black line).Figure 11
Action control of the fuzzy controller.Figure 12
Intercepted solar radiation.Figure 13
Inside humidity.Moreover, we can note from Figure11 that the control signal presents high switching; it is because of some errors in measurements that affect the response of the controller; but, in spite of this, it still leads the process to the desired trajectory, which proves its robustness and efficiency.


## 4.1. Materials and Measurements

The process is a greenhouse with the following geometrical characteristics: 1.5 m height, total width of 1 m, and total length of 1.5 m.The inputs of the system are obtained from several sensors located inside and outside the greenhouse as follows.(i)
The inside and outside temperatures were measured by two LM35 transistors; both have an AD 620 amplifier to amplify the delivered signal.(ii)
The inside humidity was measured using a humidity sensor (SY-HS-230BT).(iii)
The solar radiation was obtained by a pyranometer type LP-PYRA 03 module.In order to control the ventilation rate, the greenhouse was equipped with two fans driven by two three-phase motors. These engines have a power supply delivered by a frequency converter (microdrive FC 51 Danfoss). The rotational speed of the fans is proportional to the frequency of the three-phase power supply and so the ventilation rate inside the greenhouse depends on the frequency of the power supply. Thus, the control input will be the frequency that will be computed by the fuzzy controller.The communication between MATLAB, the sensors, and the frequency converter is carried by the data acquisition module (KUSB-3100).


## 4.2. Fuzzy Identification

The structure of the discrete TS fuzzy system that will be approximated will have a local model for every rulej having the following form [20]:(27)R(j):Tj(k+1)=αj1(k)T(k)+αj2(k)u(k)+αj3(k)To(k)+αj4(k)Rs(k)+αj5(k)Rhi(k)+αj6(k)forj=1,…,N,
where u is the frequency delivered by the frequency converter to the motors equipped with fans to control the ventilation rate inside the greenhouse.The fuzzy identification is performed in two steps. The first one is the determination of the appropriate membership functions of each input to set the antecedent part of the TS fuzzy system, which is done using fuzzy clustering based on C-means algorithm. The membership functions were fixed andwould not change during the control application. For that purpose, on the day of 25-02-2011, we gave a random action control to collect the data needed under the sample time of 15 seconds; after that we set the membership functions in Figures 4, 5, 6, 7, and 8.Figure 4
Membership function of the inside temperature.Figure 5
Membership function of the delivered frequency.Figure 6
Membership function of the outside temperature.Figure 7
Membership function of the inside humidity.Figure 8
Membership function of the intercepted solar radiation.The second step is the estimation of (27) that represents the consequent part of the TS fuzzy system. For a time-variant system like a greenhouse, the more convenient way is an online fuzzy identification using the following recursive equations based on the recursive least square algorithm [20]:(28)Kj(k)=Pj(k-1)ze(k-1)λ/βj(z(k-1))+zeT(k-1)Pj(k-1)ze(k-1),θ̂j(k)=θ̂j(k-1)+Kj(k)(y(k)-zeT(k-1)θ̂j(k-1)),Pj(k)=1λ(Pj(k-1)-Kj(k)zeT(k-1)Pj(k-1)),
where λ is a forgetting factor, zeT=(TuToRsRhi1) is the regression vector of inputs, θ̂jT=(α̂j1α̂j2α̂j3α̂j4α̂j5α̂j6) is the estimated vector of parameters, and y is the output representing the inside temperature T.Pj(0) and θ̂j(0) are obtained from an offline fuzzy identification based on an ordinary weighted least square method [20], where the data used are those collected during the day before the recursive identification and also served to build the membership functions.The estimated vectorθ̂ will be used in an adaptive fuzzy control scheme (see Figure 9) to bring the inside temperature T to the desired temperature r using the fuzzy controller output u.Figure 9
Adaptive fuzzy control scheme.


## 4.3. Fuzzy Control

The system is a first-order one so the local controller for each rule will be a PI controller. From the previous section, the transfer function of the PDC controller can be deduced as follows:(29)P(q-1)=∑j=1Nβj(z)(dj0+dj1q-1)∑j=1Nβj(z)(1-q-1)bj1,withbj1=αj2.
Consider the following polynomial having its roots inside the unit circle of the complex plane:(30)F(q-1)=1-1.0164q-1+0.5102q-2,
where its poles, having a complex form, are the following:(31)z1=0.5082+0.5020i,z2=0.5082-0.5020i.
The chosen polynomial is identified with the denominator of every closed-loop subsystem which leads to the following equality:(32)1+(dj0-1-αj1)q-1+(αj1+dj1)q-2=1-1.0164q-1+0.5102q-2,∀j=1,…,N.
The solution of (32) gives us the numerator's parameters of the PDC controller as follows: (33)dj0=-0.0164+αj1,dj1=0.5102-αj1,∀j=1,…,N.
The transfer function of the PDC controller (30) leads to the following equalities:(34)∑j=1Nβj(1-q-1)αj2u(q-1)=∑j=1Nβj(dj0+dj1q-1)E(q-1),(35)∑j=1Nβj(αj2u(q-1)-αj2q-1u(q-1))=∑j=1Nβj(dj0E(q-1)+dj1q-1E(q-1)),(36)∑j=1Nβjαj2u(q-1)=∑j=1Nβj(αj2q-1u(q-1)+dj0E(q-1)+dj1q-1E(q-1)).
From (36) a recurrent form of the control law can be deduced, which will be(37)u(k)=1∑j=1Nβj(k)αj2(k)×∑j=1Nβj(k)(αj2(k)u(k-1)+dj0(k)e(k)+dj1(k)e(k-1)),
where e represents the error between the reference output and the system output,
(38)u(k)=∑j=1Nβj(k)(α̃j2(k)u(k-1)+d̃j0(k)e(k)+d̃j1(k)e(k-1)),
whereα̃j2=αj2(k)/∑j=1Nβj(k)αj2(k),d̃j0=dj0(k)/∑j=1Nβj(k)αj2(k),d̃j1=dj1(k)/∑j=1Nβj(k)αj2(k).With this expression it is also possible to take into account the saturation of the control input since we deal with local classical controllers with recurrent form [21], and so the new expression of the control law will be(39)us(k)=∑j=1Nβj(k)(α̃j2(k)us(k-1)+d̃j0(k)e(k)+d̃j1(k)e(k-1))
With
(40)us={50Hzifu≥50Hz,uif0<u<50Hz,0ifu≤0
represent the saturated control input.


## 4.4. Practical Results

The desired outputr is 27°C. The application was conducted with a sample time equal to 15 seconds on the day of 24 March 2011. It started at 8:39 AM and ended at 6:23 PM. The results of the proposed method are shown in Figures 10 and 11. Figures 12 and 13 represent the measured disturbance inputs of the system during the application. In the beginning and in the end of the application, where the values of the intercepted solar radiation and outside temperature are very low (corresponding to the morning and later in the evening before the sun set), the controller does not take action. This is because the value of the inside temperature is below the desired temperature.Figure 10
Practical results: outside temperature (blue line) and inside temperature (black line).Figure 11
Action control of the fuzzy controller.Figure 12
Intercepted solar radiation.Figure 13
Inside humidity.Moreover, we can note from Figure11 that the control signal presents high switching; it is because of some errors in measurements that affect the response of the controller; but, in spite of this, it still leads the process to the desired trajectory, which proves its robustness and efficiency.


## 5. Conclusion

TS fuzzy systems can assimilate the evolution of complex and nonlinear process and act  consequently to provide acceptable results in terms of modeling and control. However, a stability problem arises because of the great number of rules used to model such systems. This paper offers a solution using classical controllers like a PI, a PID, or an RST for each rule and a pole placement strategy. The fuzzy controller presented combines the multicontrol approach and the performance of the classical controllers. The result is a robust nonlinear control action that can also deal with time-variant systems.



---

*Source: 537491-2011-07-19.xml*