Any methods to solve optimal control indirect method such as Pontryagin Minimum Principle?

Hi,

I would like to solve an optimal control problem using Pontryagin Minimum Principle. There are few numerical techniques with MATLAB examples using sym toolbox, bvp4c and ODE45 using shooting method. I am trying to implement using ODE45 solver by following steps:

1. Initialize states, co-state and control
2. ODE45 solver in forward time to find states.
3. ODE45 solver in backward time to find co-states
4. Find hamiltonian and performance index. If difference between previous and present hamiltonian is below given threshold the stop iteration or else go to 5.
5. Update the control with a selected step size.

I am facing many difficulties in implementation such as;

1. Difficulty in selecting initial control and co-state.
2. The state is not bounded within a range as it must be. (must be between 0 and 1, and final state must be a desired value)

I have attached the reference paper which I am trying to implement. If anyone knows any MATLAB toolbox to solve the same, kindly let me know.

Answer 1

[M]

Hi,

here are some related questions to start with:

https://mathworks.com/matlabcentral/answers/440335-how-to-implement-optimal-control-using-pontryagin-s-minimum-principle-for-hess-in-ev-application-as

https://mathworks.com/matlabcentral/answers/438910-how-to-solve-in-matlab-2018b

https://mathworks.com/matlabcentral/answers/392949-can-i-solve-equations-backward-in-time-pontryagin-maximum-principle-optimal-control-problem