......
Pinj=sdpvar(N,1);
Wgij=binvar(N,Ng,);
Wlkj=binvar(N,Ns);
Ptcw=sdpvar(N,Nm,);
Ptcl=sdpvar(N,Nm);
for s=1:Nm
for j=1:N
for i=1:Ng
B1(j,i,s)=ps(1,s)*(ai(1,i)*Psij(j,i,s)^2+bi(1,i)*Psij(j,i,s)+ci(1,i)*Wgij(j,i)); error:Incorrect use of sdpvar / subsasgn ( line 248 ) Only full matrices are allowed to use n-dimensional indexes.
end
for k=1:Ns
B2(j,k,s)=ps(1,s)*(Wlkj(j,k,s)*a);
end
B3(j,s)=ps(1,s)*(Ptcw(j,s)*kw+Ptcl(j,s)*kl)
end
C1(1,s)=sum(sum(B1));
C2(1,s)=sum(sum(B2));
end
C=sum(C1)+sum(C2)+sum(sum(B3));
......