Optimise a reference that cuts my curve into 2 equal sections
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello,
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.
0 commentaires
Réponse acceptée
Matt J
le 12 Fév 2022
Modifié(e) : Matt J
le 12 Fév 2022
%t= time, X=consumption, x=unknown midline
tc=t-t(1);
Xc=cumtrapz(t,X);
x=optimvar('x');
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
if lb>ub
disp 'Problem is infeasible'
else
xunc=trapz(t,X)/(t(end)-t(1)); %unconstrained solution
x=min(max(xunc,lb),ub); %constrained solution
end
3 commentaires
Plus de réponses (2)
William Rose
le 12 Fév 2022
If you mean that the area between the red line and the blue curve abve the line equals the area between the red line and the blue curve below the line, then the height of the red line is simply the mean value of the blue line data.
5 commentaires
dpb
le 12 Fév 2022
Attach the data as a .mat file.
Just looking at the graph, I'd venture it isn't possible to meet the constraints that the integral of the area be <95 when the magnitudes of the integrand are in the thousands. Just by a very crude approximation of the leftmost area (brown?) as roughly a rectangle of height ~(2800-1800) and a duration of 12, the integral would be roughly 1000*12 --> 12,000. Unless there's a very small scaling factor to be applied, it "just ain't a-gonna' happen".
Voir également
Catégories
En savoir plus sur Subplots dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!