hello i wrote my code.

clc

clearvars

close all

%% Initial values

t_k = 2; % number of transmitter antenna

r_k = 2; % number of receiver antenna

N = 10^4; % size of data input

N_Realization = 300; % number of channel realization in SISO channel

data_in = randi([0 1], t_k, N); % random data input

SINR_value_db = 6;

Harvestedenergy_value_dbm = -20;

H_MIMO = zeros(t_k,r_k); % defineing channel matrix

%% channel SISO

p_k = rand (1,2); % power of transmitter

p=sqrt (p_k);

k = 0; % SNR loop counter

noise_var=-90

noise_circuit=-60

for SNR_range_variable = SINR_value_db % SNR loop

for w = 1:N_Realization

k = k + 1; % adding in SNR loop counter

% for w = 1:N_Realization % Realization channel loop

H_MIMO = randn (t_k,r_k); % random SISO channel

data_out_ch_SISO =p * H_MIMO * data_in; % pass from SISO channel *

% data_out_ch_AWGN = H_AWGN_KK * data_in; % pass from AWGN channel *

% % reciver

data_out_ch_nois_MIMO(:,:,k) = awgn (data_out_ch_SISO, SNR_range_variable, noise_var); % adding noise MIMO

%data_out_ch_nois_AWGN = awgn (data_out_ch_AWGN, SNR_range_variable, -90); % adding noise AWGN

end

end

%% SINR & EH

for k=1:2;

TETA = randi([0 1], 1 ,2)

mi = randi([0 1], 1 ,2)

norm_h_mimo=norm( [H_MIMO(k,k)],1)

y= p_k * norm([H_MIMO],1)^2 ;

SINR(k) =(TETA(1,k) * p_k(1,k) * norm_h_mimo^2 ) / TETA(1,k) * sum(y,2) + TETA(1,k) * (noise_var ^ 2) + (noise_circuit ^2) ;

POWER_ratio_ID=TETA(1,k)

POWER_ratio_EH=1-TETA(1,k)

EH(k) = mi(1,k) * POWER_ratio_EH * ( sum(y,2) + (noise_var ^ 2));

end

%% optimization

```
in_snr=6;
```

e=6;

cvx_begin

%%%%for k=1:2;

variable p_k(k)

variable TETA(k)

% dual variable SNR

minimize(sum(p_k))

subject to

SINR >= SINR_value_db

EH >= -20

% minimize teta(k)

%%%% for j=1:k+1;

% k=k+1

%%% h(j,k)=inv_pos (1/sqrt(2))*(randn+i+randn)

% if j==k

% h(j,k)=h(j,k+1);

% else go to 25

% a=

%%%%(j,k)=0

% a = mtimes( j, k, ‘rdivide’ );

%%%%inv_pos -((1/in_snr)+1)* p_k(1,k) * norm_h_mimo^2+ y + noise_var +(-60/TETA(1,k))<=0;

%%%% (6/(1-teta(k)))-(a-90)<=0

%%%%1e-6 <= theta <= 1-1e-6

%%%%p >= 1e-6

%%%%end

%%%%end

% end

cvx_end

i have this result

are they true?

why i have ‘Unbounded’ status or -inf optval?