CC BY
23
ТЕОРЕМА 1. Предположим, что в точке с е V" первый дифференциал аналитической функции f(x) на компакте V" не равен тождественно нулю. Обозначим
Л/*
v,(/,c)= min ord^ —~{с) isi<n öXj
и положим
- ß*=ß(/,C)=min{ß|v1(/,c)<i2(/)+ß,ßeZ0}.
Тогда на компакте /С^с.р13 ^ = c + 7tp V" ряд Тейлора Tcf изометрически
эквивалентен своей линейной части.
ТЕОРЕМА 2 (модификация леммы Гензеля). В обозначениях теоремы 1 уравнение
Tcf(x) = О
разрешимо в компакте К^с,рр j тогда и только тогда, когда разрешимо в этом же компакте линейное уравнение
/(c)+Zf-(cX*,-c)= 0. 1=1 °Xi
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 1. Серр Ж.11. Алгебры Ли и группы Ли. М.: Мир, 1969.
УДК 517.938; 519.711.3 Е. В. Дивнсенко, В. В. Мозжилкин
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ МЕТОДОМ ПРЯМЫХ
1. Математические модели механических систем описываются краевой задачей для вектора А'(г,г)
. д2Х пдХ „ д2Х дгХ дАХ.
dt2 dt v " ' dz2 ' &3 ' 024
*<z,0)=/(z), Щ-ot
X(0,t) = q(t), ~
dz
лг(1,о = у(/), Щ-
dz
»=o
= "(0.
z=0
= w(f).
z=l 27
(1)
Здесь г - пространственная переменная, Г - время. Некоторые компоненты вектора X могут зависет ь только от переменной Л
Важной задачей исследования свойств решений системы (1) является определение предельных циклов и периодических решений. Одним из подходов к решению данной задачи является анализ с помощью рядов Фурье [ 1 ]. Суть этого подхода такова. Вектор X представляется отрезком ряда Фурье по переменной /. Коэффициенты ряда определяются из решения краевой задачи системы ОДУ по переменной г. Общее решение представляет собой линейную комбинацию нескольких линейно независимых частных решений однородной системы ОДУ и частного решения неоднородного уравнения в виде линейной функции. В результате подстановок строится нелинейное векторное уравнение для частоты колебаний и коэффициентов отрезка ряда Фурье, которое решается численно.
Такой подход позволяет определить периодические решения. Однако, если подобное решение не удается построить, это не гарантирует его отсутствия. Этот метод достаточно жёстко навязывает структуру решения исходной краевой задачи (1), которая не может быть реализована в общем случае. Если решение нелинейной алгебраической системы не найдено, или оно не удовлетворяет физическим требованиям, формально этот факт свидетельствует только о том, что конкретный решатель не справился с задачей. Данный метод сильно усложняется при применении его к исследованию произвольных движений механических систем. Поэтому целесообразно строить приближенное решение (1) таким образом, чтобы не навязывать поведение решения во времени. Одним из возможных подходов является метод прямых.
2. Разобьём счётную область г е [0,1], / > 0 на п полос прямыми
г = г{= Ш, ¿ = 0,.., п, И = Мп.
Аппроксимируем производные по пространственной переменной г конечно-разностными соотношениями
8ГХ
дгг
ч
:=хА >=-Р
Центральная разностная аппроксимация производной Х-4) (2 < / < п - 2) с р,ц = 2, коэффициентами
~А4' А4' Л4' Л4'/«4
имеет порядок аппроксимации 504А4_у(£,г), г,_2 < ^ < г,>2-
Построим соответствующие разностные аппроксимации в окрестности границ 2 = 0, 2 = 1, учитывающие граничные условия.
Разностная аппроксимация четвёртой производной по г в точке имеет вид
дХ
*Г= X! ] +
/=-1 С2
где
(а,Ь) =
13 I 4 I
2=0 2 111
_ 43 1 I
14 1
тр"
Разностная аппроксимация четвёртой производной по 2 в точке имеет вид
у=-5 №
(3)
2=1
Её коэффициенты
1 1
2 1 14 1
13 1
43 1 14 I
18 А4' 5 А«'/,«' 9 А4' 2 /,'' 14 /И' 5 А4' 3 /,' Порядок аппроксимации соотношений (2), (3) равен 2184Л >"(£,, 0> для (2) £е[0,2А], а для (3) $е[1,1-2А].
Для аппроксимации вторых и третьих производных от X по 2 воспользуемся следующими соотношениями.
Граница 2=0. Разностная аппроксимация имеет вид
д'Х
дгг
2=0
= ±а^Х1+ЬдХ
Д-1Я аппроксимации второй производной с четвёртым порядком справедлива формула с = 3 и коэффициентами
415 I 1 1
" 72 Л2' А2' * Л2'9 А2' 8 /г' 6 Л
Её погрешность аппроксимации равна - 48Л £,е[0>46].
Граница 2=1. Разностная аппроксимация имеет вид дгХ.
= > а ' х . .. +с
&
8гГ
z—ih 1=~Р
2=1
Для аппроксимации второй производной с четвёртым порядком получается формула с р = 4 и коэффициентами
(я,5) =
118 1 1 I 415 1 25
8 А2'9 А2' " Л2' /г' 72 А2' 6 А
Её пофешность аппроксимации равна -48Л £ е[1-4Л,1].
Для аппроксимации третьей производной с четвёртым порядком получается формула с р = 5 и коэффициентами
107 1 77 1 343 1 45 I " 4 А3' 2 А" 16 А'' 4 А2
Её погрешность аппроксимации равна - 1644А4у(г,£), £ е [1 - 4АД].
131
2 Л3' 16 А3' А''
В результате подстановки разностных формул в краевую задачу (1) получим задачу Коши для системы ОДУ вида
АХ1+ВХ1=Г(г1,1,Х1,...,Хп_1), ^(0) = /„ Х,=*„ ¿ = 1,2,..м-\. (4)
Эта задача Коши решается любым численным методом для ОДУ.
Система (4) может быть исследована аналитически. Легко определить её особые точки и поведение решения в их окрестности. Это позволяет качественно оценить поведение динамической системы, не ограничиваясь задачами нахождения предельных циклов и периодических решений.
Применение методов теории катастроф [2] к анализу системы (4) позволит получить информацию о характере переходных процессов в механической системе.
Данная методика была реализована программна в системе Ма1ЕаЬ и была успешно применена к исследованию динамики упругого манипулятора.
БИКЛИОГРЛФИЧКСКИЙ СПИСОК
1. Андрейченко К.П., Андрейченко Д.К. Математическое моделирование динамических систем. Саратов, 2000.
2. Постои Т., Стюарт И. Теория катастроф и её приложения. М.: Мир, 1980.
УДК 519.85
А. С. Дудова
УСЛОВИЯ ЗВЁЗДНОСТИ ЛЕБЕГОВА МНОЖЕСТВА ДИФФЕРЕНЦИРУЕМОЙ ПО НАПРАВЛЕНИЯМ ФУНКЦИИ'
I. Давно известное в действительном и комплексном анализе понятие звёздного множества ныне активно используется в рамках абстрактного выпуклого анализа [1]. Напомним, что множество Асй" называется звёздным относительно точки д: 6 Л, если
/+а(х-х')бА, Ухе А, ае[0,1]. Это - широкий класс, множеств, включающий в себя выпуклые множества, конусы, объединение выпуклых множеств с непустым пересечением. Операции сложения, умножения на число сохраняют звёздность. Кроме того, если у некоторых звёздных множеств есть общие точки звёздности, то есть относительно которых они являются звёздными, то пересечение и объединение таких множеств также будет звёздным. Это говорит о том, что с та-
' Работа выполнена при финансовой поддержке фанта Президента РФ на поддержку ведущих научных школ на выполнение научных исследований (проект ИШ-1295.2003.1).
