野值剔除的程序

function outers = fenBuTu(originData)%  采用分布图法剔除粗大误差点
                                     %  输出参数为cell矩阵
                                     %  2004.5
[m,n]=size(originData);
outers=cell(n,2);
for k=1:n % 对各列分别处理    
Data=originData(:,k);    
sortData=sort(Data);   
rownum=size(sortData,1);% 算中位数    
if mod(rownum,2) % 奇数个样本点        
Xm=sortData((rownum+1)/2); % 算四分位数        
lcol=sortData(1:round((rownum+1)/2));        
ucol=sortData(round((rownum+1)/2):end);        
rownum2=size(lcol,1);        
if mod(rownum2,2)             
F1=lcol(round((rownum2+1)/2));            
F0=ucol(round((rownum2+1)/2));        
else            
F1=(lcol(round(rownum2/2+1))+lcol(round(rownum2/2)))/2;            
F0=(ucol(round(rownum2/2+1))+ucol(round(rownum2/2)))/2;        
end
else        
Xm=(Data(round(rownum/2+1))+Data(round(rownum/2)))/2;   % 算四分位数   lcol=sortData(1:round((rownum+1)/2));        
ucol=sortData(round((rownum+1)/2):end);        
rownum2=size(lcol,1);        
if mod(rownum2,2)             
F1=lcol(round((rownum2+1)/2));            
F0=ucol(round((rownum2+1)/2));        
else            
F1=(lcol(round(rownum2/2+1))+lcol(round(rownum2/2)))/2;