function data=trim(data,outval)
%data=trim(data,outval)
% 祛除壞數據，包括NaN, Inf, 和異常大小數據
% data: 列狀數據，每列來自一個總體,返回凈化后的數據
% outval: 系數因子，離均值超過outval倍標準差被判為異常大小，默認為4

% L.J.HU 8-17-1999

if nargin<2,outval=4;end
outliers=(isnan(data)|abs(data)==inf);
[n,m]=size(data);mu=mean(data);sigma=std(data);
outliers=outliers+(abs(data-ones(n,1)*mu)>outval*ones(n,1)*sigma);
if m>1,
   data(any(outliers'),:)=[];
else
  data(find(outliers'),:)=[];
end