matlab 检查数据的周期函数

function T=periodtest（testmat）
%FUNCTION T=PERIODTEST（TESTMAT）
%TESTMAT 是输入的单变量数据，函数用于统计数据最可能存在的周期
%T返回的是周期，函数缺省画出周期图像
wolfer=testmat;
%year=sunspot（:1）;wolfer=sunspot（:2）;
%plot（yearwolfer）title（‘Sunspot Data‘）pause
Y=fft（wolfer）; %Y（1）=[];
n=length（Y）;
power=abs（Y（1:n/2））.^2;
nyquist=1/2;
freq=（1:n/2）/（n/2）*nyquist;
figureplot（freqpower）
xlabel（‘data‘）;
title（‘periodvalue‘）;
period=1./freq;
plot（periodpow