IDL编程学习之逐天进行8小时平均

15195775117

附件chao.txt是：日期-时间-臭氧浓度
（压缩包中有原始数据、代码、结果txt）

代码：
pro chao
;读出年月日  
a=read_ascii('C:\Users\YGDY\Desktop\chao.txt',DELIMITER='/',data_start=1)
  x=a.field1
  year=x[0,*]
  month=x[1,*]
  day=x[2,*]
;读出小时=0和23的索引，两者一头一尾索引对应起来，等于按天切割
  a=read_ascii('C:\Users\YGDY\Desktop\chao.txt',data_start=1)
  x=a.field1
  p1=where(x[1,*] eq 0)
  p2=where(x[1,*] eq 23)
  n=n_elements(p1)
  for i=0,n-1 do begin
    title='滑动平局'+strtrim(string(long(year[p1])),2)+$
      '年'+strtrim(string(long(month[p1])),2)+'月'+$
      strtrim(string(long(day[p1])),2)+'日'
    txt=txt(title,smooth8(transpose(x[2,p1:p2])))
  endfor
end
;对某向量求8位滑动平均（可以重复使用哟）
function smooth8,x
n=n_elements(x)
re=-999
for i=0,n-8 do begin
  re=[re,mean(x[i:i+7])]
endfor
return,transpose(re[1:-1])
end
function txt,name,x
  ;/append关键字：续写！
  openw,lun,'C:\Users\YGDY\Desktop\chao\'+name+'.txt',/get_lun;,/append
  printf,lun,x
  close,/all
end