n=100;%n是離散樣本點的個數
k=1.5;%k是閾值
a=rand（1n）*n;
%a=[2.02467.61226.91965.79265.63713.65183.18461.62230.8088.6951];

Q1=22.58779486;%上四分位線上下四分位數可以用excle中QUARTILE函數
Q2=54.38822909;%中線
Q3=77.38073935;%下四分位線
IQR=Q3-Q1;

down=Q1-k*IQR;%上下截斷
up=Q3+k*IQR;

b=a（find（ac=a（find（a>Q3））;
d=a（find（a>Q1&a%數組連接起來bdc
aa=[bdc];
hold?on
plot（[0n+3][Q1Q1]‘--‘）;%加3是讓整幅圖畫面加大，可以根據需要+4……＋5……
plot（[0n+3][Q2Q2]‘--‘）;
plot（[0n+3][Q3Q3]‘-