Script not running for different variable size.

I am running this code for Total_time = 1000 and it is working fine but when I increase this variable to 100000, it is not showing any results. can anyone help me for this.

Code:

close all;
clear;
clc;
% Defining variables
A = 20; % Total number of transmitters
p = 0.01;
Length_of_packet = 1000; % Maximum packet length in bits
Bitrate = 100; % bit rate for each transmitter in bits/sec
Trasmission_Time = Length_of_packet/Bitrate; 
Total_Time = 1000; % Total time for simulation in seconds
Zero_slot_overlapping = 0; % Calulation overlapping Timeslots
One_slot_overlapping = 0; % Calulation overlapping Timeslots
Two_slot_overlapping = 0; % Calulation overlapping Timeslots
Three_slot_overlapping = 0; % Calulation overlapping Timeslots
Four_slot_overlapping = 0; % Calulation overlapping Timeslots
Five_slot_overlapping = 0; % Calulation overlapping Timeslots
Six_slot_overlapping = 0; % Calulation overlapping Timeslots
Seven_slot_overlapping = 0; % Calulation overlapping Timeslots
Eight_slot_overlapping = 0; % Calulation overlapping Timeslots
Nine_slot_overlapping = 0; % Calulation overlapping Timeslots
Ten_slot_overlapping = 0; % Calulation overlapping Timeslots

User = zeros(A,Total_Time);
for k = 1:A
    for i = 1:(Total_Time-10)
        probability_of_transmission = rand(1);
        if probability_of_transmission > p
            User(k,i) = 0;
        else 
            User(k,i) = 1;
        end
    end
end 

% for k = 1:A
%     X = sum(User(k,:));
% end

for k = 1:A
    idx = find(User(k,:)==1);
    [a,b] = size(idx);
    for i = 1:b
        for j = 1:Trasmission_Time
            User(k, (idx(i)+j-1)) = 1;
            User(k, (idx(i)+j+10)) = 0; 
        end
    end
end

X = sum(User(User == 1));
Total_no_of_packets_transmitted = (round(X,-2))/10

Number_of_overlapping_packets_per_column = sum(User);

% % Create all possible row combinations
% a1 = 1:numel(User(1,:));
% a2 = 1:numel(User(:,1));
% a3 = combvec(a1,a2)';
% 
% % Assign columns of a3 to R1 and R2 to make the code readable
% R1 = a3(:,1);
% R2 = a3(:,2);
% 
% % finding the overlapps and counting them
% z = sum(and(User(R1,:),User(R2,:)),2);
% 
% % Table the data
% overlapping = table(R1,R2,z)

%Creating row pairs to compare 
row_pairs = nchoosek(1:size(User,1),2);
%disp(row_pairs);

number_of_pairs = size(row_pairs,1);
number_of_overlaps = zeros(number_of_pairs,1);

%Finding overlapping and counting them
for ii = 1:number_of_pairs
    number_of_overlaps(ii,1) = nnz(User(row_pairs(ii,1),:) & User(row_pairs(ii,2),:));     
    if number_of_overlaps(ii,1) == 0
        Zero_slot_overlapping = Zero_slot_overlapping + 1;
    elseif  number_of_overlaps(ii,1) == 1
        One_slot_overlapping = One_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 2
        Two_slot_overlapping = Two_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 3
        Three_slot_overlapping = Three_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 4
        Four_slot_overlapping = Four_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 5
        Five_slot_overlapping = Five_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 6
        Six_slot_overlapping = Six_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 7
        Seven_slot_overlapping = Seven_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 8
        Eight_slot_overlapping = Eight_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 9
        Nine_slot_overlapping = Nine_slot_overlapping + 1;
    elseif number_of_overlaps(ii,1) == 10
        Ten_slot_overlapping = Ten_slot_overlapping + 1;
        %     elseif n_overlaps(ii,1) == 11
        %          Four_slot_overlapping = Four_slot_overlapping + 2;
        %          Three_slot_overlapping = Three_slot_overlapping + 1;
    end
