Найти заявление в MATLAB и несколько условий

Question

Следующий код встраивается в функции

rsp = find(response_times >= current1 & response_times < current2 & response_times ~= current2); 

Здесь я ищу для индексов ответов, которые происходят между current1 и current2, где current1 и current2 времена, такие как 16.22 и 16.32 соответственно, а время отклика может быть равно current1, но не current2.

По большей части это работает, как и предполагалось, однако, каждый раз он вытягивает индекс значения, равного current2.

Кто-нибудь знает, почему это может быть так, или как я могу улучшить эту строку кода, чтобы исправить ее.

Вот пример массива этот код работает на:

response_times = [ 8.73000000000000 
    11.4300000000000 13.4800000000000  
    14.7900000000000 16.3200000000000  
    18.0400000000000 20.3800000000000  
    20.9900000000000 21.3400000000000  
    24.2800000000000 24.6800000000000 ]; 

Answer 1

Это на самом деле не (точно) равна current2. Обычно вам не следует сравнивать числа с плавающей запятой для равенства. Читайте this article для получения дополнительной информации, но суть проблемы в том, что большинство значений не могут быть представлены точно с представлением с плавающей запятой (т. Е. IEEE 754), поэтому нецелесообразно проверять на равенство. Неравенства в порядке. На самом деле есть аккуратный mini-site, посвященный основному объяснению проблемы.

Для переполнения стека matlab версия объяснения см. Это Q & A, озаглавленный "Why is 24.0000 not equal to 24.0000 in MATLAB?". Это довольно интересно читать!

Чтобы убедиться, что ваше неравенство действительно работает, попросите его вычислить abs(response_times-current2), и вы должны обнаружить, что в этих случаях значение не равно нулю, а нечто вроде 1.35e-15. Если вы хотите отклонить эти «слишком близкие» значения, включите тест, например (current2-response_times)>tol, и установите tol на что-то достаточно большое, чтобы отклонить эти очки.