matlab使用hampel滤波,去除异常值

 

原文链接：http://tecdat.cn/?p=7181

 

此示例显示了Hampel用于检测和删除异常值的过程的 实现。 

产生一个包含24个样本的随机信号x。 重置随机数生成器以获得可重复的结果。

rng default

lx = 24;
x = randn(1,lx);

围绕x的每个元素生成观察窗口。 在样本的任一边取k = 2个邻居。 产生的移动窗口的长度为2×2 + 1 = 5个样本。

k = 2;

iLo = (1:lx)-k;
iHi = (1:lx)+k;

截断窗口，以便函数在到达信号边缘时计算较小段的中值。

iLo(iLo<1) = 1;
iHi(iHi>lx) = lx;

记录每个周围窗口的中位数。 找到每个元素相对于窗口中位数的绝对偏差的中位数。

for j = 1:lx
    w = x(iLo(j):iHi(j));
    medj = median(w);
    mmed(j) = medj;
    mmad(j) = median(abs(w-medj));
end

缩放中位数绝对偏差

1G2erf−1(1/2)≈1.4826

 

以获得正态分布标准偏差的估计值。

sd = mmad/(erfinv(1/2)*sqrt(2));

查找与中位数相差超过nd = 2个标准偏差的样本。 将这些离群值替换为其周围窗口的中间值。 这是Hampel算法的本质。

 
yu = x;
yu(ki) = mmed(ki);

使用hampel 计算滤波后的信号并注释异常值。 覆盖在此示例中计算的过滤值。

 plot(yu,'o','HandleVisibility','off')
 

 

​

 

 如果您有任何疑问，请在下面发表评论。 

​

  

大数据部落 -中国专业的第三方数据服务提供商，提供定制化的一站式数据挖掘和统计分析咨询服务

统计分析和数据挖掘咨询服务：y0.cn/teradat（咨询服务请联系官网客服）

​QQ：3025393450

 

​QQ交流群：186388004 

【服务场景】  

科研项目; 公司项目外包;线上线下一对一培训;数据爬虫采集;学术研究;报告撰写;市场调查。

【大数据部落】提供定制化的一站式数据挖掘和统计分析咨询

欢迎选修我们的R语言数据分析挖掘必知必会课程！

 

转载于:https://www.cnblogs.com/tecdat/p/11578235.html