end

%disp([row_pairs number_of_overlaps]);
B = ([Zero_slot_overlapping; One_slot_overlapping; Two_slot_overlapping; Three_slot_overlapping; Four_slot_overlapping; Five_slot_overlapping; Six_slot_overlapping; Seven_slot_overlapping; Eight_slot_overlapping; Nine_slot_overlapping; Ten_slot_overlapping] );
T = num2cell(B)
[Zero_slot_overlapping, One_slot_overlapping, Two_slot_overlapping, Three_slot_overlapping, Four_slot_overlapping, Five_slot_overlapping, Six_slot_overlapping, Seven_slot_overlapping, Eight_slot_overlapping, Nine_slot_overlapping, Ten_slot_overlapping] = deal(T{:});

[Benjamin Thompson]

Have you tried numbers between 1000 and 100000? Is the processor just busy for a very long time or running out of resources or memory?

[Prashant Bhagat]

No, The script is running properly, only the problem is every time I run with 1000000, output is 0 0 0 0 0 0. That's the problem.

[Prashant Bhagat]

Thanks for you time.

[Jan]

You can simplify your code. Replace e.g.

User = zeros(A,Total_Time);
for k = 1:A
    for i = 1:(Total_Time-10)
        probability_of_transmission = rand(1);
        if probability_of_transmission > p
            User(k,i) = 0;
        else 
            User(k,i) = 1;
        end
    end
end

by

User = zeros(A, Total_Time);
User(:, 1:Total_Time - 10) = rand(A, TotalTime - 10) > p;

The code will be much shorter and easier, if you do not hide an index in the names of the variables like in "<Zero>_slot_overlapping" but use a vector and a real index (starting at 1 instead of 0).

The result for 1000000 might be correct: If the random signals are not overlapping.

Answer 1

[Jan]

The empty result looks correct: The larger the time, the fewer overlapping slots you find. For Total_Time=1000 there are about up to 40 overlapping slots, for Total_time=2000 about up to 8 and for 3000 less than 4. For 10'000 the code finds overlapping slots in rare cases only.

Answer 2

[Walter Roberson]

The problem is that you stop counting at 10 overlaps. As the time increases, the number of overlaps increases, until you get to the point where all of the overlaps have counts more than 10, leaving the low bins with no entries.

[Walter Roberson]

BinWidth = 10;

% Defining variables
A = 20; % Total number of transmitters
p = 0.01;
Length_of_packet = 1000; % Maximum packet length in bits
Bitrate = 100; % bit rate for each transmitter in bits/sec
Trasmission_Time = Length_of_packet/Bitrate; 
Total_Time = 100000; % Total time for simulation in seconds

slot_overlapping = 0;

User = zeros(A,Total_Time);
for k = 1:A
    for i = 1:(Total_Time-10)
        probability_of_transmission = rand(1);
        if probability_of_transmission > p
            User(k,i) = 0;
        else 
            User(k,i) = 1;
        end
    end
end 

% for k = 1:A
%     X = sum(User(k,:));
% end

for k = 1:A
    idx = find(User(k,:)==1);
    [a,b] = size(idx);
    for i = 1:b
        for j = 1:Trasmission_Time
            User(k, (idx(i)+j-1)) = 1;
            User(k, (idx(i)+j+10)) = 0; 
        end
    end
end

X = sum(User(User == 1));
Total_no_of_packets_transmitted = (round(X,-2))/10

Total_no_of_packets_transmitted = 19010

Number_of_overlapping_packets_per_column = sum(User);

% % Create all possible row combinations
% a1 = 1:numel(User(1,:));
% a2 = 1:numel(User(:,1));
% a3 = combvec(a1,a2)';
% 
% % Assign columns of a3 to R1 and R2 to make the code readable
% R1 = a3(:,1);
% R2 = a3(:,2);
% 
% % finding the overlapps and counting them
% z = sum(and(User(R1,:),User(R2,:)),2);
% 
% % Table the data
% overlapping = table(R1,R2,z)

