c语言实现pcm解码

static int pcm[8000][8] = {0};   //储存8000个已编码值
void Pcm_jiema()
{
	int n,k,c,f,g,m;
	float e,h,l,v,o,p;
	int i;
	int b[8]={0};                       //储存每个已编码值
	int z[8]={0,16,32,64,128,256,512,1024};               //每一段起始值
    int y[8]={1,1,2,4,8,16,32,64};                                  //对应每一段内最小间隔
	FILE *fp;
	fp = fopen("jiema.txt","w");                               //用来存解码值
	for(i=0;i<8000;i++)
	{
		for(k=0;k<8;k++)
		{
			b[k]=pcm[i][k];	
		}
	    if (b[0]==1)                                               //判断符号
		{
			c=1;
		}
	    else 
		{
			c=-1;
		}
		f=b[1]*4+b[2]*2+b[3];                               //计算位于第几段
	    g=z[f]+(b[4]*8+b[5]*4+b[6]*2+b[7]*1)*y[f];                       //计算有多少量化间隔
		e=y[f]/2.00;                                                                //对应半个最小间隔
	    h=g+e;                                                                //相加，减少误差
		l=h*c;                                                                   //恢复符号
		v=l/2048.00;                                                         //恢复原始振幅
		if (v>=0)                                                              //记录符号
	      {
		     m=1;
	      }
	    else
	      {
		     m=-1;
	      }
		v=fabs(v);
		if(v>=0&&v<0.25)                                          //反十三折线计算
		{
          o=v/16;
		}
        else if(v>=0.25&&v<0.375)
		{
           o=(v-0.125)/8;
		}
         else if(v>=0.375&&v<0.5)
		 {
            o=(v-0.25)/4;
		 }
         else if(v>=0.5&&v<0.625)
		{
            o=(v-0.375)/2;
	    }
         else if(v>=0.625&&v<0.75)
		{
            o=v-0.5;  
		}
         else if(v>=0.75&&v<0.875) 
		{
            o=(v-0.625)*2;
		}
         else if(v>=0.875&&v<=1)
		{
            o=(v-0.75)*4;
	    }
		p=o*m;                                                                              //恢复符号，解码结束
		fprintf(fp,"%f\n",p);
	}
	fclose(fp);
}

已编码值：

解码结果：

和编码前比较：

作图：