4.8、matlab信号生成与预处理--剔除异常值：hampel()函数

1、信号生成与预处理--剔除异常值简介

在信号生成和预处理过程中，有时候需要剔除异常值（outliers）以确保信号数据的准确性和可靠性。MATLAB提供了一些方法来识别和去除异常值，以下是一些常用的方法：

1. 箱线图检测异常值： 使用箱线图（boxplot）可以可视化数据的分布，并且明显显示出异常值。可以通过计算上下四分位数和异常值的定义来识别和标记异常值。

2. Z-Score方法： Z-Score方法是一种统计方法，可以衡量每个数据点与数据集的平均值的标准偏差之间的距离。通常情况下，如果Z-Score大于某个阈值（如3），则可以将数据点视为异常值。

3. Grubbs检验： Grubbs检验是一种统计测试，用于鉴别数据集中的异常值。它基于假设数据是正态分布的前提进行计算。

4. Hampel滤波器： Hampel滤波器是一种非参数方法，结合了中值滤波和标准差方法，可以有效地识别和剔除异常值。

一旦识别出异常值，可以采取以下方法来处理异常值：

• 替换异常值为数据集的均值、中位数或者上下限值
• 删除异常值
• 使用插值方法进行填充异常值

在MATLAB中，可以使用各种统计工具箱函数来实现异常值的识别和处理，以确保信号数据的准确性和可靠性。

2、语法

说明：对输入向量x应用Hampel滤波器来检测和去除异常值。

1）y = hampel(x)  参数：x：输入信号   y:预处理的输出信号

对于x的每个样本，函数计算由样本及其周围的六个样本组成的窗口的中位数，每边三个。它还使用中位数绝对偏差估计每个样本关于其窗口中位数的标准偏差。如果样本与中位数相差超过三个标准差，则将其替换为中位数。如果x是一个矩阵，那么该函数将x的每一列视为一个独立的通道。

2）[y,j] = hampel(x)  参数：x：输入信号  y:预处理的输出信号  j:离群指数(0/1表示，1离群)

返回一个逻辑矩阵，该矩阵在标识为离群值的所有点的位置为真

3）[y,j,xmedian,xsigma] = hampel(x) 

参数： x：输入信号  y:预处理的输出信号  j:离群指数(0/1表示，1离群)  xmedian&