这个477是如何来的

Question

c1=[0.04 0.04 0.04 0.35 0.35 0.26 1.2 1.2 1 0 0.0001 1.6 1.6 3.5 0 0 0 0 0.0001 0]';%1*20
x_L1=[0 0 0 0 0 0 0 0 0  0 0 0 0 0 0 0 0 16*2 0 0]';%20*1
x_U1=[25 25 55 30 30 65 1 1 1 2*20/0.9 2*20/0.9 1 1 1 1 1 1 2*80 2*20/0.9 200]';%20*1
PD=[30 31 32 35 32 50 75 100 120 130 150 150 145 135 135 120 140 110 105 80 60 50 30 30];
%生成A阵
A1=[0 0 0 1 1 1 0 0 0 -1 0.9 0 0 0 0 0 0 0 0 -1;
   0 0 0 0 0 0 0 0 0 0.9 -1 0 0 0 0 0 0 -1 0 0; %17*20 %行数加一
    0 0 0 1 0 0 -30 0 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 -1 0 0 5 0 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 1 0 0 -30 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 -1 0 0 5 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 1 0 0 -65  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 -1 0 0 10  0 0 0 0 0 0 0 0 0 0 0;
   1 0 0 1 0 0 -30 0 0  0 0 0 0 0 0 0 0 0 0 0;
   -1 0 0 -1 0 0 5 0 0  0 0 0 0 0 0 0 0 0 0 0;
   0 1 0 0 1 0 0 -30 0  0 0 0 0 0 0 0 0 0 0 0;
   0 -1 0 0 -1 0 0 5 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 1 0 0 1 0 0 -65  0 0 0 0 0 0 0 0 0 0 0;
   0 0 -1 0 0 -1 0 0 10  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 1 0 0 0 0 -1 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 1 0 0 0 0 -1 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 1 0 0 0 0 -1 0 0 0 0 0 0;
   ];

A2=[0 0 0 1 1 1 0 0 0 -1 0.9 0 0 0 0 0 0 0 0 -1;%20*20
    0 0 0 0 0 0 0 0 0 0.9 -1 0 0 0 0 0 0 -1 0 0;
    0 0 0 1 0 0 -30 0 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 -1 0 0 5 0 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 1 0 0 -30 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 -1 0 0 5 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 1 0 0 -65  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 -1 0 0 10  0 0 0 0 0 0 0 0 0 0 0;
   1 0 0 1 0 0 -30 0 0  0 0 0 0 0 0 0 0 0 0 0;
   -1 0 0 -1 0 0 5 0 0  0 0 0 0 0 0 0 0 0 0 0;
   0 1 0 0 1 0 0 -30 0  0 0 0 0 0 0 0 0 0 0 0;%x=[r1 r2 r3 p1 p2 p3 i1 i2 i3 pin pout su1 su2 su3 sd1 sd2 sd3 soc];1*18
   0 -1 0 0 -1 0 0 5 0  0 0 0 0 0 0 0 0 0 0 0;
   0 0 1 0 0 1 0 0 -65  0 0 0 0 0 0 0 0 0 0 0;
   0 0 -1 0 0 -1 0 0 10  0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0;
   0 0 0 0 0 0 -1 0 0 0 0 1 0 0 -1 0 0 0 0 0;
   0 0 0 0 0 0 0 -1 0 0 0 0 1 0 0 -1 0 0 0 0;  %停开机0-1变量
   0 0 0 0 0 0 0 0 -1 0 0 0 0 1 0 0 -1 0 0 0;  %i(t-1)没考虑
   ];
B1=[A1 zeros(17,20*23)];
B2=[zeros(20,20) A2 zeros(20,20*22)];
B3=[zeros(20,20*2) A2 zeros(20,20*21)];
B4=[zeros(20,20*3) A2 zeros(20,20*20)];
B5=[zeros(20,20*4) A2 zeros(20,20*19)];
B6=[zeros(20,20*5) A2 zeros(20,20*18)];
B7=[zeros(20,20*6) A2 zeros(20,20*17)];
B8=[zeros(20,20*7) A2 zeros(20,20*16)];
B9=[zeros(20,20*8) A2 zeros(20,20*15)];
B10=[zeros(20,20*9) A2 zeros(20,20*14)];
B11=[zeros(20,20*10) A2 zeros(20,20*13)];
B12=[zeros(20,20*11) A2 zeros(20,20*12)];
B13=[zeros(20,20*12) A2 zeros(20,20*11)];
B14=[zeros(20,20*13) A2 zeros(20,20*10)];
B15=[zeros(20,20*14) A2 zeros(20,20*9)];
B16=[zeros(20,20*15) A2 zeros(20,20*8)];
B17=[zeros(20,20*16) A2 zeros(20,20*7)];
B18=[zeros(20,20*17) A2 zeros(20,20*6)];
B19=[zeros(20,20*18) A2 zeros(20,20*5)];
B20=[zeros(20,20*19) A2 zeros(20,20*4)];
B21=[zeros(20,20*20) A2 zeros(20,20*3)];
B22=[zeros(20,20*21) A2 zeros(20,20*2)];
B23=[zeros(20,20*22) A2 zeros(20,20*1)];
B24=[zeros(20,20*23) A2];
A=[B1;B2;B3;B4;B5;B6;B7;B8;B9;B10;B11;B12;B13;B14;B15;B16;B17;B18;B19;B20;B21;B22;B23;B24];
[m,n] = size(A);

for i=1:23                  %Soc(t-1)+Pin(t)-Pout(t)=Soc(t);
    A(17+20*(i-1)+2,20*i-2)=1;
end
for k=1:23
    i=17+(k-1)*20;%su-sd-i(t)+i(t-1)=0;
    j=(k-1)*20;
    A(i+18,j+7)=1;
    A(i+19,j+8)=1;
     A(i+20,j+9)=1;
end
i=0;
j=0;
for i=1:24
    j=j+Npv(i)+Nw(i)-1;
end
pronum=j;
% S=zeros(823,730);
S=[A,zeros(477,j)];
p=0;
q=0;%新加的元素坐标
k=[];  
for p=1:24
    k(p)=0;
    for q=1:p
        k(p)=k(p)+Npv(q)+Nw(q)-1;%减过1
    end
end
p=1;
% V=[];
    for q=1:pronum
            if(q>k(p))
              p=p+1;  
            end
        S(477+q,1+(p-1)*20)=-1/L;%备用设置较为复杂,储能还未考虑
        S(477+q,2+(p-1)*20)=-1/L;
        S(477+q,3+(p-1)*20)=-1/L;
        S(477+q,20*p-1)=-0.9*1/L;
        S(477+q,24*20+q)=1;      
    end

Answer 1

matlab为什么给个python标签？话说你这么贴代码其他人看不懂的吧