function [mean_S_t,mean_delt_S_t, S_cor_t] = Cal_S_ave(x,t,sigma)
% sigma = std(x);
mean_S_t = 0;
mean_delt_S_t = 0;
max_S_m_r_t = 0;
min_S_m_r_t = 1e10;
for m = 2:1:5
    for j = 1:1:4
        r = j*sigma/2;
        S_m_N_r_t = Cal_S(m,x,r,t);
        if j==1
            max_S_m_r_t = S_m_N_r_t;
            min_S_m_r_t = S_m_N_r_t;
        else
            if S_m_N_r_t>max_S_m_r_t
                max_S_m_r_t = S_m_N_r_t;
            end
            if S_m_N_r_t<min_S_m_r_t
                min_S_m_r_t = S_m_N_r_t;
            end
        end
        mean_S_t = mean_S_t+S_m_N_r_t;
    end
    mean_delt_S_t = mean_delt_S_t + (max_S_m_r_t-min_S_m_r_t);
end
mean_S_t = mean_S_t/16;
mean_delt_S_t = mean_delt_S_t/4;
S_cor_t = mean_delt_S_t+abs(mean_S_t);