拒绝用FIND!直接找到矩阵中的位置…

%程序第一部分:做时间序列(亮点是直接列出每月日数)
%程序第二部分:做日数据的月均值(亮点是不用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的发家史,一个希望移植游戏的工程师自己动手搞了套方便的内核,一个希望教学顺利的教授编了套相似的但是开源的系统,一个极客在前人的基础上自己摸索出一套基础操作系统然后全世界的工程师帮他调试更新。这是一个从无到有,充分发挥人类主观能动性的故事,这个故事告诉我们,有想法就去做,利用目前有的条件,创造世上没有的东西,这就是价值~
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值