如何在數字數組在MATLAB找出突破點

Question

下午好，大家好，我有這樣的新的問題，我希望你能再次幫助我：

我有一個載體，它可以在下一環節發現：

https://drive.google.com/file/d/0B4WGV21GqSL5Y09GU240N3F1YkU/edit?usp=sharing

矢量繪圖看起來像這樣：

正如你所看到的，圖中有一些部分數據的行爲幾乎是線性的。這就是我說的：

我需要的是找到基於在數據某些部分的線性的突破點。而且你可能會問自己，當部分數據不是線性的時候會發生什麼，那麼算法就不會採用這個部分。

我希望你能幫助我，謝謝。

Answer 1

你想要做的就是分段線性時間序列分割。

有很多方法來解決這個問題，它們的複雜性和準確性不同。

這是最簡單的一種，稱爲滑動窗口分割：

function [breaks vals] = segment(data, max_error) 
    breaks = []; 
    vals = []; 
    left = 1; 
    for right = 2:length(data) 
     err = linear_regresion(data(left:right)); 
     if max(abs(err)) > max_error 
      breaks(end+1) = right-1; 
      vals(end+1) = data(right-1); 
      left = right; 
     end 
    end 
end 

function err = linear_regresion(data) 
    n = length(data); 
    x = (1:n)' - (n+1)/2; 
    y = data - mean(data); 
    k = sum(x.*y) ./ sum(x.^2); 
    err = y - k*x; 
end 

linear_regresion這裏是simple linear regression algorithm的實現。

在我的例子中，我使用最大絕對誤差作爲停止標準，但您可以用任何其他擬合函數代替它，例如， mean squared error。

這裏是max_error = 0.04細分數據的例子：

你可以找到更多關於這個和其他分割技術在this survey paper。