Псевдо генератор случаных чисел - C++ - Ответ 2515516

@sasha0123 · Регистрация: 03.12.2010

Author24 — интернет-сервис помощи студентам

Доброго времени суток!
Вот есть такая прорамка

C++

#include <math.h>
#include <stdio.h>
#include <conio.h>
#include <iostream>
using namespace std;
#define n 4
#define x01 0.21878104
#define x02 30121
#define ConstL 805
#define g 16
#define mer 1193
#define ts 1000000
 
double a1;
int teor;
int j;
int a2;
double a3;
unsigned int MT[624];  // стан генератора (тип даних int повинен бути 32-бітним, беззнаковим)
unsigned int index = 0;
const double PI = 3.14159265;
 
double SerKvad();
double Multik();
void InitialiseGenerator ( int seed );
void GenerateNumbers();
double ExtractNumber();
void Gistograma(double x[4][ts], int m, int gist[ts]);
double Vidxulennja(int gist[ts], int m);
double RivnomGist(double x[4][ts],int m);
float PopadVKolo (double x[4][ts]);
double Vidxulennja2(int one[33]);
double Stohast(double x[4][ts]);
int KilkistOdunu4ok(double x[4][ts]);
double Oper(int l);
double Resultat(double x[4][ts], int T);
double Korel(double x[4][ts]);
int Period(double x[4][ts]);
 
double x[4][ts];
 
 
int main()
{    int m;
    cout<<"--------------------------------------------------------------------"<<endl ;
    cout<<"Metod seredunnux kvadrativ:"<<endl;
    cout<<"--------------------------------------------------------------------" <<endl ;
    a1=x01;
    cout<<"Zgenerovano "<<ts<<" 4usel..."<<endl;
    j=0;
    for(int i=0;i<ts;i++)
    {
        x[j][i]=SerKvad(); cout<<"";
    }
    m=10;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 10 intervaliv\ndX = "<<RivnomGist(x,m)<<endl;
    m=100;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 100 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    m=1000;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 1000 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    cout<<"Metod otsinku kilkosti popadan v odunu4ne kolo!\ndX="<<PopadVKolo(x)<<endl;
    cout<<"Perevirka Stohastu4nosti rozpodily:\ndX="<<Stohast(x)<<endl;
    cout<<"p= "<<Korel(x)<<endl;
    double p=Period(x);
    if(p==0)cout<<"Dvzhuna periody:\nL>1000000"<<endl;
    else
    cout<<"Dvzhuna periody:\nL="<<p<<endl;
    cout<<"--------------------------------------------------------------------" <<endl   ;
    cout<<"--------------------------------------------------------------------" <<endl   ;
    cout<<"Myltuplikatuvnuj metod:"<<endl;
    cout<<"--------------------------------------------------------------------" <<endl   ;
    a2=x02;
    cout<<"Zgenerovano "<<ts<<" 4usel..."<<endl;
    j=1;
    for(int i=0;i<ts;i++)
    {
        x[j][i]=Multik(); cout<<"";
    }
    m=10;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 10 intervaliv\ndX = "<<RivnomGist(x,m)<<endl;
    m=100;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 100 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    m=1000;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 1000 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    cout<<"Metod otsinku kilkosti popadan v odunu4ne kolo!\ndX="<<PopadVKolo(x)<<endl;
    cout<<"Perevirka Stohastu4nosti rozpodily:\ndX="<<Stohast(x)<<endl;
    cout<<"p= "<<Korel(x)<<endl;
    p=Period(x);
    if(p==0)cout<<"Dvzhuna periody:\nL>1000000"<<endl;
    else
    cout<<"Dvzhuna periody:\nL="<<p<<endl;
    cout<<"--------------------------------------------------------------------"<<endl;
    int t;
    t=mer;
    cout<<"\n--------------------------------------------------------------------"<<endl   ;
    cout<<"Vuxr Mersona:"<<endl;
    cout<<"--------------------------------------------------------------------" <<endl  ;
    cout<<"Zgenerovano "<<ts<<" 4usel..."<<endl;
    j=2;
    InitialiseGenerator(t);
    for(int i=0;i<ts;i++)
    {
        x[j][i]=ExtractNumber();
    }
    m=10;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 10 intervaliv\ndX = "<<RivnomGist(x,m)<<endl;
    m=100;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 100 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    m=1000;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 1000 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    cout<<"Metod otsinku kilkosti popadan v odunu4ne kolo!\ndX="<<PopadVKolo(x)<<endl;
    cout<<"Perevirka Stohastu4nosti rozpodily:\ndX="<<Stohast(x)<<endl;
    cout<<"p= "<<Korel(x)<<endl;
    p=Period(x);
    if(p==0)cout<<"Dvzhuna periody:\nL>1000000"<<endl;
    else
    cout<<"Dvzhuna periody:\nL="<<p<<endl;
 
    cout<<"--------------------------------------------------------------------" <<endl   ;
    cout<<"--------------------------------------------------------------------" <<endl   ;
    cout<<"Zvu4ajnuj generator:"<<endl;
    cout<<"--------------------------------------------------------------------" <<endl   ;
    a2=x02;
    cout<<"Zgenerovano "<<ts<<" 4usel..."<<endl;
    j=3;
    for(int i=0;i<ts;i++)
    {
        x[j][i]=rand()/double(0x8000); cout<<"";
    }
    m=10;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 10 intervaliv\ndX = "<<RivnomGist(x,m)<<endl;
    m=100;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 100 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    m=1000;
    cout<<"Otsinka za dopomogojy gistogramu pru podili na 1000 intervaliv\ndX  = "<<RivnomGist(x,m)<<endl;
    cout<<"Metod otsinku kilkosti popadan v odunu4ne kolo!\ndX="<<PopadVKolo(x)<<endl;
    cout<<"Perevirka Stohastu4nosti rozpodily:\ndX="<<Stohast(x)<<endl;
    cout<<"p= "<<Korel(x)<<endl;
    p=Period(x);
    if(p==0)cout<<"Dvzhuna periody:\nL>1000000"<<endl;
    else
    cout<<"Dvzhuna periody:\nL="<<p<<endl;
system("pause");
return 0;
}
 
