Перевести код с python на c++. Экспоненциальное распределение

@xamelione25 · Регистрация: 04.02.2013

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

Помогите перевести код с python на c++ дополнительно добавив вывод значений в excel или txt из строки 25. Мне необязательно чтобы было построение графика и его отображение. Мне главное нужно чтобы был осуществлен вывод значений координат графика в excel или txt из строки 25.

print(*list(zip(bins, lambd * np.exp(-lambd * bins))), sep='\n')

Python
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
import numpy as np
if __name__ == '__main__':
    import matplotlib.pyplot as plt
    lambd1 = 0.1
    N = 10000
    def expa(lambd, N):
        import random
        import math
        xi = []
        yi = []
        sum = 0
        sig = 0
        for i in range(N):
            xi.append(random.random())
            yi.append(-1/lambd * math.log(xi[i]))
            sum += yi[i]
        my = sum / N
        for i in range(N):
            sig = (yi[i] - my) ** 2
        sig = sig / (N - 1)
        return yi
    arr = expa(lambd1, N)
    def test_exponential(lambd, arr):
        count, bins, ignored = plt.hist(arr, 30, density=True)
        print(*list(zip(bins, lambd * np.exp(-lambd * bins))), sep='\n')
        plt.plot(bins, lambd * np.exp(-lambd*bins), linewidth=1, color='y')
        i = 0
        plt.grid()
        plt.show()
    s = test_exponential(lambd1, arr)

@Catstail · 22.02.2021, 17:55

Ха, а я сделал на плюсах код, который и гистограмму строит:

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
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <ctime>
 
using namespace std;
 
int main()
{
    srand(time(0));
 
    int n=6000;
 
    double *y =new double[n],*x = new double[n];    
    double ymax,ymin,dy,hmin,hmax,lambda=0.5;
    int i,j,k;
    double h[30];
    int    z[30];
    
    for (i=0; i<n; i++)
    {
        x[i]=((double)rand() / (double)RAND_MAX);
        y[i]=(-1.0/lambda)*log(x[i]);
    }
 
    ymax=ymin=y[0];
    
    for (i=1; i<n; i++)
    {
        if (y[i]>ymax) ymax=y[i];
        if (y[i]<ymin) ymin=y[i];
    }
 
    dy=(ymax-ymin)/30;
    
    for (i=0; i<30; i++) h[i]=0;
    
    for (i=0; i<n; i++)
    {
        k=(int)((y[i]-ymin)/dy);
        h[k]++;
    }
 
    //for (i=0; i<30; i++) 
    //    cout << ymin+i*dy << " " << h[i] << endl;
 
    hmax=hmin=h[0];
    for (i=1; i<30; i++)
    {
        if (h[i]>hmax) hmax=h[i];
        if (h[i]<hmin) hmin=h[i];
    }
    
    for (i=0; i<30; i++) z[i]=40*(hmax-h[i])/(hmax-hmin);
    
    for (i=0; i<40; i++)
    {
        for(j=0; j<30; j++)
           if (z[j]<i)
              cout << " * ";
           else
              cout << "   ";
              
        cout << endl;
    }
    delete [] y;
    delete [] x;
    return 0;
}

@xamelione25 · 22.02.2021, 18:08 **[ТС]**

Catstail, почему так сильно значения игриков взлетели...ониже мизерные должны получаться???

Добавлено через 5 минут
Catstail, интересно то что звёздочки без значений запускаются...а с координатами ошибку выбивает

@Catstail · 22.02.2021, 18:14

xamelione25, ты знаешь, как строят гистограмму? Делят интервал от xmin до xmax на 30 (например) отрезков, а потом считают, сколько значений попало в каждый из них. Это большие целые числа. Что такое твои "игрики" - я не понял.

@u235 · 22.02.2021, 18:23

А нельзя ли при переводе кода задействовать GNU Scientific Library?
Там и экспоненциальное распределение и гистограмма есть...