%Creating row pairs to compare 
row_pairs = nchoosek(1:size(User,1),2);
%disp(row_pairs);

number_of_pairs = size(row_pairs,1);
number_of_overlaps = zeros(number_of_pairs,1);

%Finding overlapping and counting them
for ii = 1:number_of_pairs
    number_of_overlaps(ii,1) = nnz(User(row_pairs(ii,1),:) & User(row_pairs(ii,2),:));     
    slot = floor(number_of_overlaps(ii,1)/BinWidth) + 1;
    if slot > length(slot_overlapping)
        slot_overlapping(slot) = 1;
    else
        slot_overlapping(slot) = slot_overlapping(slot) + 1;
    end
end

x = (0:length(slot_overlapping)-1) * BinWidth;
bar(x, slot_overlapping)

[Prashant Bhagat]

Thanks you for this input.

But i want to calculate the consequtive overlapping slots from 1 to 10 so when there is 11 consequtive overlapping it will simply count it one consequitve 10 slots overlapping plus one single slot overlapping.

Can you do this?

Thanks in advance.

[Walter Roberson]

% Defining variables
A = 20; % Total number of transmitters
p = 0.01;
Length_of_packet = 1000; % Maximum packet length in bits
Bitrate = 100; % bit rate for each transmitter in bits/sec
Trasmission_Time = Length_of_packet/Bitrate; 
Total_Time = 100000; % Total time for simulation in seconds


User = zeros(A,Total_Time);
for k = 1:A
    for i = 1:(Total_Time-10)
        probability_of_transmission = rand(1);
        if probability_of_transmission > p
            User(k,i) = 0;
        else 
            User(k,i) = 1;
        end
    end
end 


for k = 1:A
    idx = find(User(k,:)==1);
    [a,b] = size(idx);
    for i = 1:b
        for j = 1:Trasmission_Time
            User(k, (idx(i)+j-1)) = 1;
            User(k, (idx(i)+j+10)) = 0; 
        end
    end
end

X = sum(User(User == 1));
Total_no_of_packets_transmitted = (round(X,-2))/10

Total_no_of_packets_transmitted = 19020

Number_of_overlapping_packets_per_column = sum(User);

% % Create all possible row combinations
% a1 = 1:numel(User(1,:));
% a2 = 1:numel(User(:,1));
% a3 = combvec(a1,a2)';
% 
% % Assign columns of a3 to R1 and R2 to make the code readable
% R1 = a3(:,1);
% R2 = a3(:,2);
% 
% % finding the overlapps and counting them
% z = sum(and(User(R1,:),User(R2,:)),2);
% 
% % Table the data
% overlapping = table(R1,R2,z)

%Creating row pairs to compare 
row_pairs = nchoosek(1:size(User,1),2);
%disp(row_pairs);

number_of_pairs = size(row_pairs,1);
number_of_overlaps = zeros(number_of_pairs,1);

%Finding overlapping and counting them
slot_overlapping = zeros(1,11);
maxslot = -inf;

for ii = 1:number_of_pairs
    number_of_overlaps(ii,1) = nnz(User(row_pairs(ii,1),:) & User(row_pairs(ii,2),:));     
    slot = number_of_overlaps(ii,1);
    maxslot = max(maxslot, slot);
    full10 = floor(slot/10);
    leftover = mod(slot,10);
    slot_overlapping(11) = slot_overlapping(11) + full10;
    if slot == 0 || leftover ~= 0
        slot_overlapping(leftover+1) = slot_overlapping(leftover+1) + 1;
    end
end

x = (0:10);
bar(x, slot_overlapping)

maxslot

maxslot = 1172

This version starts over counting when it gets to 10, so if it sees a count of (say) 137 then it would count to 10 13 times (increasing the count at the 10 slot by 13) and increase the count at the 7 slot as well.

The 10 count is huge compared to the other counts, so you can barely see them on the plot.

[Prashant Bhagat]

Thanks. I will try and let you know about it.