欧美极品高清xxxxhd,国产日产欧美最新,无码AV国产东京热AV无码,国产精品人与动性XXX,国产传媒亚洲综合一区二区,四库影院永久国产精品,毛片免费免费高清视频,福利所导航夜趣136

標題: 利用MATLAB識別1-9語音的源碼 [打印本頁]
作者: matlabxmut    時間: 2018-5-22 22:16
標題: 利用MATLAB識別1-9語音的源碼
自己做的，識別1-9的語音。


全部資料51hei下載地址：
MATLAB語音識別1-9.zip (397.37 KB, 下載次數: 31)


源碼：
  1. function trimmedX = myVAD(x)

  3. Ini = 0.1;          %Initial silence duration in seconds
  4. Ts = 0.01;          %Frame width in seconds
  5. Tsh = 0.005;        %Frame shift in seconds
  6. Fs = 16000;         %Sampling Frequency
  7. counter1 = 0;   
  8. counter2 = 0;
  9. counter3 = 0;
  10. counter4 = 0;
  11. ZCRCountf = 0;      %Stores forward count of crossing rate > IZCT
  12. ZCRCountb = 0;      %As above, for backward count
  13. ZTh = 40;           %Zero crossing comparison rate for threshold
  14. w_sam = fix(Ts*Fs);                   %No of Samples/window
  15. o_sam = fix(Tsh*Fs);                  %No of samples/overlap
  16. lengthX = length(x);
  17. segs = fix((lengthX-w_sam)/o_sam)+1;  %Number of segments in speech signal
  18. sil = fix((Ini-Ts)/Tsh)+1;            %Number of segments in silent period
  19. win = hamming(w_sam);

  21. Limit = o_sam*(segs-1)+1;             %Start index of last segment

  23. FrmIndex = 1:o_sam:Limit;             %Vector containing starting index
  24.                                       %for each segment
  25. ZCR_Vector = zeros(1,segs);           %Vector to hold zero crossing rate
  26.                                       %for all segments
  27.                                     
  28. %Below code computes and returns zero crossing rates for all segments in
  29. %speech sample  
  30. for t = 1:segs
  31.     ZCRCounter = 0;
  32.     nextIndex = (t-1)*o_sam+1;
  33.     for r = nextIndex+1:(nextIndex+w_sam-1)
  34.         if (x(r) >= 0) && (x(r-1) >= 0)
  35.          
  36.         elseif (x(r) >= 0) && (x(r-1) < 0)
  37.          ZCRCounter = ZCRCounter + 1;
  38.         elseif (x(r) < 0) && (x(r-1) < 0)
  39.          
  40.         elseif (x(r) < 0) && (x(r-1) >= 0)
  41.          ZCRCounter = ZCRCounter + 1;
  42.         end
  43.     end
  44.     ZCR_Vector(t) = ZCRCounter;
  45. end

  47. %Below code computes and returns frame energy for all segments in speech
  48. %sample
  49. Erg_Vector = zeros(1,segs);
  50. for u = 1:segs
  51.     nextIndex = (u-1)*o_sam+1;
  52.     Energy = x(nextIndex:nextIndex+w_sam-1).*win;
  53.     Erg_Vector(u) = sum(abs(Energy));
  54. end

  56. IMN = mean(Erg_Vector(1:sil));  %Mean silence energy (noise energy)
  57. IMX = max(Erg_Vector);          %Maximum energy for entire utterance
  58. I1 = 0.03 * (IMX-IMN) + IMN;    %I1 & I2 are Initial thresholds
  59. I2 = 4 * IMN;
  60. ITL = min(I1,I2);               %Lower energy threshold
  61. ITU = 5 * ITL;                  %Upper energy threshold
  62. IZC = mean(ZCR_Vector(1:sil));  %mean zero crossing rate for silence region
  63. stdev = std(ZCR_Vector(1:sil)); %standard deviation of crossing rate for
  64.                                 %silence region
  65. IZCT = min(ZTh,IZC+2*stdev);    %Zero crossing rate threshold
  66. indexi = zeros(1,lengthX);      %Four single-row vectors are created
  67. indexj = indexi;                %in these lines to facilitate computation below
  68. indexk = indexi;
  69. indexl = indexi;

  71. %Search forward for frame with energy greater than ITU
  72. for i = 1:length(Erg_Vector)
  73.     if (Erg_Vector(i) > ITU)
  74.         counter1 = counter1 + 1;
  75.         indexi(counter1) = i;
  76.     end
  77. end
  78. ITUs = indexi(1);

  80. %Search further forward for frame with energy greater than ITL
  81. for j = ITUs:-1:1
  82.     if (Erg_Vector(j) < ITL)
  83.         counter2 = counter2 + 1;
  84.         indexj(counter2) = j;
  85.     end
  86. end
  87. start = indexj(1)+1;

  89. Erg_Vectorf = fliplr(Erg_Vector);%Flips round the energy vector
  90. %Search forward for frame with energy greater than ITU
  91. %This is equivalent to searching backward from last sample for energy > ITU
  92. for k = 1:length(Erg_Vectorf)
  93.     if (Erg_Vectorf(k) > ITU)
  94.         counter3 = counter3 + 1;
  95.         indexk(counter3) = k;
  96.     end
  97. end
  98. ITUf = indexk(1);

  100. %Search further forward for frame with energy greater than ITL
  101. for l = ITUf:-1:1
  102.     if (Erg_Vectorf(l) < ITL)
  103.         counter4 = counter4 + 1;
  104.         indexl(counter4) = l;
  105.     end
  106. end

  108. finish = length(Erg_Vector)-indexl(1)+1;%Tentative finish index
  109.    
  110. %Search back from start index for crossing rates higher than IZCT
  111.    
  112. BackSearch = min(start,25);
  113. for m = start:-1:start-BackSearch+1
  114.     rate = ZCR_Vector(m);
  115.     if rate > IZCT
  116.         ZCRCountb = ZCRCountb + 1;
  117.         realstart = m;
  118.     end
  119. end
  120. if ZCRCountb > 3
  121.     start = realstart;          %If IZCT is exceeded in more than 3 frames
  122.                                 %set start to last index where IZCT is exceeded
  123. end

  125. %Search forward from finish index for crossing rates higher than IZCT
  126. FwdSearch = min(length(Erg_Vector)-finish,25);
  127. for n = finish+1:finish+FwdSearch
  128.     rate = ZCR_Vector(n);
  129.     if rate > IZCT
  130.         ZCRCountf = ZCRCountf + 1;
  131.         realfinish = n;
  132.     end
  133. end
  134. if ZCRCountf > 3
  135.     finish = realfinish;        %If IZCT is exceeded in more than 3 frames
  136.                                 %set finish to last index where IZCT is
  137.                                 %exceeded
  138. end

  140. x_start = FrmIndex(start);      %actual sample index for frame 'start'
  141. x_finish = FrmIndex(finish-1);  %actual sample index for frame 'finish'
  142. trimmedX = x(x_start:x_finish); %Trim speech sample by start and finish
  143.                                 %indices
復制代碼

作者: 792403245    時間: 2018-11-25 20:13
為啥我識別不出來....



歡迎光臨 (http://m.raoushi.com/bbs/) Powered by Discuz! X3.1