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序列的長度n

%計算經驗頻率，并算出相應的位置坐標，并繪出經驗點
p=[[1:n]/（n+1）];????????%求經驗頻率p
xp=norminv（p01）;??????%求出經驗頻率p對應的位置坐標xp
%橫坐標：x=norminv（pmusigma）標準正態分布時，mu=0sigma=1）
plot（xpX‘.‘）
hold?on?

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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????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