Simple Tuning Rules for Feedforward Compensators Applied to Greenhouse Daytime Temperature Control Using Natural Ventilation

Ana Paola Montoya-Ríos 1 , Francisco García-Mañas 2 , José Luis Guzmán 2,* and Francisco Rodríguez 2 1 Faculty of Agricultural Sciences, National University of Colombia, Medellín Campus, Medellín CO050034, Antioquia, Colombia; apmontoy@unal.edu.co 2 Department of Informatics, University of Almería, CIESOL, ceiA3, E04120 Almería, Spain; francisco.gm@ual.es (F.G.-M.); frrodrig@ual.es (F.R.) * Correspondence: joseluis.guzman@ual.es; Tel.: +34-950-214-133


Introduction
Climate control for greenhouses is critical for optimal crop production. Since a greenhouse is a closed environment, climate regulation can ensure an appropriate growth of the crop. The air temperature inside the greenhouse considerably affects to the physiological processes of the plants and, thus, it is important to regulate it in an adequate range depending on the crop variety and the geographical location of the greenhouse. This accelerates its growth and increases the production quality while economizing energy consumption [1].
Controlling the greenhouse climate requires different types of systems to cause a particular effect on the inside climatic variables. The predominant systems are ventilation and heating, but advanced greenhouses can use modern systems such as carbon dioxide injection or artificial lights. Even so, automatic control applied to greenhouses is a complicated problem. The external weather conditions should be considered as disturbances in the formulation of the control problem since they constantly cause undesirable changes in the inside climatic variables.
In comparison with the prior cited studies, in this work, a simple but effective control solution focused on disturbances rejection was pursued to regulate the greenhouse air temperature using natural ventilation. Feedforward (FF) control is a basic strategy to reduce the effect of measurable disturbances on the process output and it is commonly implemented in combination with a PID controller in order to improve the performance of the feedback loop. In Reference [22], a non-linear FF control strategy was calculated and evaluated in simulation based on a non-linear greenhouse model. In contrast, in this manuscript, a different and simpler approach is proposed to calculate linear FF compensators. The main contribution of this work is the use of simple modelling and control approaches that can be easily implemented in any real greenhouse. First, system identification methods are applied to obtain reduced linear models of the greenhouse micro-climate dynamics. Then, the resulting models are used to design the FF compensators by utilizing the innovate tuning rules proposed in Reference [23]. These tuning rules are based on linear models of the process and are focused on improving the control system performance by minimizing the Integral Absolute Error (IAE). Notice that it is the first time that these tuning rules are experimentally tested.
The document is organized as follows-in Section 2 materials and methods are described, where a detailed explanation of the experimental greenhouse, system identification methods and developed control strategies can be found. Section 3 is dedicated to discuss the most relevant results concerning the identified models, simulation results, comparison of different control schemes, and the real FF control tests performed in the experimental greenhouse. Finally, in Section 4, some conclusions and future works are mentioned.

Experimental Greenhouse
Control tests for this work were performed in an experimental greenhouse ( Figure 1) situated at the Cajamar Foundation Experimental Centre "Las Palmerillas" (Almería, Spain). This traditional "parral-type" greenhouse has a surface of 877 m 2 (37.80 m × 23.20 m) and it is equipped with modern actuators to control the micro-climate under the cover. Some examples of the available actuators are humidification and dehumidification systems, pipe heating system, carbon dioxide enrichment system and natural ventilation system.
In this work, the developed control strategies are focused on utilizing the natural ventilation system. This system consists of five zenithal windows (8.36 m × 0.73 m) placed at the roof of the cover and two lateral windows (32.75 m × 1.90 m) positioned along the north and the south sides of the greenhouse. Three AC motors (0.37 kW and 2 × 0.18 kW) are utilized to automate the vents opening depending on the desired opening percentage (from 0% to 100%) calculated by the climate controllers. The desired opening percentage is achieved by sending an electrical current to each motor during a specific time interval.
Several climate monitoring systems are distributed inside the greenhouse in order to precisely obtain measurements from a large amount of sensors. In addition, a weather station is installed outside the greenhouse to register external climate variables. In particular, a probe (HC2S3, Campbell Scientific Ltd., Shepshed, UK) is used to measure the air temperature inside and outside the greenhouse. For the external solar radiation, a pyranometer (LP02, Hukseflux, Delft, The Netherlands) is employed, and measurements of the external wind velocity are provided by an anemometer (A100L2, Vector Instruments, Rhyl, UK). These sensors are connected to a series of acquisition devices (Compact FieldPoints, National Instruments, Austin, TX, USA) and the measurements are transmitted by an Ethernet network to a supervisory control and data acquisition (SCADA) system. Experimental data for this work were obtained from the mentioned monitoring systems with a sampling period of 30 s for each climatic variable.

