以下為matlab實現程式:
clc;
clear;
%利用最大似然估計,檢測一元離群點
%氣溫測試資料,單位為攝氏度,這裡假設資料服從正太分布
tdata=[24,24.1,28.9,28.9,29.0,29.1,29.1,29.2,29.2,29.3,29.4,29,28,29.2,29.3,29,29.5,29.3,29,29,29,29,29,29.4,29.1];
len=length(tdata);
%利用最大似然估計計算均值和方差
phat=mle(tdata);
mu=phat(1);
sigma=phat(2);
%定義離群點集
outliers=;
%落在 mu+/-3sigma區間範圍之外的資料即為離群點
for i=1:len
a=abs(tdata(i)-mu)/sigma;
if a>3
outliers(length(outliers)+1)=tdata(i);
endend
證明隨機數符合正太分布的方法
因為工作需要,自己用c語言實現了matlab中的normrnd函式 隨機正太分布函式,就是生成的隨機數符合正太分布 其實相關的函式實現網上有很多,我使用了box muller方法。需要的可以網上找,這裡不做描述 本文主要是驗證自己實現的c函式有沒有真正的實現這個功能。方法就是大量資料描圖 1.用c函...
Python實現非正太分布的異常值檢測方式
工作中,我們經常會遇到資料異常,比如說瀏覽量突增猛降,交易量突增猛降,但是這些資料又不是符合正太分布的,如果用幾倍西格瑪就不合適,那麼我們如何來判斷這些變化是否在合理的範圍呢?小白查閱一些資料後,發現可以用箱形圖,具體描述如下 箱形圖 英文 box www.cppcns.complot 又稱為盒須圖...
一元錢的分法
問題 將一元錢分成1分,2分,5分有幾種分發法?思想 定義乙個n,用i,j,k分別表示1分,2分,5分的個數,i 100 2 j 5 k 寫3個for迴圈,在最內迴圈中用if語句判斷是否i j 2 k 5 100,如果是n就 1 最後把n列印出來,n即是一共多少種分法。1.三種面額都存在的情況 i,...