unbounded problem #989

BYS543 · 2021-09-11T14:16:36Z

BYS543
Sep 11, 2021

Hi Johan,

   Con = [........];
   Obj = a1+a2;
   optimize(Con,-Obj);

If i write this, and solve the maximization problem with the objective funtion Obj, which is equal to the minimization problem with the objective function -Obj . And it can be solved（although it tells the numerical issues）. However, if i write

   Con = [........];
   Con = [Con,Obj<=a1+a2];
   optimize(Con,-Obj);

it tells me the model is unbounded(also the numerical issues).
I guess the unbounded reason is because the numerical issues?

johanlofberg · 2021-09-11T14:27:27Z

johanlofberg
Sep 11, 2021
Maintainer

Impossible to answer without reproducible code

2 replies

BYS543 Sep 11, 2021
Author

Thank you, Johan.
The codes which can be solved is as follows

par1 = [617.500000000000;617.500000000000;308.750000000000;617.500000000000;617.500000000000;308.750000000000;427.500000000000;332.500000000000;617.500000000000;152;308.750000000000;152;308.750000000000;308.750000000000;308.750000000000;308.750000000000;152;152;152;76;76;76;76;76;76;152;152;152;152;76;76;76;76;76;76;308.750000000000;308.750000000000;76;76;76;152;152];
par2 = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;18.7224532224532,22.1060291060291,25.8279625779626,28.9859667359667,29.7077962577962,28.4220374220374,27.7453222453222,24.0233887733888,21.6548856548857,18.1585239085239,16.5795218295218,18.0457380457380,19.1735966735967,20.8653846153846,23.4594594594595,26.1663201663202,27.7453222453222,27.1813929313929,26.6174636174636,25.3768191268191,23.0083160083160,20.5270270270270,18.1585239085239,14.7749480249480;2.07068607068607,2.44490644490644,2.85654885654886,3.20582120582120,3.28565488565488,3.14345114345114,3.06860706860707,2.65696465696466,2.39501039501039,2.00831600831601,1.83367983367983,1.99584199584200,2.12058212058212,2.30769230769231,2.59459459459459,2.89397089397089,3.06860706860707,3.00623700623701,2.94386694386694,2.80665280665281,2.54469854469854,2.27027027027027,2.00831600831601,1.63409563409563;6.55717255717256,7.74220374220374,9.04573804573804,10.1517671517671,10.4045738045738,9.95426195426195,9.71725571725571,8.41372141372141,7.58419958419958,6.35966735966736,5.80665280665281,6.32016632016632,6.71517671517671,7.30769230769231,8.21621621621621,9.16424116424116,9.71725571725571,9.51975051975052,9.32224532224532,8.88773388773389,8.05821205821206,7.18918918918919,6.35966735966736,5.17463617463617;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;19.6715176715177,23.2266112266112,27.1372141372141,30.4553014553015,31.2137214137214,29.8627858627859,29.1517671517671,25.2411642411642,22.7525987525988,19.0790020790021,17.4199584199584,18.9604989604990,20.1455301455301,21.9230769230769,24.6486486486486,27.4927234927235,29.1517671517671,28.5592515592516,27.9667359667360,26.6632016632017,24.1746361746362,21.5675675675676,19.0790020790021,15.5239085239085;25.8835758835759,30.5613305613306,35.7068607068607,40.0727650727651,41.0706860706861,39.2931392931393,38.3575883575884,33.2120582120582,29.9376299376299,25.1039501039501,22.9209979209979,24.9480249480249,26.5072765072765,28.8461538461538,32.4324324324324,36.1746361746362,38.3575883575884,37.5779625779626,36.7983367983368,35.0831600831601,31.8087318087318,28.3783783783784,25.1039501039501,20.4261954261954;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;5.00415800415800,5.90852390852391,6.90332640332640,7.74740124740125,7.94033264033264,7.59667359667360,7.41580041580041,6.42099792099792,5.78794178794179,4.85343035343035,4.43139293139293,4.82328482328482,5.12474012474012,5.57692307692308,6.27027027027027,6.99376299376299,7.41580041580041,7.26507276507276,7.11434511434511,6.78274428274428,6.14968814968815,5.48648648648649,4.85343035343035,3.94906444906445;2.76091476091476,3.25987525987526,3.80873180873181,4.27442827442827,4.38087318087318,4.19126819126819,4.09147609147609,3.54261954261954,3.19334719334719,2.67775467775468,2.44490644490644,2.66112266112266,2.82744282744283,3.07692307692308,3.45945945945946,3.85862785862786,4.09147609147609,4.00831600831601,3.92515592515592,3.74220374220374,3.39293139293139,3.02702702702703,2.67775467775468,2.17879417879418;9.66320166320166,11.4095634095634,13.3305613305613,14.9604989604990,15.3330561330561,14.6694386694387,14.3201663201663,12.3991683991684,11.1767151767152,9.37214137214137,8.55717255717256,9.31392931392931,9.89604989604989,10.7692307692308,12.1081081081081,13.5051975051975,14.3201663201663,14.0291060291060,13.7380457380457,13.0977130977131,11.8752598752599,10.5945945945946,9.37214137214137,7.62577962577962;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;5.34927234927235,6.31600831600832,7.37941787941788,8.28170478170478,8.48794178794179,8.12058212058212,7.92723492723493,6.86382536382536,6.18711018711019,5.18814968814969,4.73700623700624,5.15592515592516,5.47817047817048,5.96153846153846,6.70270270270270,7.47609147609148,7.92723492723493,7.76611226611227,7.60498960498960,7.25051975051975,6.57380457380457,5.86486486486486,5.18814968814969,4.22141372141372;7.07484407484407,8.35343035343035,9.75987525987526,10.9532224532225,11.2259875259875,10.7401247401247,10.4844074844075,9.07796257796258,8.18295218295218,6.86174636174636,6.26507276507276,6.81912681912682,7.24532224532224,7.88461538461538,8.86486486486486,9.88773388773389,10.4844074844075,10.2713097713098,10.0582120582121,9.58939708939709,8.69438669438669,7.75675675675676,6.86174636174636,5.58316008316008;3.01975051975052,3.56548856548857,4.16580041580042,4.67515592515592,4.79158004158004,4.58419958419958,4.47505197505198,3.87474012474012,3.49272349272349,2.92879417879418,2.67411642411642,2.91060291060291,3.09251559251559,3.36538461538462,3.78378378378378,4.22037422037422,4.47505197505198,4.38409563409563,4.29313929313929,4.09303534303534,3.71101871101871,3.31081081081081,2.92879417879418,2.38305613305613;7.76507276507276,9.16839916839917,10.7120582120582,12.0218295218295,12.3212058212058,11.7879417879418,11.5072765072765,9.96361746361746,8.98128898128898,7.53118503118503,6.87629937629937,7.48440748440748,7.95218295218295,8.65384615384615,9.72972972972973,10.8523908523909,11.5072765072765,11.2733887733888,11.0395010395010,10.5249480249480,9.54261954261954,8.51351351351351,7.53118503118503,6.12785862785863;2.76091476091476,3.25987525987526,3.80873180873181,4.27442827442827,4.38087318087318,4.19126819126819,4.09147609147609,3.54261954261954,3.19334719334719,2.67775467775468,2.44490644490644,2.66112266112266,2.82744282744283,3.07692307692308,3.45945945945946,3.85862785862786,4.09147609147609,4.00831600831601,3.92515592515592,3.74220374220374,3.39293139293139,3.02702702702703,2.67775467775468,2.17879417879418;8.19646569646570,9.67775467775468,11.3071725571726,12.6897089397089,13.0057172557173,12.4428274428274,12.1465696465696,10.5171517671518,9.48024948024948,7.94958419958420,7.25831600831601,7.90020790020790,8.39397089397089,9.13461538461538,10.2702702702703,11.4553014553015,12.1465696465696,11.8996881496881,11.6528066528067,11.1096673596674,10.0727650727651,8.98648648648648,7.94958419958420,6.46829521829522;1.89812889812890,2.24116424116424,2.61850311850312,2.93866943866944,3.01185031185031,2.88149688149688,2.81288981288981,2.43555093555094,2.19542619542620,1.84095634095634,1.68087318087318,1.82952182952183,1.94386694386694,2.11538461538462,2.37837837837838,2.65280665280665,2.81288981288981,2.75571725571726,2.69854469854470,2.57276507276507,2.33264033264033,2.08108108108108,1.84095634095634,1.49792099792100;15.0987525987526,17.8274428274428,20.8290020790021,23.3757796257796,23.9579002079002,22.9209979209979,22.3752598752599,19.3737006237006,17.4636174636175,14.6439708939709,13.3705821205821,14.5530145530146,15.4625779625780,16.8269230769231,18.9189189189189,21.1018711018711,22.3752598752599,21.9204781704782,21.4656964656965,20.4651767151767,18.5550935550936,16.5540540540541,14.6439708939709,11.9152806652807;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;2.76091476091476,3.25987525987526,3.80873180873181,4.27442827442827,4.38087318087318,4.19126819126819,4.09147609147609,3.54261954261954,3.19334719334719,2.67775467775468,2.44490644490644,2.66112266112266,2.82744282744283,3.07692307692308,3.45945945945946,3.85862785862786,4.09147609147609,4.00831600831601,3.92515592515592,3.74220374220374,3.39293139293139,3.02702702702703,2.67775467775468,2.17879417879418;7.50623700623700,8.86278586278586,10.3549896049896,11.6211018711019,11.9104989604990,11.3950103950104,11.1237006237006,9.63149688149688,8.68191268191268,7.28014553014553,6.64708939708940,7.23492723492723,7.68711018711019,8.36538461538461,9.40540540540540,10.4906444906445,11.1237006237006,10.8976091476091,10.6715176715177,10.1741164241164,9.22453222453222,8.22972972972973,7.28014553014553,5.92359667359667;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;3.01975051975052,3.56548856548857,4.16580041580042,4.67515592515592,4.79158004158004,4.58419958419958,4.47505197505198,3.87474012474012,3.49272349272349,2.92879417879418,2.67411642411642,2.91060291060291,3.09251559251559,3.36538461538462,3.78378378378378,4.22037422037422,4.47505197505198,4.38409563409563,4.29313929313929,4.09303534303534,3.71101871101871,3.31081081081081,2.92879417879418,2.38305613305613;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;2.07068607068607,2.44490644490644,2.85654885654886,3.20582120582120,3.28565488565488,3.14345114345114,3.06860706860707,2.65696465696466,2.39501039501039,2.00831600831601,1.83367983367983,1.99584199584200,2.12058212058212,2.30769230769231,2.59459459459459,2.89397089397089,3.06860706860707,3.00623700623701,2.94386694386694,2.80665280665281,2.54469854469854,2.27027027027027,2.00831600831601,1.63409563409563;9.14553014553014,10.7983367983368,12.6164241164241,14.1590436590437,14.5116424116424,13.8835758835759,13.5530145530146,11.7349272349272,10.5779625779626,8.87006237006237,8.09875259875260,8.81496881496882,9.36590436590436,10.1923076923077,11.4594594594595,12.7817047817048,13.5530145530146,13.2775467775468,13.0020790020790,12.3960498960499,11.2390852390852,10.0270270270270,8.87006237006237,7.21725571725572];
par3 = [1;1;2;3;2;2;4;5;6;6;6;6;9;9;4;12;12;12;12;14;16;15;18;19;10;10;10;10;21;15;22;23;24;25;25;27;27;27;27;29;8;6];
par4 = [2;3;4;4;5;6;6;7;7;8;9;10;11;10;12;13;14;15;16;15;17;18;19;20;20;17;21;22;22;23;24;24;25;26;27;28;28;29;30;30;28;28];
par5 = [0.0300000000000000;0.0200000000000000;0.0200000000000000;0;0.0200000000000000;0.0200000000000000;0;0.0100000000000000;0.0100000000000000;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0.0200000000000000;0.0100000000000000];
par6 = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;134.000000000001,164.000000000000,180,180,180,180,180,180,180,151.000000000000,136.999999999997,150.000000000000,160.000000000001,174.999999999999,180,180,180,180,180,180,180,172.000000000000,160.999999999998,130.999999999999;0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1.77635683940025e-15,-1.77635683940025e-15,1.77635683940025e-15,0,0,0,0,0,0,0;10,10,27.0000000000015,35,35,35,35,33.0000000000012,12.0000000000007,10,10,10,10,10,27.9999999999973,35,35,35,35,33.0000000000008,24.0000000000132,10.0000000000000,0,0;10,10,10,30,30,24.9999999999995,19.0000000000002,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-2.66453525910038e-15,-1.99840144432528e-15,-4.44089209850063e-15,-1.33226762955019e-15,0,-3.44613226843649e-13,0,0,0;12.0000000000000,12.0000000000000,12.0000000000000,12.0000000000001,18.3999999999997,12.0000000000000,11.9999999999999,0,0,0,0,0,0,0,0,17.0000000000008,31.0000000000002,25.9999999999996,21.0000000000006,12,0,0,0,0];
par7 = [1;2;22;27;23;13];
options = sdpsettings('solver','gurobi');
a1 = sdpvar(30,24);
a2 = sdpvar(30,24);
a3 = sdpvar(42,24);
a4 = sdpvar(42,24);
a5 = sdpvar(42,24);
a6 = sdpvar(42,24);
a7 = sdpvar(6,24);
a8 = sdpvar(6,24);
a9 = sdpvar(30,24);
a10 = sdpvar(30,24);   
o1 = sdpvar(30,24);
o3 = sdpvar(42,24);
o4 = sdpvar(6,24);
o5 = sdpvar(30,24);
h = sdpvar(1,1);
Con = [];
for t=1:24
    for i=1:42
        Con = [Con,-a1(par3(i),t)+a1(par4(i),t)+a2(par3(i),t)-a2(par4(i),t)+a3(i,t)/100-a4(i,t)/100+a5(i,t)-a6(i,t) == 0];               
        o3(i,t)=-a5(i,t)*par1(i)-a6(i,t)*par1(i);            
    end
    for g=1:6
        Con = [Con,a1(par7(g,1),t)-a2(par7(g,1),t)+a7(g,t)-a8(g,t)<=0];
        o4(g,t)=par6(g,t)*a7(g,t)-par6(g,t)*a8(g,t);                
    end
    for i = 1:30
        from_index = find(par3(i) == i)';
        to_index = find(par4(i) == i)';
        if isempty(from_index)
            from_branch = 0;
        else
            from_branch = -ones(1,length(from_index))*(a3(from_index,t)./par5(from_index))+ones(1,length(from_index))*(a4(from_index,t)./par5(from_index)); 
        end     
        if isempty(to_index)
            to_branch = 0;
        else
            to_branch = ones(1,length(to_index))*(a3(to_index,t)./par5(from_index))-ones(1,length(to_index))*(a4(to_index,t)./par5(from_index));
        end               
        Con = [Con,from_branch+to_branch+a9(i,t)-a10(i,t) == 0];
        Con = [Con,-a1(i,t)+a2(i,t)<=1000];
        Con = [Con,a1(i,t)-a2(i,t)<=1000];    
        o1(i,t)=par2(i,t)*a1(i,t)-par2(i,t)*a2(i,t);
        o5(i,t)=-(1.05*a9(i,t)+1.05*a10(i,t));                
    end
