my code:

cvx_solver sedumi

cvx_begin

variable Bc(G,1)

minimize (a*b*sum(rel_entr(lamda_g+n*balloct*Bc+guiyihuaxishu_col.*(P-sum(Po)-Pc)./(Ug-1),fenzi_col. Pc+lamda_g+nballoct*Bc+guiyihuaxishu_col.

*(P-sum(Po)-Pc)./(Ug-1))/n…*

-(lamda_g/n+guiyihuaxishu_col.(P-sum(Po)-Pc)./((Ug-1)

-(lamda_g/n+guiyihuaxishu_col.

*n)).*(sum(Bc)-B))

*log(1-(fenzi_col.*

(c/2).(sum(Bc + Bo)-B )^2 + lamda2*Pc).*(P-sum(Po)-Pc)./(Ug-1))))/log(2)+ …*inv_pos(fenzi_col.*Bc+guiyihuaxishu_col.*Pc+lamda_g+n*balloct(c/2).

subject to

Bc>=B/G/9;

sum(Bc + Bo)-B <=0;

cvx_end

my optimation question: max the data rate