System Identification Methodology
The micro-climate inside a greenhouse is a highly complex process that involves multiple variables, phenomena, time scales and non-linear dynamics. According to referent studies in the literature [1], this micro-climate can be modelled by a set of differential equations expressing energy transfers and mass balances. These type of complex models, based on first principles, are usually obtained to accurately simulate the dynamics of the process. However, these complex models are difficult to be mathematically managed in order to calculate control strategies. For this reason, system identification methods were utilized in this work to obtain black-box simplified models for the dynamics of the process occurring in a greenhouse.
System identification techniques allows to represent complex systems by using linear models based on data collected from input-output experiments performed in the real process. In the discussed case of a greenhouse, the inside micro-climate variables are strongly affected by the outside environmental conditions (disturbances) and the state of the actuators. For this work, the inside air temperature of the greenhouse was the desirable variable to be controlled by only using the natural ventilation system. Therefore, a simple model was required to capture the effect of the main disturbances and the natural ventilation on the inside air temperature. Notice that low-order models are required to design the control approach. However, due to the complexity of the greenhouse dynamics, first, a high-order linear model is obtained to capture the process behaviour, and afterwards, the model is reduced to a low-order model for control purposes.
As a previous step, some tests were carried out in the experimental greenhouse to register data from various days for the system identification procedure. In each test, the natural ventilation system was programmed to offer a changing signal with different ventilation opening percentages during the day. Hence, the effect of the natural ventilation on the inside air temperature could be properly reproduced by the model to be identified.
In the automatic control field, one of the most used structures are the autoregressive models [24], which coefficients can be easily determined through an identification procedure. For this work, an auto-regressive with exogenous input (ARX) model was selected to obtain the high-order empirical model for the inside air temperature of the greenhouse. A simple graphical description for an ARX model is shown in Figure 2 for a Single-Input and Single-Output (SISO) system. In particular, the implemented ARX model for this work is a Multiple-Input and Single-Output (MISO) model (see Section 3.1) that can be mathematically represented in discrete-time form by the following expression: where y is the output, u i are the inputs, e is the white noise, nu is the total number of inputs and nk i is the number of samples corresponding to each input time delay. A(z) and B i (z) are polynomials with the following structure: where na and nb are the order for polynomial A and B, respectively, and a and b are the coefficients to be determined in the identification procedure with real data of the process. Once the high-order ARX model was obtained, a model reduction stage was accomplished for control purposes as commented above. In this model reduction stage, the relationship for each input with the output was modelled as a first-order model described by a transfer function expressed in Laplace domain for continuous time as follows [25]: where k i is the static gain, τ i is the time constant and L i is the time delay (dead time), all referred to the ith input. Thanks to transfer function models, the development of automatic control strategies can be easily achieved as explained in the following Section 2.3.

Automatic Control Strategies
The control problem faced in this work consists in keeping the diurnal inside air temperature (controlled variable) close to a desired reference value (setpoint) by regulating the opening of the windows (manipulated variable) of the natural ventilation. The natural ventilation system is based on the air exchange occurring between the interior and the exterior of the greenhouse. The inside air temperature of the greenhouse is generally higher than the outside air temperature (due to the greenhouse effect). Every time that the vents are opened, the cooler and more dense outdoor air enters the greenhouse causing the hotter and less dense inside air to rise and exit through the roof vents. However, the following difficulties should be taken into account in order to design an appropriate control strategy [1]:

