%?用3σ法則給出報警線
function?alarm=alarm_sigma（xmn）
%?x為數據信號，m為樣本總數，n為一開始選取的少量正常樣本數
%?example:alarm_sigma（x1600100）
alarm=zeros（m1）;
for?i=1:n
????alarm（i）=mean（x（1:n））+3*std（x（1:n））;
end
for?j=n+1:m
????if?x（j）<=alarm（j-1）
????????alarm（j）=mean（x（1:j）+3*std（x（1:j）））;
????else
????????alarm（j）=alarm（j-1）;
????end
end

plot（x）;?//畫出原始時間序列圖
hold?on;
plot（alarm?‘r-‘）;?//添加自適應預警線?

?????????
????????
????????
????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????495??2017-01-22?23:36??alarm_sigma.m

-----------?---------??----------?-----??----

??????????????????495????????????????????1