Перевести код с 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 · 24.02.2021, 16:39

Вот плотность вероятности нормального распределения:

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
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <ctime>
 
using namespace std;
 
int main()
{
    srand(time(0));
 
    int n=6000;
 
    double *y =new double[n];    
    double x1,x2,ymax,ymin,dy,hmin,hmax,sigma=3.0;
    int i,j,k;
    double h[30];
    int    z[30];
    
    // Нормальное распределение 
    // Гистограмма плотности вероятности
    
    cout << endl << "Нормальное распределение. Плотность вероятности" << endl;
    
    for (i=0; i<n; i++)
    {
        x1=((double)rand() / (double)RAND_MAX);
        x2=((double)rand() / (double)RAND_MAX);
        y[i]=sigma*cos(2*3.1415926535*x1)*sqrt(-2.0*log(x2));  
    }
 
    ymax=ymin=y[0];
    
    for (i=1; i<n; i++)
    {
        if (y[i]>ymax) ymax=y[i];
        if (y[i]<ymin) ymin=y[i];
    }
 
//cout << ymin << " "  << ymax << endl;
 
    dy=(ymax-ymin)/29;
    
    for (i=0; i<30; i++) h[i]=0;
    
    for (i=0; i<n; i++)
    {
        k=(int)((ymax-y[i])/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;
 
    return 0;
}

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

Code
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
                                                 **  **                                                                 
                                                 **  **                                                                 
                                                 **  **                                                                 
                                                 **  **                                                                 
                                                 **  **  **                                                             
                                                 **  **  **  **                                                         
                                             **  **  **  **  **                                                         
                                             **  **  **  **  **                                                         
                                             **  **  **  **  **                                                         
                                             **  **  **  **  **                                                         
                                             **  **  **  **  **                                                         
                                             **  **  **  **  **                                                         
                                             **  **  **  **  **  **                                                     
                                             **  **  **  **  **  **                                                     
                                         **  **  **  **  **  **  **                                                     
                                         **  **  **  **  **  **  **                                                     
                                         **  **  **  **  **  **  **                                                     
                                         **  **  **  **  **  **  **                                                     
                                         **  **  **  **  **  **  **                                                     
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **                                                 
                                     **  **  **  **  **  **  **  **  **  **                                             
                                 **  **  **  **  **  **  **  **  **  **  **                                             
                                 **  **  **  **  **  **  **  **  **  **  **                                             
                                 **  **  **  **  **  **  **  **  **  **  **                                             
                                 **  **  **  **  **  **  **  **  **  **  **                                             
                             **  **  **  **  **  **  **  **  **  **  **  **  **                                         
                             **  **  **  **  **  **  **  **  **  **  **  **  **                                         
                             **  **  **  **  **  **  **  **  **  **  **  **  **  **                                     
                             **  **  **  **  **  **  **  **  **  **  **  **  **  **                                     
                             **  **  **  **  **  **  **  **  **  **  **  **  **  **                                     
                         **  **  **  **  **  **  **  **  **  **  **  **  **  **  **                                     
                     **  **  **  **  **  **  **  **  **  **  **  **  **  **  **  **  **                                 
                     **  **  **  **  **  **  **  **  **  **  **  **  **  **  **  **  **  **

@xamelione25 · 24.02.2021, 19:27 **[ТС]**

Catstail, спасибо)) только подскажите а почему оно только "через раз на 10ый" запускается без этой ошибки на 49 строке.
Вызвано исключение по адресу 0x00007FF68DC83257 в MPS2.exe: 0xC0000005: нарушение прав доступа при чтении по адресу 0x0000007408CFF6A0.

@Catstail · 25.02.2021, 07:03

xamelione25, посмотрю

Новые блоги и статьи Все статьи Все блоги /
Установка нескольких штампов электронной подписи в строго определенных местах файла docx ВладимирСамохин 19.07.2026 (В!) Работа с Электронной подписью - это неотъемлемая часть современного документооборота. Но что делать, если нужно поставить несколько штампов электронной подписи в строго определенных местах. . .	сукцессия 35. Научная статья о проделанной работе anaschu 19.07.2026 Написал в формате латекс и пдф	Вангую, что это не пройдёт модерацию, и на неделе я запущу свой сервер. Hrethgir 19.07.2026 Эта публикация сейчас в песочнице и ждёт приглашения. https:/ / habr. com/ ru/ sandbox/ 295048/ начало и оглавление - Как «пернатого» заставить осваивать новые горизонты опыта через масштабирование. . .	сукцессия 33. открытые вопросы от клауде anaschu 19.07.2026 "Что накопилось за эту часть А — тринадцать правок, из которых шесть пришли из ваших вопросов и каждая оказалась реальной ошибкой, а не калибровкой: односторонний симбиоз, отсутствующий листопад,. . .
32 сукцессия anaschu 19.07.2026 сукцессия 28‑мерное ядро стабилизировано Коллеги, фиксирую разбор инженерных правок и их изоморфную проекцию на экономику, меметику и половой отбор. Модель теперь не «подкручивает» сходимость —. . .	сукцессия 31: модель микоризы - это модель ещё нескольких явлений, социальных и экономических anaschu 18.07.2026 Теория «Всего»: апдейт v1. 1. 2 — 28‑мерное ядро стабилизировано Коллеги, фиксирую разбор инженерных правок и их изоморфную проекцию на экономику, меметику и половой отбор. Модель теперь не. . .	сукцессия 30. Массив проверяющих друг друга моделей anaschu 18.07.2026 Архитектура сети взаимопроверяющих моделей микоризной сукцессии (v2. 0) Развитие тензорного ОДУ-ядра и создание кросс-платформенного калибровочного полигона Уважаемые коллеги! В продолжение. . .	Грибы - это женщины, деревья - это мужчины. Анти инь янь для союза мужчины и женщины. anaschu 18.07.2026 ГЛАВНЫЙ НАУЧНО-ФИЛОСОФСКИЙ ВЫВОД: Сексуально-Репродуктивный Капитализм против Государства Моногамии Коллеги, мы вышли на финишную прямую 20-мерного ОДУ-моделирования вековой сукцессии (ветка. . .

@xamelione25 -4 / 5 / 2 Регистрация: 04.02.2013 Сообщений: 1,861
	24.02.2021, 19:27 [ТС]
	Catstail, спасибо)) только подскажите а почему оно только "через раз на 10ый" запускается без этой ошибки на 49 строке. Вызвано исключение по адресу 0x00007FF68DC83257 в MPS2.exe: 0xC0000005: нарушение прав доступа при чтении по адресу 0x0000007408CFF6A0. Миниатюры 0

@Catstail Супер-модератор 38210 / 21143 / 4313 Регистрация: 12.02.2012 Сообщений: 34,756 Записей в блоге: 14
	25.02.2021, 07:03
	xamelione25, посмотрю 0

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

Решение