Реализуется алгоритм численного решения граничных задач для обыкновенных дифференциальных уравнений, основанный на методе коллокации и представлении решения в виде разложения по полиномам Чебышева. Предлагается вместо традиционного подхода - слияния всех условий (дифференциальных и граничных) в одну систему линейных алгебраических уравнений (СЛАУ) - перейти к методике решения задачи в несколько отдельных этапов. Вначале выделяются спектральные коэффициенты, определяющие «общее» решение исходной задачи. Трудоемкость приведения матрицы СЛАУ к диагональной форме (в случае систем ОДУ с постоянными коэффициентами) на этом этапе эквивалентна сложности умножения чебышевской матрицы коэффициентов на вектор правой части системы. На втором этапе учет граничных условий выделяет «частное» искомое решение, однозначно доопределяя недостающие коэффициенты искомого разложения. Предложенный метод может использоваться для моделирования задач классической механики. Вклад авторов: все авторы сделали эквивалентный вклад в подготовку публикации. Авторы заявляют об отсутствии конфликта интересов.
Источник: VESTNIK TOMSKOGO GOSUDARSTVENNOGO UNIVERSITETA-UPRAVLENIE VYCHISLITELNAJA TEHNIKA I INFORMATIKA
Tomsk State University Journal of Control and Computer Science