开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 屏幕保护 > 查看源码
vad.m
资源名称：vad_vqs_lpcc.rar　[点击查看]
上传用户：luvkyhj
上传日期：2013-06-25
资源大小：3k
文件大小：2k
源码类别：

屏幕保护

开发平台：

Matlab

vad.m：源码内容
  1. function [x1,x2]=vad(x)%端点检测子程序
  2. %幅度归一化到(－1，1)
  4. x=double(x);
  5. x=x/max(abs(x));
  7. FrameLen=160;
  8. FrameInc=80;
  10. amp1=10;
  11. amp2=2;
  12. zcr1=10;
  13. zcr2=5;
  14. ener1=0.1;
  15. ener2=0.05;
  17. maxsilence=3;
  18. minlen=15;
  19. status=0;
  20. count=0;
  21. silence=0;
  24. tmp1=enframe(x(1:length(x)-1),FrameLen,FrameInc);
  25. tmp2=enframe(x(2:length(x)),FrameLen,FrameInc);
  26. signs=(tmp1.*tmp2)<0;
  27. diffs=abs(tmp1-tmp2)>0.02;
  28. zcr=sum(signs.*diffs,2);
  31. amp=sum(abs(enframe(filter([1 -0.9375],1,x),FrameLen,FrameInc)),2);
  33. amp1=min(amp1,max(amp)/4);
  34. amp2=min(amp2,max(amp)/8);
  36. ener=sum((enframe(filter([1 -0.9375],1,x),FrameLen,FrameInc)).^2,2);
  38. ener1=min(ener1,max(ener)/8);
  40. x1=0;
  41. x2=0;
  42. for n=1:length(zcr)
  43.     goto=0;
  44.     switch status
  45.         case{0,1}
  46.             if amp(n)>amp1
  47.                 x1=max(n-count-1,1);
  48.                 status=2;
  49.                 silence=0;
  50.                 count=count+1;
  51.             elseif amp(n)>amp2 | zcr(n)>zcr2|ener(n)>ener2
  52.                 status=1;
  53.                 count=count+1;
  54.             else
  55.                 status=0;
  56.                 count=0;
  57.             end
  58.         case 2,
  59.             if amp(n)>amp2 | zcr(n)>zcr2|ener(n)>ener2
  60.                 count=count+1;
  61.             else
  62.                 silence=silence+1;
  63.                 if silence<maxsilence
  64.                     count=count+1;
  65.                 elseif count<minlen
  66.                     status=0;
  67.                     silence=0;
  68.                     count=0;
  69.                 else
  70.                     status=3;
  71.                 end
  72.             end
  73.         case 3,
  74.             break;
  75.     end
  76. end
  77. count=count-silence/2;
  78. x2=x1+count-1;
  80. subplot(313)
  81. plot(x);
  82. axis([1 length(x)/2 -1 1])
  83. ylabel('speech');
  84. xlabel('time');
  85. line([x1*FrameInc x1*FrameInc],[-1 1],'Color','red');
  86. line([x2*FrameInc x2*FrameInc],[-1 1],'Color','red');
  88. %subplot(412)
  89. %plot(amp)
  90. %axis([1 length(amp) 0 max(amp)])
  91. %ylabel('Amplitude');
  92. %xlabel('帧');
  93. %line([x1 x1],[min(amp),max(amp)],'Color','red');
  94. %line([x2 x2],[min(amp),max(amp)],'Color','red');
  96. %subplot(413)
  97. %plot(ener)
  98. %axis([1 length(ener) 0 max(ener)])
  99. %ylabel('Energy');
  100. %xlabel('帧');
  101. %line([x1 x1],[min(ener),max(ener)],'Color','red');
  102. %line([x2 x2],[min(ener),max(ener)],'Color','red');
  104. %subplot(414)
  105. %plot(zcr)
  106. %axis([1 length(zcr) 0 max(zcr)])
  107. %ylabel('ZCR');
  108. %xlabel('帧');
  109. %line([x1 x1],[min(zcr),max(zcr)],'Color','red');
  110. %line([x2 x2],[min(zcr),max(zcr)],'Color','red');