13.07.2015 Views

Assessment and Future Directions of Nonlinear Model Predictive ...

Assessment and Future Directions of Nonlinear Model Predictive ...

Assessment and Future Directions of Nonlinear Model Predictive ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Controlling Distributed Hyperbolic Plants 439û ∗ =ˆα(t) R(t) Lr(t) − T 0 (t)(17)Here, ˆx <strong>and</strong> ˆα are obtained using the state estimator (9), the adaptation law (13),u(¯t) is a sequence <strong>of</strong> step functions with amplitude u i (i =1,...,N u ) <strong>and</strong> durationTN u.Thevariable¯t represents time during the minimization horizon ¯t ∈ [0,H[. Theinitial condition ˆα(0) for the estimate <strong>of</strong> α is provided by the designer.Once the minimization result u(¯t) is obtained, according to a receding horizonscheme u 1 is applied to the plant at t + δ <strong>and</strong> the whole process is restarted,δ being an interval <strong>of</strong> time which is at least the time needed to compute thesolution. It is assumed that δ is much smaller than the sampling period. Thereare classes <strong>of</strong> plants, such as switched nonlinear systems [7], for which the choice<strong>of</strong> δ has a bearing on the stability properties <strong>of</strong> the predictive controller. In thecase at h<strong>and</strong>, no such problems were found.The minimization is always feasible since ũ = 0 (corresponding to u = u ∗ )preservesthe closed loop stability while satisfying the constraint (8) with V 0 = V rhc .4 Simulation ResultsSimulation results <strong>of</strong> the proposed RHC have been performed in a detailed model<strong>of</strong> the solar field obtained from first physical principles, according to [3], <strong>and</strong>calibrated with plant data. Experimental sequences for R(t) <strong>and</strong>T 0 (t) areused.The reduced model (3) uses 3 interior collocation points z =[0.113 0.500 0.887]<strong>and</strong> leads to a matrix −A/L whose eigenvalues have all strictly negative part.In order to configure the controller, the following parameter choices have beenmade: γ =1× 10 5 , K 0 = [15 15 15 15], ρ =1.5 × 10 −10 , H = 180 s <strong>and</strong> N u = 26.Figures (1) through (3) show the results. Fig. (1) shows the time evolution <strong>of</strong> thereference (easily recognizable by being a sequence <strong>of</strong> steps) <strong>and</strong> <strong>of</strong> the temperature3002502001501005000 1 2 3 4 5 6 7Time [h]Fig. 1. Solar field with RHC. Outlet oil temperature <strong>and</strong> reference <strong>and</strong> temperatureestimates at the collocation points [ o C].

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!