@Catstail · 22.02.2021, 18:27

xamelione25, да, нашли у меня небрежность. Чтобы этой ошибки не было, нужно в стр. 34 делить на 29:

C++
1
dy=(ymax-ymin)/29;

@xamelione25 · 22.02.2021, 18:36 **[ТС]**

Catstail, ну если это то же самое что и в питоне....
то почему там мизерные значения получаются а здесь на с++большие

Python
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
(0.0021717001162772694, 0.09997828535680724)
(3.309962896938033, 0.07182078338066031)
(6.617754093759789, 0.051593452588257684)
(9.925545290581544, 0.03706286989197027)
(13.2333364874033, 0.02662462494207552)
(16.541127684225057, 0.019126167384565324)
(19.848918881046814, 0.01373954673984202)
(23.156710077868567, 0.009869993335342433)
(26.464501274690324, 0.007090246154715884)
(29.77229247151208, 0.00509337634043288)
(33.080083668333835, 0.0036588973047186123)
(36.38787486515559, 0.0026284194592499587)
(39.69566606197735, 0.0018881614536856062)
(43.0034572587991, 0.0013563868821004285)
(46.311248455620856, 0.0009743792673782962)
(49.619039652442616, 0.0006999588164893277)
(52.92683084926437, 0.0005028250920192491)
(56.23462204608612, 0.00036121135587985207)
(59.542413242907884, 0.00025948117086321983)
(62.85020443972964, 0.00018640188614375483)
(66.15799563655139, 0.00013390437172131008)
(69.46578683337314, 9.619205651304792e-05)
(72.7735780301949, 6.910089355011588e-05)
(76.08136922701667, 4.963958212887101e-05)
(79.38916042383842, 3.5659280037266635e-05)
(82.69695162066017, 2.5616336766796273e-05)
(86.00474281748193, 1.8401849635330386e-05)
(89.31253401430368, 1.321922307174836e-05)
(92.62032521112543, 9.496211635440128e-06)
(95.9281164079472, 6.821734903452347e-06)
(99.23590760476895, 4.900487571202193e-06)

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
0.00030519 2736
3.23497 2016
6.46964 1396
9.7043 1117
12.939 759
16.1736 540
19.4083 371
22.643 273
25.8776 210
29.1123 163
32.347 123
35.5816 92
38.8163 52
42.051 33
45.2856 36
48.5203 16
51.755 12
54.9896 12
58.2243 11
61.459 8
64.6936 8
67.9283 6
71.163 3
74.3976 0
77.6323 1
80.867 2
84.1016 1
87.3363 2
90.571 0
93.8056 0
 
 *
 *
 *
 *
 *
 *
 *
 *
 *
 *
 *  *
 *  *
 *  *
 *  *
 *  *
 *  *
 *  *
 *  *
 *  *
 *  *  *
 *  *  *
 *  *  *
 *  *  *
 *  *  *  *
 *  *  *  *
 *  *  *  *
 *  *  *  *
 *  *  *  *
 *  *  *  *  *
 *  *  *  *  *
 *  *  *  *  *
 *  *  *  *  *
 *  *  *  *  *  *
 *  *  *  *  *  *
 *  *  *  *  *  *  *
 *  *  *  *  *  *  *
 *  *  *  *  *  *  *  *  *
 *  *  *  *  *  *  *  *  *  *
 *  *  *  *  *  *  *  *  *  *  *  *

@Catstail · 22.02.2021, 19:32

xamelione25, я не знаю, что в Питоне (это результат работы Питоновских программ). А в C++ я строю гистограмму. Как она строится - я написал выше: Делят интервал от xmin до xmax на 30 (например) отрезков, а потом считают, сколько значений попало в каждый из них. Это большие целые числа.

@xamelione25 · 23.02.2021, 00:52 **[ТС]**

Catstail, подскажите а почему по данной функции график вообще не такой какой должен быть...
хотя я функцию прописал вроде правильно...
....
не понимаю....