•
The greenhouse micro-climate is strongly affected by disturbances, both measurable and non-measurable. Thus, the designed controller should consider the effect of the outside weather conditions.

•
The motors of the natural ventilation system present two limitations: (i) actuator saturation, due to a limited opening range from 0% to 100%; and (ii) resolution, since the windows opening is performed in steps of 10%.
These issues have been considered in the proposed control approach as described in the following sections.

PID Control
The Proportional-Integral-Derivative (PID) controller is the most extended feedback control technique in most industrial processes and it is the one selected for this work because of its simplicity and high flexibility [25]. The PID controller allows to calculate a control action depending on the existing error between the measured process output and the desired setpoint. A basic PID controller can be mathematically expressed as follows: where u is the control signal and e is the control error. The main parameters of a PID controller are the proportional gain, k p , the integral time, T i , and the derivative time, T d . The values for these parameters can be calculated based on different tuning methods depending on the desired control response and the process dynamics [25]. A common formulation for a PID controller as a transfer function in Laplace domain is presented in Equation (6) and a basic feedback control scheme based on a PID controller is shown in Figure 3. In this work, a PI (Proportional-Integral) controller was selected and the derivative term was omitted mainly because of the noisy signals in the real facilities. An additional modification was made to the basic structure of a PI controller in order to take into account the saturation of the ventilation motors. An anti-windup mechanism based on the back-calculation technique was incorporated to the PI controller to prevent the undesirable "windup" effect of the integral term [25]. When the actuator is saturated, the anti-windup mechanism shown in Figure 4 is designed to reset the integrator value to zero by adding the signal e s multiplied by the tacking constant 1/T t . The parameter T t is an adjustable term to regulate how rapidly the integrator is reset. The lambda (λ) method was selected as the tuning approach to calculate the parameters for the PI controller and to impose a desired behaviour for the feedback control response. According to this method, the PI controller parameters are calculated using the following expressions [25,26]: where τ cl is the desired time constant for the closed-loop response, and k, τ and L are the parameters of the first-order model used to describe the process behaviour.

Feedforward Control
A classical feedforward control approach has been selected in this work to account for the measurable disturbance rejection problem. The basic control scheme is shown in Figure 5. The idea consists in adding a new block called feedforward compensator, G f f , that tries to reduce the effect of the disturbances before they affect to the process output [27]. The classical way to calculate the feedforward compensator is given by the following expression: where G(s) is the process model describing the dynamics between the process output, y, and the control signal, u; and G d (s) is the model that describes the dynamics between the process output and the disturbance signal, d.
As observed from Equation (9), the perfect realization of the feedforward compenstator is only feasible when there is no causality issues. Typically, there exist inversion problems because of the time delay of G d is smaller than in G. So, when perfect compensation is not possible, it is common in practice to omit the non-realizable part and to implement the rest of the compensator, or just to implement a static feedforward by using only the static gains of G d and G. These solutions allow to implement the feedforward control block and to eliminate the disturbance effect in steady-state, but the process output performance is considerably deteriorated.
Recently, this issue about feedforward compensator design has been analyzed in the literature, where simple tuning rules to account for the inversion problems have been proposed [23,28,29]. In this work, the tuning rules proposed in Reference [23] have been used. The tuning rules take into account the parameters of the process model as well as the action of the feedback controller during the design stage of the FF compensator. As a result, a guideline summary is proposed in order to minimize the IAE of the control approach while obtaining a disturbance rejection without overshoot.
The tuning rules are based on the assumptions that G and G d are described by first-order-plus-dead-time (FOPDT) models and the feedforward compensator is a lead-lag filter as follows: Then, the feedforward compensator is calculated according to the following guideline: 1. Set: 2. Calculate τ p as: 3. Calculate the compensator gain, k f f , considering the PI controller parameters (k p and T i ): Results for the design of the FF compensators with the mentioned simple rules are presented and discussed in Sections 3.3 and 3.4.

