[½»Á÷] ¡¾ÇóÖú¡¿ÓÃmatlab Ä£Äâ¸ß˹ÄâºÏ ÒÑÓÐ2È˲ÎÓë

function gussfit
%Õý̬·Ö²¼ÄâºÏ
%%PSL@CSU
%%QQ:547423688
%%Email:anyuezhiji@qq.com
%%Edit   @ 2010.4.19
x=[1,2,3,4,5,6,7];
y=[1,3,5,10,6,4,2];
%ÕâÑù
% R=[];
% for i=1:length(x)
%     R=[R repmat(x(i),1,y(i))];
% end
% [muhat,sigmahat,muci,sigmaci] = normfit(R)
%»ò
alpha=0.05;
censoring=[];
[muhat,sigmahat,muci,sigmaci] = normfit(x,alpha,censoring,y);
disp(['ÔÚ',num2str(100*(1-alpha)),'%µÄÖÃÐŶÈÏ£º']);
disp(['·åÖµmuΪ  : ',num2str(muhat,'%10.5f')]);
disp(['ÖÃÐÅÇø¼äΪ:  [',num2str(muci','%10.5f'),']']);
disp(['sigmaΪ  : ',num2str(sigmahat,'%10.5f')]);
disp(['ÖÃÐÅÇø¼äΪ:  [',num2str(sigmaci','%10.5f'),']']);
[mu1,sigma1]= normpdffit(x,y/sum(y),muhat,sigmahat);
disp(['¼ÌÐøµ÷ÓÃ×Ô¶¨Ò庯ÊýÄâºÏ£º']);
disp(['·åÖµmuΪ  : ',num2str(mu1,'%10.5f')]);
disp(['sigmaΪ  : ',num2str(sigma1,'%10.5f')]);
figure;
%»æÖÆÖ±·½Í¼
bar(x,y/sum(y));
hold on;
%»æÖƶÔÓ¦µã
plot(x,y/sum(y),'ro');
%¶¨Òå»æͼÇø¼ä
t=min(x)max(x)-min(x))/50:max(x);
%ÄâºÏ·Ö²¼ÇúÏßÊý¾Ý
nyf=normpdf(t,muhat,sigmahat);
plot(t,nyf,'r-');
mynyf=normpdf(t,mu1,sigma1);
plot(t,mynyf,'k-.','LineWidth',2);
legend('hist','hist value','normfit','mynormfit');
function [muhat,sigmahat]= normpdffit(x,y,mu0,sigma0)
%Õý̬·Ö²¼pdf²ÎÊý¹À¼Æ
p=[mu0 sigma0];
opt=optimset('TolX',0.001,'Display','off');
[x,fval]=fminsearch(@mynormfit,p,opt,x,y);
muhat=x(1);
sigmahat=x(2);
function Variance= mynormfit(p,x,y)
%¼ÆËã¶ÔÓÚÈκÎÒ»×é²ÎÊýp(1)£¬p(2)£¬¸ø³öµ±Ç°Êý¾Ýϵķ½²îÀ´¡£
Variance=var(y-normpdf(x,p(1),p(2)));

ÔÚ95%µÄÖÃÐŶÈÏ£º
·åÖµmuΪ  : 4.19355
ÖÃÐÅÇø¼äΪ:  [3.65430   4.73280]
sigmaΪ  : 1.47013
ÖÃÐÅÇø¼äΪ:  [1.17480   1.96509]
¼ÌÐøµ÷ÓÃ×Ô¶¨Ò庯ÊýÄâºÏ£º
·åÖµmuΪ  : 4.17697
sigmaΪ  : 1.45444