//metod seredunnux kvadrativ
double SerKvad()
{
    float x1,x2,y[2*n];
    int x[4*n];
    int i,j;
 
    x1=a1; //po4atkove zna4ennja
    x1=x1*x1;
    for(i=0;i<4*n;i++)
    {
        x1=x1*10;
        x[i]=x1;
        x1=x1-x[i];
    }
    for(i=n;i<3*n;i++)
    {
        y[i-n]=x[i]/pow(10,(float)i-3);
    }
    for(i=0;i<2*n;i++)
    x2=x2+y[i];
    a1=x2;
    return(x2);
}
 
//multuplikatuvnuj metod
double Multik()
{
    int x0,l;
    float x1;
    int i,j;
 
    x0=a2;
    l=ConstL;
 
    //generatsija 4usla
    x1=(x0*l)%65536;
    a2=x1;
    x1=x1/65536;
    return(x1);
}
//----------------------------------------------------------------------
//      Вихр Мерсона
//----------------------------------------------------------------------
// Початкова ініціалізація – заповнення масиву стану генератора
void InitialiseGenerator ( int seed ) {
  MT[0] = seed;
  for(int i=1;i<624;i++) {
    MT[i] = 1812433253 * (MT[i-1] ^ (MT[i-1]>>30)) + i;
  }
}
// Вибрати чергове ПВЧ  - ціле, 32 бітне, беззнакове, в інтервалі [0, 2^32)
double ExtractNumber() {
  if (index >= 624) {
    GenerateNumbers();
    index = 0;
  }
  unsigned int y = MT[index];
  y = y ^ (y>>11);
  y = y ^ ((y<<7) & (2636928640));
  y = y ^ ((y<<15) & (4022730752));
  y = y ^ (y>>18);
  index = index + 1;
  double j;
  j=y/4294967296.0;
  return (j);
}
// Змінити стан – збурення масиву MT – генерація наступних 624 чисел
void GenerateNumbers() {
  for (int i=0;i<624;i++) {
    int z=(MT[i])/2147483647;
    long int x=(MT[(i+1)%624])%2147483647;
   unsigned long int y = z + x;
    if (y%2==0) {
      MT[i] = MT[(i+397)%624] ^ (y>>1);
    } else { // if y odd
      MT[i] = MT[(i+397)%624]^(y>>1) ^ 2567483615;
    }
  }
}
 
//Metod gistogramu
void Gistograma(double x[3][ts], int m, int gist[1000])
{
    float xp,xk,per;
    xp=0.0;
    per=1.0/m;
    for(int k=0;k<m;k++)
    {
    gist[k]=0;
    xk=xp+per;
        for(int i=0;i<ts;i++)
        if(x[j][i]>=xp && x[j][i]<xk)
           gist[k]++;
    if(m==10) cout<<"Kilkist elementiv v intervali ["<<xp<<";"<<xk<<") -> "<<gist[k]<<endl;
    xp=xp+per;
    }
}
 
