该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
function p20150213
data=[3.46980194700000 8.87299024100000 5.45267387500000 ...
7.03878845400000 2.80902866300000 1.30523846700000 ...
1.93331379800000 1.68563321400000 3.53578284000000 ...
2.85851813600000 1.90010907000000 0.512363610000000 ...
1.70180435900000 2.44555779200000 5.07233999800000 ...
1.90042332700000 2.84204547400000 5.17175281800000 ...
11.7145459300000 4.72564189100000 6.04745635800000 ...
4.03151241200000 6.98932446100000 5.98141742700000 ...
3.56941330200000 2.42917961900000 4.85771550700000 ...
2.41223663300000 1.18983746400000 4.74211703100000 ...
2.89138595900000 2.56091915900000 3.63523954300000 ...
5.84919730000000 3.25543438300000 4.42831612000000 ...
6.21293288200000 3.60205852500000 0.809509249000000 ...
2.18447775300000 2.31312356000000 2.28030315900000 ...
7.50142937400000 1.52016695300000 2.95776002600000 ...
3.51898884200000 1.25590895300000 2.90813182700000];
d=sort(data); %给data排序,这里可以不排。
x=0:0.01:12;
s=gamfit(d); %估计gamma参数
figure(1);clf
hist(data,50);
hold on
s=10*gampdf(x,s(1),s(2)) % 乘以10拉伸图像,更直观一些。
plot(x,s,'r','linewidth',1.5)
set(gca,'yticklabel',{})
end