Преобразование дифференциального уравнения второго порядка к системе двух уравнений первого порядка. Написание функции для вычисления правой части системы и функции, определяющей граничные условия. Формирование начального приближения и вызов солвера.
Аннотация к работе
Решение краевой задачи на примере обыкновенных дифференциальных уравнений второго порядка с использованием MATLAB функция граничный солвер дифференциальныйРешение этой задачи в MATLAB состоит из следующих этапов: 1) преобразование дифференциального уравнения второго порядка к системе двух уравнений первого порядка; Функция должна иметь два входных аргумента: переменную x, по которой производится дифференцирование, и вектор, размер которого равен числу неизвестных функций системы. Функция, описывающая граничные условия, зависит от двух аргументов - векторов ya и yb. 4) Для задания начальной сетки и приближения предназначена функция bvpinit, обращение к которой в самом простом случае выглядит следующим образом: bvpinit (начальная сетка, начальное приближение к решению). Формирование равномерной сетки целесообразно производить функцией linspace: linspace (a, b, n) возвращающей вектор х из n равноотстоящих узлов между a и b, включая границы.