Optimise a reference that cuts my curve into 2 equal sections
1 view (last 30 days)
I am collecting data from excel, then basic idea is that i want to determine a constant value line that cuts my curve into 2 eqaual sections (area above the line and the curve = area under).
The objective is to determine the value of the red line.
Matt J on 12 Feb 2022
Edited: Matt J on 12 Feb 2022
%t= time, X=consumption, x=unknown midline
sol=solve( optimproblem('Objective',x,'Constraints',Xc-x*tc<=90,'ObjectiveSense','minimize') );
lb=sol.x; %lower bound
sol=solve( optimproblem('Objective',x,'Constraints',10<=Xc-x*tc,'ObjectiveSense','maximize') );
ub=sol.x; %upper bound
disp 'Problem is infeasible'
xunc=trapz(t,X)/(t(end)-t(1)); %unconstrained solution
x=min(max(xunc,lb),ub); %constrained solution