О дискретизации линейных дифференциальных уравнений

Рассмотрены некоторые вопросы получения дискретного описания дифференциальной системы (ДС) на равномерной сетке. Рассматриваются ДС в виде системы $n$ линейных обыкновенных дифференциальных уравнений первого порядка с постоянными коэффициентами или одно уравнение $n$-го порядка для наблюдаемого функционала состояния ДС. Изучаемые вопросы дискретизации важны для задач вариационной идентификации и аппроксимации динамических процессов моделями этого типа в конечном интервале. Дано сравнение аналитического равномерного (на основе теоремы Гамильтона–Кэли) и локальных методов дискретизации: на основе разделенных разностей и с помощью интерполяции выборок из $n+1$ отсчетов многочленами Тейлора степени $n$. Получена общая формула локальной дискретизации, прозволяющая сравнивать ее разностный и интерполяционные методы. Показано с использованием свойств обратных матриц Вандермонда, что в полученной общей формуле локальной дискретизации ее интерполяционному методу соответствуют $(n+1)$-матрицы Тейлора (из коэффициентов многочленов Тейлора), а разностному — $(n+1)$-матрицы Паскаля (из чисел треугольников Паскаля).
Показано, что невырожденность матрицы наблюдаемости ДС на сетке есть необходимое и достаточное условие как для аналитической дискретизируемости, так и для приведения дискретной системы (описания ДС сетке) к каноническому фробениусовскому виду. Он эквивалентен одному обыкновенному разностному уравнению для наблюдаемой переменной с постоянными коэффициентами. Это уравнение есть основа известного вариационного метода идентификации. Показано, что интерполяционный метод локальной дискретизации есть первое (линейное) приближение формулы равномерной аналитической дискретизации. Показано, что нулевое приближение ее не зависит от коэфффициентов ДС и есть вектор коэффициентов $n$-й разности. Показано также, что нулевое приближение матрицы наблюдаемости ДС и матрицы наблюдамости полиномиальной системы $y^{(n)}=0$ на сетке есть $n$-матрица Тейлора.