到底我的目标函数和约束条件哪里出了问题?会出现如图的情况
function [g,ceq]=ST1(x) %并网
FWT=[235.540000000000,233,227.880000000000,229.420000000000,227.600000000000,219.530000000000,219.960000000000,221.670000000000,230.850000000000,220.790000000000,215.890000000000,213.690000000000,207.990000000000,200.430000000000,207.070000000000,190.990000000000,193.370000000000,183.890000000000,172.720000000000,172.590000000000,152.090000000000,150.070000000000,139.640000000000,123.150000000000,121.040000000000,109.410000000000,107.970000000000,106.590000000000,95.7400000000000,89.6400000000000,87.2500000000000,89.0800000000000,79.2300000000000,76.8800000000000,74.4100000000000,72.6300000000000,65.1300000000000,65.0100000000000,66.4900000000000,66.6200000000000,66.1700000000000,65.7800000000000,63.1200000000000,63.2400000000000,64.8600000000000,67.1200000000000,64.6800000000000,62.3900000000000];
FPV=[0,0,0,0,3.83000000000000,10.6300000000000,17.2500000000000,25.0600000000000,30.0200000000000,34.7500000000000,50.6800000000000,78.9200000000000,70.8300000000000,106.050000000000,109.380000000000,111.670000000000,128.940000000000,150.830000000000,163.920000000000,170.150000000000,190.970000000000,192.600000000000,183.900000000000,206,198.330000000000,215.900000000000,213.470000000000,229.300000000000,218.650000000000,241.660000000000,228.480000000000,239.060000000000,228.830000000000,232,215.090000000000,219.110000000000,220.760000000000,218.910000000000,198.670000000000,189.940000000000,174.530000000000,155.490000000000,164.740000000000,146.640000000000,121.220000000000,102.630000000000,107.980000000000,70.2600000000000];
FLOAD=[160.500000000000,130.400000000000,80.4000000000000,74.7000000000000,61.4000000000000,70.7700000000000,80.1000000000000,65.8000000000000,63.3000000000000,60.1500000000000,75.6600000000000,62.1500000000000,74.1500000000000,70.3700000000000,76.6600000000000,63.4300000000000,70.4400000000000,62.1500000000000,80.1100000000000,70.1500000000000,92.1500000000000,115.110000000000,108.250000000000,120.230000000000,119.570000000000,160.550000000000,190.780000000000,200.660000000000,240.560000000000,276.780000000000,275.110000000000,269.240000000000,263.660000000000,260.660000000000,255.250000000000,256.200000000000,250.780000000000,245.450000000000,240.330000000000,231.440000000000,215.400000000000,200.750000000000,195.340000000000,185.220000000000,186.560000000000,190.300000000000,194,198];
CBT=0;
for i=1:48
CBT=CBT+x(i);
end
ceq(1)=CBT;
for i=1:48 %蓄电池充放电功率
g(i+48)=-110-x(i);
g(i+336)=x(i)-400;
end
SOC=0;
for i=1:48 %蓄电池每个时间段都要在一定范围内
SOC=SOC+x(i);
g(i+96)=100-SOC;
g(i+384)=SOC-500;
end
for i=1:48 %功率平衡约束
ceq(i+48)=FWT(i)+FPV(i)+x(i)-FLOAD(i);
end
下面是目标函数
function y = GOAL1(x) %并网
FWT=[235.540000000000,233,227.880000000000,229.420000000000,227.600000000000,219.530000000000,219.960000000000,221.670000000000,230.850000000000,220.790000000000,215.890000000000,213.690000000000,207.990000000000,200.430000000000,207.070000000000,190.990000000000,193.370000000000,183.890000000000,172.720000000000,172.590000000000,152.090000000000,150.070000000000,139.640000000000,123.150000000000,121.040000000000,109.410000000000,107.970000000000,106.590000000000,95.7400000000000,89.6400000000000,87.2500000000000,89.0800000000000,79.2300000000000,76.8800000000000,74.4100000000000,72.6300000000000,65.1300000000000,65.0100000000000,66.4900000000000,66.6200000000000,66.1700000000000,65.7800000000000,63.1200000000000,63.2400000000000,64.8600000000000,67.1200000000000,64.6800000000000,62.3900000000000];
FPV=[0,0,0,0,3.83000000000000,10.6300000000000,17.2500000000000,25.0600000000000,30.0200000000000,34.7500000000000,50.6800000000000,78.9200000000000,70.8300000000000,106.050000000000,109.380000000000,111.670000000000,128.940000000000,150.830000000000,163.920000000000,170.150000000000,190.970000000000,192.600000000000,183.900000000000,206,198.330000000000,215.900000000000,213.470000000000,229.300000000000,218.650000000000,241.660000000000,228.480000000000,239.060000000000,228.830000000000,232,215.090000000000,219.110000000000,220.760000000000,218.910000000000,198.670000000000,189.940000000000,174.530000000000,155.490000000000,164.740000000000,146.640000000000,121.220000000000,102.630000000000,107.980000000000,70.2600000000000];
FLOAD=[160.500000000000,130.400000000000,80.4000000000000,74.7000000000000,61.4000000000000,70.7700000000000,80.1000000000000,65.8000000000000,63.3000000000000,60.1500000000000,75.6600000000000,62.1500000000000,74.1500000000000,70.3700000000000,76.6600000000000,63.4300000000000,70.4400000000000,62.1500000000000,80.1100000000000,70.1500000000000,92.1500000000000,115.110000000000,108.250000000000,120.230000000000,119.570000000000,160.550000000000,190.780000000000,200.660000000000,240.560000000000,276.780000000000,275.110000000000,269.240000000000,263.660000000000,260.660000000000,255.250000000000,256.200000000000,250.780000000000,245.450000000000,240.330000000000,231.440000000000,215.400000000000,200.750000000000,195.340000000000,185.220000000000,186.560000000000,190.300000000000,194,198];
FBT=0;
for i=1:48 %1~48为蓄电池,正为向微电网放电,负为充电
FBT=FBT+0.055*0.5*2.71828*abs(x(i));
end
y=28307+FBT