Построение математической модели объекта в пространственном состоянии и синтез формирующего фильтра
СОДЕРЖАНИЕ
Задание
1 Построение математической модели объекта управления в пространстве
состояний
2 Построение сигнального графа и нахождение передаточной функции
системы
2.1 Построение
сигнального графа
2.2 Построение
структурной схемы
2.3 Нахождение передаточной функции ОУ, используя формулу Мейсона 8
2.4 Определение временных и частотных характеристик по передаточной
функции
3. Формирование
передаточной функции
3.1 Определение спектральной плотности по заданной корреляционной
функции
3.2 Получение
передаточной функции
3.3 Расчет
качества системы в
Вывод 20
Список
использованных источников
21
ЗАДАНИЕ
Таблица 1-Исходные данные
| Параметры элементов эквивалентной схемы объекта управления | Выходная
переменная | |||||||||||
| R1
Ом |
R2
Ом |
R3
Ом |
R4
Ом |
R5
Ом |
R6
Ом |
L1
Гн |
L2
Гн |
L3
Гн |
L4
Гн |
C1
мкФ |
C2
мкФ |
i4 |
| 343 | 435 | 151 | 204 | - | 319 | 20 | - | 48 | 25 | - | 21421 | |
i2
i1
i3
i4
C2
L3
Рисунок 1-Эквивалентная схема объекта управления
Корреляционная
функция:
1 ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ОБЪЕКТА
УПРАВЛЕНИЯ В ПРОСТРАНСТВЕ
Составим уравнения для контуров по второму закону Кирхгофа
,
,
.
Продифференцируем последнее уравнение,
чтобы избавиться от интеграла:
Введем фиктивные переменные, равные элементам, взятым из полученных
уравнений, но на 1 или более порядков ниже:
Найдем производные по времени от фиктивных переменных:
Выразим
токи через фиктивные переменные
Полученные выражения токов подставим в систему производных по времени от фиктивных переменных и система дополняется выражением для выходной величины
x1=7,15x2-24,7x1
x2=17,15x1-773,89x2+5,79x4+e(
x3=-1,87x4
x4=x3-8,52x4+3,02x2+e(t)
i4=0,4x4
По полученной системе уравнений записывается математическая модель объекта регулирования в нормальной форме Коши:
- уравнение наблюдения;
-уравнение выходной величины объекта.
Матрицы
будут иметь следующий вид:
Получаем
математическую модель в пространстве
состояний:
2 ПОСТРОЕНИЕ СИГНАЛЬНОГО ГРАФА И НАХОЖДЕНИЕ
ПЕРЕДАТОЧНОЙ ФУНКЦИИ СИСТЕМЫ
2.1 Построение сигнального графа
a11
x1 x1 x3 1/p x3
a21 1/p a12 a43 a34
e(t) b2 1/p a42 1/p c4 i4
x2 x2 x4 x4
a22 a44
b4
Рисунок
2-Сигнальный граф
2.2 Построение структурной схемы
x1 x1
e(t)
x2 x2
x4 x4
Рисунок 3-Структурная схема ОУ
2.3 Нахождение передаточной функции ОУ,
используя формулу Мейсона
Формула Мейсона имеет следующий вид:
где k –количество возможных путей от входа к выходу;
∆ -определитель графа;
Рк -коэффициент передачи k пути от входа к выходу;
∆k-определитель всех касающихся контуров при удалении k-го пути.
, то есть определитель равен, единица
минус сумма коэффициентов передачи всех отдельных контуров, плюс сумма
всех возможных произведений из двух не касающихся контуров, минус сумма
всех возможных комбинаций из трёх не касающихся контуров.
Определяем количество и записываем уравнения
всех k–путей. В данном случае имеются
k = 2 пути от входа к выходу
Выявляем все N – замкнутые контуры и записываем их уравнения. В системе имеется N = 6 замкнутых контура. Коэффициенты передачи для контуров имеют вид:
, , , , ,
, , , , ,
Записываем выражение для определителя. Определитель включает 6
контуров, пары
не касающихся контуров и комбинации из
трёх не касающихся контуров:
Выражение для
∆i записывается как выражение для
∆, но разрываются контура через которые
проходит пути Рi:
Запишем общее выражение передаточной
функции и произведём её преобразование:
Определим устойчивость системы по теореме Ляпунова.
Находим корни
характеристического уравнения
Система устойчива, так как все корни характеристического уравнения лежат
в левой полуплоскости.
2.4 Определение временных и частотных характеристик по передаточной
Функции
Найдем переходную функцию, она равна обратному преобразованию
Лапласа от , т.е. .
Рисунок 4-График переходного процесса
Определим прямые оценки качества:
hmax(t) =0.044
hуст(t) =0.006
время согласования t1 =0.003
время tH =0.4
время переходного процесса tn =0.006
Определим динамическую ошибку
системы или перерегулирование
Найдём весовую функцию, которая равна обратному преобразованию
Лапласа от W(p), то есть
Строим амплитудно-частотную
характеристику, заменив p
jw:
Рисунок 6-График АЧХ
Определим косвенные оценки качества:
- максимальная амплитуда 0.0478
- амплитуда при нулевой частоте 0.025
- резонансная частота 1.51
- показатель колебательности 1.912
- полоса пропускания частот (.
Строим фазо-частотную
характеристику:
Рисунок
7-График ФЧХ
3 Формирование передаточной
функции формирующего
фильтра
3.1 Определение спектральной плотности по заданной корреляционной
функции
Дана корреляционная
функция:
Спектральная
плотность случайного сигнала равна:
Рисунок
8- График спектральной плотности
3.2 Получение
передаточной функции формирующего фильтра
Получаем два
квадрата модуля частотной характеристики:
Находим корни
знаменателя:
Находим корни числителя:
Строим корни на комплексной плоскости:
Рисунок 9- Корни на комплексной плоскости
Из
корней верхней полуплоскости
Так как сомножитель знаменателя образуется из решения уравнения
, то его
можно заменить
непосредственно
этим уравнением.
Получаем:
Произведем замену :
3.3 Расчет качества системы в формирующем
фильтре
Найдем передаточную
функцию в виде:
Определим устойчивость системы по теореме Ляпунова.
Находим корни характеристического уравнения
=0
Два корня характеристического уравнения положительны, находятся в правой полуплоскости, следовательно, система неустойчива.
Найдем переходную
функцию:
Рисунок 10- График переходного процесса
в формирующем фильтре
Строим амплитудно-частотную
характеристику, заменив в передаточной
функции p
jw:
Рисунок
11- График АЧХ в формирующем фильтре
В данной работе была получена
математическая модель в
СПИСОК ИСПОЛЬЗОВАННЫХ
ИСТОЧНИКОВ
1.
Ефремова Т.А., Власов В.В. Методические
указания к выполнению
2. Справочное пособие по теории систем автоматического регулирования и управления под ред. Е.А.Санковского – Минск: Высшая школа, 1973.
3. Брофеев Ю.И. Импульсная техника. -М.: Высшая школа, 1984.
4.
Р.Ли Оптимальные оценки, определение
характеристик и управление. -М.: Наука,
1966.