Software
MATLAB R (MathWorks, MA, USA) and Simulink R were used to develop and simulate the different control strategies presented in this work. Specifically, the System Identification Toolbox TM was used to obtain the ARX model and the low-order models expressed as transfer functions. Additionally, the PI controller and the feedforward compensators were programmed in MATLAB R code in discrete-time form and embedded in a SCADA system for real tests in the experimental greenhouse.

ARX Model
The developed ARX model intended as a MISO model is shown in Figure 6. For the system identification procedure, no distinction was made among the model inputs with the purpose of obtaining a simple prediction model. However, the outside weather variables should be considered as measurable disturbances. The three main disturbances affecting the inside air temperature were selected based on previous studies [1]: external solar radiation, external air temperature and external wind velocity. The ARX model was obtained with the system identification procedure explained in Section 2.2. A series of different ventilation tests were performed in the experimental greenhouse during spring 2020 in order to record a sufficient dataset for a proper model identification. In each test, diverse operational conditions were registered focusing on different ventilation opening profiles so that the data could represent the variety of the real greenhouse dynamics. Figure 7 shows an example for the inputs of the ARX model identification. The ventilation opening graphics reveals two different behaviours. During the night, the vents are usually closed and only some periodical openings are required if an active control for inside air humidity is being performed (depending on the crop necessities). During daytime, the ventilation opening varies depending mainly on the imposed limits for the inside air temperature control. Due to the weather conditions of the south-eastern of Spain and the dimensions of the experimental greenhouse, the vents are normally fully opened when the maximum level for external solar radiation occurs. The selected data for identification present a typical range of variation for climatic variables during the spring period at the location of the greenhouse.
So, the stored input data was used with the System Identification Toolbox TM in MATLAB R and several ARX models were obtained, selecting the one with best goodness of fitting. Figure 8 shows an example of the validation results for one day. As observed, the resulting linear model captures the main dynamics of the real inside air temperature. An statistical evaluation of the error between the real temperature and the predicted temperature by the model was accomplished. The mean absolute error for the model is 0.49 • C, the standard deviation is 0.63 • C, and the maximum absolute error is 2.82 • C. These values for the error are comparable to the ones obtained in previous studies [1]. Furthermore, simulations with data from other days offered a similar response for the ARX model. Hence, it can be concluded that the ARX model is validated to predict the inside air temperature of the greenhouse for control purposes.  The identified ARX model is a high-order model (10th-order) and its corresponding set of polynomials are presented in Equations (17)- (21).
B 2 (z) = 5.803 · 10 −5 z −1 − 3.879 · 10 −5 z −2 + 3.14 · 10 −5 z −3 + 3.535 · 10 −5 z −4 − 2.451 · 10 −5 z −5 + 3.468 · 10 −5 z −6 + 2.481 · 10 −5 z −7 − 7.229 · 10 −5 z −8 + 5.14 · 10 −5 z −9 − 3.998 · 10 −5 z −10 Despite the fact that the calculated ARX model is a high-order model, the validity of these types of models is constrained to particular conditions, which means that the model is effective for prediction when the input data are similar in operational conditions to the data that were used for model identification. In consequence, different ARX models should be identified to be utilized on different seasons of a year or during a crop cycle. Nevertheless, as shown in this work, an ARX model can be a simple and a reasonably reliable tool to obtain low-order models for the input-output relationships of a complex dynamic model, mainly with the intention of designing automatic control strategies.

