用matlab实现一个简单的离群点挖掘(与时序有关)

      在视频分析中，需要用到离群点挖掘，下面给出一个简单的离群点挖掘, matlab代码如下：

 

clearclcorgX = [1 1 0 3 4 5 2 70 100  100 100 6 60 3 2 2 1 1 0]k = 3;th = 1;outlier = [];x = newArray(orgX, k)len = length(x);average = mean(x);standard = std(x);for i = 1 : len    if x(i) > 0 && x(i) - average > th * standard     % 此处没有用abs函数，所以，有些负数可能没有离群        outlier = [outlier orgX(i)];    endendoutlier

function y = newArray(x, k);n = length(x);for i = 1 : k    y(i) = 0;    y(n - i + 1) = 0;endfor i = k + 1 : n - k    y(i) = x(i) - ( min(x(i - k : i - 1)) + min(x(i + 1 : i + k)))/2;  end

        结果为：

orgX =

  Columns 1 through 18

     1     1     0     3     4     5     2    70   100   100   100     6    60     3     2     2     1     1

  Column 19

     0

x =

  Columns 1 through 10

         0         0         0    2.0000    3.0000    4.0000  -34.5000   19.0000   96.0000   96.0000

  Columns 11 through 19

   63.5000  -45.0000   56.0000   -0.5000         0    1.0000         0         0         0

outlier =

   100   100   100    60