Sir,sorry, I didn’t make it clear.I’m very sorry to take up your time.

I run the diagnostic program you mentioned the other day, as follows.

For diagnostic purposes, replace

```
variable p(6);
variable m(6,6)
...
minimize(obj)
```

with

```
p = ones(6,1);
m = ones(6,6);
...
my_objective = obj
whos my_objective
```

I can run it and get an image with no data but coordinates.

警告: Empty cvx model; no action taken.

In cvx_end (line 22)

In ceshi2 (line 74)

However, running the original program later will still prompt me that there is no scalar. What could be the reason for this?

错误使用 cvx_end (line 267)

Your objective function is not a scalar.

出错 ceshi2 (line 74)

cvx_end

For the convenience of explanation, the original code is attached：

syms V;%非负控制参数

syms t;%时隙

n=0.9;%基站间传输能量的损耗率

G= 18;%

r=1000000;%用户设备传输速率,bit,1Mb/s

c=0.2;%冷却时间的约束，最多占时隙的百分比

d=200;%m,

q= 100;%噪声功率电平10DB的平方

for V=0:20:100

for t=0:1:5000

L=normrnd(0,5.8);%生成一个正态分布随机数

L=(6.14+20*log10(200)+L)**18;*

if t==0 %冷却时间和能量队列的初始值

F=[0,0,0,0,0,0];

E=[0,0,0,0,0,0];

B=[0,0,0,0,0,0];

end

a= unifrnd(0,20,1,1);%dBm

b=unifrnd(20,40,1,1);%dBm

e =unifrnd (a,b,1,6);%dBm，

cvx_clear

cvx_begin

variable p(6) nonnegative;%非负

variable m(6,6) nonnegative;%非负

Obj = 0;

expression m1(6)

expression m2(6)

for j=1:6

for x=1:6

if x~=j

m1(j)=m1(j)+m(j,x);

end

end

for x=1:6

if x~=j

m2(j)=m2(j)+m(x,j);

end

end

F0=F(j);

E0=E(j);

B0=B(j);

e0=e(j);

p0=p(j);

Obj = Obj + (F0+V)(c/r)(3*r-(log(1+(p0*L)/300)/log(2))+3*B0)+(E0+V)**p0+E0*e0+(V(1-n)-E0)*m1+n*E0*m2;*

end

minimize(Obj)

subject to

for j=1:1:6

for x=1:6

if x~=j

m1(j)=m1(j)+m(j,x);

end

end

for x=1:6

if x~=j

m2(j)=m2(j)+m(x,j);

end

end

p(j) <=E(j)+nm2(j)-m1(j);

end

cvx_end

E1=0;

F1=0;

m3=[0,0,0,0,0,0];

m4=[0,0,0,0,0,0];

for j=1:6

for x=1:6

if x~=j

m3(j)=m3(j)+m(j,x);

end

end

for x=1:6

if x~=j

m4(j)=m4(j)+m(x,j);

end

end

F(j)=F(j)+1-log2(3*q+p(j)**L)/(3*r)-B(j)/r+log2(3q)/(3*r)-c;*

E(j)=E(j)-p(j)+e(j)-m3(j)+nm4(j) ;

B(j)=B(j)-r*(log2(3*q+p(j)**L)/(3*r)+B(j)/r-log2(3q)/(3*r))+log2(1+p(j)*L/q)/3;*

E1=E1+E(j);

F1=F1+3F(j);

end

end

figure(1);

plot(V,E1);

figure(2);

plot(V,F1);