Low-Order Models
Now, the high-order ARX model described in the previous section is reduced to low-order models to facilitate the control design stage. As a result, four different low-order models will be obtained such as shown in Figure 9. On the one hand, the relationship between the natural ventilation and the inside air temperature of the greenhouse is given by the process transfer function, G u . On the other hand, three additional transfer functions, G di , are obtained for each of the disturbances. The procedure to reduce the high-order ARX model consists in utilizing an unit step signal as input for the ventilation opening while the disturbances are set to zero. Then, the transient response for the inside air temperature predicted by the ARX model is analyzed applying the reaction curve method to calculate the characteristic parameters of a transfer function. The same actions are applicable to find a model for a measurable disturbance, in this case, the unit step signal is for one of the disturbances while setting to zero the ventilation opening and the remaining disturbances. The structure for the transfer functions was imposed as a first-order model with or without dead time, depending on the different transient responses. This procedure was smoothly accomplished with the System Identification Toolbox TM in MATLAB R . In Figure 10, a graphical comparison is presented to confirm that the transfer function models accurately capture the dynamics predicted by the ARX model. The identified first-order transfer functions are formulated in Equations (22)- (25), where the values of the parameters are rounded whereas a higher precision was preserved for MATLAB R calculations. These models are then utilized in the next section to design the PI controller and to calculate the FF compensators following the methods described in Section 2.3.

Design and Simulation of Control Strategies
The proposed control strategy in this work is a combination of a PI controller with feedforward compensators. The implementation scheme can be observed in Figure 11. This control scheme tackles the problem of temperature control for a greenhouse in a simple way, but taking into account the particularities mentioned in Section 2.3 regarding the presence of disturbances and the limitations of the actuators for natural ventilation. . Additionally, for the anti-windup mechanism, the T t coefficient (see Figure 4) was established as T t = √ T i [25]. Once the PI controller was tuned, the FF compensators were calculated. A summary for the designed FF transfer functions is presented in Table 1. Firstly, classical FF compensators, G * f f i , were determined with Equation (9). As can be observed, perfect realization was not possible since a causality problem appears due to the dead time of the process transfer function (see Equation (22)), which is greater than the dead times of the disturbances transfer functions. Thus, the time delay terms were excluded. When the simple tuning rules from Equations (13)- (15) were applied to determine the FF compensators, the same problem occurs with the time delays satisfying L d < L. However, the main advantage of these innovative rules is that the time delays can be taken into account to calculate the pole and the gain of the FF first-order transfer functions (if proven 0 < L − L d < 1.7 τ d ). Hence, slight differences in values for the parameters of the modern FF compensators compared against the classical approach can be observed in Table 1.
Simulations were executed in Simulink R to evaluate the performance of the control scheme before testing it in the experimental greenhouse. A Simulink R model was built to reproduce the structure shown in Figure 11, in which the ARX model was utilized as the predictor for the inside air temperature of the greenhouse. The sample time for the Simulink R model execution was fixed in 30 s, the same as in the real control system. Additionally, the limitation of the vents motors was taken into account with a saturation block and a function to imitate the real resolution by approximating the control signal to step intervals of 10%. For each simulation test, three variations of the control scheme were individually executed with data from the same day in order to compare the graphical differences between the following control strategies-(i) PI control, (ii) PI control plus classical FF compensators, and (iii) PI control plus FF compensators calculated with the simple tuning rules. The results for the simulation with data from 14 March 2020 are shown in Figure 12. In each graph, the time gap when the automatic control is active is indicated by two dashed vertical lines. The graph for the inside air temperature shows an acceptable control response closed to the desired temperature setpoint. Analyzing the graph for the error signal, it is clearly noticeable that, as expected, the PI controller without FF compensators presents an oscillatory response when rejecting the disturbances. When the FF compensators are included, the control response is smoother and the controlled variable is maintained closer to the setpoint. In regard to the control signals, saturation occurs for the PI controller during a longer time, but its control signal is less aggressive compare to the FF control. Nevertheless, the control signal for the FF control was expected to be more varying since the values of the three disturbances have to be cancelled in each time instant.
A different simulation was executed with data from 1 May 2020 and the results are presented in Figure 13. In this case, a higher value for the setpoint was established in order to test the control behaviour at a different operational point. Also, the presence of significant changes on external disturbances is reduced compared to the previous simulation. In the graph for the inside air temperature, it can be concluded the effectiveness of the control strategies, considering that the controlled variable is practically regulated at the imposed setpoint. In the error gap graphic, it is confirmed that the FF control is superior to the PI control, presenting an undisturbed response. Finally, in this case, the control signal for FF control is less aggressive since the changes in the disturbances occur with smaller magnitude.
In Table 2, the IAE values are presented for each control strategy in order to offer a quantitative comparison for the previous simulation results. The IAE values were calculated from the first time instant when the process output reached the setpoint until the last time instant when the output distanced from the setpoint. In both days, the IAE is significantly lower for the FF control strategies compared to the PI control, as expected. Regarding the FF control strategies, the best performance is achieved by the simple tuning rules for FF compensators, presenting the lowest values for IAE in both simulations. The IAE minimization was one of the objectives of applying the simple tuning rules for FF control. Although the IAE differences compared to the classical FF compensators may seem not important for a single day, the accumulated IAE when considering a full crop season would be relevant. Velocity (m/s)

