Фрактальное Броуновское движения ( оптимизация кода)

@krupa · Регистрация: 13.12.2010

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

В общем у меня функция генерации фрактального броуновского движения.
Делается все очень просто. Расскажу алгоритм
длина ряда составляет 2^N+1. и формируется он таким образом:
Формируется три точки нормальным законом распределения с 0 мат. ожиданием и каким-то СКО.
Например:
0.5 0.2 0.8
После этого проходит процесс интерполяции - мы усредняем значения между точками ( я не знаю как это правильно назвать)
т.е. добавляем новые элементы между элементами, вставляя их среднее значение.
x_i+1/2=(x_i+x₊₁)/2
Т.е. получаем
0.5 0.35 0.2 0.5 0.8
После этого пересчитываем СКО по формуле

Python
1
sigma = pow(0.5,H)*sigma

где H - показатель Хёрста ( вообще принимает значение от 0 до 1 )
и после этого к каждому элементу нашего ряда ( на данном этапе состоящем из 5 узлов, добавляем случайную величину сгенерированную

Python
1
random.normal(0,sigma)

и процесс повторяется: интерполяция, пересчёт СКО, добавление случайной величины.

вот то что я написал:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
def interpolation(y):
    i=0
    new_y=[]
    for i in range(len(y)-1):
        new_y.append(y[i])
        new_y.append((y[i+1]+y[i])/2)
    new_y.append(y[i])
    return new_y
 
def offset(y,sigma):
    y = [number+random.normal(0,sigma) for number in y]
    return y
H=0.7
m=0
sigma=1
y=random.normal(m,sigma,size=3)
i=0
 
while i<10:
    y=interpolation(y)
    sigma = pow(0.5,H)*sigma
    y=offset(y,sigma)
    i+=1

можно ли это как-то оптимизировать? я не прошу сделать, я прошу натолкнуть на какие-то функции, какие-то принципы работы Python.

@dondublon · 09.04.2013, 09:13

Питон, как скриптовый язык, штука вообще медленная.
Операция append медленная - память надо выделять сразу, когда знаешь, сколько элементов будет.
Можно использовать нативный код, попробуй numpy и его массивы.

Новые блоги и статьи Все статьи Все блоги /
Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .	Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .
Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .	модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	09.04.2013, 09:13
	Питон, как скриптовый язык, штука вообще медленная. Операция append медленная - память надо выделять сразу, когда знаешь, сколько элементов будет. Можно использовать нативный код, попробуй numpy и его массивы. 0