CC BY
14
№ 9 (78)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
сентябрь, 2020 г,
К РЕШЕНИЮ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ВИДА у' = f(x,y) МЕТОДОМ ХОРД
Аблялимов Олег Сергеевич
канд. техн. наук, старший научный сотрудник, и.о. профессора кафедры «Локомотивы и локомотивное хозяйство», Ташкентский государственный транспортный университет,
Узбекистан, г. Ташкент E-mail: o. ablyalimov@gmail. com
TO SOLUTION OF DIFFERENTIAL EQUATIONS FORM y' = f(x,y) BY THE CHORD METHOD
Oleg Ablyalimov
Doctor of philosophy, chief worker, acting professor of the chair «Loсomotives and locomotive economy»
Tashkent state ^атроП university, Uzbekistan, Tashkent
АННОТАЦИЯ
Приводится обоснование предложенного принципа решения обыкновенных дифференциальных уравнений методом хорд и классификация возможных разновидностей (приёмов) численного и графоаналитического способов решения таких уравнений.
ABSTRACT
The substantiation of the proposed principle for decision ordinary differential equations by the chord method and the classification possible varieties (methods) of numerical and graphic - analytical methods for solution such equations.
Ключевые слова: уравнение, численный метод, решение, поезд, задача, переменная, шаг, производная, полином, ряд, начальные условия.
Keywords: equation, numerical method, solution, train, problem, variable, step, derivative, polynomial, series, initial conditions.
Настоящие исследования являются продолжением работы [1], в которой были рассмотрены некоторые основные (базовые) численные методы решения обыкновенных дифференциальных уравнений вида у' = f(x,y) и посвящены обоснованию метода хорд, предложенного автором [4], сущность которого заключается в следующем.
С определённой точностью можно, производя необходимые построения и приняв соответствующие зависимости линейными в исходных интервалах ha и Ауа, найти средние значения переменных хс и ус в искомых (текущих) интервалах hn и Ауп. Нахождение конечных значений хп+г и уп+г на шаге интегрирования, затем производится или по производной для средних значений или пересчётом соответствующих координат.
Дано обыкновенное дифференциальное уравнение
y'=f(x,y) (1)
и известны начальные условия хп ,уп , за которые принимаются и полученные в процессе решения результаты на соответствующем шаге интегрирования.
Полагаем, что на заданном отрезке хп<х<хк решение существует и соблюдается условие Липшица
[3].
По начальным условиям имеем
У'= Г (хп, Уп) (2)
что позволяет при принятом некотором исходном интервале ка найти величину
Ауа = УпК (3)
а затем проверить полученное значение по условию
1АУа1 > ^Уа\ < Ы (4)
где Ау0 - наибольшая возможная величина исходного приращения, определяемая при значении У' = 0;
ап - принятая предельная величина исходного приращения. Принятая величина ап обеспечивает выбор таких величин ка и Ауа, в пределах которых принятие соответствующих зависимостей линейными обеспечивает заданную точность расчётов без
Библиографическое описание: Аблялимов О.С. К решению дифференциального уравнения вида y'=f(x,y) методом хорд // Universum: технические науки : электрон. научн. журн. 2020. № 9(78). URL: https://7universum.com/ru/tech/archive/item/10711
№ 9 (78)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
сентябрь, 2020 г,
их усложнения. Последнее приводит также к установлению ограничивающих связей между исходными и искомыми интервалами, то есть
hn < ha \Ауп\ < \Ауа\ < \an\
(5)
(6)
В пределах На значение Ауп может быть найдено по тангенсу угла наклона хорды кривой у(х), который определяется производной для средних значений переменных хс и ус в пределах искомых интервалов
&Уп = f(xcfyc)ha
(7)
и будет определяться характером заданного уравнения (1), причём хорды кривой y(x) могут располагаться:
а) по линии к - d (рис. 1), что будет при у' = у' = const., то есть когда в пределах Ауа производная не меняется;
б) по линии к - а, если у' = const. и значения у' уменьшается в пределах ha и Ауа. Так как угол наклона хорды k - а стал меньше, то в пределах hn = ha текущее приращение Ауп < Ауа и может быть 0 < Ауп < Ауа в зависимости от характера (1), поэтому
f(xctyc) =
ЬУп
(8)
. па
и очевидно получим хс = хп + — , которое останется неизменным для данного исходного шага, а значение ус может быть найдено на основе дополнительных построений, о чём будет сказано ниже;
в) по линии к - ai, если у' = const. и значение у' увеличивается в пределах рассматриваемого интервала, следовательно, угол а растёт. Учитывая (5) и (6)
будем иметь Ауп = Ауа, а 0 <кп<ка, так как значение Нп может изменяться в зависимости от характера уравнения (1), поэтому
(9)
откуда
гДе Ус =Уп +
кУа Ус Кхс.Ус)
h = АУа П Кхс,ус)
(10)
ЬУа
а значение хг необходимо
определить на основе дополнительных построений, о чем будет сказано ниже;
г) по к - e, когда у' = 0 и уп+1 =уп = const. Для этого крайнего случая Ауп = 0 и ус = уп.
Из произведённого анализа видно, что значения hn и Ауп будут зависеть от закономерностей изменения у' в пределах выбранных интервалов.
Изменения у' заданного уравнения (1) выявляется путём исследования зависимости у' = (р1(у) или х' = (р2(х) в системе у'у0 или х'хО. Так как результаты исследований в зависимостях у' = (р1(у) и х' = <р2(х) аналогичны, приведём их для случая у' = Фг(У).
В пределах половины величины ha и Ауа изменения приращения производной будет (рис. 1).
Ау = ^+д-^ = у'-у;1 (ii)
ду 2
гДе У а
'a=f(xn+hf,yn +
h
1
1
п
2
h
п
Рисунок 1. Возможные случаи зависимостей у'(у) и приращений переменных х,у в пределах исходных
интервалов Ауа и .
Возможные значения и характер изменения Ау' в i Л
пределах - Ау1:
а) Ау' = 0 и у' = Уа, что соответствует случаю, когда угол наклона хорды кривой у(х) не меняется (линия к - поэтому приращение Ау~ = Ауп = у'пЬп,
№ 9 (78)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
сентябрь, 2020 г,
то есть искомое приращение находят по производной для средних значений хс и ус в исходных интервалах
АУа и К;
б) Ау' = 0, то есть у^ > yh, что соответствует случаям увеличения угла наклона хорды кривой y(x), следовательно, на основании (5) и (6), а также см. рис. 1, Ауп — Ауа — const., а hn<ha и будет зависимой переменной, поэтому
h
¿Уа Кхс,Ус)
откуда
К =
¿Уа К*с,Ус)
(12)
(13)
где хс определится по соответствующей фор. Ауа
муле а ус = Уп + —
Рассмотренный случай соответствует таким дифференциальным уравнениям, для которых, в исход-
ду'Ауа dy'*ha _
ных интервалах + дх2 > 0;
в) Ау' < 0 и Уа <Уп т.е. угол наклона хорды кривой y(x) уменьшается в рассматриваемых интерва-
ду'Ауа , dy'*ha _ _
лах, что может быть при —--+ —-< 0.
^ ду2 дх 2
В таких случаях на основании (5) и (6), а также см. рис. 1, hn = ha = const., a Ауп < Ауа будет зависимой переменной, поэтому
откуда
где Хп = х„ +
¿^ = f(xc,yc),
¿Уп = f(xc,yc)ha
(14)
(15)
ha 2 '
а ус — определится по соответствующей формуле, полученной по дополнительным построениям;
i
г) Ау — —уп и тогда ~Ауа = 0, а уа — 0, то есть угол наклона хорды стал равен нулю. Ясно, что в этом случае Ауп — Ауа — 0 и уп+± — уп.
Линию расположения средних точек (с±,с2 текущих интервалов в зависимости от изменения у' в пределах исходных интервалов, можно выявить путём дополнительных построений в системе у'у0 по
найденным двум крайним точкам d(уп +---,yh) и е
(Уп, 0). 2
Точка d - соответствует случаю, когда Ауп — Ауа, Ау' = 0 и Уп=Уа = const., что позволяет определить положение хорды k - d кривой y(x) тангенсом
угла наклона её, то есть у^ — f(xn,yn) — f (хп +
f-yn+ir)-
Таким образом, точка d будет средней точкой искомого интервала Ауп в этом случае, так как у^ —
' _ АУа_ АУп Уа hn h„ .
Координаты этой точки всегда могут быть найдены по начальным условиям.
Вторая точка e характеризует второй крайний случай: изменения переменной у нет и производная у'е = 0. В этом случае, по условиям для точки е, имеем уп+х = уп и Ауа = 0, в силу выражения(б)
Ауп = 0, что следует также из у'е = ^ = 0.
"а
Следовательно, средняя точка этого случая ус = уп и находится в точке е (уп, 0).
В случаях, когда для уравнения (1) зависимость у' = ф±(у) будет в зоне между линиями п - dи h - e, которую в пределах А— принимаем линейной. Положение средних точек будет определяться точками пересечения линии п - а, отражающей изменения у' = (у) с прямой d - e, принятой за линию изменения положения средних точек искомых интервалов при Уа = Уп.
Написав уравнения прямых п - a и d - e можно найти точки пересечения их, по которым и будут найдены средние значения переменной ус (ординаты средних точек).
Аналогичные построения можно сделать для случаев У' > Уп, что приведёт к необходимости нахождения точек пересечения для определения абсцисс, то есть значений хс.
Ордината средней точки текущего интервала определяется следующим порядком.
у—уо у'
Для прямой п - a имели-= —т, откуда
Уп—Уо Уп
?=ТГУ"
Уп~Уо
(16)
где у0 - отрезок, отсекаемый на оси у при продолжении прямой п - a (корень уравнения прямой п - a).
Из соотношения
¿Уа
Уо-Уп
2(Уп~Уа) Уп
находим
_ Ауа*у'п
У°-2<у'п-у>)+У» Для прямой d - e имеем — L., откуда
' _ 2(У~Уп)Уп У ¿Уа
Подставив (17) в (16) получим
(17)
(18)
У-Уп-
у' =
¿Уа*Уп 2(Уп-Уа)..'
"vsn sa;
¿Уа*Уп Уп
2(у1-у'п)
(19)
Приравняем (19) и (18) и обозначив —
¿Уа*Уп 2(У'-У'„)
ш, получим
У-Уп-m ' _2{у-уп) ,
-т Уп ¿Уп Уп
(19')
Решаем относительно у. Умножим все члены на m и раскрыв скобки получим у Ауа — уп • Ауа — тАуа + 2ут — 2у„т = 0 или
1
n
№ 9 (78)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
сентябрь, 2020 г,
у(Ауа - 2т) = уп(Ауа + 2т) + тАуа> откуда
искомая ордината
уп(Ауа + 2т)+тАуа
У =
Ауа+2т
или заменив m получим
Аналогично, для случаев у 'а > у', абсцисс средней точки текущего интервала Нп
Хп(Уа-ук)+Уа(Хп+^г)
(23)
Уп
У =
2(Уп-у'о)) 2(у'-уЬ)
ДУа +
2АУа*Уп 2(У'-У')
л, * д,7 ■ УпАУаУп ■ АУа*Уп
и окончательно ордината средней точки текущего интервала Ауп, приняв у = ус
УпЫ-уЯ+укЬп^)
У =
Уп-Уа
(20)
Конечные значения переменных на шаге интегрирования в рассматриваемых случаях
Уп+1 =Уп+Кхс>Ус)
. ha
где хс = хп+— a
хп+1 = хп + ha
(21)
(22)
Конечные значения переменных
Y = у ■ АУ" хп+1 хп +к )
, Ауа
где Ус =Уп+— и
Уп+1 =Уп + АУа
(24)
(25)
На основании принципа решения, использованного в описанном методе хорд, практически можно иметь несколько возможных разновидностей (приёмов) численного метода решения обыкновенных дифференциальных уравнений, которые реализуются двумя способами - численным и графо - аналитическим. К численным способам относят варианты уточнения метода Эйлера и расчётов по у' и х', а последние разделяют на приёмы определения уп+1 пересчётом и по производной ус и хс. Графоаналитическим способом являются варианты расчётов и построений по z(xy).
Использование метода хорд к решению уравнений вида у' = f(x,y) было впервые показано в [4], а для решения дифференциального уравнения движения поезда в [5,6] и использовано авторами [2,7 и другие] при составлении ряда программ по нормированию расхода топлива на тягу поездов, анализу и оценке эффективности перевозочной работы локомотивов в условиях эксплуатации.
2у'а-у'п
Список литературы:
1. Аблялимов О. С. Уравнение движения поезда и некоторые методы его решения [Текст] / О. С. Аблялимов // Universum: технические науки: электрон. научн. журн. 2020. № 9 (78). URL: https://7universum.com/ru/tech/archive/item/10713 (дата обращения: 08.09.2020).
2. Аблялимов О. С. О методах исследования перевозочной работы локомотивов [Текст] / О. С. Аблялимов // Науч. - техн. конференция с участием зарубежных учёных «Ресурсосберегающие технологии на железнодорожном транспорте» / Ташкентский ин-т. инж. ж-д транспорта. - Ташкент, 2011. - С. 79 - 85.
3. Березин И. С. Методы вычислений [Текст] / И. С. Березин, Н. П. Жидков. - Изд. 2-е. - М.: ГИФМЛ, 1962. -Т. 1. - 464 с.
4. Толкачёв А. В. Решение дифференциальных уравнений методом хорд [Текст] / А. В. Толкачёв // Сб. Вычислительная и прикладная математика. ИК с ВЦ АН УзССР, Ташкент, 1972.
5. Толкачёв А. В. О методах решения уравнения движения поезда при расчётах на ЭЦВМ [Текст] / А. В. Толкачёв // Тр. ТашИИТ, вып. 88 / Ташкентский ин-т. инж. ж-д трансп. - Ташкент, 1972. - С. 79 - 87.
6. Толкачёв А. В. О численном методе решения уравнения движения поезда [Текст] / А. В. Толкачёв // «Вестник ВНИИЖТ» / Всесоюзный науч-иссл. ин-т. ж-д транспорта. - М.: Трансжелдориздат, 1972, № 7. - С. 53 - 59.
7. Толкачёв А. В. Исследование возможностей снижения расхода топлива тепловозами на тягу поездов за счёт оптимальных режимов вождения поездов на САЗ ж.д. [Текст] / А. В. Толкачёв, О. С. Аблялимов // Научно -технический отчёт ТашИИТ, ВНИТЦ № 79028894, Инв. № Б836191 / Ташкентский ин-т. инж. ж-д трансп. -Москва, 1979. - 67 с.
