63713 (Управление динамической системой)

Федеральное агентство по образованию

Государственное образовательное учреждение высшего профессионального образования

Сибирский государственный технологический университет

Факультет автоматизации и информационных технологий

Кафедра системотехники

Курсовая работа

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

Задачи

1. Подобрать аналитические выражения для функций M_g, M_c.

2. Найти равновесное состояние системы.

3. Численно решить систему уравнений моментов и управляющего устройства при начальных условиях и полученных выражениях , . Решение вести до установления значений и . Проверить совпадение и . Построить графики.

4. Линеаризовать уравнение моментов в окрестности точки равновесия, численно рассчитать линеаризованную систему. Построить графики.

5. Замкнуть систему. Представить разомкнутую систему в векторно-матричной форме для непрерывного и дискретного времени. Представить замкнутую систему в векторно-матричной форме для непрерывного и дискретного времени.

Оценить управляемость системы. Составить характеристическое уравнение системы. На основе критерия Рауса - Гурвица определить

1. значение коэффициента k = k₀, соответствующее пределу устойчивости линеаризованной системы

2. Найти корни характеристического уравнения системы и исследовать перемещение корней на комплексной плоскости при варьировании коэффициента усиления k. Построить траекторию движения корней.

3. Построить переходный процесс в системе. Уравнение решить аналитически, выполнив спектральное разложение матрицы А и использовав собственные числа и собственные вектора матрицы А.

4. Используя преобразование Лапласа, получить передаточные функции системы по каналу ux1 . На основе z-преобразования аналогичным образом получить дискретную передаточную функцию системы.

5. Выписать выражения для амплитудно-фазовой, амплитудной, фазовой, вещественной и мнимой частотных характеристик для системы. Для значения a = 0.9 построить годограф АФЧХ и графики характеристик A(), (), Re(), Im().

6. Оценить устойчивость системы по критерию Найквиста (по АФЧХ системы) и устойчивость системы по критерию Михайлова, построив для этой цели годограф Михайлова. Определить запас устойчивости системы.

Реферат

В пояснительной записке содержится 22 листа текстовой части, 19 рисунков и 1 источник данных.

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

Целью работы является исследование поведения управляемой динамической системы, описанной системой дифференциальных уравнений. На основе исходных данных находятся равновесное состояние системы, вид линеаризованной системы; исследуется соответствующая ей замкнутая система и описывается переходный процесс в этой системе; определяются частотные характеристики системы и устойчивость.

Содержание

Введение

1. Исходные данные

2. Нахождение аналитического вида функций Mc(ω) и Mg(ω,μ)

3. Нахождение равновесного состояния системы

4. Численное нахождение функций ω(t) и μ(t) равновесного состояния

5. Линеаризация и численное решение разомкнутой системы

6. Замкнутая система

7. Оценка управляемости системы

8. Оценка устойчивости системы

9 Построение переходного процесса

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

11. Амплитудная, фазовая, вещественная, мнимая и амплитудно-фазовая частотные характеристики

12. Оценка устойчивости системы по критерию Найквиста, по критерию Михайлова

Заключение

Библиографический список

Введение

Теория управления – это наука, изучающая процессы в системах управления с информационной точки зрения, обычно абстрагируясь от физической природы объектов и управляющих устройств. Процессы в автоматических системах управления изучает теория автоматического управления.

Важнейшие принципы построения систем автоматического управления:

• принцип обратной связи;

• принцип оптимальности;

• принцип адаптивности;

• принцип робастности.

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

1 Исходные данные

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

а) Уравнение моментов:

(1)

б) Уравнение управляющего устройства:

t - время, сек; J - момент инерции движущихся частей, приведенный к валу двигателя, кг*м / сек²; - угловая скорость двигателя, 1/сек; M_g, M_c - момент движущих сил и сил сопротивления, кг*м; - управляющее воздействие; u - задающее воздействие; , - параметры управляющего устройства

