cal_s.m

用 CC法计算延迟时间和潜入位数的matlab程序

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;