Задача оптимального измерения с учетом резонансов: алгоритм программы и вычислительный эксперимент

В статье описана программа, реализующая алгоритм численного метода решения задачи оптимального измерения с учетом резонансов — задачи восстановления динамически искаженного сигнала с учетом инерционности измерительного устройства и его механических резонансов, решаемой с использованием методов теории оптимального управления. Основной идеей алгоритма численного решения является представление компонент измерения тригонометрическими полиномами, которое позволяет свести задачу оптимального управления к задаче выпуклого программирования относительно неизвестных коэффициентов многочленов. Использование стандартных методов, например, градиентных, при решении задачи выпуклого программирования, в силу сложности функционала качества, приводит к неудовлетворительным результатам. Поэтому предлагается иной, более простой метод, который вместе с тем более трудоемок. В статье представлен ряд решений, позволяющих повысить скорость вычислений, блок-схема основной процедуры программы, написанной на языке С++. Для конкретной модели датчика приводятся результаты вычислительного эксперимента.