%程序第一部分:做时间序列(亮点是直接列出每月日数)
%程序第二部分:做日数据的月均值(亮点是不用if和find而直接反出数据位置,最大程度优化程序)
%function sum=m_sum(data,y1,y2)
size_d=size(data);�ta为没有日期的数据,列为站点或者矩阵点的数据,y1=1961;y2=2013;
DATA(size_d(1),size_d(2)+3)=0;
DATA(:,4:end)=data;
%直接为timemake程序
c1=[31 29 31 30 31 30 31 31 30 31 30 31];
c2=[31 28 31 30 31 30 31 31 30 31 30 31];
time=[];
for i=y1:y2
%
n=0;
if(((mod(i,4)==0&mod(i,100)~=0))|mod(i,400)==0)
c=c1;
else
c=c2;
end
for j=1:12
n=0;
for k=1:c(j)
n=n+1;
time=[time;i j n];
end
end
end
DATA(:,1:3)=time;
%计算月均值
ans1=[];
for y=y1:y2
for mon=1:12
ans1=[ans1;mean(DATA(time(:,1)==y&time(:,2)==mon,:))];%<---就是这个,一步解决所有问题,感觉超棒
end
end
今天看到了linux的发家史,一个希望移植游戏的工程师自己动手搞了套方便的内核,一个希望教学顺利的教授编了套相似的但是开源的系统,一个极客在前人的基础上自己摸索出一套基础操作系统然后全世界的工程师帮他调试更新。这是一个从无到有,充分发挥人类主观能动性的故事,这个故事告诉我们,有想法就去做,利用目前有的条件,创造世上没有的东西,这就是价值~
%程序第二部分:做日数据的月均值(亮点是不用if和find而直接反出数据位置,最大程度优化程序)
%function sum=m_sum(data,y1,y2)
size_d=size(data);�ta为没有日期的数据,列为站点或者矩阵点的数据,y1=1961;y2=2013;
DATA(size_d(1),size_d(2)+3)=0;
DATA(:,4:end)=data;
%直接为timemake程序
c1=[31 29 31 30 31 30 31 31 30 31 30 31];
c2=[31 28 31 30 31 30 31 31 30 31 30 31];
time=[];
for i=y1:y2
%
end
DATA(:,1:3)=time;
end