Octave 5.2.0两向量匹配无法找到多个元素的问题求助

我需要用for循环比较两个不同维度的向量：从CSV文件读取的energy_a，以及手动定义的xminpeaks、xmaxpeaks。要求使用Octave 5.2.0版本，目标是生成长度为28的indexmin和indexmax向量，但当前代码仅能生成长度为20的向量，部分值（例如2.8050）始终无法匹配到。

考虑到浮点数精度问题，我尝试将所有数值乘以10000转为整数后再匹配，但修正变量名笔误后仍失败，生成的向量长度为21，依旧找不到28050（对应2.8050）这类值，还出现匹配错误值的情况。

现有代码 #

第一段代码（直接浮点数比较） #

% Load data file
fid = fopen("Fluo_spectre_aragonite.csv",'r');
indata = textscan (fid , '%f %f %f %f %f %f %f %f %f %f %f %f %f',
'Delimiter', ';',
'HeaderLines' ,1);
fclose (fid);

energy_a = indata{1, 1};

%define the vectors
xminpeaks=[1.4400, 1.8600, 2.8050, 3.1050, 3.5100, 3.9000, 4.4400, 4.6650, 5.1150, 5.6850, 6.0300, 6.7350, 7.1700, 7.6500, 8.2050, 9.1800, 9.9150, 12.5400, 13.7850, 15.4800, 16.2900, 17.0550, 17.7150, 18.9600, 21.3450, 22.4950, 24.1650, 25.2600];

xmaxpeaks=[1.5750, 2.0400, 3.1050, 3.4050, 3.9000, 4.1100, 4.6650, 5.1150, 5.6850, 6.0300, 6.7350, 7.1700, 7.6500, 8.2050, 9.1800, 9.9150, 10.6800, 12.8400, 14.4600, 16.2900, 17.0550, 17.7150, 18.9600, 21.3450, 22.4950, 24.1650, 25.2600, 25.7400];

indexmin=[];
indexmax=[];

%finding the values
for i = 1:length(xminpeaks)
  index=1;
  for j = 1:length(energy_a) 
    if energy_a(j)== xminpeaks(i)
      indexmin = [indexmin index];
    endif
    if energy_a(j) == xmaxpeaks(i)
      indexmax = [indexmax index];
    endif 
    index+=1;
  endfor
endfor

第二段代码（转整数比较，存在笔误） #

xminpeaks=[1.4400, 1.8600, 2.8050, 3.1050, 3.5100, 3.9000, 4.4400, 4.6650, 5.1150, 5.6850, 6.0300, 6.7350, 7.1700, 7.6500, 8.2050, 9.1800, 9.9150, 12.5400, 13.7850, 15.4800, 16.2900, 17.0550, 17.7150, 18.9600, 21.3450, 22.4950, 24.1650, 25.2600];

xmaxpeaks=[1.5750, 2.0400, 3.1050, 3.4050, 3.9000, 4.1100, 4.6650, 5.1150, 5.6850, 6.0300, 6.7350, 7.1700, 7.6500, 8.2050, 9.1800, 9.9150, 10.6800, 12.8400, 14.4600, 16.2900, 17.0550, 17.7150, 18.9600, 21.3450, 22.4950, 24.1650, 25.2600, 25.7400];

xminpeaks10000=[xminpeaks*10000];
xmaxpeaks10000=[xmaxpeaks*10000];
energy_a10000=[energy_a*10000];

indexmin=[];
indexmax=[];

%on cherche les index de chaque intervalles de valeurs englobant les pics
for i = 1:length(xminpeaks)
  index=1;
  for j = 1:length(energy_a)
    % 注：原代码存在笔误，xminpeaks1000应为xminpeaks10000，xmaxpeaks1000应为xmaxpeaks10000
    if energy_a10000(j)== xminpeaks1000(i)
      indexmin = [indexmin index];
    endif
    if energy_a10000(j) == xmaxpeaks1000(i)
      indexmax = [indexmax index];
    endif 
    index+=1;
  endfor
endfor

问题根源与解决方法 #

1. 浮点数精度本质问题 #

部分十进制小数（如2.8050）无法被二进制浮点数精确表示，即使乘以10000后，实际存储值可能是28049.999999999996而非理想的28050，导致直接相等判断==失败。

2. 改用近似匹配（核心解决方法） #

不要用==直接比较，改用误差容忍范围内的近似匹配，设置一个极小的阈值（如1e-6），判断两个值的绝对差是否小于该阈值。同时移除多余的index变量，直接使用循环变量j作为索引：

% 替换原循环内的判断逻辑
for i = 1:length(xminpeaks)
  for j = 1:length(energy_a) 
    if abs(energy_a(j) - xminpeaks(i)) < 1e-6
      indexmin = [indexmin j];
    endif
    if abs(energy_a(j) - xmaxpeaks(i)) < 1e-6
      indexmax = [indexmax j];
    endif 
  endfor
endfor

3. 验证数值实际存储值 #

可以打印energy_a中疑似匹配值的精确存储，确认差异：

% 查找energy_a中接近2.8050的数值
target = 2.8050;
candidates = energy_a(abs(energy_a - target) < 1e-3);
disp('接近2.8050的数值：');
disp(candidates);
disp('精确浮点表示：');
disp(sprintf('%.15f', candidates));

4. 向量化优化（可选） #

Octave支持向量操作，用ismembertol函数（Octave 4.4及以上支持）可以直接实现近似匹配，替代嵌套循环，效率更高：

% 生成长度为28的indexmin
[~, indexmin] = ismembertol(xminpeaks, energy_a, 1e-6);
% 生成长度为28的indexmax
[~, indexmax] = ismembertol(xmaxpeaks, energy_a, 1e-6);

内容的提问来源于stack exchange，提问作者Brclsbtn