Построение математической модели управления в пространстве состояний

Построение сигнального графа и структурной схемы системы управления. Расчет передаточной функции системы по формуле Мейсона. Анализ устойчивости по критерию Ляпунова. Синтез формирующего фильтра. Оценка качества эквивалентной схемы по переходной функции.
Краткое сожержание материала:

Размещено на

Размещено на

СОДЕРЖАНИЕ

ИСХОДНЫЕ ДАННЫЕ

1. Построение математической модели ОУ в пространстве состояния

2. Построение сигнального графа и структурной схемы системы

3. Нахождение передаточной функции системы по формуле Мейсона

4. Анализ устойчивости системы по критерию Ляпунова

5. Определение прямых оценок качества системы

6. Синтез формирующего фильтра

7. Оценка качества эквивалентной схемы по переходной функции

ВЫВОД

Список используемой литературы

ИСХОДНЫЕ ДАННЫЕ

Рисунок 1. Эквивалентная схема объекта управления

R₁ = 345 Ом

R₂ = 375 Ом

R₃ = 127 Ом

R₄ = 244 Ом

R₆ = 250 Ом

L₃ = 47 Гн

L₄ = 25 Гн

C₁ = 40511 мкФ

С₂ = 23007 мкФ

i₂=?

1. ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ОУ В ПРОСТРАНСТВЕ СОСТОЯНИЯ

e(t) i₄

Рисунок 2. Структурная схема объекта управления

Cоставим уравнения по 2-му закону Кирхгоффа для контуров:

(1)

В системе уравнений (1) следует избавиться от всех интегралов, продифференцировав в данном случае первое и четвертое уравнения этой системы:

Используя метод условного интегрирования, вводим фиктивные переменные, равные элементам, взятым из уравнений на 1 или более порядков ниже.

(2)

Находим производные по времени от фиктивных переменных и, применяя предыдущие уравнение, выражаем зависимостями от токов и фиктивных переменных.

(3)

Из систем уравнений (1) - (3) найдем выражения для токов через фиктивные переменные:

(4)

Подставим числовые данные в систему (4) и найдем i₁, i₂, i₃, i₄:

Полученные значения токов подставим в систему (3) и найдем , , , , дополнив систему выражением для выходной величины:

По полученной системе уравнений и уравнению для выходной величины объекта регулирования запишем математическую модель в нормальной форме Коши:

- уравнение наблюдения;

- уравнение выходной величины объекта,

где A, B, C, D - матрицы, Х - матрица внутренних переменных, U - матрица входных переменных (в данном случае ЭДС).

Матрицы будут иметь вид:



Получаем математическую модель в пространстве состояний:

2. ПОСТРОЕНИЕ СИГНАЛЬНОГО ГРАФА СИСТЕМЫ И СТРУКТУРНОЙ СХЕМЫ

Размещено на

Размещено на

Рис. 3. Граф системы

Рис. 4. Структурная схема системы

3. НАХОЖДЕНИЕ ПЕРЕДАТОЧНОЙ ФУНКЦИИ СИСТЕМЫ ПО ФОРМУЛЕ МЕЙСОНА

Найдем передаточную функцию объекта управления непосредственно из графа системы по формуле Мейсона, которая имеет вид:

; ,

где k - количество прямых возможных путей от входа к выходу;

? - определитель графа;

P_k - коэффициент передачи k-го пути от входа к выходу;

?_k - определитель всех касающихся контуров при удалении k-го пути;

- сумма коэффициентов передачи всех отдельных контуров;

- сумма всех возможных произведений из двух не касающихся контуров;

- сумма всех возможных комбинаций из трех не касающихся контуров.

Запишем уравнения всех прямых путей от входа к выходу:

Определим все замкнутые контуры и запишем их уравнения:



Запишем выражение для определителя системы:

Запишем определители путей:

Запишем и посчитаем выражение передаточной функции, подставив все найденные величины:

4. АНАЛИЗ УСТОЙЧИВОСТИ СИСТЕМЫ ПО КРИТЕРИЮ ЛЯПУНОВА

Проведём оценку устойчивости САУ по критерию Ляпунова, для чего найдем корни характеристического уравнения:

Так как все 4 корня характеристического уравнения находятся в левой части, то согласно критерию устойчивости по Ляпунову система является устойчивой.

5. ОПРЕДЕЛЕНИЕ ПРЯМЫХ ОЦЕНОК КАЧЕСТВА СИСТЕМЫ

Определение временных и частотных характеристик по передаточной функции.

Найдем переходную функцию, она равна обратному преобразованию Лапласа от , т.е.

Размещено на

Размещено на

Рис. 5. График переходного процесса

Найдем весовую функцию системы и построим ее график.

Размещено на

Размещено на

Рис. 6. График весовой функции

Определим частотные характеристики системы. Приведем передаточную функцию системы к частотной форме записи, заменив оператор Лапласа р на jщ:

Амплитудо-частотная характеристика находится по формуле:



Рис. 7. График амплитудо-частотной характеристики

Найдем фазо-частотную характеристику и построим ее график.

Рис. 8. График фазо-частотной характеристики

По переходной функции определим прямые оценки качества системы:

Перерегулирование

По графику видно, что система является неустойчивой, поэтому прямые оценки качества провести в этом случае невозможно.

По графику АЧХ определим косвенные оценки качества ОУ:

Максимальное значение амплитуды:

A_max(щ)=0,5

Резонансная частота:

щ_р=10 , Гц

Амплитуда сигнала при нулевой частоте

А(0)=0

Показатель колебательности:

М=A_max/A(0)=0

Полоса пропускания

6. СИНТЕЗ ФОРМИРУЮЩЕГО ФИЛЬТРА

По корреляционной функции K_Х(ф) определим спектральную плотность S_Х(щ) для белого шума (S_V(щ)=const), который подается на вход формирующего фильтра.

Корреляционная функция задана выражением:

Найдем спектральную плотность случайного сигнала и построим ее график:



Размещено на

Размещено на

Рис. 9. Спектральная плотность

По заданным статическим характеристикам S_x и S_V определим передаточную функцию формирующего фильтра Ш(р), для чего воспользуемся уравнением:

В рамках данной курсовой работы спектральная плотность S_V(щ) принимается равной 1. Получаем, что квадрат модуля частотной характеристики определяется следующим образом:

Найдем корни числителя и знаменателя .

Найдем корни числителя:

Определим корни знаменателя:

Строим корни на комплексной плоскости:

Размещено на

Размещено на

Рис. 10. Корни на комплексной плоскости

Из корней верхней полуплоскости формируем выражение для Ш(jщ):

Так как сомножитель знаменателя щ-64j образуется из решения уравнения 64+jw, то его можно заменить непосредственно этим уравнением:



Запишем передаточную функцию формирующего фильтра, заменив jщ на р:

сигнальный передаточный ляпунов фильтр

Представим ОУ в виде:

Размещено на

Размещено на

Рисунок 11. Эквивалентная схема с фильтром

Для этой системы определим общую передаточную функцию:

Определим устойчивость полученной системы по критерию Раусса для эквивалентной схемы.

Характеристическое уравнение имеет вид:

Вычисляем третью строку таблицы Раусса:

R₃ = a₀/a₁=1/205=0,00488

C₃₁ = a₂ - R₃a₃=17988 - 0,00488·672923= 14704,13576

C₃₂ = a₄ - R₃a₅=6654506 - 0,00488·16868944=6572185,55328

C₃₃ = a₆ - R₃a₇=4814540 - 0,00488·444344= 4812371,60128

Вычисляем четвёртую строку таблицы Раусса:

R₄ = a₁/С₃₁=205/14704,13576=0,01394

C₄₁ = a₃ - R₄C₃₂=672923- 0,01394·6572185,55328= 581306,73339

C₄₂ = a₅ - R₄C₃₃=16868944- 0,01394·481237...