Форум программистов, компьютерный форум CyberForum.ru

Сгенерировать последовательность 60 случайных чисел с экспоненциальным законом распределения с параметром - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ определить головной класс http://www.cyberforum.ru/cpp-beginners/thread239044.html
Как сослаться в классе B на класс A? Пишет IntelliSense: нестатическая ссылка не член должна указываться относительно заданного объекта Как обойти? class A { int x; class B {
C++ Найти сумму елементов главной и вспомог. диагоналей. Дана задача: Найти сумму элементов главной и вспомогательной диагоналей и вывести на экран сообщение какая сумма больше. и формула есть i-n-j-1 вот что покачто готово: #include<iostream.h> void main() { http://www.cyberforum.ru/cpp-beginners/thread238989.html
C++ введённые 2 строки делятся на 3, в 3-й строке веводятся оперделённые части 1-й и 2-й строк
Помогите склеить такую штуку: 1. введённые 2 строки (S1 и S2) делятся на 3 части strlen (S1)/3; strlen (S2)/3; 2. в третей строке нужно чтобы выводило 1-ю часть S1, 2-ю часть S2 и 3-ю часть S1. пример. S1= asdtyhnhu S2= mfjlencjy S3=asdlennhu
Нарисовать любую геометрическую фигуру C++
Здравствуйте! Помогите пожалуйста, как использовать в Turbo C графику с помощью библиотеки #include <graphics.h>? Хотя-бы вывести что нибудь (точку, линию окружность...)
C++ последовательность 1 3 5 7 ... 8 6 4 2 0 http://www.cyberforum.ru/cpp-beginners/thread238952.html
нужно было написать алгоритм 1 3 5 7 ... 8 6 4 2 0. Причём половину ячеек массива занимают нечётные числа в порядке возрастания, половину - чётные в порядке убывания. Написал код, компилируется, да не так, как нужно.Подскажите, пожалуйста!Возможно наличие дурацких ошибок:) #include "TXLib.h" int main() { txCreateWindow (500, 500); const int SIZE = 22; int x = {1}; int...
C++ что такое trunc в pascal и на какой оператор его заменить в С++?? ... REPEAT m1:=int(w1/h); m:=trunc(m1); for j:=1 to m do begin ... подробнее

Показать сообщение отдельно
alex_x_x
бжни
 Аватар для alex_x_x
2441 / 1646 / 84
Регистрация: 14.05.2009
Сообщений: 7,163
06.02.2011, 20:43     Сгенерировать последовательность 60 случайных чисел с экспоненциальным законом распределения с параметром
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
#include <cmath>
#include <cstdlib>
#include <vector>
#include <algorithm>
#include <iostream>
#include <iterator>
 
struct s_exp{
  s_exp( const double param ) : m_dParam(param) {}
  double operator()(){
    double dA = 1.0*rand()/RAND_MAX;
    return -1.0/m_dParam * log( 1 - dA );
  }
private:
  double m_dParam;
};
 
int main(){
  std::vector<double> v( 60 );
  std::generate( v.begin(), v.end(), s_exp(0.8) );
  std::sort( v.begin(), v.end() );
  std::copy( v.begin(), v.end(), std::ostream_iterator<double>( std::cout, " " ) );
}
Добавлено через 36 минут
ну чтото вроде

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
#include <cmath>
#include <cstdlib>
#include <vector>
#include <algorithm>
#include <iostream>
#include <iterator>
#include <functional>
#include <numeric>
#include <ctime>
 
const int N = 60, N_STEP = 10;
const double PARAM = 0.8;
 
struct s_exp{
  s_exp( const double param ) : m_dParam(param) {}
  double operator()(){
    double dA = 1.0*rand()/RAND_MAX;
    return -1.0/m_dParam * log( 1 - dA );
  }
private:
  double m_dParam;
};
 
struct s_transf{
  s_transf( const double dFirst ) : m_dPrev( dFirst )  {}
  double operator()( const double dEl ){
    double dRes = dEl - m_dPrev;
    m_dPrev = dEl;
    return dRes;
  } 
private:
  double m_dPrev;
};
 
struct s_disp{
  s_disp( const double dAvg ) : m_dAvg( dAvg ) {}
  double operator()( const double x, const double y ){
    return x + (m_dAvg-y)*(m_dAvg-y);
  }
private:
  double m_dAvg;
};
 
int main(){
  srand(time(0));
  std::vector<double> v( N );
  std::vector<double> v_res( N-1 );
  std::generate( v.begin(), v.end(), s_exp(PARAM) );
  std::sort( v.begin(), v.end() );
  std::copy( v.begin(), v.end(), std::ostream_iterator<double>( std::cout, " " ) ); 
  std::cout << std::endl << "<<>>" << std::endl;
  std::transform( ++v.begin(), v.end(), v_res.begin(), s_transf( *v.begin() ) );
  std::sort( v_res.begin(), v_res.end() );
  std::copy( v_res.begin(), v_res.end(), std::ostream_iterator<double>( std::cout, " " ) ); 
  std::cout << std::endl; 
  double dAvg = 1.0 * std::accumulate( v_res.begin(), v_res.end(), 0 ) / v_res.size();
  double dDisp = 1.0 * std::accumulate( v_res.begin(), v_res.end(), 0, s_disp(dAvg) ) / v_res.size();
  std::cout << "avg = " << dAvg << std::endl << "disp = " << dDisp << std::endl;
  
  double dStep = ( v_res.back() - v_res.front() ) / (N_STEP+1);
  std::cout << "step = " << dStep << std::endl;
  std::vector<double>::iterator it = v_res.begin();
  for( size_t i=0;i<=N_STEP;++i ){
    for( ;*it<=dStep*(i+1);++it ){
      std::cout << *it << " ";
    }
    std::cout << std::endl;
  }
  std::cout << v_res.back();
}
 
Текущее время: 05:55. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru