双累积曲线的绘制

双累积曲线是检验两个参数间关系一致性及其变化的常用方法，它可用于水文气象要素一致性的检验，在水文学上，通过累加的方式去增加展现出序列的规律性。具体的理论知识可以参考穆兴民老师的关于双累积曲线的相关论文，本文主要记录自己初在做双累积曲线时候，初步学习Matlab，写自己在的一点心得。不好勿喷。纯属爱好

注：流域径流数据中，第一列为年份，第二列为年径流量，第三列为气温，第四列为降水量。

clc,clear

A=xlsread('E:\原始数据\气象水文数据\流域径流1959-2018.xlsx');

data_flow = A(:,2);

data_tem = A(:,3);

data_prec = A(:,4);

flow1=data_flow(1);

tem1=data_tem(1);

prec=data_prec(1);

%

Sum_flow=zeros(58,1);

Sum_tem=zeros(58,1);

Sum_prec=zeros(58,1);

for i=1:58;

  Sum_flow(i)= sum(data_flow(1:i));

  Sum_tem(i)= sum(data_tem(1:i));

  Sum_prec(i)= sum(data_prec(1:i));

end

years=(1959:2016)';

xlswrite('cumu.xlsx',[years,Sum_flow,Sum_tem,Sum_prec])