测试MATLAB的广义极值分布函数GEV

nRows   = 1000;  nColumn= 250; 
data = trnd(8, nRows , nColumn);    %1000行250列，自由度为5的student t的随机数
maxSeries = max(data );                   % 250 列最大值，也就是所谓的极值
paras=gevfit(maxSeries);                  %根据极值参数，参数为形状参数k≠0，尺度参数sigma(σ)，位置参数mu(?)
k=paras(1); sigma=paras(2); mu=paras(3); 
histogram(maxSeries,2:20);            %绘制极大值x的频数分布图,指定x轴指定范围在2-20年内,否则无法匹配高度
x= linspace(2,20,100);                      %产生1000个[2,20]的数据列
%1.根据指定的x序列值，计算GEV分布的概率值（GEV的3参数已知）
y_pdf= nColumn*gevpdf(x, k,sigma,mu );
line(x,y_pdf);                                      %这是概率分布曲线，如下图所示