电力系统经济运行及市场竞价Matlab.zip

function ans =jiaoyanjieguo(jizucanshu,D,A,geidingfuhe,time)%初步结果进行验证,使各机组功率均在范围内, %输出为11行1列的数组,数据包括各机组功率和总耗量 h=0;%标志位,当有越限时,h值为越限机组中a值最大的机组号 l=0;%所有越限机组中最大的a值 m=0;%a值最大的机组的边界值 for i=1:10%对第一到第十的机组验证,若有越限,找出a值最大的机组,把机组号赋值给h,最大的斜率a赋值给l,边界值赋值给m if D(i,1)<jizucanshu(i,5)%判断是否小于最小值 if jizucanshu(i,2)>l l=jizucanshu(i,2); h=i; m=jizucanshu(i,5); end else if D(i,1)>jizucanshu(i,6)%判断是否大于最大值 if jizucanshu(i,2)>l l=jizucanshu(i,2); h=i; m=jizucanshu(i,6); end end end end if h>0 %有越限,更改A矩阵中a值最大的越限机组座对应的行参数 A(h,11)=0; A(h,h)=-jizucanshu(h,3)/m; D=chubujieguo(jizucanshu,geidingfuhe,A,time); %重新求得各机组的功率 jiaoyanjieguo(jizucanshu,D,A,geidingfuhe,time);%再次进行验证 else %上述循环后至无越限,G存放总耗量,最后连同各机组功率一起存放至E中,即E为输出结果 E=zeros(11,1);%存放输出结果 G=zeros(1,1);%存放总耗量 for i=1:10 G(1,1)=G(1,1)+jizucanshu(i,2)*D(i,1)^2+jizucanshu(i,3)*D(i,1)+jizucanshu(i,4);%求总耗量 end E=[D;G];%输出结果 ans=E; end end