How to specify a nonlinear mpc controller for continuous time delay differential equation state function?

2023 3 14
1 个回答
回答已采纳
更新时间：2023 3 20
14 次查看（30 天）

0 个投票

I have a system of continuous-time nonlinear delay differential equations with two states and two explicit time delays. I would like to build a nonlinear mpc controller for this system, but I am not sure how to write a state function with time delays.
The system:
Where are the state varaibles, are constants, u is the manipulated variable, and are the continuous-time nonlinear state equations.
Is this possible to do using nlmpc ?

请先登录，再回答此问题。

 采纳的回答

1 个投票

You can basically add states to help model the delays. So your new discretized state vector would be [x(k) y(k) x(k-1) y(k-1) ... x(k-tau1/Ts) y(t-tau2/Ts)]. Then you can use the respective state function with nlmpc.

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Kayvon
Kayvon 2023-3-14
Thanks this worked great. I was able to get it to work using the following setup:
For
nlobj = nlmpc(6,6,1);
nlobj.Model.StateFcn = @stateFunc;
function xdot = stateFunc(x,u)
xdot = [f1(x(1),x(2),x(5),x(6),u(1)); % x(t+1) = f1(x,y,x(t-2),y(t-2),u)
f2(x(1),x(2),x(5),x(6),u(1)); % y(t+1) = f2(x,y,x(t-2),y(t-2),u)
x(1)-x(3); % x(t) = x(t-1)
x(2)-x(4); % y(t) = y(t-1)
x(3)-x(5); % x(t-1) = x(t-2)
x(4)-x(6)] % y(t-1) = y(t-2)
end
I'm guessing when matlab discretizes the system it computes not so that's why I have etc instead of just passing etc.
However, having a delay longer than 1 or 2 timesteps would introduce a lot of new variables and make the computation very slow. Is there a way around that?
Can't think of any other way that would work. Maybe another option would be to fit a data-driven model? If # of states becomes high with an analytical model, you may still be able to capture the dynamics with the original number of states if you have data. Model Predictive Control Toolbox supports neural ODEs since R2022b so you can give it a try
https://www.youtube.com/watch?v=q93jlLNwrVU
Kayvon
Kayvon 2023-3-15
Ok makes sense. I will look into that. Thanks for your help!
Kayvon
Kayvon 2023-3-20
In this case I was able to rescale the time variable in the delay differential equations so that the needed time delays were small which made the computation tractable. But in other cases I like the idea of fitting a data-driven model and I'm going to look into that as well anyways.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Linear Plant Specification 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by