实验数据过滤数组

Question

我在过滤带有实验数据的数组时遇到问题 - 数据是长数组（双精度）。这个数组的索引只是时间（秒），值可以是负数或正数。问题是有时实验数据是完全错误的（（ - 这是由于 到我们设备中的错误。此错误值与其他数组值非常不同。例如 前 100 秒的值等于 +75 +（小增量值）或 - 40 -（小增量值），但是在 3 秒期间的值是 +270（这是我上面说过的错误）。如何从数组中清除这些值？平均值计数可能不是最好的方法，因为如果我从数组中获取每个值来求和，那么这个错误数据会破坏整个统计数据（由于它的“权重”很大）。我已经绘制了这个数据 - 图片很清晰但是.. 如果这个错误是由白噪声引起的，那么我可以实现例如卡尔曼滤波器，但首先没有用于计算值的函数，其次这个错误不是白噪声。

Answer 1

假设我们大多数时候都有正确的价值观。如果信号是永久的（delta 很小），我们可以使用中值过滤数据：

delta = 1e-3;
a=[75 * ones(1,20) + rand(1,20) * delta, 270 * ones(1,3) + rand(1,3) * delta];
m = median(a);
filtered_a = a(m - 2 * delta < a & a < m + 2 * delta);