У меня стоит задача сформировать звуковой сигнал (меандр) из 8-ми гармоник, каждая из которой может иметь или не иметь сдвиг по фазе (на пи).
#data - массив с байтами, в формате ['01010101', '11110000']
ADS = len(data) * SPP #Вычисляю размер массива под сигнал
signal = []
for i in range(ADS):
signal.append(0)
# Формируем по одному периоду сигнала на каждый байт
for b in range(len(data)):
d = 0
for N in range(1, 17, 2):
for i in range(SPP):
index = i + (b * SPP)
signal[index] += harm(N, FREQ, i, SPS, int(data[b][d]))
#s.next()
d += 1
# Функция для вычисления отдельной гармоники
def harm(N, freq, cs, sps, ph):
arg = freq*N*2*np.pi*(cs/sps)
if ph:
arg += np.pi
return AMPLITUDE*(4/np.pi)*np.sin(arg)/N
Хотел бы узнать, могу ли я как-то оптимизировать данный участок кода. Я думал использовать map, но не придумал, как заменить i (номер текущего отсчёта). Заранее спасибо!
SPS
?harm
посмотреть. Ну и так навскидку можно попробовать декорироватьharm
декораторомfunctools.lru_cache
. Но тут будет ещё завесить от того, насколько велико значениеSPP
. Если не поможет, тогда копать в сторонуNumba
, либо разбираться в вычислениях и как-то их векторизовать черезNumpy
, но это придётся понять, что вообще делает этот алгоритм.