[/quote] А нельзя ли как нибудь по-простому, численно, методом Рунге-Кутта порядка эдак шестого? А то у меня есть система неслабых дифференциальных уравнений и хотелось процесс их решения немного ускорить. Я так понимаю, что основная трудность в том, чтобы реализовать плавующюю точку и тригонометрические функции? [/quote] В тех самых цифровых диф.анализаторах все было с фиксированной запятой, идея была в том, что каждое уравнение реализовалось на накапливающем сумматоре - процессорном элементе. И между этими сумматорами пересылались сигналы в виде приращений - переполнений сумматоров. Ну прям как сигналы между процессами в VHDL. И синусы-косинусы получались сами по себе - как результаты решения уравнений. Назывались они уравнения Шеннона.
А насчет интегрирования уравнений. Не лезет красиво в ПЛИС плавающая запятая. Лучше фиксированная, но длинная. Я эффектно использовал арифметику рациональных дробей- все лучше, чем с плавающей запятой. Хорошая книжка Д.М.Смит Математическое и цифровое моделирование для инженеров и исследователей. Машиностроение-1980. Опять-таки, по книжке, берется аналоговый решатель диф.уравнений и заменяется на цифровой эквивалент. Все выливается в схему типа рекурсивного фильтра, основанную на теореме Котельникова.
|