P=[1 0.101;
2 0.103;
3 0.105;
4 0.11;
5 0.106;
6 0.099;
7 1.88;
8 0.109;
9 0.111;
10 0.102];
P_X=P(:,1);
P_Y=P(:,2);
y = P_Y-mean(P_Y);% 求解残差
u = std(P_Y);% 求解标准差
id= find(abs(y) > 1.2*u); % 发现粗大误差序号,采用的是1.2σ准则
P_Y(id) = [];% 剔除坏点
P_X(id) = []; % 剔除坏点对应的x
P=[P_X P_Y]; %数据重组
剔除坏点后的P如下图所示: