close

Вход

Забыли?

вход по аккаунту

?

Построение общей функции Ляпунова для семейства механических систем с одной степенью свободы.

код для вставкиСкачать
УДК 531.36
Вестник СПбГУ. Сер. 10, 2013, вып. 4
И. Е. Мурзинов
ПОСТРОЕНИЕ ОБЩЕЙ ФУНКЦИИ ЛЯПУНОВА
ДЛЯ СЕМЕЙСТВА МЕХАНИЧЕСКИХ СИСТЕМ
С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ∗)
Исследование устойчивости систем c переключениями является одной из актуальных задач современной теории управления [1–3]. Система с переключениями представляет собой гибридную динамическую систему, состоящую из семейства подсистем
и закона переключения, определяющего в каждый момент времени, какая подсистема
активна [1, 2]. Системы такого рода широко применяются в задачах управления механическими, энергетическими, электроэнергетическими системами, при управлении
технологическими процессами, а также в ряде других областей [1, 2]. Во многих случаях необходимо проверить, является ли система устойчивой либо диссипативной при
произвольных переключениях. Данные задачи возникают естественным образом, если
закон переключения неизвестен либо слишком сложен для того, чтобы использовать
его точное представление.
Основной подход к решению задач на исследование устойчивости или диссипативности системы при произвольных переключениях состоит в построении общей для всех
подсистем функции Ляпунова, удовлетворяющей теореме Ляпунова об асимптотической устойчивости [4] либо теореме Йосидзавы о равномерной диссипативности [5].
Этот подход успешно применяется (см. [1, 2, 6] и цитируемую там литературу). Тем
не менее условия существования общей функции Ляпунова не полностью изучены даже для семейства линейных автономных систем [6–8]. Для двумерных и трехмерных
линейных систем получены необходимые и достаточные условия существования общей
квадратичной функции Ляпунова [3, 7], однако для систем большей размерности это
удалось сделать только при дополнительных требованиях, например требовании коммутативности матричных коэффициентов [2, 9].
Еще сильнее данная задача усложняется для механических систем. Их движение
описывается уравнениями второго порядка, что приводит к возникновению некоторых
особых свойств. В таких системах переключения могут происходить в действующем
силовом поле. Причины переключений могут быть как внутренними, например использование компьютера или микропроцессора в контуре управления [10], так и внешними,
если движение механической системы происходит в среде с изменяющимися параметрами. Даже если действующие силы являются линейными, при таких переключениях
условия коммутативности будут заведомо нарушены [11], поэтому соответствующие результаты о существовании общей квадратичной функции Ляпунова неприменимы к механическим системам. Такие особые свойства механических систем приводят к необходимости выделять их в отдельный класс систем с переключениями, который обладает
некоторыми теоретическими особенностями и несомненной практической значимостью.
В настоящей статье исследована нелинейная гибридная механическая система
Мурзинов Илья Евгеньевич – студент, 199034, Санкт-Петербургский государственный университет;
e-mail: murz42@gmail.com.
∗) Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 13-01-00376) и Санкт-Петербургского государственного университета (НИР
№ 9.38.674.2013).
c И. Е. Мурзинов, 2013
49
с одной степенью свободы, которая описывается дифференциальным уравнением второго порядка, содержащим два параметра (коэффициенты демпфирования и жесткости).
Предполагается, что эти параметры могут переключаться с одних значений на другие. Для соответствующего семейства подсистем получены достаточные условия существования общей функции Ляпунова специального вида. Их выполнение обеспечивает
асимптотическую устойчивость или диссипативность соответствующей гибридной системы при произвольных переключениях.
Постановка задачи. Рассмотрим семейство нелинейных механических систем с одной степенью свободы
ẍ + ai ẋxα + bi xβ = 0,
i = 1, . . . , N.
(1)
Уравнения такого вида называются уравнениями Льенара [12]. Здесь x(t) – состояние
системы, ai и bi – постоянные положительные коэффициенты, α > 0 – рациональное
число с четным числителем и нечетным знаменателем, β > 1 – рациональное число
с нечетными числителем и знаменателем. Каждая подсистема имеет асимптотически
устойчивое положение равновесия x = ẋ = 0 [13].
Зададим кусочно-постоянную функцию σ = σ(t) : [0, +∞)−→Q = {1, . . . , N }, определяющую закон переключения. В данной статье предполагается, что на каждом ограниченном временном интервале закон переключения имеет конечное число точек разрыва, которые называются моментами переключения, и принимает постоянное значение между любыми двумя последовательными моментами переключения. Такие законы
переключения будем называть допустимыми. Семейство (1) вместе с законом переключения описывают гибридную систему
ẍ + aσ ẋxα + bσ xβ = 0.
(2)
Задача 1. Получить условия на коэффициенты ai и bi , при которых нулевое решение уравнения (2) асимптотически устойчиво для любого допустимого закона переключения.
Задача 2. Получить условия на коэффициенты ai и bi , при которых решения уравнения (2) предельно ограничены для любого допустимого закона переключения. В этом
случае система называется диссипативной [5].
При α = 0, β = 1 и α = 0, β > 1 необходимые и достаточные условия существования общей функции Ляпунова специального вида для семейства (1) получены
соответственно в [14] и [15].
В качестве функции Ляпунова выбираем функцию
V =
xβ+1
y2
+ c1
+ c2 xμ y + c3 xy γ ,
2
β+1
(3)
где c1 > 0, c2 > 0, c3 < 0; μ и γ – рациональные числа с нечетными числителями
и знаменателями, μ 1, γ 1.
Исследование асимптотической устойчивости. Производная функции Ляпунова (3) в силу i-го уравнения семейства (1) имеет вид
V̇ = (c1 − bi )xβ y + c2 μxμ−1 y 2 + c3 y γ+1 − ai xα y 2 − ai c2 xα+μ y −
− ai c3 γxα+1 y γ − bi c2 xμ+β − c3 bi γxβ+1 y γ−1 = W (x, y).
50
Найдем условия на коэффициенты ai и bi и показатели степеней α и β, при которых
можно выбрать параметры c1 , c2 , c3 , μ, γ таким образом, чтобы функция Ляпунова (3)
являлась положительно-определенной, а ее производная в силу подсистем из семейства (1) – отрицательно-определенной. Тогда функция Ляпунова (3) будет удовлетворять теореме Ляпунова об асимптотической устойчивости [5].
Для нахождения таких условий применим следующие леммы.
n
Лемма 1 [4]. Рассмотрим функцию f (x) =
cs xρs s , где cs > 0, ρs – положиs=1
тельные рациональные числа с четными числителями и нечетными знаменателями. Если для функции r(t, x), заданной при x < H и t 0, справедлива оценка
|r(t, x)| L|x1 |ϕ1 |x2 |ϕ2 ...|xn |ϕn , где H > 0, L > 0, ϕi 0, i = 1, . . . , n, то, для того
чтобы при любом значении коэффициента L функция f (x)+r(t, x) была положительно определена, необходимо и достаточно выполнения неравенства
n
ϕs
s=1
ρs
> 1.
Лемма 2 [4]. Рассмотрим функцию f (x) = c1 xρ11 + c2 xρ22 + c3 xδ11 xδ22 , где c1 , c2 , c3 –
положительные коэффициенты, ρi , δi – положительные рациональные числа с четными числителями и нечетными знаменателями. Будем считать, что ρδ11 + ρδ22 < 1.
Пусть для функции r(t, x), заданной при t 0, x < H, справедлива оценка |r(t, x)| L|x1 |ϕ1 |x2 |ϕ2 , H > 0, L > 0, ϕi 0, i = 1, 2. Тогда, для того чтобы при любом значении коэффициента L функция f (x) + r(t, x) была положительно определена, необходимо и достаточно выполнения неравенств
ϕ1 +
ρ1 − δ 1
ϕ2 > ρ1 ,
δ2
ρ2 − δ 2
ϕ1 + ϕ2 > ρ2 .
δ1
При помощи леммы 1 находим следующие необходимые условия положительной
определенности функции Ляпунова (3):
β+1
,
2
2β
.
γ
β+1
μ
Предположим cначала, что
α
2
+
< 1.
μ+β
γ+1
Используя лемму 2, получим, что для отрицательной определенности функции W (x, y)
необходимо, чтобы выполнялись условия
μ = β − α,
β
,
γ β−α
μ α + 1,
γ
γ
γ
2β−2α−1
β−α ,
2α+1
α+1 ,
β+2α+1
β+1 .
51
Если β 2α + 1, то эти условия могут быть записаны в виде
μ = β − α,
γ
(4)
2β−2α−1
β−α .
Заметим, что при выполнении условий (4) функция Ляпунова (3) положительно определена.
Если β < 2α+1, то применение леммы 2 приводит к противоречию в системе условий
отрицательной определенности функции W (x, y). Поэтому в данном случае пользуемся
леммой 1, с помощью которой получаем следующие необходимые условия положительной определенности функции Ляпунова (3) и отрицательной определенности функции
W (x, y):
2β
,
γ = β+1
μ=
β+1
2 .
Применение лемм 1 и 2 позволяет не учитывать слагаемые более высоких порядков
при исследовании отрицательной определенности функции W (x, y), что существенно
упрощает задачу в случае β 2α + 1. Однако при β < 2α + 1 анализ оставшихся слагаемых представляет собой довольно сложную проблему, поэтому в настоящей статье
рассмотрены только случаи β > 2α+1 и β = 2α+1. Стоит заметить, что при β = 2α+1
системы
"
ẋ = y,
ẏ = −ai xα y − bi xβ , i = 1, . . . , N,
соответствующие уравнениям семейства (1), являются обобщенно-однородными [13].
Случай β > 2α + 1. В этом случае μ = β − α, и функцию W (x, y) можно записать
в виде
W = xα (c1 − bi − ai c2 )xβ−α y − ai y 2 − bi c2 x2β−2α + c2 (β − α)xβ−α−1 y 2 −
− ai c3 γxα+1 y γ − c3 bi γxβ+1 y γ−1 + c3 y γ+1 .
Если γ > 2β−2α−1
β−α , то для отрицательной определенности функции W (x, y) достаточно
отрицательной определенности функции
g(x, y) = (c1 − bi − ai c2 )xβ−α y − ai y 2 − bi c2 x2β−2α ,
являющейся квадратичной формой относительно переменных xβ−α и y.
В соответствии с критерием Сильвестра функция g(x, y) отрицательно определена
тогда и только тогда, когда выполнены неравенства
4bi ai c2 > (c1 − bi − ai c2 )2 ,
i = 1, . . . , N,
'
ai bi c2 > |c1 − bi − ai c2 | ,
i = 1, . . . , N.
или
2
Следовательно,
√
√
' 2
' 2
a i c 2 − b i < c1 <
a i c2 + b i ,
52
i = 1, . . . , N.
Для существования требуемого числа c1 необходимо и достаточно, чтобы имели место
неравенства
√
' √
'
i, j = 1, . . . , N.
a i c2 − b i < a j c2 + b j ,
Получаем условия
'
√ √
√ √
c2 ai + aj > √bi − 'bj ,
√
√
√
c2
ai − aj < b i + b j ,
i, j = 1, . . . , N.
Для всех i, j = 1, . . . , N положим
'
bi − bj
√
√ ,
ai + aj
√
A=
⎧
⎪
⎨+∞,
max
i,j=1,...,N
'
√
B=
bi + bj
⎪
⎩ min √
√ ,
ai >aj
ai − aj
если ai = aj ,
(5)
i, j = 1, . . . , N,
если не все ai равны друг другу.
(6)
Таким образом, справедлива следующая теорема.
Теорема 1. Пусть β > 2α + 1. Тогда для существования общей функции Ляпунова вида (3) для семейства (1), удовлетворяющей требованиям теоремы Ляпунова об
асимптотической устойчивости, достаточно, чтобы выполнялось неравенство A<B,
где A и B определяются формулами (5) и (6).
Следствие 1. Пусть β > 2α + 1. Тогда для асимптотической устойчивости системы с переключениями (2) достаточно, чтобы выполнялось неравенство A<B, где
A и B определяются формулами (5) и (6).
З а м е ч а н и е 1. Теорема 1 дает конструктивный алгоритм построения функции
Ляпунова.
Следствие 2. Если переключения происходят только в скоростных силах ( bs =
b = const > 0), то положение равновесия x = ẋ = 0 системы (1) является асимптотически устойчивым при любом допустимом законе переключения и любых положительных значениях коэффициентов ai и b.
Следствие 3. Если переключения происходят только в потенциальных силах
( as = a = const > 0), то положение равновесия x = ẋ = 0 системы (1) является
асимптотически устойчивым при любом допустимом законе переключения и любых
положительных значениях коэффициентов a и bi .
Следствие 4. Если система (1) содержит две подсистемы, то положение равновесия x = ẋ = 0 системы (1) является асимптотически устойчивым при любом
допустимом законе переключения и любых положительных значениях коэффициентов ai и bi .
З а м е ч а н и е 2. Следствия 3 и 4 не справедливы для линейных систем ( α = 0,
β = 1) [6].
Случай β = 2α + 1. В данном случае для отрицательной определенности функдостаточно отрицательной определенности
ции W (x, y) при μ = β − α, γ > 2β−2α−1
β−α
функции
g(x, y) = (c1 − bi − c2 ai )xα+1 y + c2 (α + 1 − ai )y 2 − c2 bi x2α+2 .
53
Эта функция является квадратичной формой относительно переменных xα и y. Пользуясь критерием Сильвестра, получаем систему неравенств, выполнение которых обеспечивает положительную определенность функции Ляпунова и отрицательную определенность ее производной:
2 c2 bi ai − c2 (α + 1) > |c1 − bi − c2 ai |,
i = 1, . . . , N,
ai
,
i = 1, . . . , N,
c2 <
α+1
2
c1 > c2 (α + 1).
Из этой системы легко определить условия существования параметра c1
'
'
bi − c2 ai − c2 (α + 1) < bj + c2 aj − c2 (α + 1) ,
i, j = 1, . . . , N.
(7)
Для нахождения условий разрешимости этих неравенств применим результаты, полученные в [14]. В указанной работе была изучена устойчивость линейной системы
(α = 0, β = 1) и на аналогичном этапе проведено исследование отрицательной определенности производной функции Ляпунова, которое свелось к вопросу существования
решения системы (7). √
√
√
bi −
√
bj
bi −
bj
Положим Θij = √α+1 , ξij = √α+1 . При всех i, j = 1, . . . , N определим числа
min
2 max
ρij , ρmax
следующим образом:
ij , ρij
а) при bi < bj положим
"
*
ai
aj
,
ρmax
=
min
;
ij
α+1 α+1
б) при bi > bj и
√
ai aj (α + 1)ξij положим
ρmin
ij
"
= max
aj
ai
,
α+1 α+1
*
+ 1,
= −1;
ρmax
ij
+
,
√
aj
ai
в) при bi > bj , ai aj < (α + 1)ξij и min α+1
, α+1
|ai − aj | (α + 1)ξij положим
ρmin
ij =
=
ρmax
ij
2
2
(α + 1)2 (ai + aj ) − 2 ai aj − (α + 1)ξij
ξij
2
(ai − aj )2 + 4(α + 1)ξij
,
2
2
(α + 1)2 (ai + aj ) + 2 ai aj − (α + 1)ξij
ξij
;
2
(ai − aj )2 + 4(α + 1)ξij
√
aj
ai
г) при bi > bj , ai aj < (α + 1)ξij и min{ α+1
, α+1
}|ai − aj | > (α + 1)ξij положим
ρmin
ij =
54
2
2
(α + 1)2 (ai + aj ) − 2 ai aj − (α + 1)ξij
ξij
2
(ai − aj )2 + 4(α + 1)ξij
,
"
ρmax
ij
= min
ai
aj
,
α+1 α+1
д) при ai < aj положим
ρ2ijmax =
е) при ai aj и
*
;
ai
;
α+1
'
aj (ai − aj ) > (α + 1)Θij положим
ρ2ijmax =
aj
;
α+1
'
aj (ai − aj ) (α + 1)Θij положим
Θ2ij (α + 1)2 (ai + aj ) + 2 ai aj − (α + 1)Θ2ij
ρ2ijmax =
.
(ai − aj )2 + 4(α + 1)Θ2ij
ж) при ai aj и
Пусть теперь
A = max ρmin
ij ,
(8)
2 max
B = min ρmax
.
ij , ρij
(9)
i,j
i,j
При A < B система (7) имеет решение, что обеспечивает существование искомых значений параметров c1 и c2 функции Ляпунова (3). Таким образом, справедлива следующая
теорема.
Теорема 2. Пусть β = 2α + 1. Тогда для существования общей функции Ляпунова вида (3) для семейства (1), удовлетворяющей требованиям теоремы Ляпунова об асимптотической устойчивости, достаточно, чтобы выполнялось неравенство
A < B, где A и B определяются формулами (8) и (9).
Следствие 5. Пусть β = 2α + 1. Тогда для асимптотической устойчивости системы с переключениями (2) достаточно, чтобы выполнялось неравенство A < B, где
A и B определяются формулами (8) и (9).
Пример. Рассмотрим семейство, состоящее из трех уравнений:
ẍ + x2 ẋ + 2x7 = 0,
ẍ + 3x2 ẋ + x7 = 0,
ẍ + 2x2 ẋ + kx7 = 0,
где k – положительный коэффициент. Функцию Ляпунова будем строить в виде
V (x, y) =
x8
y2
+ c1
+ c2 x4 y + c3 xy γ .
2
8
(10)
Здесь γ > 9/5, c1 > 0, c2 > 0, c3 < 0. Пользуясь теоремой 1, выпишем достаточное
условие существования для этого семейства общей функции Ляпунова типа (10):
"
*
*
"
1 k−2 k−1
3
max
,
,
, k + 1, k + 2 .
< min
4
3
5
2
Решая это неравенство, получим условие k < 13/2.
55
Исследование диссипативности. Для исследования системы (2) на равномерную
диссипативность обратимся к теореме Йосидзавы [5]. Функцию Ляпунова опять строим
в виде (3). Производная этой функции – W (x, y). Для того чтобы найти условия на параметры функции Ляпунова, при которых существуют Δ1 > 0, Δ2 > 0, такие, что
сама функция положительна при (x, y) > Δ1 , а ее производная – отрицательна при
(x, y) > Δ2 , применим лемму (см. [16], с. 180). Для изучаемого случая она может
быть записана таким образом.
Лемма 3. Рассмотрим функцию f (x) = c1 xρ11 + c2 xρ22 + c3 xδ11 xδ22 , где c1 , c2 , c3 –
положительные коэффициенты; ρi , δi – положительные рациональные числа с четδ1 δ2
ными числителями и нечетными знаменателями. Будем считать, что
+
> 1.
ρ1 ρ2
Пусть для функции r(t, x), заданной при t 0, x > H, справедлива оценка |r(t, x)| L|x1 |ϕ1 |x2 |ϕ2 , H > 0, L > 0, ϕi 0, i = 1, 2. Тогда, для того чтобы при любом значении коэффициента L существовало положительное число Δ, такое, что функция
f (x) + r(t, x) < 0 при x > Δ и t 0, необходимо и достаточно выполнения неравенств
ρ1 − δ 1
ρ2 − δ 2
ϕ2 < ρ1 ,
ϕ1 + ϕ2 < ρ2 .
(11)
ϕ1 +
δ2
δ1
З а м е ч а н и е 3. Если неравенства (11) в лемме 3 нестрогие, то такие условия
являются необходимыми.
Используя свойства обобщенно-однородных функций [5], выпишем необходимые
условия существования положительного Δ1 , такого, что функция Ляпунова V (x, y) > 0
при (x, y) > Δ1 . Имеем
β+1
μ
,
2
(12)
2β
.
γ
β+1
Из первого условия леммы 3 с учетом замечания 3 для слагаемых (c1 − bi )xβ y
и ai c2 xα+μ y в выражении для W (x, y) получаем соответственно неравенства μ β − α
и μ β − α. Следовательно, должно иметь место равенство μ = β − α. Рассматривая
оставшиеся слагаемые с учетом этого соотношения, определяем необходимые условия
существования положительного Δ2 , такого, что W (x, y) < 0 при (x, y) > Δ2 .
Если β 2α + 1, то
μ = β − α,
γ
2β−2α−1
β−α .
Поскольку μ = β − α и μ 1, с помощью леммы 3 возможно исследовать только
случай β α + 1. Заметим, что при таких значениях параметров μ и γ условия (12)
выполнены.
Повторяя рассуждения, проведенные при изучении асимптотической устойчивости,
приходим к следующим результатам.
Теорема 3. Пусть α + 1 β < 2α + 1. Тогда для существования общей функции
Ляпунова вида (3) для семейства (1), удовлетворяющей требованиям теоремы Йосидзавы о равномерной диссипативности, достаточно, чтобы выполнялось неравенство
A < B, где A и B определяются формулами (5) и (6).
Следствие 6. Пусть α + 1 β < 2α + 1. Тогда для равномерной диссипативности
системы с переключениями (2) достаточно, чтобы выполнялось неравенство A < B,
где A и B определяются формулами (5) и (6).
56
Теорема 4. Пусть β = 2α + 1. Тогда для существования общей функции Ляпунова
вида (3) для семейства (1), удовлетворяющей требованиям теоремы Йосидзавы о равномерной диссипативности достаточно, чтобы выполнялось неравенство A < B, где
A и B определяются формулами (8) и (9).
Следствие 7. Пусть β = 2α+1. Тогда для равномерной диссипативности системы
с переключениями (2) достаточно, чтобы выполнялось неравенство A < B, где A и B
определяются формулами (8) и (9).
Литература
1. DeCarlo R., Branicky M., Pettersson S., Lennartson B. Perspectives and results on the stability and
stabilisability of hybrid systems // Proc. IEEE. 2000. Vol. 88. P. 1069–1082.
2. Hai Lin, Antsaklis P. J. Stability and stabilizability of switched linear systems: a survey of recent
results // IEEE Trans. Automat. Control. 2009. Vol. 54, N 2. P. 308–322.
3. Shorten R., Wirth F., Mason O., Wulff K., King C. Stability Criteria for Switched and Hybrid
Systems // SIAM Rev. 2007. Vol. 49, N 4. P. 545–592.
4. Александров А. Ю., Платонов А. В. Устойчивость движений сложных систем. СПб.: НИИ химии
С.-Петерб. ун-та, 2002. 79 с.
5. Демидович Б. П. Лекции по математической теории устойчивости. М.: Наука, 1967. 469 с.
6. Liberzon D., Morse A. S. Basic Problems in Stability and Design of Switched Systems // IEEE
Control Syst. Magazin. 1999. Vol. 59, N 15. P. 59–70.
7. Пакшин П. В., Поздяев В. В. Критерий существования общей квадратичной функции Ляпунова
множества линейных систем второго порядка // Изв. РАН. Теория и системы управления. 2005. Вып. 4.
С. 22–27.
8. Narendra K. S., Balakrishnan J. A common Lyapunov function for stable LTI systems with commuting
A-matrices // IEEE Trans. Automat. Control. 1994. Vol. 39, N 12. P. 2469–2471.
9. Unsolved Problems in Mathematical Systems and Control Theory / eds. by V. D. Blondel
& A. Megretski. Princeton: Princeton University Press, 2004. 334 p.
10. Kosov A. A., Vassilyev S. N., Zherlov A. K. Logic-based controllers for hybrid systems // Intern.
J. Hybrid Syst. 2004. Vol. 4, N 4. P. 271–299.
11. Aleksandrov A. Yu., Chen Y., Kosov A. A., Zhang L. Stability of Hybrid Mechanical Systems with
Switching Linear Force Fields // Nonlinear Dynamics and Systems Theory. 2011. Vol. 11, N 1. P. 53–64.
12. Зубов В. И. Математические методы исследования систем автоматического регулирования. Л.:
Судостроение, 1959. 324 с.
13. Rouche N., Habets P., Laloy M. Stability Theory by Liapunov’s Direct Method. New York
e. a.: Springer, 1977. 300 p.
14. Васильев С. Н., Косов А. А. Общие функции Ляпунова и вектор-функции сравнения Матросова в анализе гибридных систем // Труды X Междунар. Четаевской конференции. Т. 2. Секция 2:
Устойчивость. 12–16 июня. Казань, 2012. С. 162–176.
15. Aleksandrov A. Yu., Murzinov I. E. On the Existence of a Common Lyapunov Function for a Family
of Nonlinear Mechanical Systems with One Degree of Freedom // Nonlinear Dynamics and Systems Theory.
2012. Vol. 12, N 2. P. 137–143.
16. Александров А. Ю., Платонов А. В. Метод сравнения и устойчивость движений нелинейных
систем. СПб.: Изд-во С.-Петерб. ун-та, 2012. 263 с.
Статья рекомендована к печати проф. А. П. Жабко.
Статья поступила в редакцию 30 мая 2013 г.
Документ
Категория
Без категории
Просмотров
6
Размер файла
302 Кб
Теги
построение, степенью, система, одной, функции, свобода, механической, ляпунова, общее, семейство
1/--страниц
Пожаловаться на содержимое документа