一开始理解有歧义，以为2008Q1季度，就是求第一个值的均值和标准差，所以均值永远是第一个值，标准差
没有计算，一直为初值0
改后的程序：

data ex;
input
date $ x;
cards;
2008Q1 1
2008Q1 1
2008Q2 2
2008Q2 2
2008Q3 3
2008Q4 4
2008Q4 3
2009Q1 2
2009Q2 1
2009Q3 3
2009Q3 2
;
run;

data ex1;
set ex;
year=input(scan(date,1,'Q'),4.);
q=input(scan(date,2,'Q'),1.);
run;

proc iml;
use ex1;
read all var{year q x} into g;
q=j(nrow(g),2,0);     
/*给q矩阵赋初值为0，用来存储mean，std*/
do i=1 to nrow(g);   
  /*从矩阵第一行循环到最后一行*/
   if i=1 then do;   
   /*如果是第一行*/
   n=((g[,1]=g[i,1])&(g[,2]=g[i,2]))[+];  
   /*这里计算year=2008且Q=1的记录条数n*/
   q[i,1]=g[1:n,3][:];              
   /*求前n行的均值作为2008Q1的均值*/
   q[i,2]=((g[1:n,3]-g[1:n,3][:])##2)[:]##0.5;
   /*这里其实是求前n行的std，iml中的std函数不是很好用，这里其实是按标准差公式算的，即每个记录减去均值，求平方和的均值，最后开根号*/
   end;
   else if all(g[i,1:2]=g[i-1,1:2]) then do;   
   /*如果不是第一行，且该行year的值等于上一行year的值，Q的值也相等，那该行的均值标准差为空值*/
   q[i,1]=.;                                    
   q[i,2]=.;
   end;
  else do;                                    
  /*如果都不属于上述情况，即每年，每个季度的第一个值，就求它上面所有行的均值以及标准差*/
   q[i,1]=g[1:i-1,3][:];
   q[i,2]=((g[1:i-1,3]-g[1:i-1,3][:])##2)[:]##0.5;   
   /*之前的程序其实求的方差，这边添加##0.5,表示矩阵每个元素开0.5次方，即开根，也可以用sqrt函数，效果一样*/
end;
end;
create ans from q[colname={'mean','std'}]; /*将q矩阵导入到ans数据集中，列名称分别为mean 和 std*/
append from q;
quit;
data ans;
merge ex ans;
run;

您好！首先真的非常感激您大半夜的回复，请注意休息喔！

这里第一行的问题是解决了，但是我发现了另一个问题，就是您说的“ /*如果都不属于上述情况，即每年，每个季度的第一个值，就求它上面所有行的均值以及标准那*/，就没能体现移动平均的效果了，而是对它之前的所有行进行了计算。我希望得到的是  比如对于 2009q3,mean 和std 是根据 2008Q3——2009Q2或者是2008Q4——2009Q3（两种定义，均只包含4个季度，我提出问题的时候用的是第一种定义，但是现在更希望使用第二种定义方法），而您上面的程序是对2008Q1~2009Q2求的平均和std。又如2012Q2的mean和std，想要的是根据2011Q3-2012Q2来计算的。
不管如何，真心感谢您的帮助！

我觉得2008Q1有计算值的话，逻辑和随后的不一致。而且计算值和2008Q2总是相同的。京剧

亲，我可以接受稍微改变下移动平均的定义，比如2008q2的mean std是由2008Q1 和2008Q2的数据来计算的,又如2012Q2的mean和std，想要的是根据2011Q3-2012Q2来计算的。
每季度的记录超过1W,而且具体多少条是不定的，所以我觉得不能采用你说的那个办法。

其实要点不是1w还是两万。我的最大的n可以是十万，一百万。如果10w，10×4×5 =200w。即使是2千万，程序的效率也没有问题。
如果是如你的定义。程序改动即可

用新的定义，好像我的程序有点问题。京剧其实我的第一个程序，2008Q2计算用的是2008Q1的数据。2008Q3用的是2008Q1+Q2数据依次类推。.  

非常感谢您！这下99。99%对了。
标准差公式是减去均值算平方和然后除以自由度（n-1）,在小样本的时候还是有差别，我的样本很大，没什么影响。如果精益求精的话，您可以再修改下
不过太感谢您了，我会把您的回答设置成最佳答案。
再次感谢！

好吧，无偏估计的话，是要除以n-1

将计算标准差的公式改为q[i,2]=(((g[locate,3]-g[locate,3][:])##2)[+]/(ncol(locate)-1))##0.5;