function S_m_N_r_t = Cal_S(m,x,r,t)
N = length(x);
S_m_N_r_t = 0;
l = floor(N/t);    %子序列長度
for s = 1:1:t
    X = x(s:t:s+(l-1)*t);  %保證子序列長度一樣
    Cs_m_N_r_t = correlative_integral(m,X,r,t);
    Cs_1_N_r_t = correlative_integral(1,X,r,t);
    S_m_N_r_t = S_m_N_r_t+(Cs_m_N_r_t-Cs_1_N_r_t^m);
end
S_m_N_r_t = 1/t*S_m_N_r_t;