function out_data=huadong(indata,k)
%indata为输入数据,两列以上,第一列为时间,第二列开始为观测数据
%k为滑动窗口的宽度
[m,n]=size(indata);
tmp_k=floor(k/2);
out_data=nan(m,n);
for i=1:m
out_data(i,1)=indata(i,1);
end
for j=2:n
for i=1:m-2*tmp_k
tmp_sum=0;
for t=i:i+k-1
tmp_sum=tmp_sum+indata(t,2);
end
tmp_sum=tmp_sum/k;
out_data(i+tmp_k,j)=tmp_sum;
end
end
end