[filename,filepath]=uigetfile('*.*','请选择文件');
%计算水文极端,年最大,年连续5日最大,连续干日,连续极端径流低值日数
%wqianghy@163.com
data=dlmread(datapath);
data(1)=[];
datalength=length(data);
numyear=datalength/365;
AM1X=zeros(numyear,1); %求年最大值
for i=1:numyear
tmin=365*(i-1)+1;
tmax=365*i;
AM1X(i)=max(data(tmin:tmax));
end
AM5X=zeros(numyear,1); %求连续5日最大
data5x=zeros(datalength,1);
for i=1:datalength
if i<5
data5x(i)=mean(data(1:5));
else
tmi=i-5+1;
tmx=i;
data5x(i)=mean(data(tmi:tmx));
end
end
for i=1:numyear
tmin=365*(i-1)+1;
tmax=365*i
%计算水文极端,年最大,年连续5日最大,连续干日,连续极端径流低值日数
%wqianghy@163.com
%南京大学流域水文过程模拟实验室
datapath=strcat(filepath,filename);data=dlmread(datapath);
data(1)=[];
datalength=length(data);
numyear=datalength/365;
AM1X=zeros(numyear,1); %求年最大值
for i=1:numyear
tmin=365*(i-1)+1;
tmax=365*i;
AM1X(i)=max(data(tmin:tmax));
end
AM5X=zeros(numyear,1); %求连续5日最大
data5x=zeros(datalength,1);
for i=1:datalength
if i<5
data5x(i)=mean(data(1:5));
else
tmi=i-5+1;
tmx=i;
data5x(i)=mean(data(tmi:tmx));
end
end
for i=1:numyear
tmin=365*(i-1)+1;
tmax=365*i