Цифровой генератор синусоиды, Нужен алгоритм Опции

[aal]

Требуется сформировать 2 синуса (естественно в цифре). Можно пойти табличным методом - но для 8 бит и 1к хватит, а для 18 бит и метра может мало оказатся. Есть другой метод http://www.gaw.ru/html.cgi/txt/doc/...v/max2000_5.htm .
Но интересные эфекты: даже при инт64 вылазеют гармоники на уровне -70 дб + к этому и частота не точно заданная.
Попробывал так:
X(n)=k*X(n-1)-X(n-2), где k=2*cos(2*Pi*F/Fd);
для 0 градусов X(-1)=0, X(-2)=-A*sin(2*Pi*F/Fd);
для 90 градусов X(-1)=A, X(-2)=A*sin(pi/2+2*Pi*F/Fd)
90 гр. частота на сотые доли процентов выше и уползает вперёд по отношению к 0 гр.
Сейчас хочу через приблежённое вычеслени по Тейлору попробывать. Но явно по вычислениям более накладно будет.
Может ещё есть способ с минимум вычислений сгенерить синус?
(Всё должно в итоге оказатся в FPGA Спартан3).

[Shamil]

Можно пойти табличным методом - но для 8 бит и 1к хватит, а для 18 бит и метра может мало оказатся.

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

Для примера: Нам требовалось иметь синусы 425 Гц, 2100 Гц, 2600 Гц на выходе
ИКМ кодека (8000 выборок/сек по 8 бит).

Сначала думали придеться ставить ПЗУ, но поигравшись с подбором
кратных частот нашли замены: 421,053 Гц (19 выборок), 2105,26 (19 выборок),
2600 (40 выборок).

Причем Synplify очень хорошо оптимизирует таблицы, и вышеприведенные синусы
заняли в ПЛИС Actel Ex256 не больше 10 %.