Here, I also paste my code:
clc
clear all
close all
K = 4;
E = 3;
wk=[-100 0 50 180; -100 100 0 180];
we=[-180 20 160;0 170 -100];
dh=100;
V_max=50;
T = 5;
delta_t=0.5;
N = T/delta_t;
luo_0 = 10^(-60/10);
delta = 1e-6;
epsluo_u = 1e-14;
siga = 10000;
k_oula = 0.58;
f_max = 510^9;
fTD_max = 310^9;
C_k = 10^3;
k_TD=1e-27;
Pk_max = 0.00110^(30/10);
Pu_J_max = 0.00110^(30/10);
H_J=60;
H_H=100;
kesi = 0.5;
qr_J=zeros(2,N);
qr_H=zeros(2,N);
sitaslot=linspace(0,2pi,N);
r_H=20;
r_J=20;
qr_J = [r_Jsin(sitaslot);r_Jcos(sitaslot)];
qr_H = [r_Hsin(sitaslot+pi);r_Hcos(sitaslot+pi)];
p_kr = 0.01Pk_maxones(K,N);
p_jr = 0.01Pu_J_maxones(1,N);
alpha_km_r = zeros(KK,N);
D_khr = zeros(K,N);
D_jer = zeros(E,N);
D_jhr = zeros(1,N);
D_ke = zeros(K,E);
for n=1:N
for k=1:K
D_khr(k,n) = H_H^2 + norm(qr_H(:,n)-wk(:,k)).^2;
end
for e=1:E
D_jer(e,n) = H_J^2 + norm(qr_J(:,n)-we(:,e)).^2;
end
D_jhr(n) = (H_J-H_H)^2 + norm(qr_J(:,n)-qr_H(:,n)).^2;
end
belta_k_r = D_khr;
for k=1:K
for e=1:E
D_ke(k,e) = norm(wk(:,k)-we(:,e))^2;
end
end
D_mhr = D_khr;
for k=1:K
for n=1:N
alpha_km_r((k-1)K+1 : Kk, n) = (D_khr(:,n) > D_khr(k,n));
end
end
alpha_mk_r = 1 - alpha_km_r;
alpha_ke1 = zeros(KK,N);
alpha_ke2 = zeros(KK,N);
alpha_ke3 = zeros(KK,N);
g_ke = zeros(KE,N);
for e=1:E
for k=1:K
if e==1
for n=1:N
alpha_ke1((k-1)K+1 : Kk, n) = (D_ke(:,e) > D_ke(k,e));
end
end
if e==2
for n=1:N
alpha_ke2((k-1)K+1 : Kk, n) = (D_ke(:,e) > D_ke(k,e));
end
end
if e==3
for n=1:N
alpha_ke3((k-1)K+1 : Kk, n) = (D_ke(:,e) > D_ke(k,e));
end
end
end
end
for k=1:K
for e=1:E
for n=1:N
g_ke(k + (e-1)K, n) = luo_0D_ke(k,e)^(-1.5)*kesi;
end
end
end
h_khr = luo_0./D_khr;
h_mhr = h_khr;
h_jhr = luo_0./D_jhr;
h_jer = luo_0./D_jer;
pai_1kr = zeros(K,N);
temp1 = p_kr.*h_mhr;
for k=1:K
if k==1
temp2=alpha_km_r(1:K,:).*temp1;
pai_1kr(1,
= sum(temp2)+ h_jhr.p_jr + epsluo_u;
end
if k==2
temp3=alpha_km_r(K+1:2K,:).temp1;
pai_1kr(2,
= sum(temp3)+ h_jhr.p_jr + epsluo_u;
end
if k==3
temp4=alpha_km_r(2K+1:3K,:).temp1;
pai_1kr(3,
= sum(temp4) + h_jhr.p_jr + epsluo_u;
end
if k==4
temp5=alpha_km_r(3K+1:4K,:).*temp1;
pai_1kr(4,
= sum(temp5)+ h_jhr.*p_jr + epsluo_u;
end
end
gamma_khr = p_kr.*h_khr./(pai_1kr);
t_kr = gamma_khr;
bar_pai_1kr = D_khr;
pai_2kr = t_kr.*pai_1kr;
bar_pai_2mr = p_kr.*h_khr;
bar_pai_3jh_r = p_jr.*h_jhr;
yita_1mr=D_khr;
yita_2mr=D_jhr;
bar_t_ke1_r = zeros(K,N);
bar_t_ke2_r = zeros(K,N);
bar_t_ke3_r = zeros(K,N);
Kesi_ke1_r = zeros(K,N);
Kesi_ke2_r = zeros(K,N);
Kesi_ke3_r = zeros(K,N);
bar_Kesi_je_r = zeros(E,N);
for k=1:K
bar_t_ke1_r(k,
= p_kr(k,:).*g_ke(k,:)./(sum( p_kr.*alpha_ke1((k-1)K+1:kK,:).*g_ke(1:K,
) + p_jr.*h_jer(1,
+ epsluo_u); %for e1
Kesi_ke1_r(k,
= sum( p_kr.*alpha_ke1((k-1)K+1:kK,:).*g_ke(1:K,
) + p_jr.*h_jer(1,
+ epsluo_u;
bar_t_ke2_r(k,
= p_kr(k,:).*g_ke(k+K,:)./(sum( p_kr.*alpha_ke2((k-1)K+1:kK,:).g_ke(K+1:2K,
) + p_jr.*h_jer(2,
+ epsluo_u);%for e2
Kesi_ke2_r(k,
= sum( p_kr.*alpha_ke2((k-1)K+1:kK,:).g_ke(K+1:2K,
) + p_jr.*h_jer(2,
+ epsluo_u;
bar_t_ke3_r(k,
= p_kr(k,:).g_ke(k+2K,:)./(sum( p_kr.alpha_ke3((k-1)K+1:kK,:).g_ke(2K+1:3K,
) + p_jr.*h_jer(3,
+ epsluo_u);%for e3
Kesi_ke3_r(k,
= sum( p_kr.alpha_ke3((k-1)K+1:kK,:).g_ke(2K+1:3K,
) + p_jr.*h_jer(3,
+ epsluo_u;
end
for e=1:E
bar_Kesi_je_r(e,
= p_jr.h_jer(e,:);
end
tau_je_r = D_jer;
alpha_km_wave = zeros(KK,N);
lamada1_km = 1ones(KK,N);
lamada2_km = 1ones(KK,N);
lamada3_km = 1ones(KK,N);
for k=1:K
alpha_km_wave= (alpha_km_r+alpha_km_r.^2+siga*lamada1_km+siga*lamada2_km.*alpha_km_r)./(1 + alpha_km_r.^2); % step 1: solve lpha_k_wave by fixed other variables
end
cvx_solver Mosek
cvx_begin
variable Sk
variable sita_1k_opt(K,N)
variable sita_2k_opt(K,N)
variable sita_3k_opt(K,N)
variable tk_opt(K,N)
variable pj_opt(1,N)
variable pk_opt(K,N)
variable bar_tke_opt(K*E,N)
variable alpha_km_opt(K*K,N)
variable fk_opt(K,N)
variable qH(2,N)
variable qJ(2,N)
variable pai_1k_opt(K,N)
variable pai_2k_opt(K,N)
variable bar_pai_1k_opt(K,N)
variable bar_pai_2m_opt(K,N)
variable bar_pai_3jh_opt(1,N)
variable yita_1m_opt(K,N)
variable yita_2m_opt(1,N)
variable Rk_loc_opt(K,N)
variable Kesi_ke_opt(K*E,N)
variable tau_ke_opt(E,N)
variable bar_Kesi_je_opt(E,N)
variable belta_k_opt(K,N)
expression alpha_mk_opt(K*K,N)
alpha_mk_opt = 1-alpha_km_opt;
maximize (Sk)
subject to
Sk >= 0;
1/N*sum(sita_1k_opt + Rk_loc_opt,2) >= Sk;
sita_1k_opt <= sita_2k_opt - sita_3k_opt;
-sita_2k_opt >= rel_entr_quad(1,1 + tk_opt)/log(2);
sita_3k_opt >= log2(1+bar_t_ke1_r) + 1/log(2)./(1+bar_t_ke1_r).*(bar_tke_opt((1-1)*K+1:1*K,:)-bar_t_ke1_r);
sita_3k_opt >= log2(1+bar_t_ke2_r) + 1/log(2)./(1+bar_t_ke2_r).*(bar_tke_opt((2-1)*K+1:2*K,:)-bar_t_ke2_r);
sita_3k_opt >= log2(1+bar_t_ke3_r) + 1/log(2)./(1+bar_t_ke3_r).*(bar_tke_opt((3-1)*K+1:3*K,:)-bar_t_ke3_r);
for n=1:N
for k=1:K
if k==1 % k = 1
i=[2 3 4];
norm([sum((alpha_km_r(i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(i,n)-bar_pai_2m_opt(i,n))/4 - (alpha_km_r(i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 - epsluo_u/2 - 1/2 + pai_1k_opt(k,n)/2, ...
(alpha_km_opt(2,n)+bar_pai_2m_opt(2,n))/2, (alpha_km_opt(3,n)+bar_pai_2m_opt(3,n))/2, (alpha_km_opt(4,n)+bar_pai_2m_opt(4,n))/2 ]) <= sum((alpha_km_r(i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(i,n)-bar_pai_2m_opt(i,n))/4 ...
- (alpha_km_r(i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 + pai_1k_opt(k,n)/2- epsluo_u/2 + 1/2;
end
if k==2
i=[1 3 4];
norm([sum((alpha_km_r(K+i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(K+i,n)-bar_pai_2m_opt(i,n))/4 - (alpha_km_r(K+i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 - epsluo_u/2 - 1/2 + pai_1k_opt(k,n)/2, ...
(alpha_km_opt(5,n)+bar_pai_2m_opt(1,n))/2, (alpha_km_opt(7,n)+bar_pai_2m_opt(3,n))/2, (alpha_km_opt(8,n)+bar_pai_2m_opt(4,n))/2 ]) <= sum((alpha_km_r(K+i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(K+i,n)-bar_pai_2m_opt(i,n))/4 ...
- (alpha_km_r(K+i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 + pai_1k_opt(k,n)/2 - epsluo_u/2 + 1/2 ;
end
if k==3
i=[1 2 4];
norm([sum((alpha_km_r(2*K+i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(2*K+i,n)-bar_pai_2m_opt(i,n))/4 - (alpha_km_r(2*K+i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 - epsluo_u/2 - 1/2 + pai_1k_opt(k,n)/2, ...
(alpha_km_opt(9,n)+bar_pai_2m_opt(1,n))/2, (alpha_km_opt(10,n)+bar_pai_2m_opt(2,n))/2, (alpha_km_opt(12,n)+bar_pai_2m_opt(4,n))/2 ]) <= sum((alpha_km_r(2*K+i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(2*K+i,n)-bar_pai_2m_opt(i,n))/4 ...
- (alpha_km_r(2*K+i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 + pai_1k_opt(k,n)/2 - epsluo_u/2 + 1/2;
end
if k==4
i=[1 2 3];
norm([sum((alpha_km_r(3*K+i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(3*K+i,n)-bar_pai_2m_opt(i,n))/4 - (alpha_km_r(3*K+i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 - epsluo_u/2 - 1/2 + pai_1k_opt(k,n)/2, ...
(alpha_km_opt(13,n)+bar_pai_2m_opt(1,n))/2, (alpha_km_opt(14,n)+bar_pai_2m_opt(2,n))/2, (alpha_km_opt(15,n)+bar_pai_2m_opt(3,n))/2 ]) <= sum((alpha_km_r(3*K+i,n)-bar_pai_2mr(i,n)).*(alpha_km_opt(3*K+i,n)-bar_pai_2m_opt(i,n))/4 ...
- (alpha_km_r(3*K+i,n)-bar_pai_2mr(i,n)).^2/8)-bar_pai_3jh_opt(n)/2 + pai_1k_opt(k,n)/2 - epsluo_u/2 + 1/2 ;
end
end
end
%
for n=1:N
for k=1:K
norm([(tk_opt(k,n)+pai_1k_opt(k,n))/2,(t_kr(k,n)-pai_1kr(k,n))(tk_opt(k,n)-pai_1k_opt(k,n))/4 - (t_kr(k,n)-pai_1kr(k,n))^2/8 + pai_2k_opt(k,n)/2 - 1/2 ]) <=…
(t_kr(k,n)-pai_1kr(k,n))(tk_opt(k,n)-pai_1k_opt(k,n))/4 - (t_kr(k,n)-pai_1kr(k,n))^2/8 + pai_2k_opt(k,n)/2 + 1/2;
norm([(pai_2k_opt(k,n)+bar_pai_1k_opt(k,n))/2,(pai_2kr(k,n)-bar_pai_1kr(k,n))*(pai_2k_opt(k,n)-bar_pai_1k_opt(k,n))/4 - (pai_2kr(k,n)-bar_pai_1kr(k,n))^2/8 + pk_opt(k,n)*luo_0/2 - 1/2 ]) <=...
(pai_2kr(k,n)-bar_pai_1kr(k,n))*(pai_2k_opt(k,n)-bar_pai_1k_opt(k,n))/4 - (pai_2kr(k,n)-bar_pai_1kr(k,n))^2/8 + pk_opt(k,n)*luo_0/2 + 1/2;
norm([(bar_pai_2m_opt(k,n)-yita_1m_opt(k,n))/2,(bar_pai_2mr(k,n)+yita_1mr(k,n))*(bar_pai_2m_opt(k,n)+yita_1m_opt(k,n))/4 - (bar_pai_2mr(k,n)+yita_1mr(k,n))^2/8 - pk_opt(k,n)*luo_0/2 - 1/2 ]) <=...
(bar_pai_2mr(k,n)+yita_1mr(k,n))*(bar_pai_2m_opt(k,n)+yita_1m_opt(k,n))/4 - (bar_pai_2mr(k,n)+yita_1mr(k,n))^2/8 - pk_opt(k,n)*luo_0/2 + 1/2;
end
end
%
%
for n=1:N
for k=1:K
for e=1:E
if e==1
norm([(Kesi_ke_opt(k,n)-bar_tke_opt(k,n))/2,(Kesi_ke1_r(k,n)+bar_t_ke1_r(k,n))*(Kesi_ke_opt(k,n)+bar_tke_opt(k,n))/4 - (Kesi_ke1_r(k,n)+bar_t_ke1_r(k,n))^2/8 - pk_opt(k,n)g_ke(k,n)/2 - 1/2 ]) <=…
(Kesi_ke1_r(k,n)+bar_t_ke1_r(k,n))(Kesi_ke_opt(k,n)+bar_tke_opt(k,n))/4 - (Kesi_ke1_r(k,n)+bar_t_ke1_r(k,n))^2/8 - pk_opt(k,n)*g_ke(k,n)/2 + 1/2;
Kesi_ke_opt((e-1)*K+k,n) <= sum(alpha_ke1((k-1)*K+1:k*K,n).*g_ke(1:K,n)) + bar_Kesi_je_opt(e,n) + epsluo_u;
end
if e==2
norm([(Kesi_ke_opt(K+k,n)-bar_tke_opt(K+k,n))/2,(Kesi_ke2_r(k,n)+bar_t_ke2_r(k,n))*(Kesi_ke_opt(K+k,n)+bar_tke_opt(K+k,n))/4 - (Kesi_ke2_r(k,n)+bar_t_ke2_r(k,n))^2/8 - pk_opt(k,n)*g_ke(K+k,n)/2 - 1/2 ]) <=...
(Kesi_ke2_r(k,n)+bar_t_ke2_r(k,n))*(Kesi_ke_opt(K+k,n)+bar_tke_opt(K+k,n))/4 - (Kesi_ke2_r(k,n)+bar_t_ke2_r(k,n))^2/8 - pk_opt(k,n)*g_ke(K+k,n)/2 + 1/2;
Kesi_ke_opt((e-1)*K+k,n) <= sum(alpha_ke2((k-1)*K+1:k*K,n).*g_ke(1+K:2*K,n)) + bar_Kesi_je_opt(e,n) + epsluo_u;
end
if e==3
norm([(Kesi_ke_opt(2*K+k,n)-bar_tke_opt(2*K+k,n))/2,(Kesi_ke3_r(k,n)+bar_t_ke3_r(k,n))*(Kesi_ke_opt(2*K+k,n)+bar_tke_opt(2*K+k,n))/4 - (Kesi_ke3_r(k,n)+bar_t_ke3_r(k,n))^2/8 - pk_opt(k,n)*g_ke(2*K+k,n)/2 - 1/2 ]) <=...
(Kesi_ke3_r(k,n)+bar_t_ke3_r(k,n))*(Kesi_ke_opt(2*K+k,n)+bar_tke_opt(2*K+k,n))/4 - (Kesi_ke3_r(k,n)+bar_t_ke3_r(k,n))^2/8 - pk_opt(k,n)*g_ke(2*K+k,n)/2 + 1/2;
Kesi_ke_opt((e-1)*K+k,n) <= sum(alpha_ke3((k-1)*K+1:k*K,n).*g_ke(1+2*K:3*K,n)) + bar_Kesi_je_opt(e,n) + epsluo_u;
end
end
end
end
for n=1:N
norm([(bar_pai_3jh_opt(n) - yita_2m_opt(n))/2,(bar_pai_3jh_r(n)+yita_2mr(n))*(bar_pai_3jh_opt(n)+yita_2m_opt(n))/4 - (bar_pai_3jh_r(n)+yita_2mr(n))^2/8 - pj_opt(n)*luo_0/2 - 1/2 ]) <=...
(bar_pai_3jh_r(n)+yita_2mr(n))*(bar_pai_3jh_opt(n)+yita_2m_opt(n))/4 - (bar_pai_3jh_r(n)+yita_2mr(n))^2/8 - pj_opt(n)*luo_0/2 + 1/2;
for e=1:E
norm([(bar_Kesi_je_opt(e,n)+tau_ke_opt(e,n))/2,(bar_Kesi_je_r(e,n)-tau_je_r(e,n))*(bar_Kesi_je_opt(e,n)-tau_ke_opt(e,n))/4 - (bar_Kesi_je_r(e,n)-tau_je_r(e,n))^2/8 + pj_opt(n)*luo_0/2 - 1/2 ]) <=...
(bar_Kesi_je_r(e,n)-tau_je_r(e,n))*(bar_Kesi_je_opt(e,n)-tau_ke_opt(e,n))/4 - (bar_Kesi_je_r(e,n)-tau_je_r(e,n))^2/8 + pj_opt(n)*luo_0/2 + 1/2;
end
end
for n=1:N
for k=1:K
norm([H_H,(qH(:,n)-wk(:,k))', (bar_pai_1k_opt(k,n)-1)/2 ]) <= (bar_pai_1k_opt(k,n)+1)/2;
yita_1m_opt(k,n) <= H_H^2 + norm(qr_H(:,n)-wk(:,k))^2 + 2*(qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n));
end
end
for n=1:N
for k=1:K
if k==1
Tem=[2 3 4];
for index=1:3
ii=Tem(index);
norm([(belta_k_opt(k,n)+alpha_km_opt((k-1)*K+ii,n))/2,(belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 - 1/2 ]) <= belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n)*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 + 1/2;
end
end
if k==2
Tem=[1 3 4];
for index=1:3
ii=Tem(index);
norm([(belta_k_opt(k,n)+alpha_km_opt((k-1)*K+ii,n))/2,(belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 - 1/2 ]) <= belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n)*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 + 1/2;
end
end
if k==3
Tem=[1 2 4];
for index=1:3
ii=Tem(index);
norm([(belta_k_opt(k,n)+alpha_km_opt((k-1)*K+ii,n))/2,(belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 - 1/2 ]) <= belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n)*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 + 1/2;
end
end
if k==4
Tem=[1 2 3];
for index=1:3
ii=Tem(index);
norm([(belta_k_opt(k,n)+alpha_km_opt((k-1)*K+ii,n))/2,(belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 - 1/2 ]) <= belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n)*(belta_k_opt(k,n)-alpha_km_opt((k-1)*K+ii,n))/4 - (belta_k_r(k,n)-alpha_km_r((k-1)*K+ii,n))^2/8 + (qr_H(:,n)-wk(:,k))'*(qH(:,n)-qr_H(:,n)) +...
(norm(qr_H(:,n)-wk(:,k))^2)/2 + H_H^2/2 + 1/2;
end
end
end
end
for n=1:N
for k=1:K
norm([H_H,(qH(:,n)-wk(:,k))', (belta_k_opt(k,n)-1)/2 ]) <= (belta_k_opt(k,n)+1)/2;
end
yita_2m_opt(n) <= (H_H-H_J)^2 + norm(qr_H(:,n)-qr_J(:,n))^2 + 2*(qr_H(:,n)-qr_J(:,n))'*(qH(:,n)-qr_H(:,n)) - 2*(qr_H(:,n)-qr_J(:,n))'*(qJ(:,n)-qr_J(:,n));
for e=1:E
norm([H_J,(qJ(:,n)-we(:,e))', (tau_ke_opt(e,n)-1)/2 ]) <= (tau_ke_opt(e,n)+1)/2;
end
end
for n=1:N-1
norm(qH(:,n+1)-qH(:,n)) <= V_max*delta_t;
norm(qJ(:,n+1)-qJ(:,n)) <= V_max*delta_t;
end
log2(1+t_kr)+1/log(2)./(1+t_kr).*(tk_opt-t_kr) <= fk_opt/C_k/1e6;
Rk_loc_opt <= fTD_max/C_k/1e6;
0<=fk_opt;
fk_opt<=f_max;
0<=pj_opt;
pj_opt<=Pu_J_max;
0<=pk_opt;
pk_opt<=Pk_max;
0 <= alpha_km_opt;
alpha_km_opt <= 1;
qH(:,1)==qH(:,end);
qJ(:,1)==qJ(:,end);
cvx_end