Операции над большими числами, Операции над большими числами Опции

[NickSmith]

Для программирования DDS синтезатора, мне надо вычислять код частоты..
Исходные:
Есть константа в виде 5050. Далее мне нужно эту константу разделить на 1000 затем полученный результат разделить на 1 умноженную на 2^28
У полученного результата то что будет после запятой округлить в большую сторону..
Как все это выполнить на асме?? Я нашел только стандартные процедуры деления, которые описаны в апнотах. Но у них нет деления для нецелых числе. Точней только для 1 байтных чисел есть..

Подскажите, как бы все это реализовать??

[NickSmith]

Имеется формула код фазы = ( a*b )/c где а =значение фазы в градусах, b = 2^12, c = 360 например: (90*4096)/360=1024. Значение фазы у меня хранится в виде 90*10^2.
Далее я 2^12/360 = 11,38. 11,38 я возвожу в 10^2. В итоге у меня получается 900*1138=1024200 и что бы мне привести это к нужному виду, я делю это на 1000..
А можно ли как то это свести к одному действию?? желательно делению...

Сообщение отредактировал NickSmith - Mar 14 2010, 10:13

[ae_]

Имеется формула код фазы = ( a*b )/c где а =значение фазы в градусах, b = 2^12, c = 360 например: (90*4096)/360=1024. Значение фазы у меня хранится в виде 90*10^2.
...

360°/4096=0.088, т.е. хранить в 12-bit регистре число градусов 0...360 не выйдет точнее, чем 0.1° (округляем 0.088 в большую сторону до целого десятичного разряда).
Отсюда и формат хранения переменной должен быть не 90*10^2, а 90*10=900, т.е. 0...3599
Далее, надо получить из 900 код 1024, а из 3599->4095. N = 0x91AD
т.е. берёте фазу в виде 0...3599, умножаете на 0x91AD (используя алгоритм беззнакового mul16x16 из appnote), результат (32-bit) сдвигаете влево на один разряд. старшие два байта = код фазы
(900*N)<<1 = 0x2002434<<1 = 0x4004868. отбрасываем младшие 16 бит, получаем 0x400 (1024)
(3599*N)<<1 = 0x7FFFF23<<1 = 0xFFFFE46. -//- получаем 0xFFF (4095)