皮尔逊三型曲线matlab编程

皮尔逊三型曲线的matlab源代码，可以供大家利用，进行水文频率计算分析

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