Disturbance 3 -External wind velocity
Control ON Control OFF Figure 13. Results for simulation with data from 1 May 2020. As observed in the results for the simulations, in order to control the inside air temperature during daytime, a reachable setpoint should be fixed for the midday hours. If the automatic control is activated at dawn, a lower setpoint would be established. As the sun rises during the diurnal period, the external solar radiation and the external air temperature cause an increase in the inside air temperature. Thus, the setpoint would be rapidly surpassed and the actuators would be saturated to maintain the vents with an opening of 100%. This problem can be solved by activating the automatic control at a desired time instant by means of a bumpless control transfer mechanism [25]. When the midday hours are reached, this mechanism will allow to switch to the automatic controller if the process was being manually controlled until that time instant, for example. The bumpless control transfer mechanism was enabled in Simulink R for the PI controller and it was also implemented in code for the real tests in the greenhouse, as explained in the following Section 3.4.

Real Tests with Feedforward Control
For the control tests in the experimental greenhouse, only the FF compensators calculated with the simple tuning rules were evaluated, being the first time ever that these rules are tested in a real process. The control scheme from Figure 11 was adapted to a discrete-time script coded in MATLAB R to be on-line executed in the SCADA system of the greenhouse. The sample time of the control loop corresponds to 30 s.
For each control test, a procedure was followed to activate the control algorithm only during the midday hours. Initially, another climate controller for the greenhouse was active and the natural ventilation was regulated mainly to maintain the inside air humidity in an optimal range for the crop. When the inside air temperature approximates the midday hours and the solar radiation is close to its maximum value, the previously described bumpless control transfer mechanism allowed to switch the automatic control to the PI controller and the FF compensators implemented for this work.
Three tests were executed on different days in the greenhouse. The tests were performed on the last weeks of spring, when the most unfavorable weather conditions for greenhouses in Almería appear, so it was an excellent opportunity to evaluate the developed control strategy under hostile daytime temperatures. These tests were not carried out on continuous days, but the results are combined and presented in Figure 14 for a better visualization. Before analyzing the results, it is important to highlight that the FF compensators were tested in a different manner on each day to analyze better their effects. On the first day, FF compensators for external solar radiation and external air temperature were active, excluding the compensator for external wind velocity. On the second day, all the FF compensators were utilized. Finally, on the third day, only the compensator for external air temperature was active.
As can be observed in Figure 14, a satisfactory response is obtained for the inside air temperature, which is controlled at different setpoints. In the error gap graphic, the presented errors are contained in an interval smaller than ±1 • C. Therefore, it can be confirmed the remarkable performance of the control strategy considering the limitations of the ventilation motors and all the simplifications that were carried out to obtain the parameters for the PI controller and the FF compensators. Figure 15 shows the best control results corresponding to day 2 in detail. The inside air temperature is correctly controlled at the desired setpoint, with no significant deviations. In this day, all the FF compensators were active and, thus, the control signal for the ventilation opening is visibly varying. All the disturbances seem to be well rejected by the FF compensators on this windy day with an external air temperature that was not as stable as compared to day 1 and 3. Velocity (m/s)

