find is not working in a loop, it skips numbers

조회 수: 7 (최근 30일)
Katerina F
Katerina F 2025년 2월 24일
댓글: Katerina F 2025년 2월 24일
Hello,
I am using a loop with a "find" to indentify the rows in the third column that have the same number (e.g. 0.75 etc), then I select which to chose based on a criterio. In some cases for example 0.82, 0.93 the function find returns a zero array instead of an array with ones at the locations of the number. Depending on the limits of the initial range ("sisto" in the code below) it skips different number. I tried to correct using the mask instead of the find but I am still having the same results. Could you please help me with this? (I attach the text file with the input matrix) Thanks.
load('HP4140AT5OR40DOS.txt');
LINE=HP4140AT5OR40DOS;
x=LINE(:,3);
sisto = 0.60:0.01:1;
s=size(sisto,2)
for i=1:s
sisto(i)
find = (x==sisto(i));
MA=LINE(find,:)
dast=MA(:,5); %the storage days
ouf=abs(40-dast)% how close they are to 30
[minValue, minIndex] = min(ouf)
point=MA(minIndex,:)%the point I select for this wind contribution
end
  댓글 수: 2
Stephen23
Stephen23 2025년 2월 24일
"find is not working in a loop, it skips numbers"
FIND is working. Your numbers are different.
"Could you please help me with this?"
Do not use EQ to compare exact equality of binary floating point numbers. Always compare using a tolerance, e.g.:
tol = 1e-5; % select this value depending on the data
abs(A-B)<tol
Steven Lord
Steven Lord 2025년 2월 24일
If you're using a release that has it, use the isapprox function.

댓글을 달려면 로그인하십시오.

채택된 답변

Torsten
Torsten 2025년 2월 24일
편집: Torsten 2025년 2월 24일
Never test two double precision numbers for equality. Because of floating point arithmetic, they will almost always be classified as different.
Replace
find = (x==sisto(i));
by
find = abs(x-sisto(i))<1e-6;
And better don't use the name "find" for the variable. You overwrite the internal MATLAB function "find":

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Loops and Conditional Statements에 대해 자세히 알아보기

태그

Community Treasure Hunt

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

Start Hunting!

Translated by