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) %比较密度函数
用于简单计算某个分位数下的密度值,数据越多,其结果越接近。