УДК 681.3.06:681.323(519.6)
Я.Е. Ромм, Г.А. Джанунц
КУСОЧНАЯ ЛИНЕАРИЗАЦИЯ ЗАДАЧИ КОШИ ДЛЯ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ*
Кусочная линеаризация задачи Коши для обыкновенных дифференциальных уравнений (ДУ) строится на основе интерполяционного полинома Ньютона. Коэффициенты полинома восстанавливаются по его корням с использованием параллельного алгоритма. На текущем подынтервале по разностным приближениям интерполируется правая часть , , помощью его первообразной. Процесс итеративно повторяется до минимизации погрешности приближения правой части.
Кусочная линеаризация; полином Ньютона; задача Коши; обыкновенные дифферен-.
Ya.E. Romm, G.A. Dzhanunts
PIECEWISE LINEARIZATION OF THE CAUCHY PROBLEM FOR THE ORDINARY DIFFERENTIAL EQUATIONS
Piecewise linearization of the Cauchy problem for the ordinary differential equations (DE) is based on Newton's interpolation polynomial. The polynomial coefficients are reconstructed from its roots, using the parallel algorithm. At the current subinterval on difference approximations the right part DE is interpolated, coefficients of the interpolation polynomial are calculated, the solution is approaching with its primitive. The process is repeated iteratively until the approximation error of the right side is minimized.
Piecewise linearization; Newton's interpolation polynomial; Cauchy problem; ordinary differential equations.
Известные методы линеаризации позволяют преобразовать задачу Коши для обыкновенных дифференциальных уравнений (ОДУ) к линейному виду. Однако проблема линеаризации сохраняется при ограничениях на точность такой линеаризации и при требовании достаточного числа производных от правой части [1]. Ниже рассматривается вопрос о построении линеаризации задачи Коши на основе разностных методов при помощи кусочно-полиномиальной интерполяции разностных значений в ограничениях общего вида. Специфика метода, посредством которого достигается искомая точность линеаризации, обсуждается применительно к
, .
Описание исходного метода кусочно-полиномиальной аппроксимации . -полиномиальной аппроксимации функций на основе интерполяционного полинома Ньютона [2]. Аппроксимация действительной функции u = u (x) от одной действительной переменной на произвольном промежутке [а, Р] выполняется сле-. , которых покрывает [а, Р]:
*
Работа поддержана грантом РФФИ по проекту № 10-07-00178a от 2010 г. 26
р-і
[а, р]=у гх,., х,.+л, (1)
, = о
для определенности принимается Р = 2к, ке {0,1,...}. Пусть априори задана граница е абсолютной погрешности аппроксимации функции. При каждом , на , -м подынтервале из (1) строится интерполяционный полином Ньютона Т ,п (г) с
х - х{ х{+1 - х,
равноотстоящими узлами, где г =--------, к =-----------расстояние между узлами.
к п
Степень полинома п выбирается одинаковой на всех подынтервалах и минимальной при условии:
ы (х)-'
*(х )-Тіп ( )|^£> х є[ хі’ хі+0’ і = 0’ Р-1.
эютона преобраз фициентами: Тіп (і)= Ріп (х), і =
Полином Ньютона преобразуется к каноническому виду с числовыми коэф-
X - X;
к где
п ________
Р,п (х) = X а,<.х1 , х е[х!, х!+1], , = 0, Р -1. (2)
(: = 0
С этой целью рассматриваемый полином записывается в виде
1 П
1 = 1 }! к = о
П /V //
Т іп (0 = и (хі0 ) + :рП( - k), 1 =■
где х!] = х! + ] к, ] = 0, п , - узлы интерполяции, А 1ы10 - конечная разность 1 -го порядка в точке х10. Вычисляются конечные разности Аы{0 = ы(х, 1)-ы(х,0), Акы, 0 =Ак-1ы, 1 -Ак-1ы, 0, к = 2, п, значения которых обозначаются
ЪЦ =АЧ0, } = Ъ п .
В этих обозначениях
1 -1
п П( - к)
Т ,п ( ) = ы (х, 0) + %Ь,. к=\ . (3)
1=1 1!
1 -1
Каждое произведение Р, (г ) = П( - к) , -
к = 0
ложением на множители, где к = 0, , -1 его нули, по которым можно восстановить коэффициенты: Р- (г) = 0 + й]1 г + 2 г2 +... + ^ г1. Ниже используются
обозначения гг = 1,1 = 0, 1 -1. Применяется отличное от формулы Виета соот-:
1 0 0 0
й.. > і і - гі-< 1 0 0
і(і-1) І =п (! = 1 0 - ^-е 0 0
о йі 0 0 0 0 - ^-е 0 1 - V
(4)
Доказательство (4), инвариантный относительно номера коэффициента алгоритм вычисления и программная реализация приводятся в [3].
Найденные согласно (4) коэффициенты , I = 0, і, можно сохранить в памяти компьютера для любой необходимой степени і, что предполагается ниже. Полином (3) переводится в (2) приведением подобных с заменой переменной:
п
-Е а
И=1
г'
(5)
где
, (6)
, =< ]'■
Как показывает эксперимент, при сравнительно общих условиях достигается точность приближения функций с границей погрешности порядка 10~20 [2].
Табличные производную и первообразную полинома (5) можно использовать для аппроксимации производной от функции и определенного интеграла:
и {X)
1Ёаіеге~1, \и(х)йх=аЕ іп(г)йг=аЕЕ
а /- =1 ■ л - ■ л л л
Р
Р-1 п
Р-1
,£+1
і=0 о
і = 0 /:=0
£ +1
Высокая точность такой аппроксимации подтверждается численным экспериментом [2].
Изложенный метод непосредственно применим для быстрого вычисления функций в правых частях ОДУ, однако ниже он видоизменяется для аппроксимации всей правой части, приближения решения на основе первообразной от аппроксимирующего полинома и для линеаризации задачи Коши.
-
Коши с использованием разностных методов. Пусть для простоты рассматривается задача Коши для дифференциального уравнения первого порядка:
у'(х ) = / (- У)) У (х 0 ) = У 0
(7)
которую требуется решить на произвольно фиксированном промежутке [а, Ь]. Предполагается, что на [а, Ь] выполнены все условия существования и единственности решения задачи Коши. С целью избежания несущественных оговорок предполагается, что длина промежутка [а, Ь] кратна описываемой ниже величине Ни. Ку-
а
сочно-полиномиальное приближение видоизменяется для решения задачи Коши с применением разностных методов следующим образом.
На [а, Ь] строится система пересекающихся подынтервалов равной длины
т-1
[а Ь]= и[ х, x-■+l^, (8)
I = 0
где для всех г из (8) выполняется
хг+1 = хг + Ни, X0 = а, хт = Ь, (9)
Х!+1 - Х{
где и - степень конструируемого интерполяционного полинома, Н =---------------- -
и
шаг интерполяции на г -м подынтервале, который совпадает с шагом разностно-.
Номер подынтервала г интерпретируется как номер г -й группы из и +1 шагов, на которых строится текущая интерполяция решения. Ввиду кратности Ь - а и
’>- а
Ни выполняется равенство Т = -
Ни
Пусть а1 = х1, Ь1 = х1+1, где х1, х1+1 из (8), (9). Приближенное решение задачи (7) на [а, Ь] сводится к последовательному приближению на подынтервалах , Ь^, г = 0, Т -1, где значение функции в начальной точке каждого подынтервала равно её значению в конечной точке предыдущего подынтервала: У (а; ) = У (Ь;-1).
В соответствии с формулой (9) и принятыми обозначениями отрезок \^а1, Ь^ разбивается на и равных частей узлами интерполяции:
хгр = а{ + pH, р = 0, и . (10)
(10) -
,
Фгр = /(р,У(р), Р = 0 u,. (11)
где значения угр, р = 1, и определяются по разностному методу, например по :
У1р = У г (-1) + Н • / (хг (-1), У г (-1)), р = ^ , (12)
По условиям (11) строится интерполяционный полином Ньютона степени и относительно независимой переменной х, который в общем случае по схеме (3)-(6) приводится к виду полинома с числовыми коэффициентами:
V
, (13)
V іп (Х) = аі 0 + Е аИ
г=1
= Е —“7^ > ЬИ =А'Фі 0.
]=/ ]!
к ,
где аі 0 =Ф і 0, а
Полином (13) интерполирует функцию правой части (7), приближая производную искомого решения. Приближение решения строится как первообразная от
(13) с постоянной, принимающей значение у і 0. Построенный полином
ЪЛх ) =
і* п
Гуіп (х)йх = Уі0 + пЕ
(=0
1+1
X - X:
Л(+1
(14)
принимается за приближение у (х) на г -м подынтервале:
у (х )~Z г (х), х £ [^ аг, .
Полином (14) вычисляется по схеме Горнера с заменой г = -
(15)
ъ;(г ) =
( ({ а л
аіп { + і (п-1,
п +1
п
г + -
п-2)
V
п-1
г + ... + -
а_10
1
г.
(16)
Наличие коэффициентов (13) влечет готовые коэффициенты для (16). Аналогичное приближение строится на \^а1+1, Ьг+1^ и т. д., до исчерпания промежутка
[а, Ь].
Приближение решения существенно уточняется путем неоднократного решения задачи на подынтервале, где значения у!р в (11) вычисляются не разност-, -
(14), . -
ная точность достигается программной оптимизацией степени интерполяционного полинома и шага разностного метода [4]. Очевидно, в качестве исходного разностного метода можно выбрать метод более высокого порядка.
,
решения задачи Коши для систем ОДУ, приводится в [4].
-
задачи Коши используется непосредственно для построения линеаризации.
Кусочно-полиномиальная линеаризация задачи Коши. Путь требуется линеаризовать систему ОДУ
<1¥- = ^ (х, У),
йх
У (х0 ) = У0
(17)
где ^ ( У) = ( (( У), ( ( У) ..., (х, У) У = (у! (х), у2 (х) , ..., ук (х)) ,
Уо = (у 01, у 02,..., у 0 N), на промежутке [а, Ь] в предположении существования и единственности решения. Для решения поставленной задачи систему (17) требует-
ся представить в виде
йУ ~
-----= Л • У + Ь ,
йх
(18)
где матрица А и вектор Ь линейно зависят от х, решение системы (18) прибли-(17).
п
— = Аі • У, Ьі = 0, і = 0, Т -1,
йх
где Т из (8), элементы А. приближают компоненты производных, деленных на
При моделировании процесса с помощью задачи Коши для ОДУ кусочнополиномиальные приближения решения можно строить параллельно для множества различных начальных значений и множества границ погрешности. Моделирующие решения воспроизводятся параллельно по множеству подынтервалов -достаточно одновременно считать из памяти хранимые коэффициенты интерполирующих решение полиномов. Линеаризация (19) влечет возможность компьютерного анализа устойчивости решений систем ОДУ путем вычисления соответственных матричных произведений [5] на основе (19) параллельно по всем подынтерва-
, (19)
этом случае априори известны.
Численное моделирование с помощью данного метода характеризуется малой , -ма метода и использования памяти компьютера.
1. Березин И.С., Жидков Н.П. Методы вычислений. - М.: Физматгиз, 1962. Т. 2. - 640 с.
2. Аксайская Л.Н. Разработка и исследование параллельных схем цифровой обработки сигналов на основе минимизации временной сложности вычисления функций: Автореф. дис. ... канд. техн. наук. - Таганрог: ТТИ ЮФУ, 2008.
3. Ромм Я.Е. Локализация и устойчивое вычисление нулей многочлена на основе сортировки // Кибернетика и системный анализ. - 2007. - № 2. - С. 161-174.
4. РоммЯ.Е., ДжанунцГ.А. Компьютерный метод разностно-анадитического решения обыкновенных дифференциальных уравнений на основе интерполяционного полинома Ньютона / ТГПИ. - Таганрог, 2009. - 40 с. Деп. в ВИНИТИ 18.06.09, № 379-В2009.
5. РоммЯ.Е. Моделирование устойчивости по Ляпунову на основе преобразований разно-
//
моделирование. - 2008. - № 12. - С. 105-118.
. . ., . . .
,
У г
йУ_
йх
0
0
(19)
0
і = 0, Т-1.
V
0
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
Ромм Яков Евсеевич
Таганрогский государственный педагогический институт.
E-mail: romm@List.ru.
347926, г. Таганрог, ул. Инициативная, 48.
Тел: 8634601899.
Кафедра информатики; заведующий кафедрой; д.т.н.; профессор.
Джанунц Гарик Апетович
E-mail: janunts@inbox.ru.
Тел.: +79185069024.
Romm Y akov Evseevich
Taganrog State Pedagogical Institute.
E-mail: romm@List.ru.
48, Initsiativnaya Street, Taganrog, 347926, Russia.
Phone: +7634601899.
The Department of Information; Head the Department; Dr. of Eng. Sc.; Professor.
Dzhanunts Garik Apetovich
E-mail: janunts@inbox.ru.
Phone: +79185069024.
The Department of Information; Postgraduate Student.
УДК 621.391:519.21
E.B. Моисеева
АЛГОРИТМ ИДЕНТИФИКАЦИИ ПРОМЫШЛЕННОГО ОБЪЕКТА ПО ЕГО ВРЕМЕННОЙ И ЧАСТОТНОЙ ХАРАКТЕРИСТИКАМ,
ДЛЯ ЦЕЛЕЙ ОБУЧЕНИЯ НА ТРЕНАЖЕРНОМ СТЕНДЕ
Представлен метод идентификации реальной переходной характеристики, полученной с помощью учебного тренажерного комплекса “Иерархическая автоматизированная система контроля и управления процессом нагрева". Идентификация является обязательным элементом и наиболее сложной стадией процесса решения прикладных задач. В общем виде задача идентификации заключается в определении оператора объекта, преобразующего входные воздействия в выходные.
Идентификация; транспортное запаздывание; набор инерционных звеньев; учебный .
E.B. Moiseeva
ALGORITHM OF IDENTIFICATION OF THE INDUSTRIAL OBJECT WITH ITS CHARACTERISTICS FOR TIME AND FREQUENCY FOR PURPOSES OF LEARNING ON TRAINING STAND
Presented a method to identify the actual transfer characteristic obtained through academic training complex, "Hierarchical self-matizirovannaya control and manage the process of heating. Iden-Katsia is a must and the most difficult stage in the process of solving practical problems. In general, the identification problem is to identify the operator of a facility that converts the input actions on the weekend.Identification; distance-velocity lag; set of inertial links; educational training stand.
Identification; distance-velocity lag; set of inertial links; educational training stand.