%% 測試MATLAB的廣義極值分佈函數GEV
blocksize = 1000;
nblocks = 250;
t = trnd(5,blocksize,nblocks);%1000行250列，自由度爲5的student t的隨機數
x = max(t); % 250 列最大值，也就是所謂的極值
paramEsts=gevfit(x) %根據極值樣本擬合最佳參數，參數1-3分別爲：形狀參數k≠0，尺度參數sigma(σ)，位置參數mu(µ)
hist(x,2:20); %繪製極大值x的頻數分佈圖,指定X範圍在2-20年內
xgrid = linspace(2,20,1000);%產生2-20數據列，數量爲1000個
%根據指定的x值和3個參數，獲取對應的概率值，再乘總數就是頻數，繪製出曲線圖
line(xgrid,nblocks*gevpdf(xgrid,paramEsts(1),paramEsts(2),paramEsts(3)) );
ps:當k<0時，GEV是PⅢ極值分佈；K>0時，是PⅡ極值分佈；