久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標題: lpc語音識別源碼 [打印本頁]
作者: titan7149    時間: 2017-4-29 14:43
標題: lpc語音識別源碼
lpc語音識別源碼，這是我根據趙力那本“語音信號處理”后面的源代碼寫的LPC美爾倒譜特征的C源代碼，基本上沒做什么修改，需要的朋友就可以不用自己再錄入了。
程序調試通過，可以直接用。
由于這段代碼是我一個程序的一部分，我轉過來可能會少一些定義什么的，所以并不是非常完整。

  1. //Calculate LPC coefficients coefficients based on autocorrelation
  2. void corref(int ip,float cor[],float alf[])

  4. {
  5. int i,j;
  6. float refi,alfsave,resid;
  7. alf[1]=(-cor[1]);
  8. alf[0]=1;
  9. resid=1-cor[1]*cor[1];
  10. for(i=2;i<=ip;i++)
  11. {
  12. refi=cor[i];
  13. for(j=1;j<i;j++)
  14. refi+=alf[j]*cor[i-j];
  15. refi/=resid;
  16. alf[i]=(-refi);
  17. for(j=1;2*j<=i;j++)
  18. {
  19. alfsave=alf[j];
  20. alf[j]=alfsave-refi*alf[i-j];
  21. if(2*j!=i)
  22. alf[i-j]-=refi*alfsave;
  23. }
  24. }
  25. }

  27. //Calculate LPC cepstrum based on alf
  28. void alfcep(int ip,float alf[],float cep[],int n)

  30. {
  31. int m,i;
  32. float ss;
  33. m=0;
  34. goto start;
  35. next:
  36. m+=1;
  37. ss=0;
  38. i=0;
  39. goto step3;
  40. step2:
  41. ss-=alf[i]*cep[m-i];
  42. step3:
  43. i+=1;
  44. if(i>ip) goto step4;
  45. if(i<m) goto step2;
  46. ss-=alf[m]*m;
  47. step4:
  48. cep[m]=ss;
  49. start:
  50. if(m<n) goto next;
  51. for(m=2;m<=n;m++)
  52. cep[m]/=m;
  53. }

  55. //Calculate mel-warped LPC cepstrum
  56. void cepmel(float cep[],float mel[])
  57. {
  58. int k,n;
  59. float u0[ORDER+2],u1[ORDER+2];
  60. float alpha=ALPHA;
  61. //Calculation mel[0]*
  62. u0[ORDER+1]=0;
  63. for(n=ORDER;n>=0;--n)
  64. u0[n]=cep[n]+alpha*u0[n+1];
  65. //mel[0]=u0[0];
  66. //Calculation mel[1]*
  67. u1[ORDER+1]=0;
  68. for(n=ORDER;n>=0;--n)
  69. u1[n]=((1-alpha*alpha)*u0[n+1]+alpha*u1[n+1]);
  70. mel[1]=u1[0];
  71. //Calculation mel[2],mel[3],......,mel[ORDER] *
  72. for(k=2;k<=ORDER;++k)
  73. {
  74. for(n=0;n<=ORDER+1;++n)
  75. u0[n]=u1[n];
  76. for(n=ORDER;n>=0;--n)
  77. u1[n]=u0[n+1]+alpha*(u1[n+1]-u0[n]);
  78. mel[k]=u1[0];
  79. }
  80. }

  82. main()
  83. {
  84. int rad[][]; //端點檢測后的語音數據
  85. float cor[ORDER+1]; //autocorrelation coefficients ORDER＝14
  86. float alf[ORDER+1]; //LPC
  87. float cep[ORDER+1]; //LPC cepstrum
  88. float mel[ORDER+1];
  89. float mel_data[length][ORDER];
  90. //mel-warped LPC cepstrum，length為端點檢測后的每個語音數據的幀長
  91. int mel_v=0;

  93. //Auto-Correlation Calculation
  94. correl(rad[count],cor,W_SIZE,ORDER);
  95. //LPC Analysis Program
  96. corref(ORDER,cor,alf);
  97. //LPC Cepstrum Coefficient Calculation
  98. alfcep(ORDER,alf,cep,ORDER);
  99. alf[0]=cep[0]=cor[0]; //cor[0] is power
  100. //MEL Cepstrum Coefficient Calculation
  101. cepmel(cep,mel);
  102. //printf("\n Frame %d :",count);

  104. for(i=0;i<ORDER;i++)
  105. {
  106. mel_data[mel_v][i]=mel[i+1];
  107. //printf("\nFrame %d: Mel[%d][%d]= %g",count,mel_v,i,mel_data[mel_v][i]);
  108. }
  109. }
復制代碼

下載：
00063058lpc.rar (1.54 KB, 下載次數: 13)






歡迎光臨 (http://www.zg4o1577.cn/bbs/) Powered by Discuz! X3.1
主站蜘蛛池模板： 伊人99| 精品一区二区三区日本 | 一级毛片在线播放 | av电影一区 | 亚洲精品久久视频 | 亚洲第一av网站 | 久久狠狠 | 欧美日韩1区2区 | 成人一区在线观看 | 久草视频网站 | 久久久久久蜜桃一区二区 | 6080yy精品一区二区三区 | 欧美日韩成人网 | 国产成都精品91一区二区三 | 一级黄色毛片 | 在线观看中文字幕亚洲 | 欧产日产国产精品国产 | 狠狠操狠狠 | 久久国产高清 | 国产xxxx搡xxxxx搡麻豆 | 国产精品久久久久久久久久久久 | 91 视频网站 | 国产三级国产精品 | 欧美在线视频一区 | www久| 色视频在线免费观看 | 一级毛片视频在线 | 黄色国产| 国产情品 | 亚洲精品视频一区二区三区 | 国产乱人伦精品一区二区 | 久久久久国产精品 | 国产精品久久久久久久久久久免费看 | 日本一区二区三区四区 | 国产亚洲一区精品 | 福利视频一二区 | 7799精品视频天天看 | 337p日本欧洲亚洲大胆鲁鲁 | 欧美456| 日韩有码在线观看 | 国产免费一区二区 |