检测时域离散信号的周期波动性

[屋卡]

如这样的一组数据生成信号
数据

引用

http://pan.baidu.com/s/1i431rc5
信号图像

要检测的

load('data.mat');
[row,col,w]=size(T);
x_sub=1:col;  %原图
plot(x_sub,T(x_sub));
T2=zeros(1,col);
for xy=11:col-10
    num=0;
    for i=xy-10:xy+10
        num=num+T(i);
    end
    T2(xy)=num/21;
end
hold on
x_sub=1:col;  %滤波
plot(x_sub,T2(x_sub));

要检测它的周期波动性，现在的想法是对信号滤波后测峰峰值，但是不知应该咋样梳理