Disturbance 3 -External wind velocity
Control ON Control OFF Figure 15. Real control test from 24 May 2020. All the feedforward compensators were active.
On day 2, the automatic control was activated at 12:35 by means of the previously described bumpless control transfer mechanism. During the first minutes, the control strategy determined that it was necessary to close the vents to reach the setpoint for the inside air temperature. Since the external solar radiation was still rising and the external air temperature presented a peak at 13:00, the FF compensators caused the ventilation to open progressively before the inside air temperature had reached the setpoint. Thanks to this change in the control signal, the inside air temperature did not increase excessively. Once the solar noon happened at 14:10, the control error remained near to zero, since the temperature inside the greenhouse was stabilized at the setpoint. However, in order to regulate the temperature near to the setpoint, the FF compensators had to deal with external air temperature variations and specially with the external wind velocity, which presented several rapid changes. For this reason, the control signal is significantly varying, and some transitions from 100% to 0% of ventilation opening occurred in time intervals of 15 min. As the external solar radiation decreases, only intermediate values (between 20% and 60%) were necessary for the ventilation opening. This tendency for the control signal was logical to avoid the decrease of the inside air temperature until the automatic control was disabled at 19:30.

Conclusions
In this paper, a control strategy based on simple tuning rules for feedforward compensators was designed and tested in a real greenhouse to control the inside air temperature by means of a natural ventilation system. The main objective of this work was to propose a simple but effective control scheme for a complex non-linear process, such as the micro-climate inside a greenhouse. The feedback controller was designed as a PI controller, which parameters were calculated by the lambda tuning method. The limitations of the ventilation actuators were taken into account with an anti-windup mechanism. For the FF control, the classical approach was compared to a series of novel tuning rules proposed in Reference [23]. These simple tuning rules have been tested in a real process for the first time in the present work.
A system identification methodology was accomplished to obtain a high-order ARX model. The main advantage of using an ARX model to simulate the dynamics of the greenhouse climate was the simplicity of the procedure to identify the model, since only a data-set of particular climatic variables was required. However, a remarkable disadvantage is the need to recalculate a new ARX model for different seasons of the year, due to the strong dependence of the greenhouse micro-climate on the meteorological variations of its location. The high-order ARX model was then reduced to low-order models expressed as first-order transfer functions. The transfer functions offered simplified models to easily calculate the PI controller parameters and the FF compensators. Some simulations were executed and the benefits of including FF compensators for disturbances rejection were evidenced when compare to PI control without FF compensation.
The designed control strategy presented positive results for the real tests performed in the experimental greenhouse. The inside air temperature was effectively controlled at different setpoints, under diverse operational conditions, and control errors were contained in an interval smaller than ±1 • C. Although disturbances rejection was satisfactory considering the limitations of the ventilation actuators, some slight imperfections appeared in the real tests due to the modeling simplifications and other non-measurable disturbances affecting the inside air temperature.
Future works may be focused on identifying more models for different seasons of the year, in order to design controllers and FF compensators for a wide range of operational conditions. Furthermore, additional FF compensators for important measurable disturbances, such as the inside soil surface temperature, could be included to improve the performance of the proposed control strategy.
Author Contributions: A.P.M.-R. was the responsible of the work, the design of the control approach and the writing of the manuscript. F.G.-M. was the responsible of the experimental results and the implementation issues. J.L.G. contributed to the elaboration of the manuscript and revision, and the control process design. F.R. was the responsible of the discussion section and the revision of the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This work has been funded by the National R+D+i Plan Project DPI2017-85007-R of the Spanish Ministry of Science, Innovation and Universities and ERDF funds.