double Vidxulennja(int gist[1000],int m)
{
    double dx=0.0;
    for(int i=0;i<m;i++)
    {
        dx=dx+(teor-gist[i])*(teor-gist[i]);
    }
    dx=dx/m;
    return(dx);
}
 
double RivnomGist(double x[3][ts],int m)
{
   int gist[1000];
   teor=ts/m;
   Gistograma(x,m,gist);
   double dx;
   dx=Vidxulennja(gist,m);
   return(dx);
}
 
float PopadVKolo (double x[3][ts])
{
    unsigned int k=0;
    for(int i=0;i<ts-1;i=i+2)
    if(x[j][i]*x[j][i]+x[j][i+1]*x[j][i+1]<1)
    k++;
    float dx;
    double teoret;
    teoret=ts/2*(PI/4);
    dx=k/teoret;
    return(dx);
}
 
int KilkistOdunu4ok(double Stoh)
{
    int one = 0;
    unsigned ii=Stoh*4294967296;
    for(int t=0;t<32;t++)
       {    if(ii&1==1)
        {
            one++;
        }
            ii=ii>>1;
        }
       return(one);
  }
double Stohast(double x[3][ts])
{
int one[33];
int h;
    memset(one, 0, sizeof(one));
   for(int i=0;i<ts;i++)
     {
         h=KilkistOdunu4ok(x[j][i]);
         one[h]++;
     }
double dX=Vidxulennja2(one);
return (dX);
}
 
double Vidxulennja2(int one[33])
{
    double teoret[33];
    double dx=0.0;
    int y=32;
    double d2,d1,d3,d4;
    d1=Oper(y);
    d2=pow(0.5,32);
    for(int i=0;i<33;i++)
    {
        int k=y-i;
        d3=Oper(i);
        d4=Oper(k);
        teoret[i]=ts*d2*(d1/(d3 * d4));
        dx+=(teoret[i]-one[i])*(teoret[i]-one[i]);
    }
    dx=dx/33;
    return(dx);
}
 
double Oper(int l)
{
double res=1;
if(l<=1) return (1);
for(int i=1;i<=l;i++) res*=i;
return(res);
}
 
double Korel(double x[3][ts])
{
    double mas[100];
    double max = Resultat(x,1);
    int tulda = 1;
    for (int i = 2; i <= 100; i++)
    {
        if (max < Resultat(x,i)) {
            max = Resultat(x,i);
            tulda = i;
        }
    }
    cout<<"Dovir4a jmovirnist: B=";
    double B;
    if(max<0) B=max*(-1);
    else B=max;
    B*=sqrt(ts);
    cout<<B<<endl;
    cout<<"Perevirka koreljatsiji rozpodily";
    cout<<"\nT = "<<tulda<< endl;
    return(max);
}
 
int Period(double x[3][ts])
{
    double zn = x[j][ts-1];
    int i=ts-2;
    do{
       if(zn==x[j][i])
       {
          return (ts-1-i); }
          else i--;
    }
 
    while (i>0);
 
   return 0;
}
 
double Resultat(double x[3][ts], int T)
{
      double Dxi, Dxit, Res;
      double dija1, d1=0, d2=0, d3=0;
 
      dija1=1.0/(ts-T);
 
      for(int i =0;i<ts-T;i++)
      d1+=x[j][i];
 
      Dxi=dija1*d1-dija1*dija1*d1*d1;
 
      for(int i=0;i<(ts-T);i++)
      d2+=x[j][i+T];
 
      Dxit= dija1*d2-dija1*dija1*d2*d2;
 
 
      d1=0;
      d2=0;
 
      for(int i=0;i<ts-T;i++)
      {
          d1+=x[j][i]*x[j][i+T];
          d2+=x[j][i];
          d3+=x[j][i+T];
      }
 
      Res=dija1*d1 - dija1*d2*d3;
      d1=Dxi*Dxit;
      d1=sqrt(d1);
      Res=Res /d1;
      return(Res);
}

компилирую ее и выдает ошыбку :
1>c:\users\администратор\documents\visual studio 2010\projects\b2\b2\b2.cpp(351): error C2668: 'sqrt' : ambiguous call to overloaded function
1> d:\programfiles\vs 2010\vc\include\math.h(589): could be 'long double sqrt(long double)'
1> d:\programfiles\vs 2010\vc\include\math.h(541): or 'float sqrt(float)'
1> d:\programfiles\vs 2010\vc\include\math.h(127): or 'double sqrt(double)'
1> while trying to match the argument list '(int)'

