ARM, RISC-V контроллеры

А можно было просто воспользоваться формулой Эйлера: exp(j*w*t) = cos(w*t)+j*sin(w*t), и вспомнить немного школьной математики: exp[j*w*(t+∆t)] = exp(j*w*t)*exp(j*w*∆t) = [cos(w*t)+j*sin(w*t)]*[cos(w*∆t)+j*sin(w*∆t)] = cos(w*t)*cos(w*∆t)-sin(w*t)*sin(w*∆t)+j*[cos(w*t)*sin(w*∆t)+sin(w*t)*cos(w*∆t)], откуда немедленно находим: cos[w*(t+∆t)] = cos(w*t)*cos(w*∆t)-sin(w*t)*sin(w*∆t) и sin[w*(t+∆t)] = cos(w*t)*sin(w*∆t)+sin(w*t)*cos(w*∆t). Считая, что: cos[w*(t+∆t)] == cos_n+1 и sin[w*(t+∆t)] == sin_n+1, находим формулы рекурсии: cos_n+1 = cos_n*cos(w*∆t)-sin_n*sin(w*∆t) и sin_n+1 = cos_n*sin(w*∆t)+sin_n*cos(w*∆t) Эти формулы, конечно, несколько сложнее, так как требуют двух констант: sin(w*∆t) и cos(w*∆t), четырех умножений и двух сложений, но зато позволяют получить одновременно синус и косинус, причем с заранее известной одинаковой(!) аплитудой не зависящей от частоты синусоиды и, кроме того, не требуют, чтобы параметры sin(w*∆t) или cos(w*∆t) были много меньше единицы. При этом, ессно, должно выполняться соотношение: cos(w*∆t)²+sin(w*∆t)² = 1..