Octave 如何有效地比较矩阵元素

如何解决Octave 如何有效地比较矩阵元素

我有一个大矩阵，我想比较相邻元素以验证值是否接近相等。

例如，在这个序列 1006,1004,999,1000,1003,6,1005,1003 ..... 值 6 并不“接近”1003 或 1005。 我想要一种有效的方法来进行比较。

这是查找范围之外的任何内容的缓慢代码。在我的旧电脑上需要 190 秒。

谢谢。

big = 1e5;
tic;
a = 0;
x = rand(100,big);
for ii = 1:100
  for jj = 1:big-1;
    y = x(ii,jj) / x(ii,jj+1);
    if (or(y < 0.999,y > 1.001))   a++;
    endif  
  endfor;
endfor;
toc

解决方法

使用矢量化代替循环：

y = x(:,1:end-1) ./ x(:,2:end);
a = nnz(y < 0.999 | y > 1.001);