Numerical precision when using colon operator

2 次查看（过去 30 天）
Xingwang Yong
Xingwang Yong 2022-4-2
评论： Stephen23 2022-4-2
I am doing simulation, which basically is iteration over time. Given the delta time dt and total time inter_pulse_space*num_pulse, the full time frame time1 can be easily calculated. I want to observe signal at time2 only. Surprisingly, time2 is not a subset of time1. Is this due to numerical precision?
clear
dt = 0.1/1e3; % 0.1 millisecond
num_pulse = 100;
inter_pulse_space = 40*dt;
time1 = dt:dt:inter_pulse_space*num_pulse;
time2 = inter_pulse_space/2 + (0:num_pulse-1)*inter_pulse_space;
ismember(time2,time1)
format long
time1(180)
time2(5)

请先登录，再回答此问题。

采纳的回答

KSSV
KSSV 2022-4-2
You cannot use ismember for floating point numbers. If you want to have values at time2, you can use interp1. Read about interp1.
If you want to comapre floating point numbers refer:
https://www.geeksforgeeks.org/problem-in-comparing-floating-point-numbers-and-how-to-compare-them-correctly/
Also read about ismembertol.
  1 个评论
Stephen23
Stephen23 2022-4-2
You certainly can use ISMEMBER for floating point numbers.
If you understand the finite precision of binary floating point numbers and expect/require exact matches only, then there is absolutely nothing stopping you from using ISMEMBER.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 NaNs 的更多信息

标签

产品


版本

R2020a

Community Treasure Hunt

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

Start Hunting!

Translated by