THETA_t = [zeros(m,d) zeros(m,d) zeros(m,d) zeros(m,d); sdpvar(m,d) zeros(m,d) zeros(m,d) zeros(m,d);...
sdpvar(m,d) sdpvar(m,d) zeros(m,d) zeros(m,d); sdpvar(m,d) sdpvar(m,d) sdpvar(m,d) sdpvar(m,d)];
ETA_t = [sdpvar; sdpvar; sdpvar; sdpvar];
constraints1 = [];
constraints1 = [constraints1, abs(ETA_t(i,:))+norm(THETA_t(i,:),1)*phi_max <= Umax]; ---- eq(17)