对366（天）×53（年）求月平均

输入样本为 366×53的矩阵，求月平均

monstruct01 = [0 31 28 31 30 31 30 31 31 30 31 30 31]; % 平年月
acc_m01 = cumsum(monstruct01);
monstruct02 = [0 31 29 31 30 31 30 31 31 30 31 30 31]; % 闰年月
acc_m02 = cumsum(monstruct02);

for year = 1960:2012   % 起始年份
    sign = leapyear(year);
    for i = 2:13
        if  sign == 0
            acc_m = acc_m01;
        else
            acc_m = acc_m02;
        end
        beg = acc_m(i-1)+1;
        ter = acc_m(i);
        montlymean(i-1,year-1959) = mean(longtan_q(beg:ter,year-1959));
    end
end