function?P_fit
X=[18500?17700?13900?13300?12800?12100?12000?11500?11200?10800?10800?10700?10600?10500?9690?8500?8220?8150?8020?8000?7850?7450?7290?6160?5960?5950?5590?5490?5340?5220?5100?4520?4240?3650?3220];????????%輸入流量的資料
X=sort（X）;??????????????%sort是把流量從小到大排列
X=fliplr（X）;????????????%把流量x從大到小排列
n=length（X）;????????????%求出流量資料X序列的長(zhǎng)度n

%計(jì)算經(jīng)驗(yàn)頻率，并算出相應(yīng)的位置坐標(biāo)，并繪出經(jīng)驗(yàn)點(diǎn)
p=[[1:n]/（n+1）];????????%求經(jīng)驗(yàn)頻率p
xp=norminv（p01）;??????%求出經(jīng)驗(yàn)頻率p對(duì)應(yīng)的位置坐標(biāo)xp
%橫坐標(biāo)：x=norminv（pmusigma）標(biāo)準(zhǔn)正態(tài)分布時(shí)，mu=0sigma=1）
plot（xpX‘.‘）
hold?on?

%繪制海森機(jī)率格紙
x=0:10:40000;???????????%縱坐標(biāo)
m=mean（X）???????????????%計(jì)算樣本X的均值
K=X/m;??????????????????%計(jì)算模比系數(shù)K
Cv=sqrt（（1/（n-1））.*sum（（K-1）.^2））;??????%計(jì)算變差系數(shù)Cv（利用的是課本51頁(yè)樣本的無(wú)偏估計(jì)）
%Cs=（sum（（K-1）.^3））/（Cv^3）/（n-3）;????????%計(jì)算偏態(tài)系數(shù)Cs
Cs=Cv;
alfa=4/（Cs^2）;
beta=（m*Cv*Cs）/2;
a0=m-2*Cv*m/Cs;
P=1-gamcdf（x-a0alfabeta）;?????%計(jì)算理論頻率P
XP=norminv（P01）;??????????????%計(jì)算理論頻率P所對(duì)應(yīng)的坐標(biāo)XP
plot（XPx‘r-.‘）
gtext（‘Cs=Cv‘‘FontSize‘8）
Cs=2*Cv;
alfa=4/（Cs^2）;
beta=（m*Cv*Cs）/2

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????156052??2012-11-19?15:42??p3\p3圖.fig

?????文件??????99151??2012-11-19?15:42??p3\P3圖形.jpg

?????文件???????2369??2012-11-19?15:37??p3\P_fit.m

?????文件???????2369??2012-11-19?15:46??p3\P_fit.txt

?????目錄??????????0??2013-07-19?15:55??p3

-----------?---------??----------?-----??----

???????????????259941????????????????????5