Add the second matrix to the code, concatenate it with ‘H’ along dimension 3, then use a for loop to address and plot them.
H1=[0 1 1 0 1 0 1 1 1 0;
0 1 1 1 0 1 0 0 1 1;
1 1 1 0 0 1 1 1 0 0;
1 0 0 1 1 0 0 1 1 1;
1 0 0 1 1 1 1 0 0 1];
H2 = randi([0 1], size(H1)); % Second Matrix
Hm = cat(3,H1,H2); % Concatenate To 3D Array
hold on
for k = 1:size(Hm,3)
H = Hm(:,:,k); % Set 'H' To Correct Matrix
dvi = sum(H,1);
dci = sum(H,2)';
[adv,bdv] = hist(dvi,unique(dvi));
[adc,bdc] = hist(dci,unique(dci));
dv = sum((adv./sum(adv)).*bdv);
dc = sum((adc./sum(adc)).*bdc);
set(gca,'XTick',[]);
set(gca,'YTick',[]);
adv(adv == 0) = NaN;
stem(bdv,adv,'linewidth',2);
% hold on;
adc(adc == 0) = NaN;
stem (bdc,adc, 'linewidth',2);
xlim([min(min(bdv),min(bdc))-1 max(max(bdv),max(bdc))+1]);
end
hold off
xlabel('Degree of nodes');
ylabel('Number of nodes');
grid on;
axis on;
legend('dv','dc','location','NorthEast');
This works with the matrices, however it does not appear to produce the result of the original code with the random matrix I created here to test the code. Since I am not certain what that is, I defer to you to solve that. It might be necessary to use 'DisplayName' to create the legend correctly.
.
