Автокореляционная функция

[DaRKSPiRT]

спасибо вам огромное! общий принцип понял, как только напишу конечный вариант сразу же выкину сюда, ждите....Mind working...

[DaRKSPiRT]

#include <iostream>
#include <fstream>
#include <vector>
#include <iterator>
#include <iostream.h>
#include <math.h>

int main()
{
 std::vector<double> buffer; //создаем буфер для входного файла;
 std::vector<double> akf;    //-||- для вычисленных значений АКФ;
 std::ifstream in("in.txt");  //Обозначаем входной файл;

   if(in) {
   typedef std::istream_iterator<int> ii;
     std::copy(ii(in), ii(), std::back_inserter(buffer));
    in.close(); //загружаем все целые числа в вектор и закрываем файл;
   buffer.insert(buffer.end(),buffer.begin(),buffer.end());
    int n=buffer.size();//Определяем размер буфера;
    int norm=0;
    for (int i=0; i<n; i++)
   {

   if(norm<buffer[i]) norm=buffer[i];//определяем нормирующее число;
   int k=floor(i/500); //Определение номера сегмента;
   int n1=k*500;       //Задание начальной точки егмента;
   int n2=(k+1)*500;   //Конечной точки сегмента;
   int l=i-n1;         //Задержка на сегменте;
   if (n2>n) n2=n-1;   //Предотвращение ошибки граничных условий;
   double sum=0;       //Создаем переменную суммы;



       for(int j=0; j<n2-n1; j++)     //Цикл вычисления АКФ на сегментах
       {                              //по 500 отсчетов с задержкой i;
       double p1=buffer[n1+j];        //первый множитель;
       double p2=buffer[n2-l+j];      //второй множитель;
       sum+=p1*p2;                    //автосумма произведения множителей;
       }                              //конец вычисления суммы;
       akf.push_back(sum);            //заполняем массив значениями АКФ;


   }

   norm=norm*norm;

    std::vector<double> rezult; //Создание вектора нормированных значений;
   for(int i=0;i<500;i++)       //на выходе будет 500 чисел;
   {                            //начало цикла нормировки;
     int h=floor(n/500);        //колличество сегментов по 500 значений;
     double rez=0;
     for(int j=0;j<h;j++) rez+=akf[j*500+i];//сумма соот-щих значений каждого сегмента;
     rezult.push_back(rez/norm); //заносим вычисления в вектор результатов;
   }



    // теперь сохраним результаты в другой файл
     std::ofstream out("out.txt");
    if(out) {
       std::copy(rezult.begin(), rezult.end(), std::ostream_iterator<int>(out, " ")); 
     out.close();
     } else std::cout << "error: out.txt" << std::endl; //Ошибка если выходной файл неверен;

   } else std::cout << "error: in.txt" << std::endl;//Ошибка если входной файл неверен;

  return 0;
 }

Это более обширная версия того что я хотел.
еще бы это все визуализировать. например как здесь:

std::ifstream in("in.txt");

задать переменную имени файла? вместо "in.txt"

[filosof_x86]

#include <iostream>
#include <fstream>
#include <string>

int main()    
{  
  std::string file_name;
  std::getline(std::cin, file_name);
  std::ifstream in(file_name.c_str());
  // ...
  return 0;  
}

[DaRKSPiRT]

if(защита.курсовой=5) beer

[filosof_x86]

Сообщение от DaRKSPiRT

if(защита.курсовой=5)

Выражение этого ифа всегда будет истиной, ибо защита.куросовой присваивается 5, а потом преобразуется в булевое, а все, что отлично от нуля – истина

[DaRKSPiRT]

пиво навсегда))