
- •Глава 1 динамические системы в задачах обработки навигационной информации
- •1.1. Математические модели динамических систем и методы их описания
- •1.1.1. Определение и классификация динамических систем
- •1.1.2. Определение выхода линейных динамических систем с помощью фундаментальной матрицы и весовой функции
- •1.1.3. Передаточная функция стационарных линейных динамических систем
- •1.1.4. Частотная характеристика стационарных линейных динамических систем
- •1.1.5. Взаимосвязь между основными функциями, используемыми при описании линейных динамических систем
- •Основные функции, используемые при описании линейных систем
- •Взаимосвязь различных функций, используемых для описания поведения линейных стационарных систем
- •1.1.6. Определение модели линейной динамической системы в пространстве состояний по заданной передаточной функции
- •1.1.7. Основные свойства линейных динамических систем
- •Задачи к разделу
- •Контрольные вопросы
- •1.2. Линейные стационарные динамические системы
- •1.2.1. Интеграторы
- •1.2.2. Фильтры Баттерворта
- •1.2.3. Модель акселерометра
- •1.2.4. Модель микромеханического гироскопа
- •1.2.5. Простейшая модель ошибок построения вертикали в инерциальной системе
- •Вклад уходов гироскопов в ошибки определения нп в автономном режиме в течение 3 мин
- •Вклад ошибок акселерометров в ошибки определения нп в автономном режиме в течение 3 мин
- •Задачи к разделу
- •Контрольные вопросы
- •1.3. Дискретизация и моделирование линейных
- •1.3.1. Дискретизация непрерывных систем
- •1.3.2. Основные методы описания lti-объектов в Matlab
- •Создание lti-объектов
- •Извлечение информации о моделях
- •1.3.3. Особенности дискретизации стационарных систем
- •Методы дискретизации, используемые в функции c2d
- •Методы вычисления матричной экспоненты
- •Контрольные вопросы
- •1.4. Задание для моделирования
- •Пример выполнения задания в Matlab
- •Заключение к главе 1
1.1.5. Взаимосвязь между основными функциями, используемыми при описании линейных динамических систем
Введенные в предыдущих разделах функции, используемые для описания поведения динамических систем, могут быть сведены в табл. 1.1.1.
Т а б л и ц а 1.1.1
Основные функции, используемые при описании линейных систем
Название функции |
Обозначение |
Назначение |
Примечание |
Фундаментальная матрица |
|
Определяет общее решение
|
Используется как для стационарных, так и нестационарных систем |
Весовая функция, (импульсная переходная функция) |
|
Устанавливает связь выхода и входа при нулевых начальных условиях:
|
|
Передаточная функция |
|
Устанавливает связь изображений выхода и входа:
|
Используется для стационарных систем |
Частотная характеристика |
|
Устанавливает связь гармонических колебаний на выходе и входе:
|
Взаимосвязь между этими функциями для случая стационарных линейных систем отражается в табл. 1.1.2.
Т а б л и ц а 1.1.2
Взаимосвязь различных функций, используемых для описания поведения линейных стационарных систем
,
Отыски-ваемая функция |
Задаваемая функция |
|||
|
|
|
|
|
|
|
* |
* |
* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В
табл. 1.1.2 использованы обозначения,
введенные в приложении 1 для прямых
(
,
)
и обратных
,
преобразований
Фурье и Лапласа.
На диагонали этой таблицы приведены выражения, которые можно трактовать как исходные определения соответствующих функций для рассматриваемого класса систем, а вне диагонали – выражения, устанавливающие взаимосвязь между ними.
Следует
также отметить, что выражение
,
определяющее тот факт, что ЧХ может быть
найдена как преобразование Фурье от
весовой функции, следует воспринимать
с учетом замечания, приведенного в конце
подраздела 1.1.4.
Звездочка * в верхней графе означает, что не существует однозначного соответствия между задаваемыми и отыскиваемыми функциями. Так, например, зная ПФ, невозможно однозначно описать систему в форме пространства состояния. Причины отсутствия такого перехода обсуждаются в следующем подразделе.
1.1.6. Определение модели линейной динамической системы в пространстве состояний по заданной передаточной функции
Ранее предполагалось, что в качестве исходного описания системы задано ее представление в форме пространства состояния. При решении прикладных задач нередко возникает необходимость получить модель системы в пространстве состояния, располагая ее передаточной функцией. Обсудим этот вопрос более подробно.
Передаточной функции (1.1.22) соответствует уравнение, записанное в операторной форме, т.е. с помощью преобразования Лапласа
(1.1.37)
Во временной области этому уравнению можно однозначно сопоставить дифференциальное уравнение с постоянными коэффициентами -го порядка
(1.1.38)
Специфика этого уравнения заключается в том, что в правой части у него содержатся q-е производные от входного воздействия. Суть обсуждаемой задачи заключается в получении системы в форме Коши с помощью передаточной функции. Задача определения уравнений состояния по ПФ сводится к известной в теории дифференциальных уравнений задаче приведения линейных уравнений n-го порядка к системе уравнений в нормальной форме Коши [1, 2]. Простейший пример решения этой задачи рассмотрен в приложении 2 (пример 6).
Заметим, что решение задачи приведения линейных уравнений n-го порядка к системе уравнений в нормальной форме Коши неоднозначно, т.е. одной ПФ соответствует множество описаний системы в пространстве состояний в виде (1.1.14), (1.1.15). Покажем это.
Пусть задана ПФ
,
соответствующая системе (1.1.14), (1.1.15). С
использованием невырожденной, не
зависящей от времени
матрицы
,
введем новую переменную
. (1.1.39)
Поскольку
и
,
а следовательно,
,
можем записать
.
Или
.
Выход системы с новой переменной связан как
.
Вычисляя ПФ для преобразованной системы, имеем
Таким образом, для системы, имеющей ПФ , имеем новое описание в виде:
; (1.1.40)
, (1.1.41)
в котором ,
,
,
.
(1.1.42)
Преобразование называется преобразованием базиса, а , как отмечается в [73, с. 455], преобразованием подобия. Очевидно, что существует бесконечное множество преобразований подобия. Поскольку ПФ систем, связанных преобразованием подобия, совпадают, то каждой ПФ также соответствует бесконечное множество описаний в пространстве состояний. Следует также иметь в виду, что одной и той же ПФ, кроме того, могут соответствовать различные описания в пространстве состояний, не связанные между собой преобразованием подобия [1, с. 50].
Пример 1.1.4. Проиллюстрируем неоднозначность представления системы с заданной ПФ в пространстве состояний на примере системы, имеющей ПФ
,
(1.1.43)
в
которой
– круговая частота.
Рассмотрим систему, задаваемую уравнениями:
(1.1.44)
.
Поскольку
в этом случае
,
,
,
то с использованием (1.1.19) нетрудно
убедиться в том, что ПФ для этой системы
имеет вид (1.1.43).
Введя
преобразование
с матрицей
и используя (1.1.41), получим следующее
представление:
(1.1.45)
,
в
котором
,
,
.
Легко убедиться, что и для этой системы ПФ также имеет вид (1.1.43).
Рассмотрим два наиболее важных варианта получения описания системы в пространстве состояния для SISO-системы, ПФ которой задана в виде (1.1.22).
Первый вариант основан на использовании так называемой управляемой канонической формы для системы в пространстве состояния, для которой матрицы, определяющие описание системы записываются как [1]:
;
(1.1.46)
,
.
(1.1.47)
Этим матрицам соответствуют уравнения:
(1.1.48)
.
(1.1.49)
Матрица (1.1.46) называется матрицей Фробениуса. Нетрудно убедиться в том, что соответствующий ей характеристический полином, определяемый из уравнения , задается следующим соотношением:
.
Иными словами, стоящие на последней
строке элементы матрицы
представляют собой коэффициенты
характеристического полинома.
Используя (1.1.20), можно показать, что для системы (1.1.48), (1.1.49) ПФ будет задаваться выражением (1.1.22).
Пример 1.1.5. Убедимся в справедливости приведенного утверждения на двухмерном примере, т.е. когда
,
,
.(1.1.50)
Используя эти
соотношения и выражения (1.1.20), (П2.28) для
ПФ и матрицы
,
запишем
,
откуда с очевидностью следует
. (1.1.51)
Особенность приведенного варианта получения описания системы в пространстве состояний заключаются в том, что входное воздействие присутствует в правой части уравнения только для последней компоненты вектора состояния, а выходной сигнал формируется в виде линейной комбинации компонент вектора состояния. Рассмотрим другой вариант получения описания системы в пространстве состояния для SISO-системы, ПФ которой задана в виде (1.1.22).
Этот вариант основан на использовании
наблюдаемой канонической формы,
для которой матрица динамики имеет вид
(1.1.46), а матрицы
и
записываются как
Элементы
матрицы
вычисляются через коэффициенты
многочленов
методом неопределенных коэффициентов.
Можно также использовать следующую
рекуррентную формулу
Этим матрицам соответствуют уравнения:
(1.1.52)
И в этом случае легко убедиться, что для такой системы передаточная функция будет задаваться выражением (1.1.22).
Особенность второго варианта представления системы в пространстве состояний с ПФ вида (1.1.22) заключается в том, что одно и тоже входное воздействие присутствует в правых частях уравнений для всех компонент вектора состояний, а выход системы соответствует его первой компоненте.
Пример 1.1.6. Убедимся в справедливости приведенного утверждения на примере системы второго порядка, т.е. покажем, что при
;
;
(1.1.53)
передаточная функция будет иметь вид (1.1.51).
Используя эти соотношения и выражения (1.1.20), (П2.28) для ПФ и матрицы , запишем
.
Откуда с очевидностью получаем (1.1.51).