[n,t]=size(stationmonthtp);
b=0;
i=1;
num=zeros(12,39);
for y=1979:1:2017
for m=1:1:12
num(m,i)=sum(stationmonthtp(stationmonthtp(:,2)==m,4));%每月的降水和存到num中
end
i=i+1;
end
yy = stationmonthtp(:,1);
mm = stationmonthtp(:,2);
pp = stationmonthtp(:,4)
for year = 1979:2017
yb = year-1978;
for month = 1:12
ind = yy == year & mm == month; 可以看看有没有选对位置
pms(month,yb) = sum(pp (ind));
end
end
不知这样是否可行,给个例子