computeWSC_NUAV.m

function [WSC_Val, DeltaComp] = computeWSC_NUAV(A, E, UAVs, dAB, gammaA, gammaJ, channelParam, boolFade )
    if boolFade==1
        [WSC_Val, DeltaComp] = computeWSC_NUAV_Rice(A, E, UAVs, dAB, gammaA, gammaJ, channelParam );
    else
        %   A:              Alice's Position        1x3
        %   B:              Bob's Position          1x3
        %   E:              Eve's Positions         nEx3
        %   UAV1:           Position of UAVs        nUAVx3
        %   Rj:             Orbit radius of UAVs    1x1
        %   hj:             UAVs height             1x1
        %   dAB:            Distance from A to B	1x3     Can be actual or estimated
        %   gammaA:         Tx SNR at A             1x1
        %   gammaJ:         Tx SNR at UAVs          1x1
        %   channelParam:   Channel and PL parameters

        nUAV = size(UAVs,1);
        B = A;
        B(1) = B(1) + dAB;

        % Channel parameters
        phi         = channelParam (1);
        omega       = channelParam (2);
        alpha       = channelParam (3);
        alpha_AG    = channelParam (4);
        ne_LOS      = channelParam (5);
        ne_NLOS     = channelParam (6);
        Rs          = channelParam (7);
        KbVal       = channelParam (8);
        KeVal       = channelParam (9);


        % Parameters regarding Eve
        dAE = transpose(sqrt( ( A(:,1) - E(:,1) ).^2 + ( A(:,2) - E(:,2) ).^2 ));
        OmegaAE = dAE.^alpha;

            % UAVs
        dJE = sqrt( ( UAVs(:,1) - E(:,1)' ).^2 + ( UAVs(:,2) - E(:,2)' ).^2  + ( UAVs(:,3) - E(:,3)' ).^2);
        Theta_JE = (180/pi) * asin(UAVs(:,3)./dJE);
        PLOS_JE = 1./(1 + phi * exp( -omega*( Theta_JE - phi ) ) );
        LJE = PLOS_JE.*(abs(dJE).^alpha_AG)*ne_LOS + (1-PLOS_JE).*(abs(dJE).^alpha_AG)*ne_NLOS;
        gJE = 1./LJE;
        intJE = gammaJ*sum(gJE,1);

        aNJ = gammaA;
        bNJ = gammaA;

        % Parameters regarding Bob
        OmegaAB = dAB.^alpha;

            % UAVs
        dJB = sqrt( ( UAVs(:,1) - B(1) ).^2 + ( UAVs(:,2) - B(2) ).^2  + ( UAVs(:,3) - B(3) ).^2);
        Theta_JB = (180/pi) * asin(UAVs(:,3)./dJB);
        PLOS_JB = 1./(1 + phi * exp( -omega*( Theta_JB - phi ) ) );
        LJB = PLOS_JB.*(abs(dJB).^alpha_AG)*ne_LOS + (1-PLOS_JB).*(abs(dJB).^alpha_AG)*ne_NLOS;
        gJB = 1./LJB;
        intJB = gammaJ*sum(gJB,1);


            % Secrecy metrics
        aJ = gammaA ./ ( 1 + intJB );
        bJ = gammaA ./ ( 1 + intJE );

        SOP_J   = 1 - (exp( -(OmegaAB./aJ ).*(2^Rs - 1) ))*( 1 ./ ( (2^Rs)*(OmegaAB./OmegaAE).*( bJ./aJ ) + 1 ) );
        SOP_NJ  = 1 - (exp( -(OmegaAB./aNJ).*(2^Rs - 1) ))*( 1 ./ ( (2^Rs)*(OmegaAB./OmegaAE).*(bNJ./aNJ) + 1 ) );

        DeltaComp = (1-SOP_J)./(1-SOP_NJ);

        coverage    =  sum( DeltaComp(:)>=1);
        efficiency  = mean( DeltaComp(:) );
        WSC_Val     = coverage*efficiency;
    end
end