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

  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[][]; //端點(diǎn)檢測(cè)后的語(yǔ)音數(shù)據(jù)
  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為端點(diǎn)檢測(cè)后的每個(gè)語(yǔ)音數(shù)據(jù)的幀長(zhǎng)
  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. }
復(fù)制代碼

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






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