Автокорреляционная функция и частота основного тона (ЧОТ)

int main(int argc, char* argv[])
{
    FILE *f,*fton;
    int size;
    short *buf;
 
 
 
 
    if (!argv[1])   {   cout<<"Enter the name of file!"<<endl;      system("pause");        return -1;  }
 
    f=fopen(argv[argc-1],"rb"); //исходный файл с голосом
    if (f==NULL){   cout<<"Open error!"<<endl;      system("pause");        return -1;  }
 
    fseek(f,0,2);
    size=ftell(f)/sizeof(short);//Кол-во отсчетов
 
    fseek(f,0,0);
    buf=(short*)malloc(size*sizeof(short));
    fread(buf,sizeof(short),size,f); //считываем отсчеты
    int bs=(int)(2*fd * 0.013f);
    short in[208];//bs
    int t;
    fton=fopen("Freq.pcm","w");
 
    for(int i=0;i<size/bs;i++)
    {
        for(int j=0;j<=bs;j++)
        {in[j]=buf[i*bs+j];}
 
        t=(int)TON(in);
        fprintf(fton,"%d \n ",t);
    }
    
 
    fclose(f);
    fclose(fton);
    return 0;
}

#include <stdio.h>
#include <iostream>
 
using namespace std;
 
int fd = 8000; //Частота дискретизации
 
double TON (const short *in) //размер массива равен двум максимальным периодам ОТ
{
    int minOT   = fd * 0.0030f; // минимальный шаг - 3 мс (330 Гц)
    int maxOT   = fd * 0.013f;  // максимальный шаг - 13 мс (77 Гц)
    int deltaOT = fd * 0.0001f; // шаг - 0,1 мс
    if (deltaOT==0) deltaOT=1;
 
        int OT=0;//ЧОТ
    int maxAC=-10000L;
 
    for (int pitch=minOT;pitch<=maxOT;pitch+=deltaOT)
    {
        double max=-10000L,val=0;
 
                for(int j=0;j<maxOT;j++)
        {
            double p=in[j];
            double q=in[j+pitch];
 
            val=p*q;       //Считаем значение автокорреляции
            if(val>max) max=val; //находим максимум при данном значении ОТ
        }
 
        if (max>maxAC)//находим максимум АКФ
        {
            maxAC=max;
            OT=pitch;
        }
    }
    return (double)fd/(OT); //возвращаем ЧОТ
}
 
int main(int argc, char* argv[])
{
    FILE *f,*fton;
    int size;
    short *buf;
 
//    if (!argv[1])   {   cout<<"Enter the name of file!"<<endl;      system("pause");        return -1;  }
 
    f=fopen("ton200.pcm", "rb"); //исходный файл с голосом
    if (f==NULL)
    {
        cout<<"Open error!"<<endl;
        system("pause");
        return -1;
    }
 
    fseek(f,0,2);
    size=ftell(f)/sizeof(short);//Кол-во отсчетов
 
    fseek(f,0,0);
    buf=(short*)malloc(size*sizeof(short));
    fread(buf,sizeof(short),size,f); //считываем отсчеты
    int bs=(int)(2*fd * 0.013f);
    short in[208];//bs
    int t;
    fton=fopen("Freq.pcm","w");
 
    for(int i=0;i<size/bs;i++)
    {
        for(int j=0;j<bs;j++)
        {in[j]=buf[i*bs+j];}
 
        t=(int)TON(in);
        fprintf(fton,"%d \n ",t);
    }
 
 
    fclose(f);
    fclose(fton);
    return 0;
}