end
Con = [Con,a1>=0,a2>=0];
Con = [Con,a3>=0,a4>=0];   
Con = [Con,a5>=0,a6>=0];      
Con = [Con,a7>=0,a8>=0]; 
Con = [Con,a9>=0,a10>=0];   
%%% constraint needs to pay attention
Con = [Con,h==sum(sum(o1+o5))+sum(sum(o3))+sum(sum(o4))];
optimize(Con, -h, options);

The code which is unbounded, is as follows

 par1 = [617.500000000000;617.500000000000;308.750000000000;617.500000000000;617.500000000000;308.750000000000;427.500000000000;332.500000000000;617.500000000000;152;308.750000000000;152;308.750000000000;308.750000000000;308.750000000000;308.750000000000;152;152;152;76;76;76;76;76;76;152;152;152;152;76;76;76;76;76;76;308.750000000000;308.750000000000;76;76;76;152;152];
par2 = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;18.7224532224532,22.1060291060291,25.8279625779626,28.9859667359667,29.7077962577962,28.4220374220374,27.7453222453222,24.0233887733888,21.6548856548857,18.1585239085239,16.5795218295218,18.0457380457380,19.1735966735967,20.8653846153846,23.4594594594595,26.1663201663202,27.7453222453222,27.1813929313929,26.6174636174636,25.3768191268191,23.0083160083160,20.5270270270270,18.1585239085239,14.7749480249480;2.07068607068607,2.44490644490644,2.85654885654886,3.20582120582120,3.28565488565488,3.14345114345114,3.06860706860707,2.65696465696466,2.39501039501039,2.00831600831601,1.83367983367983,1.99584199584200,2.12058212058212,2.30769230769231,2.59459459459459,2.89397089397089,3.06860706860707,3.00623700623701,2.94386694386694,2.80665280665281,2.54469854469854,2.27027027027027,2.00831600831601,1.63409563409563;6.55717255717256,7.74220374220374,9.04573804573804,10.1517671517671,10.4045738045738,9.95426195426195,9.71725571725571,8.41372141372141,7.58419958419958,6.35966735966736,5.80665280665281,6.32016632016632,6.71517671517671,7.30769230769231,8.21621621621621,9.16424116424116,9.71725571725571,9.51975051975052,9.32224532224532,8.88773388773389,8.05821205821206,7.18918918918919,6.35966735966736,5.17463617463617;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;19.6715176715177,23.2266112266112,27.1372141372141,30.4553014553015,31.2137214137214,29.8627858627859,29.1517671517671,25.2411642411642,22.7525987525988,19.0790020790021,17.4199584199584,18.9604989604990,20.1455301455301,21.9230769230769,24.6486486486486,27.4927234927235,29.1517671517671,28.5592515592516,27.9667359667360,26.6632016632017,24.1746361746362,21.5675675675676,19.0790020790021,15.5239085239085;25.8835758835759,30.5613305613306,35.7068607068607,40.0727650727651,41.0706860706861,39.2931392931393,38.3575883575884,33.2120582120582,29.9376299376299,25.1039501039501,22.9209979209979,24.9480249480249,26.5072765072765,28.8461538461538,32.4324324324324,36.1746361746362,38.3575883575884,37.5779625779626,36.7983367983368,35.0831600831601,31.8087318087318,28.3783783783784,25.1039501039501,20.4261954261954;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;5.00415800415800,5.90852390852391,6.90332640332640,7.74740124740125,7.94033264033264,7.59667359667360,7.41580041580041,6.42099792099792,5.78794178794179,4.85343035343035,4.43139293139293,4.82328482328482,5.12474012474012,5.57692307692308,6.27027027027027,6.99376299376299,7.41580041580041,7.26507276507276,7.11434511434511,6.78274428274428,6.14968814968815,5.48648648648649,4.85343035343035,3.94906444906445;2.76091476091476,3.25987525987526,3.80873180873181,4.27442827442827,4.38087318087318,4.19126819126819,4.09147609147609,3.54261954261954,3.19334719334719,2.67775467775468,2.44490644490644,2.66112266112266,2.82744282744283,3.07692307692308,3.45945945945946,3.85862785862786,4.09147609147609,4.00831600831601,3.92515592515592,3.74220374220374,3.39293139293139,3.02702702702703,2.67775467775468,2.17879417879418;9.66320166320166,11.4095634095634,13.3305613305613,14.9604989604990,15.3330561330561,14.6694386694387,14.3201663201663,12.3991683991684,11.1767151767152,9.37214137214137,8.55717255717256,9.31392931392931,9.89604989604989,10.7692307692308,12.1081081081081,13.5051975051975,14.3201663201663,14.0291060291060,13.7380457380457,13.0977130977131,11.8752598752599,10.5945945945946,9.37214137214137,7.62577962577962;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;5.34927234927235,6.31600831600832,7.37941787941788,8.28170478170478,8.48794178794179,8.12058212058212,7.92723492723493,6.86382536382536,6.18711018711019,5.18814968814969,4.73700623700624,5.15592515592516,5.47817047817048,5.96153846153846,6.70270270270270,7.47609147609148,7.92723492723493,7.76611226611227,7.60498960498960,7.25051975051975,6.57380457380457,5.86486486486486,5.18814968814969,4.22141372141372;7.07484407484407,8.35343035343035,9.75987525987526,10.9532224532225,11.2259875259875,10.7401247401247,10.4844074844075,9.07796257796258,8.18295218295218,6.86174636174636,6.26507276507276,6.81912681912682,7.24532224532224,7.88461538461538,8.86486486486486,9.88773388773389,10.4844074844075,10.2713097713098,10.0582120582121,9.58939708939709,8.69438669438669,7.75675675675676,6.86174636174636,5.58316008316008;3.01975051975052,3.56548856548857,4.16580041580042,4.67515592515592,4.79158004158004,4.58419958419958,4.47505197505198,3.87474012474012,3.49272349272349,2.92879417879418,2.67411642411642,2.91060291060291,3.09251559251559,3.36538461538462,3.78378378378378,4.22037422037422,4.47505197505198,4.38409563409563,4.29313929313929,4.09303534303534,3.71101871101871,3.31081081081081,2.92879417879418,2.38305613305613;7.76507276507276,9.16839916839917,10.7120582120582,12.0218295218295,12.3212058212058,11.7879417879418,11.5072765072765,9.96361746361746,8.98128898128898,7.53118503118503,6.87629937629937,7.48440748440748,7.95218295218295,8.65384615384615,9.72972972972973,10.8523908523909,11.5072765072765,11.2733887733888,11.0395010395010,10.5249480249480,9.54261954261954,8.51351351351351,7.53118503118503,6.12785862785863;2.76091476091476,3.25987525987526,3.80873180873181,4.27442827442827,4.38087318087318,4.19126819126819,4.09147609147609,3.54261954261954,3.19334719334719,2.67775467775468,2.44490644490644,2.66112266112266,2.82744282744283,3.07692307692308,3.45945945945946,3.85862785862786,4.09147609147609,4.00831600831601,3.92515592515592,3.74220374220374,3.39293139293139,3.02702702702703,2.67775467775468,2.17879417879418;8.19646569646570,9.67775467775468,11.3071725571726,12.6897089397089,13.0057172557173,12.4428274428274,12.1465696465696,10.5171517671518,9.48024948024948,7.94958419958420,7.25831600831601,7.90020790020790,8.39397089397089,9.13461538461538,10.2702702702703,11.4553014553015,12.1465696465696,11.8996881496881,11.6528066528067,11.1096673596674,10.0727650727651,8.98648648648648,7.94958419958420,6.46829521829522;1.89812889812890,2.24116424116424,2.61850311850312,2.93866943866944,3.01185031185031,2.88149688149688,2.81288981288981,2.43555093555094,2.19542619542620,1.84095634095634,1.68087318087318,1.82952182952183,1.94386694386694,2.11538461538462,2.37837837837838,2.65280665280665,2.81288981288981,2.75571725571726,2.69854469854470,2.57276507276507,2.33264033264033,2.08108108108108,1.84095634095634,1.49792099792100;15.0987525987526,17.8274428274428,20.8290020790021,23.3757796257796,23.9579002079002,22.9209979209979,22.3752598752599,19.3737006237006,17.4636174636175,14.6439708939709,13.3705821205821,14.5530145530146,15.4625779625780,16.8269230769231,18.9189189189189,21.1018711018711,22.3752598752599,21.9204781704782,21.4656964656965,20.4651767151767,18.5550935550936,16.5540540540541,14.6439708939709,11.9152806652807;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;2.76091476091476,3.25987525987526,3.80873180873181,4.27442827442827,4.38087318087318,4.19126819126819,4.09147609147609,3.54261954261954,3.19334719334719,2.67775467775468,2.44490644490644,2.66112266112266,2.82744282744283,3.07692307692308,3.45945945945946,3.85862785862786,4.09147609147609,4.00831600831601,3.92515592515592,3.74220374220374,3.39293139293139,3.02702702702703,2.67775467775468,2.17879417879418;7.50623700623700,8.86278586278586,10.3549896049896,11.6211018711019,11.9104989604990,11.3950103950104,11.1237006237006,9.63149688149688,8.68191268191268,7.28014553014553,6.64708939708940,7.23492723492723,7.68711018711019,8.36538461538461,9.40540540540540,10.4906444906445,11.1237006237006,10.8976091476091,10.6715176715177,10.1741164241164,9.22453222453222,8.22972972972973,7.28014553014553,5.92359667359667;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;3.01975051975052,3.56548856548857,4.16580041580042,4.67515592515592,4.79158004158004,4.58419958419958,4.47505197505198,3.87474012474012,3.49272349272349,2.92879417879418,2.67411642411642,2.91060291060291,3.09251559251559,3.36538461538462,3.78378378378378,4.22037422037422,4.47505197505198,4.38409563409563,4.29313929313929,4.09303534303534,3.71101871101871,3.31081081081081,2.92879417879418,2.38305613305613;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;2.07068607068607,2.44490644490644,2.85654885654886,3.20582120582120,3.28565488565488,3.14345114345114,3.06860706860707,2.65696465696466,2.39501039501039,2.00831600831601,1.83367983367983,1.99584199584200,2.12058212058212,2.30769230769231,2.59459459459459,2.89397089397089,3.06860706860707,3.00623700623701,2.94386694386694,2.80665280665281,2.54469854469854,2.27027027027027,2.00831600831601,1.63409563409563;9.14553014553014,10.7983367983368,12.6164241164241,14.1590436590437,14.5116424116424,13.8835758835759,13.5530145530146,11.7349272349272,10.5779625779626,8.87006237006237,8.09875259875260,8.81496881496882,9.36590436590436,10.1923076923077,11.4594594594595,12.7817047817048,13.5530145530146,13.2775467775468,13.0020790020790,12.3960498960499,11.2390852390852,10.0270270270270,8.87006237006237,7.21725571725572];
par3 = [1;1;2;3;2;2;4;5;6;6;6;6;9;9;4;12;12;12;12;14;16;15;18;19;10;10;10;10;21;15;22;23;24;25;25;27;27;27;27;29;8;6];
par4 = [2;3;4;4;5;6;6;7;7;8;9;10;11;10;12;13;14;15;16;15;17;18;19;20;20;17;21;22;22;23;24;24;25;26;27;28;28;29;30;30;28;28];
par5 = [0.0300000000000000;0.0200000000000000;0.0200000000000000;0;0.0200000000000000;0.0200000000000000;0;0.0100000000000000;0.0100000000000000;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0.0200000000000000;0.0100000000000000];
par6 = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;134.000000000001,164.000000000000,180,180,180,180,180,180,180,151.000000000000,136.999999999997,150.000000000000,160.000000000001,174.999999999999,180,180,180,180,180,180,180,172.000000000000,160.999999999998,130.999999999999;0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1.77635683940025e-15,-1.77635683940025e-15,1.77635683940025e-15,0,0,0,0,0,0,0;10,10,27.0000000000015,35,35,35,35,33.0000000000012,12.0000000000007,10,10,10,10,10,27.9999999999973,35,35,35,35,33.0000000000008,24.0000000000132,10.0000000000000,0,0;10,10,10,30,30,24.9999999999995,19.0000000000002,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-5.55111512312578e-16,-2.66453525910038e-15,-1.99840144432528e-15,-4.44089209850063e-15,-1.33226762955019e-15,0,-3.44613226843649e-13,0,0,0;12.0000000000000,12.0000000000000,12.0000000000000,12.0000000000001,18.3999999999997,12.0000000000000,11.9999999999999,0,0,0,0,0,0,0,0,17.0000000000008,31.0000000000002,25.9999999999996,21.0000000000006,12,0,0,0,0];
par7 = [1;2;22;27;23;13];
options = sdpsettings('solver','gurobi');
a1 = sdpvar(30,24);
a2 = sdpvar(30,24);
a3 = sdpvar(42,24);
a4 = sdpvar(42,24);
a5 = sdpvar(42,24);
a6 = sdpvar(42,24);
a7 = sdpvar(6,24);
a8 = sdpvar(6,24);
a9 = sdpvar(30,24);
a10 = sdpvar(30,24);   
o1 = sdpvar(30,24);
o3 = sdpvar(42,24);
o4 = sdpvar(6,24);
o5 = sdpvar(30,24);
h = sdpvar(1,1);
Con = [];
for t=1:24
    for i=1:42
        Con = [Con,-a1(par3(i),t)+a1(par4(i),t)+a2(par3(i),t)-a2(par4(i),t)+a3(i,t)/100-a4(i,t)/100+a5(i,t)-a6(i,t) == 0];               
        o3(i,t)=-a5(i,t)*par1(i)-a6(i,t)*par1(i);            
    end
    for g=1:6
        Con = [Con,a1(par7(g,1),t)-a2(par7(g,1),t)+a7(g,t)-a8(g,t)<=0];
        o4(g,t)=par6(g,t)*a7(g,t)-par6(g,t)*a8(g,t);                
    end
    for i = 1:30
        from_index = find(par3(i) == i)';
        to_index = find(par4(i) == i)';
        if isempty(from_index)
            from_branch = 0;
        else
            from_branch = -ones(1,length(from_index))*(a3(from_index,t)./par5(from_index))+ones(1,length(from_index))*(a4(from_index,t)./par5(from_index)); 
        end     
        if isempty(to_index)
            to_branch = 0;
        else
            to_branch = ones(1,length(to_index))*(a3(to_index,t)./par5(from_index))-ones(1,length(to_index))*(a4(to_index,t)./par5(from_index));
        end               
        Con = [Con,from_branch+to_branch+a9(i,t)-a10(i,t) == 0];
        Con = [Con,-a1(i,t)+a2(i,t)<=1000];
        Con = [Con,a1(i,t)-a2(i,t)<=1000];    
        o1(i,t)=par2(i,t)*a1(i,t)-par2(i,t)*a2(i,t);
        o5(i,t)=-(1.05*a9(i,t)+1.05*a10(i,t));                
    end
