#广义极值分布#matlab
如何选择GEV分布作为最优概率分布拟合月径流序列，并分别计算其概率密度最大处所对应流量作为各月最适宜生态流量？
要求：必须是能够绘制出拟合曲线，并把结果表示出来（类似正态分布那种的）
这个有现成的 matlab 工具包，你可以找找，直接调用
不知道你这个问题是否已经解决, 如果还没有解决的话:究竟哪一部分的值才作为尾部数据，来进行帕累托分布呢？
这时候就要定一个阈值u，超出这个阈值的数据点才作为拟合的对象，并且将其减去阈值，然后得到的值才是我们要去拟合的数据点。
%% 读取数据
load('yb.mat');
yb1 = yb;
czyb=yb1(yb1>u)-u; %准备拟合的数据点
ybs=length(czyb); %超出阈值的数据点个数
%% 参数估计(最大似然估计）
paramEsts = gpfit(czyb);
kHat = paramEsts(1); %形状系数
sigmaHat = paramEsts(2); %尺度系数
这里的拟合非常简单，其实我之前不知道matlab中有这个的内置包，我就自己写了最小二乘法，结果是一样的，但是运行速度较慢。参数估计很简单就完成了。