#毕业论文# 改了一上午了，求不出解。

for i = 1:N
    for t=1:24
    if t<ts(i) | t>td(i)
        con = con+[Pev(t,i)==0,SOC(t,i)==0];
    end
    end
end
for i = 1:N
SOC(ts(i),i)=SOC0(i);
end
for i=1:N
for t=(ts(i)+1):td(i)
SOC(t,i)=SOC(t-1,i)+(Pev(t,i)/Eev);
con = con + [SOC(td,i) >= SOCend(i)];
end
end
objective=(Cgrid')*Pev*1;
options=sdpsettings('verbose',0);
sol=solvesdp(con,objective,options);
if sol.problem==0
    solution = value(x)
else
    disp('***kkk');
   sol.info
   yalmiperror(sol.problem)
end