Script not running for different variable size.

1 次查看（过去 30 天）
Prashant Bhagat
Prashant Bhagat 2022-8-15
评论： Prashant Bhagat 2022-8-16
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{:});
  4 个评论
显示 2更早的评论
Jan
Jan 2022-8-16
编辑：Jan 2022-8-16
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.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Jan
Jan 2022-8-16
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.

更多回答（1 个）

Walter Roberson
Walter Roberson 2022-8-16
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.
  4 个评论
显示 2更早的评论
Walter Roberson
Walter Roberson 2022-8-16
Ran in:
% 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.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Chebyshev 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by