1. 程式人生 > >利用matlab計算在任一資料上的密度函式值

利用matlab計算在任一資料上的密度函式值

clear; u=random(‘Normal’,0,1,1,1000)’;%原始資料 e1=u; quant=0.1; [f2,x2] = ksdensity(e1); qqq=quantile(e1,quant); [xx,index]=sort(abs(x2(?-qqq)); A1=[x2(index(1)) x2(index(2))]; A2=[index(1) index(2)]; %利用插值法 stao=(f2(A2(2))(qqq-A1(1)) + f2(A2(1))(-qqq+A1(2)))/(A1(2)-A1(1)) normpdf(norminv(quant,0,1),0,1) %比較密度函式

資料越多,其結果越接近。