Nonlinear optimization with differential equations in Constraint

2018 5 18
2 个回答
更新时间：2020 9 26
4 次查看（30 天）

1 个投票

Hello, I have a nonlinear optimization problem like min f(X), where X is a vector of four elements. And the constraints are in the form of dX/dt=g(X,u), where u is the decision variable and it is defined by: u=h(X). I tried to use fmincon and ode23 together, however, I do not know what to do with u=h(X) as the constraint. If I want to define that as a nonlinear constraint in fmincon, how to define the function as X is not still known. Thank you.

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

回答（2 个）

Nikhil Negi
Nikhil Negi 2018-5-18

1 个投票

i can think of two ways you can try both try replacing the dX/dt = g(X,u) with dX/dt = g(X,h(X)) so eventually you have dX/dt = p(X) and then continue with how you did it.
or you can define u as a separate function and pass it as a parameter when calculating dX/dt.

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

Sareh Agheb
Sareh Agheb 2018-5-19
编辑：Sareh Agheb 2018-5-19
Thanks Nikhil for your reply. I tried to define all the equations based on X, and make X the decision variable of fmincon instead of u. But, I still cannot get an accurate result. Here is the main body of my code (the simple version):
x0 = [0.1,0,-0.0002,0];
alpha=2;
x=zeros(N,4);
for ct1 = 0:Ts:Duration
options = optimoptions('fmincon','Algorithm','sqp','Display','none');
COSTFUN = @(x) MyobjFunction(u,x,x0) ;
nonlcon = []; lb=[]; ub=[];
x = fmincon(COSTFUN,x,[],[],[],[],lb,ub,nonlcon,options);
u=-alpha*x(:,2);
end
function costtot = MyobjFunction(u,x,x0)
x00=x0;
for k=1:N
tspan=[ct1+(k-1)*Ts ct1+k*Ts];
[~, x] = ode23(@(t,x) derivative(u,x),tspan, x00);
cost(k) =0.5*J1*(x(end,2)^2)+0.5*J2*(x(end,4)^2);
costtot=costtot+cost(k);
x00=[x(end,1),x(end,2),x(end,3),x(end,4)];
end
end
function dx = derivative(u,x)
dx = zeros(4,1);
dx(1)=x(2);
dx(2)=(xref1/J1)-((E*V)/(J1*X12))*sin(x(1)-x(3))-(-alpha*x(2))/J1;
dx(3)=x(4);
dx(4)=(xref2/J2)-((E*V)/(J2*X12))*sin(x(3)-x(1));
end
Nikhil Negi
Nikhil Negi 2018-5-21
hello sareh, can you provide the mathematical equations also. that would be really helpful
Nikhil Negi
Nikhil Negi 2018-5-21
编辑：Nikhil Negi 2018-5-21
i don't think your x0 is getting updated to x00=[x(end,1),x(end,2),x(end,3),x(end,4)] should it get updated for the next window of time?? i'm not really familiar with model predictive control but from your code it seems like you want it to get updated.
Sareh Agheb
Sareh Agheb 2018-5-22
Thank you Nikhil for your help. x00 and x is getting updated in MyobjFunction, but the problem is the results of x is not returned to the main program. So, what I get for both x and u is just zero.

请先登录，再进行评论。

Sareh Agheb
Sareh Agheb 2018-5-21
编辑：Sareh Agheb 2018-5-21

1 个投票

Hello Nikhil. I am going to apply Model Predictive Control on the above optimization problem. So, the objective should be minimized over the control time horizon. U is the control input and x is the states. Since, I am using a feedback control, u=-kx2. u is a vector of N by 1, and x is a matrix of N by 4, where N is the prediction horizon.

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Nafis Irtija
Nafis Irtija 2020-9-26
Hello Sareh, were you able to solve this problem? I am trying to solve an optimization problem similar to yours, and I am having a similar problem. I am unable to find any answer.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Nonlinear Control 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by