Вернуться к обсуждению:
Псевдо генератор случаных чисел C++

Следующий ответ

Другие темы раздела
C++ Пожалуйсто помагите написать эти две задачи на с++ Пожалуйсто помагите написать эти две задачи на с++ ссесию мне бе них не закрыть а никто не понимает((((( ребят если кто знает как сделать помагите буду очень блогодарна!!!!!)))=* Дан двух мерный массив размерностью N x M ( N и M вводиться с клавиатуры). Найти в строках наименьшие элементы. Произвести вывод данного массива в виде матрицы. Дан двухмерный массив размерностью N x M (N и M... https://www.cyberforum.ru/ cpp-beginners/ thread450814.html	C++ Найти количество всех положительных элементов двумерного массива Дан двумерный массив размерностью N x M (N и М вводятся с клавиатуры). Найти количество всех положительных элементов данного массива. Произвести вывод данного массива в виде матрицы. помоги пожалуста!
C++ mingw32, SEH, std::exception, setjmp/longjmp, signal https://www.cyberforum.ru/ cpp-beginners/ thread450793.html Начал разбираться с вопросом - как понял - есть несколько вариантов: (в порядке обнаружения мною) 1) LibSEH - работает. http://www.programmingunlimited.net/siteexec/content.cgi?page=mingw-seh 2) <eh.h>, _set_se_translator - долго мучался с заголовками, а когда скомпилил - не смог слинковать. 3) <excpt.h> - видимо, недоделанная версия LibSEH - работает неадекватно (или я чего-то не...	C++ Работа с рядками https://www.cyberforum.ru/ cpp-beginners/ thread450789.html Дан текст. Групы символов, розделеные пропусками которые не содержат пропусков всередині себя, будем называть словами словами.Посчитать , сколько слов в начинаются и заканчиваются одной и той же буквой.
C++ Работа с рядками Дан текст, который составляется из букв, цифр и пропусков. Последовательность которая начинается символом и в которой нет пропусков назовём идентификатором. Вывести все идентификаторы которыесодержатся в тексте и подсчитать сколько раз каждый из них встречался.	C++ Заполнить матрицу по образцу Заполнить матрицу размером nm числами 1,2,3, … nm, оюходя её соответственно варианту: https://www.cyberforum.ru/ cpp-beginners/ thread450782.html
C++ В одномерном массиве найти наименьшую по длине убывающую последовательность В одномерном массиве найти наименьшую по длине убывающую последовательность https://www.cyberforum.ru/ cpp-beginners/ thread450780.html	C++ как можно реализовать это в коде? Здравствуйте. Если вам не трудно скажите как можно в коде реализовать это: Вариантов множество. Например, заводится массив с указателями на текущие элементы лент; затем с использованием этого просматриваются текущие элементы лент и среди них выявляется минимальный; после "вытаскивания" минимального элемента сдвигается соответствующий указатель на один элемент и процесс повторяется.
C++ Из vector в массив. Доброго всем утра коллеги! Собственно, есть вектор с шаблоном std::vector<float> vect; и теперь надо каждый из элементов вектора передать в массив for(unsigned int i=0; i<vect.size(); i++){ float arr={0, 0, 0}; arr = vect.at(i)//пишет: error C2106: '=' : left operand must be l-value https://www.cyberforum.ru/ cpp-beginners/ thread450772.html	C++ Преобразовать массив по правилу Дан одномерный массив целых чисел, состоящий из N элементов. Преобразовать массив таким образом, чтобы сначала располагались все элементы, отличающиеся от максимального не более чем на 20%, а потом - все остальные. помоги пожалуста! https://www.cyberforum.ru/ cpp-beginners/ thread450770.html
Лишние символы при чтении из файла C++ Здравствуйте! Есть текстовый файл следующего содержания: DROP FUNCTION IF EXISTS NumericOnly; CREATE FUNCTION NumericOnly (val VARCHAR(255)) RETURNS VARCHAR(255) BEGIN DECLARE idx INT DEFAULT 0; IF ISNULL(val) THEN RETURN NULL; END IF;	C++ Командная строка https://www.cyberforum.ru/ cpp-beginners/ thread450763.html Здравствуйте. Возможна тема, подобная этой уже создавалась, в любом случае не нашел и, если что - прошу прощения. Суть проблемы такова: написал элементарную программу: #include <iostream> using namespace std;

	22.02.2012, 12:05