2004B数模竞赛 电力分配预案

Capacity=xlsread('C:\Users\Administrator\Desktop\数模代码合集\数模第二次优秀论文研读\第三题\各机组段容量.xls');
Price=xlsread('C:\Users\Administrator\Desktop\数模代码合集\数模第二次优秀论文研读\第三题\各机组的段价.xls');
Speed=xlsread('C:\Users\Administrator\Desktop\数模代码合集\数模第二次优秀论文研读\第三题\各机组的爬坡速率.xls');
format short g;
[n,m]=size(Capacity); %n=8个机组 m=10段
PL=982.4;  %负荷需求
art=zeros(n,1);  %蚂蚁标识符
sum=0;  %和值
g=zeros(n,1); %g对应每个机组的段容量和
part=zeros(8,1); %完整与否
while(1)
    MIN=Price(1,art(1)+1);  %每次循环默认第一个为最小值
    for i=2:n
        if art(i)==10
            continue; %蚂蚁走到底
        end
        if Price(i,art(i)+1)<MIN
            MIN=Price(i,art(i)+1); %逐行检索最小值
        end
    end  %找到最小值MIN
    [I,J]=find(Price(:,:)==MIN); %确定MIN位置
    if size(I,1)~=1  %有多个MIN
        I=I(1); 
        J=J(1); %记录第一个MIN
        Price(I,J)=inf;  %暂时把MIN值改为inf 以免影响后面的循环
    end
    art(I)=art(I)+1;  %蚂蚁标识符对应位置+1
    sum=sum+Capacity(I,J); %求和
    if sum<PL
        g(I)=g(I)+Capacity(I,J); %记录每个机组的段容量

    elseif sum==PL
        g(I)=g(I)+Capacity(I,J);
        disp('sum=PL');
        break;  
    else
        g(I)=g(I)+(Capacity(I,J)-sum+PL); %调整sum和g确保sum=PL
        sum=PL;
        art(I)=art(I)-1; %第一遍的bug 最后一个应该减去1
        part(I)=1; %不完整最后一个
        break;  %出循环条件
    end
end
%---------------step1---------------------
Price=xlsread('C:\Users\Administrator\Desktop\数模代码合集\数模第二次优秀论文研读\第三题\各机组的段价.xls');
g0=[120 73 180 80 125 125 81.1 90]'; %初始g0
v=Speed.*60; %各机组的爬坡速率矩阵
t=1/4; %十五分钟为一个交易周期
g2=zeros(8,1);%辅助矩阵
flag=zeros(8,1);
while(1)
    for i=1:n
        if g(i)-g0(i)>=v(i)*t %爬不到这么高的情况
            g(i)=g0(i)+v(i)*t; %新的最大值g(i)
            g2(i)=g(i);%辅助矩阵来确定位置
            for j=1:10 %对每个段容量的循环
                g2(i)=g2(i)-Capacity(i,j); %依次减
                if g2(i)==0 %刚好等于0 
                    art(i)=j;
                    break;
                elseif g2(i)<0 %小于0，即最后一段只取了一部分，art(i)=j-1