там вообще парабола должна быть как на синем графике

а получается вообще не понятно что (см. то что в консоле)

Кликните здесь для просмотра всего текста

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
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <ctime>
#define _USE_MATH_DEFINES
using namespace std;
int main()
{
    srand(time(0));
    double pi = acos(-1);
    int n = 10000;
    int m2 = 10;
    int m3 = 33;
    int sigma2 = 1.5;
    int sigma3 = 3;
    double* y = new double[n], * x = new double[n];
    double* y1 = new double[n], * x1 = new double[n];
    double* y2 = new double[n], * x2 = new double[n];
    double a, b, c, ymax, ymin, y1max, y1min, y2max, y2min, dy, dy1, dy2, hmin, hmax, h1min, h1max, h2min, h2max, lambda = 0.1;
    int i, j, k, k1, k2;
    int h[30]; int h1[30]; int h2[30];
    int z[30]; int z1[30]; int z2[30];
    for (i = 0; i < n; i++)
    {
        while (1)
        {
            a = ((double)rand() / (double)RAND_MAX);
            if (a > 0) break;
        }
        x[i] = a;
        y[i] = (-1.0 / lambda) * log(x[i]);
    }
    for (i = 0; i < n; i++)
    {
        while (1)
        {
            b = ((double)rand() / (double)RAND_MAX);
            if (b > 0) break;
        }
        x1[i] = b;
        y1[i] = sigma2 * cos(2 * pi * x1[i]) * sqrt(-2*log(x1[i])) + m2;
    }
    for (i = 0; i < n; i++)
    {
        while (1)
        {
            c = ((double)rand() / (double)RAND_MAX);
            if (c > 0) break;
        }
        x2[i] = c;
        y2[i] = sigma3 * cos(2 * pi * x2[i]) * sqrt(-2 * log(x2[i])) + m3;
    }
 
    ymax = ymin = y[0];
    y1max = y1min = y1[0];
    y2max = y2min = y1[0];
    for (i = 0; i < n; i++)
    {
        if (y[i] > ymax) ymax = y[i];
        if (y[i] < ymin) ymin = y[i];
 
        if (y1[i] > y1max) y1max = y1[i];
        if (y1[i] < y1min) y1min = y1[i];
 
        if (y2[i] > y2max) y2max = y2[i];
        if (y2[i] < y2min) y2min = y2[i];
    }
 
    dy = (ymax - ymin) / 29;
    dy1 = (y1max - y1min) / 29;
    dy2 = (y2max - y2min) / 29;
    for (i = 0; i < 30; i++) h[i] = 0;
    for (i = 0; i < 30; i++) h1[i] = 0;
    for (i = 0; i < 30; i++) h2[i] = 0;
    for (i = 0; i < n; i++)
    {
        k = (int)((y[i] - ymin) / dy);
        h[k]++;
 
        k1 = (int)((y1[i] - y1min) / dy1);
        h1[k1]++;
 
        k2 = (int)((y2[i] - y2min) / dy2);
        h2[k2]++;
    }
 
    for (i = 0; i < 30; i++)
        cout << "Xi = " << ymin + i * dy << "; Yi = " << h[i] << endl;
    hmax = hmin = h[0];
    for (i = 1; i < 30; i++)
    {
        if (h[i] > hmax) hmax = h[i];
        if (h[i] < hmin) hmin = h[i];
    }
    for (i = 0; i < 30; i++) z[i] = 40 * (hmax - h[i]) / (hmax - hmin);
    for (i = 0; i < 40; i++)
    {
        for (j = 0; j < 30; j++)
            if (z[j] < i)
                cout << " * ";
            else
                cout << "   ";
        cout << endl;
    }
    delete[] y;
    delete[] x;
 
    for (i = 0; i < 30; i++)
        cout << "Xi = " << y1min + i * dy1 << "; Yi = " << h1[i] << endl;
    h1max = h1min = h1[0];
    for (i = 1; i < 30; i++)
    {
        if (h1[i] > h1max) h1max = h1[i];
        if (h1[i] < h1min) h1min = h1[i];
    }
    for (i = 0; i < 30; i++) z1[i] = 40 * (h1max - h1[i]) / (h1max - h1min);
    for (i = 0; i < 40; i++)
    {
        for (j = 0; j < 30; j++)
            if (z1[j] < i)
                cout << " * ";
            else
                cout << "   ";
        cout << endl;
    }
    delete[] y1;
    delete[] x1;
 
    for (i = 0; i < 30; i++)
        cout << "Xi = " << y2min + i * dy2 << "; Yi = " << h2[i] << endl;
    h2max = h2min = h2[0];
    for (i = 1; i < 30; i++)
    {
        if (h2[i] > h2max) h2max = h2[i];
        if (h2[i] < h2min) h2min = h2[i];
    }
    for (i = 0; i < 30; i++) z2[i] = 40 * (h2max - h2[i]) / (h2max - h2min);
    for (i = 0; i < 40; i++)
    {
        for (j = 0; j < 30; j++)
            if (z2[j] < i)
                cout << " * ";
            else
                cout << "   ";
        cout << endl;
    }
    delete[] y2;
    delete[] x2;
    return 0;
}

