I want to return the value of x , description below

%%%%%% I want to return the x value from File 2 to File 1.

%%%% File 1
% Total Capacity After Water_filling and Genetic Algorithm
clc;
clear;

M = 60; %Number of Antennas BS
N = (6:10).^2; %Number of Elements IRS
K = 1:10; %Number of Single Antennas UE
L = 3; %Number of paths
fc = 28e9; %Carrier freq
SNR = -10; %Signal to Noise Ratio
data.Pt = 1; %Total Power
data.Pn = data.Pt/(10^(SNR/10)) ; %Noise Power

C_A = zeros(numel(N),numel(K));

Max_iter = 1;

for iter = 1:Max_iter
    for n = N
        
        %%% Millimeter wave Channel
        G = M_Wave_channel_Matrix(fc, M, n, L, 'BS', 'IRS'); %NxM
        data.G = G;
        
        F = zeros(1,n);
        H = zeros(1,M);
        
        for k = K
            
            %%% channel IRS-UE
            fk = M_Wave_channel_Matrix(fc, n, 1, L, 'IRS', 'UE'); %1xN
            F(K==k,:) = fk; %KxN
            data.F = F;
            
            %%% channel BS-UE
            hk = M_Wave_channel_Matrix(fc, M, 1, L, 'BS', 'UE'); %1xM
            H(K==k,:) = hk; %KxN
            data.H = H;
            
            
            %%% Genetic Algorithm
            FT = @(x)-Genetic_A_Set(x,data);
            A = []; % No other constraints
            b = [];
            Aeq = [];
            beq = [];
            lb = zeros(1,n);
            ub = ones(1,n);
            opts = optimoptions(@ga, ...
                'PopulationSize', 150, ...
                'MaxGenerations', 200, ...
                'EliteCount', 10, ...
                'FunctionTolerance', 1e-2);
            [x,fval] = ga(FT,n,A,b,Aeq,beq,lb,ub,[],opts);
            C_ga = -fval;
            C_A(N==n,K==k) = C_A(N==n,K==k) + 1/Max_iter*(C_ga); %Algorithm Capacity
            
        end
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% File 2
% Function of Genetic Algorithm
function C_ga = Genetic_A_Set(x,data)

theta = (2*pi)*(x); %Reflecting Angle IRS
Phi = diag(exp(1i*theta)); %Phase Shift Matrix for IRS. Amplitude = 1

HT = data.H+data.F*Phi*data.G; %Channel Matrix
do_waterfill = 1; % Flag

%% Rep_Water_filling
while do_waterfill
    
    W = pinv(HT); %Precoding Matrix
    W_bar = W./sqrt(sum(abs(W.^2),1));
    D_Matrix = HT*W_bar; %Diagonal Matrix
    D_Square = (abs(diag(D_Matrix))').^2; % Channel gains
    P = waterfill(data.Pt,data.Pn./D_Square);
    if any(P==0)
        HT(P==0,:) = []; %Delete User Row
    else
        do_waterfill = 0;
    end               
end

R = P.*D_Square/data.Pn;
C_ga = sum(log2(1+R)); %Capacity

end

[Soniya Jain]

At this line, where you are calling the function of file 2,

%%% Genetic Algorithm
FT = @(x)-Genetic_A_Set(x,data);

The function is returning the value of C_ga, so do you want the value of C_ga or x or both to be returned by the function?

Also, in the function, value of x is nowhere modifying, so why do you want the the function to return the value of x?

[Matthew Worker]

Actully I made a mistake I want to return the size of R, because I want to compare it with the C_ga.

This function (FT = @(x)-Genetic_A_Set(x,data)) I'm using it for Genetic Algorithm, so I was thinking maybe it is better to create another function to return the sizes of R.

[Soniya Jain]

% Function of Genetic Algorithm
function C_ga = Genetic_A_Set(x,data)

So, at this line, replace C_ga with the variable name suppose size_R which captures the size of R in the code, like

size_R = size(R);

Note: size(R) will return 1X2 array, refer https://www.mathworks.com/help/matlab/ref/size.html for more information of size() function.

[Matthew Worker]

I changed the size(R) to length (R) so I'll return only scalar value, I think even the length of R is ok for me to compare it with C_ga, but the promblem is by replacing function C_ga = Genetic_A_Set(x,data) with

function length_R = Genetic_A_Set(x,data) is not a good solution for me because I want to return C_ga also so I can use it with the genetic algorithm, so I want to return the length of R using another function because I only need it for comparison at the end of the code.

Genetic Algorthim function will return 50 values of C_ga and I expect to get 50 values of R.

[Walter Roberson]

What is the difference between this and https://www.mathworks.com/matlabcentral/answers/867095-the-loops-in-the-first-picture-are-required-the-description-is-below?s_tid=srchtitle ? I don't like trying to solve the same question in two different discussions.

[Matthew Worker]

ok, It is completely different.

here I want to return the length of R (Length_R) from File 2, but the problem is if you read the decision below that I need to return also C_ga which I will put in FT and use for the Genetic Algorithm in file one, but Length_R I need to return it to file 1 to compare it with C_A.

[Image Analyst]

Original question

I want to return the value of x , description below

dg = sum(log2(1+R)); %Capacity

end

ok, It is completely different.

here I want to return the length of R (Length_R) from File 2, but the problem is if you read the decision below that I need to return also C_ga which I will put in FT and use for the Genetic Algorithm in file one, but Length_R I need to return it to file 1 to compare it with C_A.

[Rena Berman]

(Answers Dev) Restored edit

Answer 1

[Soniya Jain]

% Function of Genetic Algorithm
function C_ga = Genetic_A_Set(x,data)

So, instead of writing another function to find length of R, you can return the length of R in this function only,

% Function of Genetic Algorithm
function [C_ga, length_R] = Genetic_A_Set(x,data)

[Matthew Worker]

After using .. function [C_ga, length_R] = Genetic_A_Set(x,data)

The (FT) will have both values ​​and I need only the C_ga value in FT because I will use the FT in Genetic Algorithm function as following :

1- FT = @(x)-Genetic_A_Set(x,data);

2- [x,fval] = ga(FT,n,A,b,Aeq,beq,lb,ub,[],opts);

3- C_ga = -fval;

Can you tell me how exactly I can call the function twice in file 1 .. the first one to return C_ga in FT and the second one to return the length_R.

Thank you

[Soniya Jain]

Instead of calling the function, this way

FT = @(x)-Genetic_A_Set(x,data);

you can first store the return values in some variable in file1 and can use those variables further in the code, like this

[C_ga,lengthR] = Genetic_A_Set(x,data);
FT = @(x) - C_ga;

[Matthew Worker]

So I can use it like this right to return the value of lengthR

R_L =@(x) lengthR

[Soniya Jain]

Yes, i guess.

[Matthew Worker]

[Soniya Jain]

Hey, can you also share the source code?

[Matthew Worker]

I can send it to you by email.

[Soniya Jain]

If possible, send it here only.

[Matthew Worker]

I think I should return R_length by using another function. right