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.

170 H.G. Bock et al.change anymore, but in an even larger neighborhood where the linearization isstill valid, see the illustration in Figure 1. In the case that only approximations<strong>of</strong> Jacobian <strong>and</strong> Hessian are used within the QP, we still obtain‖w 1 − w ∗ (x 0 )‖≤κ ‖x ′ 0 − x 0 ‖ ,with κ being small if the quality <strong>of</strong> the approximations is good.It is interesting to note that the good prediction properties are independentfrom the class <strong>of</strong> optimization problems. Therefore, the scheme can be appliedto solve both tracking problems <strong>and</strong> problems with an economic objective. Anexample for the latter can be found in [14].3.2 St<strong>and</strong>ard Real-Time Iteration SchemeLet us now consider the full real-time scenario, where we want to solve a sequence<strong>of</strong> optimization problems P (x(t)) where x(t) is the system state that changescontinuously with time <strong>and</strong> which is used as initial value x 0 in problem (8).In the st<strong>and</strong>ard real-time iteration scheme [8, 10] we proceed as follows:Start with an initial guess (w 0 ,λ 0 ,µ 0 ), <strong>and</strong> perform the following steps fork =0, 1,...:1. Preparation: Based on the current solution guess (w k ,λ k ,µ k ), compute allfunctions <strong>and</strong> their exact Jacobians that are necessary to build the QP (13),<strong>and</strong> prepare the QP solution as far as possible without knowledge <strong>of</strong> x 0 (seeSection 4.1 for more details). This corresponds to the initialization neededfor the initial value embedding stated above.2. Feedback Response: at time t k , obtain the initial value x 0 := x(t k )fromthe real system state; solve the QP (13) to obtain the step ∆w k =(∆s x 0 k ,∆sz 0k ,∆u 0k,...), <strong>and</strong> give the approximation ũ(x(t k )) := u 0k + ∆u 0kimmediately to the real system.3. Transition: Set the next solution guess asw k+1 := w k + ∆w k ,λ k+1 := λ QPk , <strong>and</strong> µ k+1 := µ QPk .3.3 Nominal Stability <strong>of</strong> the Real-Time Iteration SchemeA central question in NMPC is nominal stability <strong>of</strong> the closed loop. For thereal-time iteration scheme, the state vector <strong>of</strong> the closed loop consists <strong>of</strong> thereal system state x(t k ) <strong>and</strong> the content (w k ,λ k ,µ k ) <strong>of</strong> the prediction horizon inthe optimizer. Due to the close connection <strong>of</strong> system <strong>and</strong> optimizer, stability <strong>of</strong>the closed loop system can only be addressed by combining concepts from both,NMPC stability theory <strong>and</strong> convergence analysis <strong>of</strong> Newton-type optimizationmethods. For the st<strong>and</strong>ard real-time iteration scheme this analysis has been carriedout in [11], <strong>and</strong> for a related scheme with shift in [12]. In these papers, pro<strong>of</strong>s<strong>of</strong> nominal stability <strong>of</strong> the closed loop are given under reasonable assumptions.The class <strong>of</strong> feedback controls for shrinking horizon problems is treated in [9].

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

Saved successfully!

Ooh no, something went wrong!