ncff_frame.m

清浊音判别方法,用matlab软件开发,如自相关和能量方法等

% NCCF 函数
% X:输入序列  N：序列长度 tao：时间延迟
function p = ncff_frame(N,tao,X)
 
   sequence = X([(tao + 1) : N]);
   sequence_n = X([1:(N-tao)]);
   accumulate = sequence .*sequence_n;
   result = sum(accumulate);
   % 计算s*s--指列
   sequence_huang = X([tao + 1 :N]);
   sequence_X = sequence_huang.*sequence_huang;
   sum_X = sum(sequence_X);
   
   % 计算s 的延迟平方
   sequence_defer = X([1 :N-tao]);
   square_sequence = sequence_defer.*sequence_defer;
   sum_X_defer = sum(square_sequence);
   p = result/sqrt(sum_X*sum_X_defer);