Функции M_g, M_c заданы таблицами 1 и 2, численные значения коэффициентов определены в таблице 3

Таблица 1 – Зависимость M_g от и

	0.0	0.2	0.4	0.6	0.8	1.0
0.00	53.50	55.62	57.54	59.86	61.98	64.10
10.60	46.28	48.63	50.99	53.35	55.71	58.06
21.20	36.48	39.08	41.67	44.27	46.86	49.46
31.80	24.11	26.95	29.78	32.61	35.44	38.27
42.40	9.17	12.24	15.31	18.38	21.45	24.52
53.00	0.00	0.00	0.00	1.58	4.89	8.19

Таблица 2 – Зависимость M_с от

	0.00	10.60	21.20	31.80	42.40	53.00
Мс	10.70	13.50	20.22	30.84	45.37	63.82

Таблица 3 – Значение параметров системы

J	m	R1	R2	C
0.06	10.03	19.40	1.03	1.03

Начальные условия: t = 0; = 0; = 0; ; u = 0.5.(3)

2 Нахождение аналитического вида функций M_c(ω) и M_g(ω,μ)

динамическая система (1)

Аналитический вид функции момента движущих сил M_c(ω) находится методом наименьших квадратов:

Аналитический вид функции момента движущих сил M_g(ω,μ) находится методом наименьших квадратов. Сначала по столбцам при различных μ вычисляется матрица ABC зависимости M_g(ω,μ) от μ. Первый столбец матрицы ABC вычисляется при μ=0 из системы:

Остальные столбцы заполняются аналогично при μ равном 0, 0.2, 0.4, 0.6, 0.8, 1.

Матрица ABC выражает зависимость функции M_g(ω,μ) от ω при различных μ. При этом функция M_g(ω,μ) имеет вид:

Строки матрицы выражают зависимость слагаемых (А(μ), В(μ) и С(μ)) функции M_g(ω,μ) от μ, соответственно 1-ая строка А(μ), 2-ая строка В(μ), 3-я строка С(μ). А(μ), В(μ) и С(μ) имеют вид:

Коэффициенты при μ вычисляются методом наименьших квадратов из матрицы ABC по строкам. Так для А(μ) по первой строке матрицы ABC из системы

Аналогично находим аналитический вид В(μ) и С(μ). Получаем:

3 Нахождение равновесного состояния системы

Найдем равновесное состояние системы при следующих условиях . Подставим эти условия в систему (1), получим систему вида:

Решая систему численно, получаем равновесное состояние системы при ω₀=34.54 и μ₀=0.5. Построим графики M_c(ω) и M_g(ω,μ) при разных μ₀. На рисунке 1 жирными сплошными линиями отмечены графики M_c(ω) и M_g(ω,μ) при μ₀=0.5

Рисунок 1 – Графики функций M_c(ω) и M_g(ω,μ)

4 Численное нахождение функций ω(t) и μ(t) равновесного состояния

Для того чтобы из системы (1) найти функции ω(t) и μ(t), необходимо понизить степень системы, то есть избавиться от производных второго порядка. Для этого введем функцию Z(t)= μ'(t), получим систему вида:

(2)

Решая систему численно, получаем табличные значения ω(t) и μ(t), по которым строим графики ω(t) (рисунок 2) и μ(t) (рисунок 3). По графикам хорошо видно, что ω(t) и μ(t) стремятся к равновесным значениям ω0=31.948 и μ0=0.5, ω(t)→ 31.948, μ(t) →0.5, что соответствует вычислениям.

Рисунок 2 – График функции ω(t)

Рисунок 3 – График функции μ(t)

5 Линеаризация и численное решение разомкнутой системы

Линеаризуем систему (2) в окрестности точки равновесия. Для этого выведем систему из равновесия, придав u, μ, ω малые приращения ∆u, ∆μ, ∆ω→0. Соответственно придается приращение Z, ∆Z→0.

(3)

Теперь разложим функции M_c(ω) и M_g(ω,μ) в ряды Тейлора по формулам: