prob_analysis.asv

matlab从入门到精通第6章书后源码对初学者很有帮助

% prob_analysis.m
% 产生随机数据
x=[trand(20,1,50) trnd(10,3,100)+5];
% 设置混合概率密度函数
pdf_norm_mix=@(x,p,mu1,mu2,sigma1,sigma2) ...
    p*normpdf(x,mu1,sigma1)+(1-p)*normpdf(x,mu2,sigma2);
% 设置参数数值
pStart=0.5;
muStart=quantile(x,[0.25 0.75]);
sigmaStart-sqrt(var(x)-0.25*diff(nuStart).^2);
start=[pStart muStart sigmaStart sigmaStart];
% 设置参数的上下限
lb=[0 -inf -inf 0   0];
ub=[1  inf  inf inf inf];
% 设置求解属性
options=statset('MaxIter',300,'MaxFunEvals',600);
paramEsts=mle(x,'pdf',pdf_norm_mix,'start',start, ...
    'lower',lb,'upper',ub,'options',options);
% 绘制基础数据的直方图
bins=-2.5:0.5:7.5;
h=bar(bins,histc(x,bins)/(length(x)*0.5),'histc');
set(h,'FaceColor',)