statistic_zm.m

来自「清浊音判别方法,用matlab软件开发,如自相关和能量方法等」· M 代码 · 共 27 行

function statistic_zm(sequence_m,sequence_z)
      length_m = length(sequence_m);
      length_z = length(sequence_z);
      
      for n = 0:0.01:1
          if n == 0
             m_sum(1) = sum(sequence_m == 0 );    
          else
              space_f = sum(sequence_m < (n - 0.01));
              space_b = sum(sequence_m <= n);
              m_sum((n/0.01) + 1) = space_b - space_f;
          end
      end
      m_sum = m_sum/length_m;
      subplot(2,1,1);
      plot((0:0.01:1),m_sum);
      title('m-sum-statistic');
      
      for k = 1 : 200
          zeros_f = sum(sequence_z < (k - 1));
          zeros_b = sum(sequence_z <= k );
          z_sum(k) = zeros_b - zeros_f;
      end
      z_sum = z_sum/length_z;
      subplot(2,1,2);
      plot(z_sum);
      title('z-su-_statistic');