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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 嵌入式/单片机编程 > DSP编程 > 查看源码
vector_quantization.c
资源名称：speaker_recognition.rar　[点击查看]
上传用户：bossps2lzz
上传日期：2022-07-07
资源大小：522k
文件大小：2k
源码类别：

DSP编程

开发平台：

C/C++

vector_quantization.c：源码内容
  1. #include <math.h>
  2. struct buffer {
  3. short data[100][256];
  4. };
  6. #define dim 8
  7. #define nn 8
  8. float codebk[8][128];
  10. void vetcor_quantization(struct buffer *data) {
  12. nn=no_tr_vec;
  13.         mtemp=1;
  14.         for (k=0;k<dim;++k) {
  15.           codebk[0][k]=0.0;
  16.           for (i=0;i<nn;++i) {
  17.             codebk[0][k]+=data[i][k];
  18.           }
  19.           codebk[0][k]/=(float) nn;
  20.         }
  22.         while (mtemp < cb_size) {
  24.           mtemp=mtemp*2;
  25.           for(i=0;i<mtemp;++i) {  /* splitting */
  26.             j=i/2;
  27.             del = 0.001;
  28.             if(i == i/2*2) del=-0.001;
  29.             for (k=0;k<dim;++k)
  30.               yy[i][k]=codebk[j][k]*(1.0+del);
  31.           }
  33.           totd2=10.0e15;
  35.           flag1 = 0;
  36.           while ( flag1 == 0) {
  37.             totd1=0.0;
  38.             for(i=0;i<mtemp;++i) {
  39.               for(k=0;k<dim;++k)
  40.                 codebk[i][k]=0.0;
  41.               bin[i]=0;
  42.             }
  43.             for(i=0;i<nn;++i){
  44.               j=0;
  45.               dist2=0.0;
  46.               for(k=0;k<dim;++k) {
  47.                 dist2+=(data[i][k]-yy[j][k])*(data[i][k]-yy[j][k]);
  48.               }
  49.               index=0;
  50.               for(j=1;j<mtemp;++j){
  51.                 dist1=0.0;
  52.                 for(k=0;k<dim;++k)
  53.                   dist1+=(data[i][k]-yy[j][k])*(data[i][k]-yy[j][k]);
  54.                 if(dist1 < dist2){
  55.                   dist2=dist1;
  56.                   index=j;
  57.                 }
  58.               }
  59.               ++bin[index];
  60.               for(k=0;k<dim;++k)
  61.                 codebk[index][k]+=data[i][k];
  62.               totd1=totd1+dist2;
  63.             }  /* end for i<nn */
  65.             for(j=0;j<mtemp;++j){
  66.               if(bin[j] > 0) {
  67.                 for(k=0;k<dim;++k){
  68.                   codebk[j][k]/=(float) bin[j];
  69.                   yy[j][k]=codebk[j][k];
  70.                 }
  71.               }
  72.             }
  74.     totd1/=(float) nn*dim;
  75.             drel=(totd2-totd1)/totd1;
  76.     k=fflush(res_file);
  77.             flag1 = 1;
  78.             if(drel > 0.001) {
  79.               flag1 = 0;
  80.               totd2=totd1;
  81.             }
  82.           }   /* end while flag1 */
  84.         }   /* end while mtemp<cb_size */
  86. dist=totd1;