@xamelione25 · 23.02.2021, 01:51 **[ТС]**

Catstail, на питоне оно вот такой вид имеет

Python
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
import numpy as np
 
if __name__ == '__main__':
    import matplotlib.pyplot as plt
    m2, sigma2 = 10, 1.5
    m3, sigma3 = 33, 3
    lambd1 = 0.1
    N = 10000
    def norm(m, sigma, N):
        import random
        import math
        xi = []
        yi = []
        for i in range(N):
            xi.append([random.random(), random.random()])
            yi.append(sigma * math.cos(2 * math.pi * xi[i][0]) * math.sqrt(-2 * math.log(xi[i][1])) + m)
        return yi
    arr2 = norm(m2, sigma2, N)
    arr3 = norm(m3, sigma3, N)
 
    def test_normalvariate(mu, sigma, arr):
        count, bins, ignored = plt.hist(arr, 30, density=True)
        i = 0
        plt.plot(bins, 1/(sigma * np.sqrt(2 * np.pi)) * np.exp( - (bins - mu)**2 / (2 * sigma**2) ), linewidth=1, color='y')
        plt.grid()
        plt.show()
 
    s = test_normalvariate(m2, sigma2, arr2)
    s = test_normalvariate(m3, sigma3, arr3)

Добавлено через 33 минуты
Catstail, смутное сомнение.... но мне даже кажется что парабола почему-то перевернутая верх-тормашками)))

@Catstail · 23.02.2021, 07:09

xamelione25, второй график я не строил.

@u235 · 23.02.2021, 07:20

xamelione25, вот ваш "причесанный" код.

Python
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
import numpy as np
import matplotlib.pyplot as plt
 
m2, sigma2 = 10, 1.5
m3, sigma3 = 33, 3
lambd1 = 0.1
N = 10000
 
def norm(m, sigma, N):
    xi = np.random.random(size=(N, 2))
    yi = sigma * np.cos(2 *np.pi * xi[:,0]) * np.sqrt(-2 * np.log(xi[:,1])) + m
    return yi
        
def test_normalvariate(mu, sigma, arr):
    count, bins, ignored = plt.hist(arr, 30, density=True)
    i = 0
    plt.plot(bins, 1/(sigma * np.sqrt(2 * np.pi)) * np.exp( - (bins - mu)**2 / (2 * sigma**2) ), linewidth=1, color='y')
    plt.grid()
    plt.show()
 