end
Con = [Con,a1>=0,a2>=0];
Con = [Con,a3>=0,a4>=0];   
Con = [Con,a5>=0,a6>=0];      
Con = [Con,a7>=0,a8>=0]; 
Con = [Con,a9>=0,a10>=0];   
%%% constraint needs to pay attention
Con = [Con,h<=sum(sum(o1+o5))+sum(sum(o3))+sum(sum(o4))];
optimize(Con, -h, options);

johanlofberg Sep 11, 2021
Maintainer

Unrecognized function or variable 'gen'.

BYS543 · 2021-09-12T00:41:14Z

BYS543
Sep 12, 2021
Author

Well, I have corrected my code, this can be run without no definition.

0 replies

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

unbounded problem #989

{{title}}

{{editor}}'s edit

{{editor}}'s edit

Replies: 2 comments 2 replies

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

Select a reply

unbounded problem #989

BYS543 Sep 11, 2021

Replies: 2 comments · 2 replies

johanlofberg Sep 11, 2021 Maintainer

BYS543 Sep 11, 2021 Author

johanlofberg Sep 11, 2021 Maintainer

BYS543 Sep 12, 2021 Author

BYS543
Sep 11, 2021

Replies: 2 comments 2 replies

johanlofberg
Sep 11, 2021
Maintainer

BYS543 Sep 11, 2021
Author

johanlofberg Sep 11, 2021
Maintainer

BYS543
Sep 12, 2021
Author