For loop with decimal increment not seen as integer or logical value
9 次查看(过去 30 天)
显示 更早的评论
Goal:
- Have for loop repeatedly create matrices of random integers
- Create new matrix where columns = variance for each row of random integer matrices
- Compare variance to increment ranging from 0:0.01:4 (0-4 in increments of 0.01) and retain each matrix of 0s and 1s generated this way
The following code produces the desired results when the increment = 1 (not 0.01)
T = 10;
W = zeros(3,10);
for trial = 1:T
mtest = randi([1,5],3,4)
W(:,trial) = var(mtest,0,2);
end
minc = zeros(3,10,5);
for mindex = 0:1:4
minc(:,:,mindex+1) = W >= mindex
end
Problem 1: The same code does not work when the increment is set to 0.01:
linc = zeros(3,10,401);
for lindex = 0:0.01:4
linc(:,:,lindex+1) = W >= lindex
end
Error = Index in position 3 is invalid. Array indices must be positive integers or logical values.
Problem 2: After much reading, I think the problem is that "lindex+1" produces indices of 1.01, 1.02...etc. Yet both of the below solutions produce the same error.
Solution 1: Add 1 to remove 0 from index; multiply by 100 such that index should proceed as 100, 101, 102...etc.
linc = zeros(3,10,401);
for lindex = 0:0.01:4
linc(:,:,(lindex+1)*100) = W >= lindex
end
Solution 2: Start index at 0.01 instead of 0 and multiply by 100 such that index should proceed as 100, 101, 102...etc.
linc = zeros(3,10,401);
for lindex = 0.01:0.01:4
linc(:,:,lindex*100) = W >= lindex
end
Why does the code work if increment = 1 but not 0.01? How can I fix this? Thank you in advance, this is a wonderful community.
3 个评论
采纳的回答
KSSV
2020-1-23
编辑:KSSV
2020-1-23
lindex = 0:0.01:4 ;
linc = zeros(3,10,401);
for i = 1:length(lindex)
linc(:,:,i) = W >= lindex(i) ;
end
Please note that, your indices should not be zero, or fractions or negative. In your original code in problem 1, you have used fractions, zero as indices, therefore you gor errors.
3 个评论
KSSV
2020-1-23
When the step size was 1, you had all integers as indices and there was no error. When the step size was 0.001, you have fractions and indices cannot be fractions.
Walter Roberson
2020-1-23
If you take values constructed with the : operator with increment 0.01 and you multiply the values by 100, then roughly 1/10th of the results will not be integers.
MATLAB does not represent numbers in decimal: it represents them in binary. In binary, 1/100 requires an infinite repeating number to represent exactly, for precisely the same reason that in decimal, 1/7 needs an infinite repeating number to exactly represent. Therefore any finite binary truncation of a multiple of 1/100 will not exactly represent the fraction, and when you multiply the inexact representation by 100 you will not always get out an integer.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!