function [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL)
RestraintCount=size(PVi,1)+size(PGi,1)+Busnum; %约束条件数
t_Bal_volt=Volt(Balance);
Volt=sparse(ones(1,Busnum));
%Volt(Balance)=t_Bal_volt;
Volt(Balance)=1;
UAngel=sparse(1,Busnum);
Init_Z=sparse(ones(1,RestraintCount));
Init_W=sparse(-1*ones(1,RestraintCount));
Init_L=sparse(ones(1,RestraintCount));
Init_U=sparse(ones(1,RestraintCount));
Init_Y=sparse(1,2*Busnum);%与学姐一致
tPU=sparse(GenU(:,2));% 发电机有功上限
tQU=sparse(PVQU(:,1));% 无功上限
tPL=sparse(GenL(:,2));% 发电机有功下限
tQL=sparse(PVQL(:,1));% 无功下限
PG(PGi)=(tPU+tPL)/2;
QG(PVi)=(tQU+tQL)/2;
end