if __name__ == '__main__':
    arr2 = norm(m2, sigma2, N)
    arr3 = norm(m3, sigma3, N)
    s = test_normalvariate(m2, sigma2, arr2)
    s = test_normalvariate(m3, sigma3, arr3)

Если нужен вывод как на С++
11я строка должна быть такой:

Python
1
yi = sigma * np.cos(2 *np.pi * xi[:,0]) * np.sqrt(-2 * np.log(xi[:,0])) + m

@xamelione25 · 23.02.2021, 12:07 **[ТС]**

u235, немного не понял ... Вы говорите ...если нужен на с++ ...а показываете на питоне)))

Добавлено через 23 минуты
u235, я это уже сделал ...см в 48 посте в коде строки 41 и 51

Добавлено через 1 час 34 минуты
Catstail,

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

второй график я не строил.

да я это понимаю ...
Просто подскажите ... там для построения второго графика достаточно только функцию поменять и все, или нужно еще дополнительно какие-то изменения вносить???

@u235 · 23.02.2021, 12:23

xamelione25, у вас на C++ выдает бимодальную гистограмму. А на Python - одномодальную. Приведенный код на C++ не эквивалентен коду на Python. Я показал, где это исправить для идентичности гистограмм. Проблема в $https://www.cyberforum.ru/cgi-bin/latex.cgi?{{x}_{i}}^{*}$
для C++ (это тот же x_i) и в Python (это два случайных x) у вас это совершенно разные величины, из-за этого и разные результаты.

@xamelione25 · 23.02.2021, 12:37 **[ТС]**

u235, но то что вы написали - это же на питоне или на с++?

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

yi = sigma * np.cos(2 *np.pi * xi[:,0]) * np.sqrt(-2 * np.log(xi[:,0])) + m

Добавлено через 47 секунд
u235, насколько я понимаю - это на питоне конструкция, а как мне тогда на с++ её прописать?

Добавлено через 11 минут

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

у вас на C++ выдает бимодальную гистограмму. А на Python - одномодальную. Приведенный код на C++ не эквивалентен коду на Python. Я показал, где это исправить для идентичности гистограмм. Проблема в
для C++ (это тот же xi) и в Python (это два случайных x) у вас это совершенно разные величины, из-за этого и разные результаты.

да ...я икс со звездочкой видел..... не придал особистого значения....

но как это все правильно на с++ написать ...ни бум бум

@u235 · 23.02.2021, 12:37

xamelione25, например задублируйте 35-39 строки, введите новые переменные типа b2, x1_ подумайте.. Это ваше задание, вам его и делать.
Разберитесь что такое икс со звездочкой в задании. У меня задания нет.

@Catstail · 23.02.2021, 14:32

xamelione25, я сегодня попозже построю и второй.

@xamelione25 · 23.02.2021, 20:06 **[ТС]**

Catstail, расскажите, а там разве изменения самой функции недостаточно - то есть там для этого какие-то иные еще манипуляции вычислительные нужно делать?

@Catstail · 24.02.2021, 15:22

xamelione25, у тебя тема - "экспоненциальное распределение". Я тебе в плюсах построил гистограмму плотности вероятности. Она такая, как и должна быть (смотри здесь). То, что ты ошибочно называешь "параболой" - это "нормальное распределение" смотри здесь. Тебе и это нужно?

@xamelione25 · 24.02.2021, 15:34 **[ТС]**

Catstail,

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

Тебе и это нужно?

Да ... мне нужно выполнить показательное и нормальное распределение ...
вот показательное уже есть, нужно еще для нормального построить ...

@xamelione25 · 24.02.2021, 15:36 **[ТС]**

Catstail, я сделал дубль кода, забил туда функции, а график какой-то неправильный строиться ... см.48пост ... из-за "x*" - не знаю как правильно его оформить

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Программа принимает математическое выражение в виде строки и выдаёт его производную в виде строки и вычисляет значение производной при заданном х Логарифм записывается как: (x-2)log(x^2+2) -. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .

