%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?以下為測試3sigma法去除粗大誤差的代碼
%?使用時，將其拷貝到命令窗口，回車執行即可
%?clear
%?clc
%?X?=?ones（1?50）;
%?X（1）?=?1.3;
%?X（2）?=?1.1;
%?X（7）?=?1.8;
%?X（10）?=?1.5;
%?X（12）?=?0.7;
%?X（13）?=?0.9;
%?X（18）?=?0.5;
%?X（20）?=?0.3;
%?[dest?index]?=?EliminateGrossErrors（X）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%?去除粗大誤差
function?[dest?index]?=?EliminateGrossErrors（src）
dest?=?src;?????%?去除粗大誤差之后的數據
index?=?[];?????%?每次去除的數據編號
%?3sigma法則
while?1
????y?=?dest?-?mean（dest）;??????????????????%?殘差
????s?=?sqrt（sum（y?.^?2）?/?（numel（y）?-?1））;?%?標準差
????id?=?find（abs（y）?>?3?*?s）;??????????????%?粗大誤差序號
????if?isempty（id）??????????????????????????%?判斷是否還含有粗大誤差
????????break;
????else
????????if?isempty（index）???????????????