Doing an example find the optimal state and optimal control based on minimizing the performance index, Fmincon error of scalar value
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Example to find the optimal state and optimal control based on minimizing the performance index. I have written the code but the error shown in fmincon is "Supplied objective function must return a scalar value." I have checked the matrices but still the error is there. I have attached te codes, kindly tell how to resolve the problem.
EXAMPLE:
Find the optimal state and optimal control based on minimizing the performance index 𝐽=∫ (𝑥(𝑡) − 1/2 𝑢(𝑡)^2 ) 𝑑𝑡 , 0 ≤ 𝑡 ≤ 1 subject to 𝑢(𝑡) = 𝑥̇(𝑡) + 𝑥(𝑡) with the condition 𝑥(0) = 0, 𝑥(1) = 1 2 (1 − 1/𝑒 )^ 2 where 𝐽𝑒𝑥𝑎𝑐𝑡 = 0.08404562020 In this example the initial approximation is 𝑥1 (𝑡) = 1 2 (1 − 1/𝑒 ) ^2
CODES:
File1;
function F = cost_function(x)
global def;
global m;
% global P_alpha_1;
s=def.k ;
C3=x(1,(s+1):2*s);
;
u=(C3*m.H) ;
F=(x-(1/2)*(u*u'));
File 2:
function F = system_of_equations(x)
global def;
global m;
global init;function F = system_of_equations(x)
global def;
global m;
global init;
global P_alpha_1;
s=def.k ;
C1=x(1,1:1);
C3=x(1,(s+1):s);
x1=('C1*P_alpha_1*m.H') + init(1);
u=('C3*m.H') ;
D_alpha1_x1= 'C1*m.H';
M=Haar_matrix(s);
HC=M(:,s);
%%control law1
% F = horzcat( D_alpha1_x1 - u , ...
F = horzcat( D_alpha1_x1 - u , ...
('C1*P_alpha_1*HC') + init(1) - (1/2*((1-exp(-1))^2)) ) ;
end
File 3:
alpha_1=1;
k=8; %no. of Haar wavelets
b=1; %Total number of days to plot
initialize(alpha_1,k,b )
global m;
global init;
global def;
global P_alpha_1;
P_alpha_1=fractional_operation_matrix(k,alpha_1,b,m.H);
s=def.k;
% C3=x(1,((2*s)+1):3*s);
% u=(C3*m.H) ;
% cost function=(1/2)*u.^2;
x0=zeros(2,2*s);
% system_of_equations(x)
A = [];
bb = [];
Aeq = [];
beq = [];
lb = [];
ub = [];
fun = @cost_function;
nonlcon=@system_of_equations;
x = fmincon(fun,x0,A,bb,Aeq,beq,lb,ub,nonlcon)
2 comentarios
Matt J
el 17 de Feb. de 2022
It would be easiest if you run the code here in the forum, rather than on your local computer. That way, we can all see what it's doing in the same environment.
Ver también
Categorías
Más información sobre Introduction to Installation and Licensing en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!