Hi Hammad,
While I cannot explicitly solve your problem, I have some general tips when it comes to Optimization Problems with MATLAB and a possible bug in your code to list out.
The following lines should appear outside the for loop,
SoCeop = 1; % initial stored electricity, assume it is full
Electricity_discharged(1)=0;
Battery_Stored(1)=Capacity_Battery*SoCeop; %Initial electricity stored inside the battery in the beggining of the year
Currently this code is run on every single iteration of the loop, while this doesn't really change anything about the output of the code, I just wanted to point it out.
And when it comes to Optimization problems in MATLAB, it is often better to use the Problem-Based approach for solving the problem instead of using solver functions directly.
For example, in your use-case it becomes quite difficult to read through all the logic being used for your code because the MATLAB code and the optimization problem are so coupled together that it becomes difficult to separate out and focus on the actual problem being optimized. Whereas the Problem-Based approach produces code that is quite similar to the equations that you would write on a sheet of paper as part of the problem statement.
The following page gives a starter guide on how to use the Problem-Based approach for solving optimization problems:
If you have computed the parameters and input values correctly as per the provided code and your problem statement and you are still not getting any feasible points, then there really is no other option but to increase the size of the potential feasibility region by loosening the restrictions on the bounds of the decision variables like increasing the battery capacity's upper bound and so on.
I have another question as well, how is the initial guess for the decision variable being used to compute the constraints for optimizing the decision variable itself?
For example:
The "Electricity_discharged" variable is used in "Aeq" and the "Electricity_charged" variable is used in "Beq" to formulate the equality constraints for the problem. However these two variables were computed indirectly with the help of the "Capacity_Battery" variable which is initialized to 2000 at the start of the code and is itself a decision variable. From a code point of view all this works, but logically speaking from the optimization problem's point of view, it does not make sense to do this. Maybe a reformulation of the problem itself is in order?
Hope I was able to provide some insights and thinking points to take this further!
