给出matlab代码，使用多年水位年极大值，通过皮尔逊三型统计方法，计算100年重现期水位

由于缺少数据，无法编写完整的代码。以下是可能的代码框架，您可以根据自己的数据进行调整。

% 假设我们有一个名为"water_levels"的向量，其中包含多年水位数据

% 计算年极大值
annual_max = zeros(1, num_years);
for i = 1:num_years
annual_max(i) = max(water_levels(year_starts(i):year_ends(i)));
end

% 计算经验分布函数（EDF）
edf = 1 - ecdf(annual_max);

% 计算皮尔逊三型参数
pearson_params = gpfit(edf);

% 计算100年重现期水位
return_period = 100;
rp_prob = 1 - 1/return_period;
rp_water_level = gpinv(rp_prob, pearson_params(1), pearson_params(2), pearson_params(3));