@xamelione25 -4 / 5 / 2 Регистрация: 04.02.2013 Сообщений: 1,860
	22.02.2021, 18:08 [ТС]
	Catstail, почему так сильно значения игриков взлетели...ониже мизерные должны получаться??? Добавлено через 5 минут Catstail, интересно то что звёздочки без значений запускаются...а с координатами ошибку выбивает 0

@Catstail Супер-модератор 38173 / 21108 / 4307 Регистрация: 12.02.2012 Сообщений: 34,709 Записей в блоге: 14
	22.02.2021, 18:14
	xamelione25, ты знаешь, как строят гистограмму? Делят интервал от xmin до xmax на 30 (например) отрезков, а потом считают, сколько значений попало в каждый из них. Это большие целые числа. Что такое твои "игрики" - я не понял. 0

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,759
	22.02.2021, 18:23
	А нельзя ли при переводе кода задействовать GNU Scientific Library? Там и экспоненциальное распределение и гистограмма есть... 0

@Catstail Супер-модератор 38173 / 21108 / 4307 Регистрация: 12.02.2012 Сообщений: 34,709 Записей в блоге: 14
	22.02.2021, 19:32
	xamelione25, я не знаю, что в Питоне (это результат работы Питоновских программ). А в C++ я строю гистограмму. Как она строится - я написал выше: Делят интервал от xmin до xmax на 30 (например) отрезков, а потом считают, сколько значений попало в каждый из них. Это большие целые числа. 1

@Catstail Супер-модератор 38173 / 21108 / 4307 Регистрация: 12.02.2012 Сообщений: 34,709 Записей в блоге: 14
	23.02.2021, 07:09
	xamelione25, второй график я не строил. 0

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,759
	23.02.2021, 12:23
	xamelione25, у вас на C++ выдает бимодальную гистограмму. А на Python - одномодальную. Приведенный код на C++ не эквивалентен коду на Python. Я показал, где это исправить для идентичности гистограмм. Проблема в $https://www.cyberforum.ru/cgi-bin/latex.cgi?{{x}_{i}}^{*}$ для C++ (это тот же x_i) и в Python (это два случайных x) у вас это совершенно разные величины, из-за этого и разные результаты. 0

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,759
	23.02.2021, 12:37
	xamelione25, например задублируйте 35-39 строки, введите новые переменные типа b2, x1_ подумайте.. Это ваше задание, вам его и делать. Разберитесь что такое икс со звездочкой в задании. У меня задания нет. 0

@Catstail Супер-модератор 38173 / 21108 / 4307 Регистрация: 12.02.2012 Сообщений: 34,709 Записей в блоге: 14
	23.02.2021, 14:32
	xamelione25, я сегодня попозже построю и второй. 1

@xamelione25 -4 / 5 / 2 Регистрация: 04.02.2013 Сообщений: 1,860
	23.02.2021, 20:06 [ТС]
	Catstail, расскажите, а там разве изменения самой функции недостаточно - то есть там для этого какие-то иные еще манипуляции вычислительные нужно делать? 0

@Catstail Супер-модератор 38173 / 21108 / 4307 Регистрация: 12.02.2012 Сообщений: 34,709 Записей в блоге: 14
	24.02.2021, 15:22
	xamelione25, у тебя тема - "экспоненциальное распределение". Я тебе в плюсах построил гистограмму плотности вероятности. Она такая, как и должна быть (смотри здесь). То, что ты ошибочно называешь "параболой" - это "нормальное распределение" смотри здесь. Тебе и это нужно? 0

@xamelione25 -4 / 5 / 2 Регистрация: 04.02.2013 Сообщений: 1,860
	24.02.2021, 15:36 [ТС]
	Catstail, я сделал дубль кода, забил туда функции, а график какой-то неправильный строиться ... см.48пост ... из-за "x*" - не знаю как правильно его оформить 0

Перевести код с python на c++. Экспоненциальное распределение

Решение