- •Анализ представленных результатов позволяет сформулировать следующие выводы.
- •Литература
- •Среднее время ожидания заявок
- •Литература
- •Н.А. Рубина, Ю.Г. Кирчин
- •Литература
- •Приложение
- •Access
- •Power Point
- •Литература
- •Word
- •Ecxel
- •КОНЦЕПТУАЛЬНЫЕ МОДЕЛИ В БАЗАХ ЗНАНИЙ
- •Введение
- •Рис. 4. Окно меню "Химический состав"
- •Рис. 5. Окно меню "Поиск"
- •Литература
- •Введение
- •Технологии проектирования ВсС
- •Традиционные подходы к проектированию ВсС. Ключевыми чертами традиционного процесса проектирования микропроцессорных вычислительных систем следует считать:
- •Архитектурные абстракции сквозного проектирования ВсС
- •Опыт использования архитектурных абстракций в проектировании ВсС
- •Заключение
- •Литература
- •ДИНАМИЧЕСКИЕ ИНСТРУМЕНТАЛЬНЫЕ КОМПОНЕНТЫ
- •Введение
- •Обзор вариантов построения тестовых систем
- •Заключение
- •Литература
- •УНИВЕРСАЛЬНАЯ МОДЕЛЬ ПОСЛЕДОВАТЕЛЬНОСТНЫХ СХЕМ В САПР
- •А.Г. Зыков, О.Ф. Немолочнов, В.И. Поляков
- •Рис. Универсальная модель последовательностной схемы
- •Рис.1. Схема пересчёта
- •ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ ОПТИЧЕСКИХ ПРИБОРОВ СО ВСТРОЕННОЙ СИСТЕМОЙ СТАБИЛИЗАЦИИ ИЗОБРАЖЕНИЯ
- •Рис. 2. Имитационная модель ОССИ. (АЛ – алгоритм имитации события; Мt – временной фактор)
- •Введение
- •Интегральные параметры.
- •Спектральные параметры.
- •(4) Интегральная яркость изображения вычисляется по формуле
- •Структурные параметры.
- •ВОССТАНОВЛЕНИЕ ЛИНИЙ ЭКСТРЕМУМОВ СЛОЖНЫХ КАРТИН ИНТЕРФЕРЕНЦИОННЫХ ПОЛОС
- •М.В. Волков
- •Введение
- •Обработка одномерных сигналов на основе модификации гистограмм
- •Обработка картины полос
- •Восстановление линий экстремумов интерференционных полос
- •Примеры обработки реальных интерферограмм
- •Заключение
- •Литература
- •О.В. Павлушко
- •Dimage 7
- •Olimpus E-10
- •Заключение
- •ОПТИЧЕСКИЕ ВОЛОКОННЫЕ УСИЛИТЕЛИ ИНФОРМАЦИОННЫХ СИГНАЛОВ
- •Н. С. Макаров
- •Введение
- •Эрбиевые волоконные усилители
- •Тулиевые волоконные усилители
- •ВКР волоконные усилители
- •Гибридные усилители
- •Стокс-антистоксовые ВКР-усилители
- •Заключение
- •Введение
- •Метод фильтрации Калмана
- •Фильтр Калмана второго порядка
- •Фильтрация акустических сигналов
- •Заключение
- •А. Акунова, А.В. Ушаков
- •Литература
- •2. Постановка задачи
- •3. Синтез алгоритма адаптации
- •1. Введение
- •3. Синтез алгоритма управления
- •Основной результат предлагаемой работы сформулирован в следующей теореме.
- •Заключение
- •Литература
- •И.В. Мирошник, А.Н. Шалаев
- •А. А. Мельников, Е. В. Рукуйжа, А. В. Ушаков
- •О.В. Слита, И.В. Мирошник
- •Литература
- •ОЦЕНКА ЗАПАСОВ УСТОЙЧИВОСТИ СИСТЕМ С ИНТЕРВАЛЬНЫМИ ПАРАМЕТРАМИ
- •Введение
- •Основные положения
- •Литература
- •КОНВЕРГЕНЦИЯ ЭЛЕКТРОННЫХ И ПЕЧАТНЫХ СРЕДСТВ ИНФОРМАЦИИ
- •СИНТЕЗ ИНФОРМАЦИОННО-ЭНЕРГЕТИЧЕСКИХ СВЯЗЕЙ В СТРУКТУРАХ ИЗОБРЕТЕНИЙ
- •А.Б. Бушуев
- •Акунов Т.А., Ушаков А.В. Связь алгебраических спектров собственных
- •значений и сингулярных чисел в задаче обеспечения стабильности
& |
γ |
T |
(D |
T |
P + PD)x , |
|
|
ˆ |
|
x |
(6) |
||||
θ = |
|
|
|
||||
2 |
|
|
|||||
|
|
|
|
|
|
ˆ |
|
где коэффициент γ > 0 влияет на скорость изменения параметра θ. |
|
||||||
Допущение. Предположим, что существует малый коэффициент |
γ > 0 , такой |
||||||
что величина θ&ˆ незначительна и для производной по матрице P выполнено условие:
ρP − P& > 0 ,
где число ρ > 0 .
Основной результат предлагаемой работы сформулирован в следующей теореме. Теорема. Пусть пара (A +βD, B) - полностью управляема для любых значений β .
Пусть допущение выполнено, тогда существует малое положительное число γ такое,
что закон управления (3) при выполнении условия (5), обеспечивает асимптотическую устойчивость положения равновесия x = 0.
Доказательство. Рассмотрим функцию Ляпунова вида:
V = x |
T |
|
|
~ |
2 |
/ γ , |
|
|
|
|
|
|
|
(7) |
||
|
Px + θ |
|
|
|
|
|
|
|
|
|||||||
~ |
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
где θ = θ− θ - параметрическая невязка. |
|
|
|
|
||||||||||||
Дифференцируя уравнение (7), получаем: |
|
|
||||||||||||||
& |
= x |
T |
|
& |
|
~~& |
T |
T |
|
+ θ(D |
T |
P |
+ PD) − 2PBBP)x . |
(8) |
||
V |
|
Px + 2 |
θθ/ γ + x |
|
( A P + PA |
|
||||||||||
Из уравнения (4) находим: |
|
|
|
|
|
|
||||||||||
A |
T |
P + |
|
|
|
|
|
ˆ |
T |
P + PD) |
|
(9) |
||||
|
PA − PBBP = −2α P − θ(D |
|
|
|||||||||||||
и, подставляя в (8), имеем следующее неравенство:
& |
≤ x |
T & |
~~& |
/ γ + x |
T |
ˆ |
T |
T |
θ(D |
T |
P + PD)x . |
(10) |
V |
Px + 2 |
θθ |
|
(−2αP − θ(D |
|
P + PD))x + x |
|
|
|
|
|
~ |
|
ˆ |
|
~& |
|
& |
& |
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
||
Учитывая, что θ = θ− θ и |
|
θ = θ − θ = −θ для выражения (10), получаем: |
|
||||||||||||||||||||||
|
|
& |
≤ x |
T & |
~ |
|
T |
(D |
T |
P + PD)x + x |
T |
|
|
|
~ |
|
T |
P + PD))x = |
|
||||||
|
|
V |
Px |
− θ x |
|
|
|
(−2αP + θ(D |
|
|
|||||||||||||||
= x |
T & |
|
T |
Px . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(11) |
||
Px − 2α x |
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
|
|
|
|
||||||
Предположим, что |
коэффициент |
выбран |
таким, |
что скорость изменения |
|||||||||||||||||||||
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
∂P & |
|
∂P |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
& |
|
ˆ |
|
|
|
|
|
|
|
|
|
||
параметра |
θ незначительна и матрица |
P = |
ˆ θ |
+ |
|
|
& |
|
|
|
|
|
|||||||||||||
∂α |
α мала (см. допущение). Тогда, |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂θ |
|
|
|
|
|
|
|
|
|||
найдется такое α, что для неравенства (11) будет выполнено: |
|
|
|
|
|
||||||||||||||||||||
V& ≤ −δxT Px ≤ 0 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(12) |
|||||
где число δ > 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
равновесия x = 0 |
|
|
||||||
Из неравенства |
(12) |
следует, |
что |
положение |
устойчиво |
по |
|||||||||||||||||||
Ляпунову, |
вектор |
x |
- ограниченный и квадратично интегрируемый, |
а сигнал |
ˆ |
||||||||||||||||||||
θ - |
|||||||||||||||||||||||||
ограничен. |
Ограниченность |
x |
и |
|
ˆ |
обеспечивает ограниченность |
правой части |
||||||||||||||||||
|
θ |
||||||||||||||||||||||||
дифференциального уравнения (2) с законом управления (3). Последнее, гарантирует выполнение следующего условия:
lim x(t) = 0 ,
t→∞
что наряду с устойчивостью по Ляпунову положения равновесия |
x = 0, обеспечивает |
|||||
также и его асимптотическую устойчивость. |
|
|
||||
Замечание 1. Следует отметить, что выбор параметра γ |
представляет |
собой |
||||
отдельную задачу, решение которой не имеет места в рамках предлагаемой работы. |
||||||
Замечание 2. Выбор параметра α, удовлетворяющего |
условию (5), |
может |
||||
производиться из следующих соображений: |
|
|
||||
α = α0 + |
|
ˆ |
|
, |
|
(13) |
|
|
|
||||
|
θ |
|
|
|||
230
где число α0 выбирается из условия: |
|
|
ˆ |
D} > 0 . |
(14) |
Re λ i {α0 I + A + θ0 |
||
Заключение
Результаты предлагаемой статьи не претендуют на универсальность и общность, а лишь расширяют класс подходов управления линейной системой в условиях неопределенности ее параметров. В работе рассматривается линейный стационарный объект, переменные состояния которого подлежат измерению. Предполагая, что ряд несложных условий выполнен, выбирается закон управления, обеспечивающий асимптотическую устойчивость по части переменных состояния системы, а именно по всем переменным объекта, но не его регулятора.
Литература
1.Фомин В.Н., Фрадков А.Л., Якубович В.А. Адаптивное управление динамическими объектами. М. Наука, 1981.
2.Фрадков А.Л. Адаптивное управление в сложных системах. М.: Наука, 1990.
3.Андриевский Б.Р., Фрадков А.Л., Избранные главы теории автоматического управления с примерами на языке MATLAB. СПб.: Наука, 1999.
4.Мирошник И.В., Никифоров В.О., Фрадков А.Л. Нелинейное и адаптивное управление сложными механическими системами. СПб.: Наука, 2000.
5.Никифоров В.О. Робастное управление линейным объектом по выходу // Автоматика и телемеханика. 1998. № 9.
6.Никифоров В.О., Фрадков А.Л. Схемы адаптивного управления с расширенной ошибкой. Обзор // Автоматика и телемеханика. 1994. №9.
7.Бобцов А.А., Лямин А.В., Сергеев К.А. Синтез закона адаптивного управления для стабилизации не точно заданных нестационарных объектов // Изв. Вузов Приборостроение. - 2001. - №3.
8.Дроздов В.Н., Мирошник И.В., Скорубский В.И. Системы автоматического управления с микроЭВМ. Л.: Машиностроение. Ленингр. отд-ие, 1989.
231
КРОНЕКЕРОВСКАЯ МАТРИЧНАЯ МОДЕЛЬ ДИНАМИЧЕСКИХ ПРОЦЕССОВ С МОДУЛЯЦИЕЙ
Л.В. Кожевникова, А.В. Ушаков
Решается задача конструирования модели динамических процессов с модуляцией. Для указанных целей используются возможности кронекеровских векторных и матричных структур, получаемую при этом избыточную размерность модели предлагается редуцировать с помощью оператора сужения.
Введение. Постановка задачи
Системы управления и следящие системы с модуляцией составляют заметную часть практики автоматического управления. Модуляторами в таких системах являются: сельсины, поворотные трансформаторы, индуктивные датчики, полудисковые модуляторы лучистой энергии и т.д. [1].
Однако теоретические исследования процессов в системах с модуляцией, иногда именуемых системами на несущей переменного тока [2], в последнее время заметно сократились, причем это произошло на фоне интенсификации исследования динамических систем с использованием возможностей метода пространства состояния
(МПС) [3, 4].
Метод пространства состояния характеризуется значительной алгебраизацией общей теории систем, в рамках которой получены хорошие наработки применительно к кронекеровским матричным структурам [5]. Обнаруживается, что возможности этих структур могут быть эффективно использованы применительно к исследованию систем с модуляцией, причем получаемые модельные представления процессов с модуляцией оказываются максимально адаптированными к инструментарию МПС. Проблемам конструирования моделей динамических процессов с модуляцией на основе используемых возможностей кронекеровских векторных и матричных структур, инвариантных относительно размерности вход – выходных отношений посвящается данная работа.
Свойства векторных и матричных кронекеровских структур
Для целей погружения в сформулированную выше проблему приведем определение векторных и матричных кронекеровских структур, а так же те их свойства, которые непосредственно сориентированы на построение матричной модели процессов
с модуляцией. |
|
|
|
|
|
|
|
|
|
X и Y , |
||||
Определение 1 |
(О.1). Кронекеровским произведением двух векторов |
|||||||||||||
X Rn Y Rm , называется вектор X Y , |
составленный из сепаратных произведений |
|||||||||||||
{X iY j ; i = |
|
; j = |
|
} |
их элементов так, что становится справедливым представление |
|||||||||
1, n |
1, m |
|||||||||||||
X Y = col{X iY j ; i = |
|
|
; j = |
|
} , где |
X Y R nm |
(1) |
|||||||
1, n |
1, m |
|||||||||||||
Примечание 1 |
(П.1). Очевидно, |
кроме |
кронекеровского произведения |
|||||||||||
X Y векторов может быть построено также произведение Y X векторов, причем, в |
||||||||||||||
общем случае эти произведения оказываются не коммутативными так, что |
|
|||||||||||||
X Y ≠Y X |
(О.2). Если размерности векторов X и Y одинаковы, |
|
||||||||||||
Определение 2 |
то на их |
|||||||||||||
кронекеровском произведении X Y может быть построено согласованное сужение |
||||||||||||||
этого произведения ( X Y ) S , задаваемого представлением |
|
|||||||||||||
(X Y )S = col{X iYi ; i = |
|
} |
|
|
(2) |
|||||||||
1, n |
|
|
||||||||||||
Примечание 2 (П.2). Согласованное сужение кронекеровского векторного |
||||||||||||||
произведения X Y |
может быть осуществлено |
с помощью оператора сужения с |
||||||||||||
матрицей S вида
232
S = diag{[0Ti−1 M 1 M 0Tn−i ]; i = |
1, n |
} |
(3) |
так, что становится справедливой запись |
|
||
( X Y ) S = S( X Y ) |
(4) |
||
В качестве одного из свойств кронекеровского произведения векторов рассмотрим правило дифференцирования кронекеровских векторных произведений по скалярному параметру, причем в основном сосредоточимся на случае, когда скалярным параметром является время.
Свойство 1 (СВ.1). Дифференцирование векторной кронекеровской структуры в виде их кронекеровского произведения осуществляется по правилам дифференцирования сложной функции, представленной в мультипликативной форме так, что
d |
• |
|
|
|
|
( X (t) Y (t)) =( X (t) Y (t)) = X&(t) Y (t) + X (t) Y&(t) |
(5) |
||
dt |
||||
d |
|
• |
|
|
( X (t) Y (t) Z(t)) = |
( X (t) Y (t) Z(t)) = |
|
||
dt |
(6) |
|||
|
|
|||
X&(t) Y (t) Z(t) + X (t) Y&(t) Z(t) + X (t) Y (t) Z&(t)
Определение 3 (О.3). Кронекеровским произведением прямоугольных матриц A Rn×m , B R p×q называется матрица ( A B) размерности (np × mq) , составленная в
силу соотношения |
|
||||
A B = col{row( Ai, j B; j = |
|
);i = |
|
} |
(7) |
1, m |
1, n |
||||
Примечание 3 (П.3): Кронекеровское произведение произвольных прямоугольных |
|||||
матриц не обладает коммутативностью так, что |
|
||||
A B ≠ B A |
(8) |
||||
Задача конструирования матричной модели динамических процессов с модуляцией в своей основе использует квадратные матрицы, коими являются матрицы состояния системы, конечномерного источника внешнего воздействия и конечномерного источника модулирующего сигнала, поэтому ниже имеется в виду
класс квадратных матриц. |
квадратных матриц A Rn×n и |
Определение 4 (О.4). Кронекеровской суммой |
|
B Rm×m называется матрица ( A B) , размерности |
(nm × nm) , составленная в силу |
соотношения |
|
A B = A I B + I A B , |
(9) |
где I A , I B - единичные матрицы, согласованные по размерности соответственно с
матрицами А и В.
Примечание 4 (П.4). Для кронекеровской суммы квадратных матриц А и В, а в общем случае и произвольного числа матриц, существует альтернативное название – преобразование Сильвестра матриц, обозначаемое Si{A, B) , тогда оказывается
справедливой запись
A B = A I B + I A B = Si{A, B} |
(10) |
Для случая 3-х квадратных матриц А, В, С кронекеровская сумма или их преобразование Сильвестра будет записано в форме
Si{A, B, C} = A B C = A I B IC + I A B IC + I A I B C |
(11) |
Отметим тут же, что как и кронекеровское произведение матриц кронекеровская сумма не коммутативна.
Кронекеровские матричные структуры обладают следующими свойствами.
233
Свойство 2 (СВ.2). Алгебраический спектр собственных значений кронекеровского произведения A B квадратных матриц A R n×n и B R m×m как матричной функции от матриц обладает тем свойством, что его элементы образованы попарными произведениями собственных значений кронекеровски перемножаемых матриц
σ{ A B} = {μ K : det( μI − A B) = 0; μ K = λ Ai λBj ; i = 1, n; j = 1, m; k = 1, mn} (12)
Свойство 3 (СВ.3). Алгебраический спектр собственных значений кронекеровской суммы A B квадратных матриц A R n×n и B R m×m как матричной функции от матриц обладает тем свойством, что его элементы образованы попарными суммами
собственных значений кронекеровски суммируемых матриц |
|
||||||
σ{A B} ={νl : det(νl I − A B) = 0;νl = λAi + λBj ; i = |
|
; j = |
|
; l = |
|
} |
(13) |
1, n |
1, m |
1, mn |
|||||
В (12) и (13) λAi и λBj собственные значения матриц А и В. |
|
||||||
Сделаем следующее примечание к свойствам (СВ.2) и (СВ.3).
Примечание 5 (П.5). Алгебраические спектры собственных значений кронекеровских произведений A B и B A в силу (12) совпадают, аналогичным свойством обладают и спектры кронекеровских сумм A B и B A.
Свойство 4 (СВ.4). Определитель кронекеровского произведения квадратных
матриц удовлетворяет соотношению |
|
det(A B) = (det A) m (det B)n , |
(14) |
где A Rn×n и B Rm×m .
Свойство 5 (СВ.5). След кронекеровской суммы квадратных матриц удовлетворяет соотношению
tr( A B) = m trA + n trB , |
(15) |
где A Rn×n и B Rm×m .
Свойство 6 (СВ.6). Ранг кронекеровского произведения квадратных матриц удовлетворяет условию
rang ( A B) = rangA rangB , |
(16) |
где A Rn×n и B Rm×m .
Для целей решения поставленной задачи полезно напомнить основные свойства кронекеровских матричных произведений произвольных матриц, которые оказываются необходимыми при преобразованиях матричных композиций, содержащих в своем составе эти произведения.
Свойство 7 (СВ.7).
(P Q)(W V ) = PW QV |
(17) |
Свойство 8 (СВ.8). |
(18) |
(P + Q) R = P R + Q R |
|
P (Q + R) = P Q + P R |
(19) |
P (Q R) = (P Q) R |
(20) |
В (17) – (20) матрицы P, Q, R, W, V имеют произвольные размерности, не противоречащие правилам перемножения и сложения матриц.
Свойство 9 (СВ.9).
P Q = (P IQ )(I P Q) |
(21) |
(P1 Q1 )(P2 Q2 )K(PR QK ) = (P1 P2 KPK ) (Q1Q2 KQK ) |
(22) |
(P Q)−1 = P−1 Q−1 |
(23) |
I (P1 P2 KPK ) = (I P1 P1 )(I P2 P2 )K(I PK PK ) |
(24) |
В выражениях (21) – (24) IQ, IP и т.д. – единичные матрицы, по размерности согласованные с соответствующими матрицами Q, P и другими аналогичными.
234
Свойство 10 (СВ.10). Оператор сужения кронекеровского произведения векторов с матрицей сужения S удовлетворяет соотношению
S(PX QZ ) = S(P Q)( X Z ) |
(25) |
Основной результат
Воспользуемся приведенными в предыдущем разделе свойствами векторных и матричных кронекеровских структур для целей построения динамической модели процессов в линейной многомерной непрерывной системе с модуляцией. При построении моделей процессов будем предполагать, что источник внешнего воздействия (ИВВ) является конечномерным и он представим автономной системой; будем предполагать, что модулирующий сигнал также является конечномерным и потому источник модулирующего сигнала (ИМС) также представим автономной системой. Таким образом полное априорное описание задачи приобретает вид
X&(t) = FX (t) + Gv(t); X (0); Y (t) = CX (t) |
(26) |
Z&(t) = ΓZ(t); Z(0); g(t) = PZ(t) |
(27) |
ZM (t) = ΓM ZM (t); ZM (0); gM (t) = PM ZM (t) |
(28) |
В модели (26) многомерной непрерывной системы Х – векторы состояния; v – вектор внешнего воздействия; Y – вектор выхода; X Rn ; v,Y Rm ; F, G, C - матрицы состояния входа и выхода соответственно; F Rn×n , CT , G Rn×m .
В модели (27) источника внешнего воздействия Z и g - вектора состояния и выхода ИВВ соответственно; Z Rl ; g Rm ; Γ, P - матрицы состояния и выхода;
Γ Rl×l ; P R m×l .
В модели (28) источника модулирующего сигнала ZM и gM - векторы состояния и выхода ИМС соответственно; Z M R K ; gM Rm ; ΓM , PM - матрицы состояния и выхода
ИМС; ΓM RK×K ; PM Rm×K .
Процесс формирования модулированного внешнего воздействия v(t) представим
в форме |
|
v(t) = col{g j (t) gMj (t); j =1, m} |
(29) |
Нетрудно видеть, что процесс модуляции внешнего воздействия в форме (29) допускает представление его в виде кронекеровского произведения векторов с последующим сужением, что может быть представлено в форме:
v(t) = S(g(t) gM (t)) |
(30) |
Если теперь учесть правила формирования g(t) и gM (t) представленных в (27) и
(28), то (30) в силу свойств кронекеровских произведений матриц может быть записано в форме
v(t) = S(g(t) gM (t)) = S(PZ(t) PM ZM (t)) = S(P PM )(Z(t) ZM (t)) (31)
Выражение (31) представляет модулированный сигнал v(t) как функцию состояния системы с вектором состояния Z(t) ZM (t) .
Сформируем динамическую систему, с вектором состояния Z(t) ZM (t) ,
опираясь на модели (27) и (28), описывающих процессы по компонентам введенного вектора являющегося кронекеровским произведением векторов, а также свойства матричных кронекеровских структур. В результате получим следующую цепочку равенств
235
