事件定义为持续或短暂的满足某一阈值条件。例如极端高温为定义为日最高气温大于95%分位数的记录,若有一天的气温条件满足上式,则算是一次极端高温事件,若连续10天都满足大于95%分位的条件,则这10天算是一次连续的极端高温事件。程序的目的就是算出一个区域内事件的发生频率、持续时间和起始日期。
%数据预处理。满足事件条件的标为1,不满足的标为0,存入列向量a中
%%连续事件个数统计,输入量a为列向量,单个站点事件,发生记为1,否则记为0
%%
%k为起始坐标,m为终止坐标
%C:连续个数,Z起始位置,S对应持续长度
function [C Z S]=lianxu1(a)
n=length(a);
i=1;
z=[]; %连续事件发生的位置
size=[];%对应z的事件持续月份
COUNT=0;%连续的个数
while(i
if(a(i+1)-a(i)==1)%事件开始
k=i;%k为起始坐标
for j=k:n-1
%数据预处理。满足事件条件的标为1,不满足的标为0,存入列向量a中
%%连续事件个数统计,输入量a为列向量,单个站点事件,发生记为1,否则记为0
%%
%k为起始坐标,m为终止坐标
%C:连续个数,Z起始位置,S对应持续长度
function [C Z S]=lianxu1(a)
n=length(a);
i=1;
z=[]; %连续事件发生的位置
size=[];%对应z的事件持续月份
COUNT=0;%连续的个数
while(i