How to make sure each column in a matrix sums exactly to 1?
1 view (last 30 days)
Show older comments
Hello,
I would need to create a matrix of random numbers between 1 and 0, where the diagonal of the matrix is composed of zeros only and where each column sums to 1. I approximately reached the result by the following code, but some of the columns add up to slightly less than 1 (around 0.97). Why doesn't Matlab manage to normalize all of them to 1? Is there a way to solve the problem?
Here's the code:
a = [0:0.0001:1];
b(1:10000) = (0.5/10000);
F=randsrc(100,100,[a;0.5 b]);
columnSums = sum(F);
N=100;
for n=1:N
F(n,n)=0;
end
for i = 1:numel(columnSums)
F(:,i) = F(:,i)./columnSums(i);
end
Thank you! Marco
0 Comments
Accepted Answer
Lukas Bystricky
on 7 Aug 2015
It's because you're calculating the column sums and then setting the diagonal entries to 0. If you calculate the sums after you modify the diagonal then what you wrote would work.
More Answers (1)
Torsten
on 7 Aug 2015
a = [0:0.0001:1];
b(1:10000) = (0.5/10000);
F=randsrc(100,100,[a;0.5 b]);
N=100;
for n=1:N
F(n,n)=0;
end
columnSums = sum(F);
for i = 1:numel(columnSums)
F(:,i) = F(:,i)./columnSums(i);
end
Best wishes
Torsten.
See Also
Categories
Find more on Calendar in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!