97040

СИНТЕЗ И АНАЛИЗ АЛГОРИТМОВ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ МЕТОДАМИ ТЕОРИИ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

Курсовая

Информатика, кибернетика и программирование

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

Русский

2015-10-13

601.5 KB

1 чел.

Министерство образования и науки Российской Федерации

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

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

им Гагарина Ю.А.

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

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

по дисциплине «Основы теории оптимального управления» на тему:

СИНТЕЗ И АНАЛИЗ

АЛГОРИТМОВ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ

МЕТОДАМИ ТЕОРИИ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

Выполнил

студент гр.УИТ-31

Маслов Ю.О.

Проверила

Торгашова О.Ю.

Саратов 2012

Задание на курсовую работу

Задача 1

Оптимальное управление в RL-цепи

Для электрической схемы содержащей источник питания e(t), активные сопротивления r, R  и индуктивность L

необходимо :

а) определить оптимальный закон изменения напряжения источника питания e(t), приводящий к изменению тока через сопротивление R и индуктивность L в схеме от заданного начального значения  i(t0)=i0  до заданного конечного значения i(t1)=i1 , такой, чтобы суммарная энергия активных потерь (затрачиваемая на нагрев) при этом изменении была минимальной;

б) определить оптимальный закон изменения тока i(t) , соответствующий оптимальному закону изменения e(t);

в) вычислить энергию активных потерь в схеме при оптимальном режиме изменения напряжения и тока ( e(t), i(t) ) и сравнить ее с энергией активных потерь, затрачиваемой на нагрев при линейном изменении тока в схеме от начального до конечного значения;

г) построить графики оптимальных и линейных изменений ЭДС и токов.

Значения параметров элементов схемы в зависимости от варианта задания приведены в таблице 1.

                                                                            Таблица 1

Номер     варианта

r

[Ом]

R

[Ом]

L

[Г ]

i(t0)

[А]

i(t1)

[А]

18

7,5

60

0,023

0,3

0,9

 Полагать t0=0 , t1=10-3 c .

Задача 2

Рис. 4

Для электрической схемы, содержащей источники ЭДС e(t) и источник тока J(t), активные сопротивления R1, R2, R3, индуктивность L и емкость С (рис. 4), необходимо:

а) определить оптимальные законы изменения напряжения источника питания e0(t) и источника тока J0(t), приводящих к изменению тока через индуктивность L в схеме от заданного начального значения iL0  до заданного конечного значения iL1, и изменению напряжения на обкладках конденсатора С в схеме от заданного начального значения uC0 до заданного конечного значения uC1 так, чтобы суммарная энергия активных потерь (затрачиваемая на нагрев) при этом изменении была минимальной;

б) определить оптимальные законы изменения тока i0L(t) и напряжения на обкладках конденсатора u0C(t), соответствующие оптимальным законам изменения e0(t), J0(t);

в) вычислить энергию активных потерь в схеме при оптимальном режиме изменения напряжений источников питания, тока и напряжения на обкладках конденсатора (e0(t), J0(t), i0L(t), u0C(t)) и сравнить ее с энергией активных потерь, затрачиваемой на нагрев при линейных изменениях тока iЛL(t) и напряжения на обкладках конденсатора uЛC(t) в схеме от начальных значений до конечных;

г) построить графики оптимальных и линейных изменений ЭДС, источника тока, тока и напряжения на обкладках конденсатора.

Значения параметров элементов схемы:

зад.

R1,

Ом

R2,

Ом

R3,

Ом

L,

Гн

С,

мкФ

iL0,

А

iL1,

А

uC0,

В

uC1,

В

18

60

80

75

0,023

5

0,3

0,9

45

90

Полагать   t0=0 , t1=10-3 c .

  1.  Задача 1

Оптимальное управление в RL-цепи

Описание объекта управления

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

и имеет вид дифференциального уравнения

,               (1)

где  x(t)=i(t) , u(t)=e(t) ,p, b – числа, равные p = – R/L, b = 1/L  

              

Конструирование функционала – критерия оптимальности.

Критерий оптимальности – квадратичный функционал

                       

где  - симметричная, неотрицательно-определенная матрица чисел, размерами ;  - симметричная, положительно-определенная матрица чисел размерами .

В данном случае вместо матриц используются числа, поэтому критерий оптимальности будет иметь другой вид:

,  (2)

Это выражение представляет собой суммарную энергию активных потерь в схеме за время  t1–t0.         

Запишем выражение для активной мощности потерь на сопротивлениях  r и R :

 или    ,

Таким образом,  q = R=60,  m = 1/r =1/7,5 =0,133,  n = 0.

Формулировка задачи как вариационной задачи на условный экстремум.

Для этого необходимо рассматривать в качестве уравнения связей уравнение системы (1), а в качестве функционала – функционал (2).

Таким образом, получаем следующую вариационную задачу:

Определить функции x(t) и u(t) доставляющие экстремум функционалу  

,         

при граничных условиях

      ,                       

и при дополнительном условии (уравнении связи)

накладываемом на функции x(t), u(t) , в классе которых ищется экстремум.

Синтез оптимального алгоритма управления.

1.4.1. Получение уравнений вариационной задачи.

Введем вспомогательную функцию (функцию Лагранжа)

                      

где  - пока неизвестная функция называемая неопределенным множителем Лагранжа.

Рассматриваемая задача называется задачей Лагранжа.

 

(t) – неопределенный множитель Лагранжа, Φ – функция Лагранжа.

Запишем уравнения Эйлера для функции  (они называются уравнениями Эйлера – Лагранжа)

                 

Решим совместно уравнения Эйлера – Лагранжа  и уравнение связи.  Это система трех уравнений для определения трех неизвестных x(t), u(t), (t).

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

              (3)

                      (4)

                    (5)

 Здесь (3), (4) – уравнения Эйлера – Лагранжа. К этим уравнениям добавлено уравнение объекта  (5) (уравнение связи).

Отыскание решения уравнений вариационной задачи.

Уравнения (3) – (5) решаются в следующем порядке:

1)  Выразим u(t) из (4):

 

Затем подставим его в  (5).

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

,          (6)

с коэффициентами

a11 = pnb/2m = -2608,                

a12 = b2/2m =7103,  

a21 = 2q – n2/2m =120,      

a22 = nb/2m – p=2608.

Получим систему уравнений:

 

              

2) Запишем систему (6) в матричной форме

,        (7)       

где

        ,

        .

3) Запишем решение уравнения (7) в соответствии с формулой Коши в виде:

,                       (8)           

где

    – вектор начальных условий, а матричная экспонента определяется по формуле Лагранжа – Сильвестра

  ,

где 1 , 2 – собственные числа матрицы А. Е - единичная матрица.

Найдем собственные числа матрицы А из условия. Получим:

,

,

.

Из системы следует, что для нахождения  и  необходимо знать начальные значения  и .

(начальное положение объекта) задано,  неизвестно.

Для определения неизвестного начального значения множителя Лагранжа (t0), входящего в (8) необходимо:

а) запишем (8) для момента времени  t1

или

   ,  

   ,                  (9)

где  e11 , е12 , е21 , е22  – элементы матрицы (числа):

  

        

б) определим (t0) из первого уравнения системы (9)

Получили, что  

Решаем уравнение (7):

 

4) Решив уравнение (7), запишем выражения для оптимальной траектории и оптимального управления:

- оптимальная траектория

- оптимальное управление

Анализ процессов в системе.

Анализ процессов при оптимальном режиме

Анализ процессов при оптимальном режиме построим графики x(t) , u(t)  на  интервале  t[t0,t1]. Этот интервал разбивается на 10 частей и вычисляются значения x(t) и u(t) в этих точках.

      t

x(t)

u(t)

0

0,15

0,287

0,0001

0,093745

1,038391

0,0002

0,063278

2,075415

0,0003

0,050216

3,683328

0,0004

0,050967

6,304424

0,0005

0,065738

10,65969

0,0006

0,098591

17,94715

0,0007

0,158565

30,17137

0,0008

0,262155

50,69491

0,0009

0,437857

85,16323

0,001

0,734001

143,0576


Анализ процессов при линейном изменении тока i(t)

Полагая, что ток изменяется линейно от заданного начального состояния до заданного конечного состояния

xЛ(t) = kt + d        ( iЛ(t) = kt + d),

(величины k , d  найдем из условия прохождения iЛ(t) и uЛ(t) через заданные начальное и конечное значения)

xЛ(0) = d=0,15  ,      xЛ(0.001) = 0.001k + 0,15 = 0,75   

                                                   k=600

                                        xЛ(t) = 600t + 0,15  

запишем на основе (1)    

выражение для закона управления uЛ(t) , обеспечивающее такое линейное изменение

.

По полученным данным построим графики процессов xЛ(t), uЛ(t).

t

xл(t)

uл(t)

0

0,15

25,65

0,0001

0,21

31,35

0,0002

0,27

37,05

0,0003

0,33

42,75

0,0004

0,39

48,45

0,0005

0,45

54,15

0,0006

0,51

59,85

0,0007

0,57

65,55

0,0008

0,63

71,25

0,0009

0,69

76,95

0,001

0,75

82,65


Сравнительная оценка процессов в схеме при оптимальном и линейном режимах.

Вычислим энергию активных потерь при оптимальном режиме, подставив в (2)  x(t) и u(t).

Вычислим энергию активных потерь при линейном режиме путем подстановки в (2) xЛ(t)  и  uЛ(t).

Сравнивая полученные величины, делаем вывод, что суммарная энергия активных потерь (затрачиваемая на нагрев) при линейном режиме больше, что говорит о целесообразности работы схемы в оптимальном режиме.

                         

Задача 2

Описание объекта управления

Математическая модель объекта получается на основе законов Кирхгофа и имеет вид линейных дифференциальных уравнений, векторно-матричный вид которых

   ,                       

где    матрицы чисел размерами  и  соответственно.

Начальное положение (состояние) объекта:.               

Конечное положение (состояние) объекта: .                 

Уравнения электрической схемы можно получить в виде:

,   (1).                                 (2)

Эти уравнения записаны в форме Коши.   

Введем обозначения:

,.

Тогда выражения (1), (2) можно записать в виде: 

     (3)

Уравнения для выходных переменных

,  .           (4)

Объединяя (3) и (4), получаем математическую модель электрической схемы:

  (5)

Где при  Ф,  Гн, [1/Ом],  Ом соответствующие матрицы получают следующие численные значения:

.

Начальное положение объекта:

Конечное положение объекта:

Конструирование функционала (критерия оптимальности).

Находится выражение для суммарной энергии активных потерь в схеме за время  t1–t0 , которое в общем виде можно представить в виде квадратичного функционала

  (6)

где  – симметричная неотрицательно-определенная матрица чисел размерами ;  – симметричная положительно-определенная матрица чисел размерами .

    

Таким образом,

, N=0

Формулировка задачи как вариационной задачи на условный экстремум.

Для этого необходимо рассматривать в качестве уравнения связей уравнение системы (1), а в качестве функционала – функционал (2).

Таким образом, получаем следующую вариационную задачу:

Определить функции x(t) и u(t) доставляющие экстремум функционалу  

,         

при граничных условиях

,

и при дополнительных условиях (уравнениях связи)

накладываемом на функции x(t), u(t) , в классе которых ищется экстремум.

Синтез алгоритма оптимального управления – решение вариационной задачи на условный экстремум.

Составляется функция Лагранжа:

.        

Записываются уравнения Эйлера-Лагранжа:

            (7)

Вычисляются составляющие соотношений (5.6):

;  

;

                 

.

Тогда уравнения Эйлера-Лагранжа принимают вид:

                            (8)

К последним уравнениям добавляется уравнение связи (уравнение объекта) получается следующая система уравнений:

                          (9)

   (10)

Найдем U:

и подставим в первые два уравнения (10):

  (11)

Введем вектор . Тогда систему (10) с учетом (11) можно представить в виде:

  (12)

где Р – блочная матрица, имеющая вид:

  (13)

Решение (12) в соответствии с формулой Коши:

  (14)

Вычислив , ее можно представить следующим образом:

  (15)

где  - функциональные матрицы размеров .

Тогда из выражений (14), (15):

  (16)

Для определения начального значения множителя Лагранжа запишем следующие соотношения:

  (17)

Из первого уравнения (17) можно определить начальное условие множителя Лагранжа:

  (18)

Теперь можно записать из выражений (16):

  (19)

Оптимальное управление определится выражением:

Таким образом найдены соотношения для оптимальной траектории и оптимального управления:

Анализ процессов в системе

Анализ процессов при оптимальном режиме

Анализ процессов при оптимальном режиме построим графики X(t) , U(t)  на  интервале  t[t0,t1]. Этот интервал разбивается на 1000 частей и вычисляются значения x(t) и u(t) в этих точках.

Рис. 1. График оптимального изменения напряжения uC(t)

Рис. 2. График оптимального изменения тока i2(t)


Рис. 3. График оптимального изменения  источника тока I(t)

Рис. 4. График оптимального изменения ЭДС e1(t)

Анализ процессов при линейном  режиме

Полагая, что ток и напряжения изменяются линейно от заданного начального состояния до заданного конечного состояния

xЛ(t) = kt + d        

(величины k , d  найдем из условия прохождения iЛ(t) и uЛ(t) через заданные начальное и конечное значения)

xЛ1(0.001) = 0.001k + 45 = 90,       k1=45000,          xЛ1(t) = 45000t + 45

xЛ2(0.001) = 0.001k + 0.3 = 0.9,     k2=600,              xЛ2(t) = 600t + 0.3

Запишем на основе (11)  

По полученным данным построим графики процессов XЛ(t), UЛ(t).

Рис. 6. График линейного изменения напряжения uC(t)

Рис. 7. График линейного изменения тока iL(t)

Рис. 8. График линейного изменения источника тока I(t)

Рис. 9. График линейного изменения ЭДС e(t)

Сравнительная оценка процессов в схеме при оптимальном и линейном режимах

Вычислим энергию активных потерь при оптимальном режиме

J=0.03511

Вычислим энергию активных потерь при линейном режиме путем подстановки в (2) xЛ(t)  и  uЛ(t).

Сравнивая полученные величины, делаем вывод, что суммарная энергия активных потерь (затрачиваемая на нагрев) при линейном режиме больше, что говорит о целесообразности работы схемы в оптимальном режиме.


Заключение

В процессе выполнения курсовой работы были  разработаны алгоритмы изменения режима работы электрической схемы, содержащей активные и реактивные элементы, которые обеспечивают минимизацию энергии активных потерь при переходе от одного режима работы схемы к другому. Также был определен вид закона изменения ЭДС источника питания (управляющего воздействия) и была проанализирована работа схемы при действии этой ЭДС.

r

e(t)

L

R

i(t)


 

А также другие работы, которые могут Вас заинтересовать

27161. Варианты формата CD 221 KB
  Однако значительная информационная ёмкость нового носителя 740 Мбайт навела специалистов на мысль использовать его в качестве элемента постоянной памяти для хранения архивных данных. Каждый кадр как уже описывалось в главе 3 содержит в себе 24 исходных информационных символа байта. В формате CDROM эти 24 символа являются обезличенными и могут нести в себе какую угодно информацию лишь бы она была преобразована в двоичную форму и организована в байты. Изза наличия этой избыточности диск CDROM имеет меньшую информационную ёмкость до...
27162. Digital Versatile Disc (DVD) 187 KB
  Digital Versatile Disc DVD 12. История появления DVD К концу 1994 года в технической прессе стали появляться сообщения о том что известный тандем SONY PHILIPS подаривший миру технологию CD готов представить на суд потребителю еще более совершенный носитель идеально подходящий для записи информации практически любого характера. В процессе работы над новым носителем несколько раз менялось его название отражая основные намерения разработчиков на том или ином этапе: MMCD MultiMediaCD; HDDVD High Density Digital Video Disc; HDCD...
27163. Система магнитооптической записи звука «Минидиск» 224.5 KB
  Звуковые характеристики Число каналов Детонации 2 или 1 отсутствуют Формат данных Частота дискретизации кГц Кодирование сжатие данных Модуляция канальный код Система защиты от ошибок 441 ATRAC EFM 814 ACIRC Оптические характеристики Длина волны излучения лазера нм Числовая апертура объектива Мощность излучения лазера при записи мВт Метод записи 780 045 25 50 главный пучок Модуляция магнитного поля 11. В общем случае магнитооптический эффект это изменение оптических свойств вещества в зависимости от его...
27164. СИСТЕМА ЦИФРОВОЙ МАГНИТНОЙ ЗАПИСИ ЗВУКА R-DAT 182.5 KB
  Описание формата RDAT Rotary Head Digital Audio Tape Recorder это система цифровой звукозаписи на магнитную ленту шириной 381 мм равную ширине ленты в обычной аналоговой компакткассете с помощью вращающихся головок. В отличие от формата CD здесь предусмотрено не только воспроизведение программ но и возможность их записи с высоким качеством. Режим I предназначен для записи и воспроизведения программ с частотой дискретизации 48 кГц при 16 разрядном линейном квантовании.
27165. Система защиты от ошибок 494.5 KB
  В магнитофонах формата RDAT так же как и в формате CD для борьбы с искажениями используется комплексная система защиты от ошибок включающая в себя два кода РидаСоломона С1[3228] и С2[3226] и двунаправленный способ перемежения данных. Проверочная матрица НР кода С1 показана на рисунке 9 а расположение символов внутри кодового слова задано векторстолбцом VP показанным на рисунке 10. Порождающий многочлен GPX этого кода имеет вид: или Поскольку кодовое слово кода С1 содержит четыре проверочных символа то этот код способен...
27166. КАНАЛЬНОЕ КОДИРОВАНИЕ 108.5 KB
  Считается что если после такого преобразования число канальных бит высокого уровня равно числу канальных бит низкого уровня то постоянная составляющая всей комбинации будет равна нулю DSV = 0. Оставшиеся 103 комбинации пришлось выбрать из тех которые имеют ненулевое значение DSV. Однако вместо одной 10разрядной комбинации каждому из этих 103 8разрядных символов поставлены в соответствие две отличающиеся друг от друга только знаком DSV. Одна из них имеет значение DSV = 2 другая DSV = 2.
27167. СЛУЖЕБНАЯ ИНФОРМАЦИЯ 172 KB
  Эта информация кроме специально отведенной для нее зоны данных субкода размещается еще и в символах кода идентификации ID которые имеют место в начале каждого блока как в зоне данных ИКМ W1 и W2 так и в зоне данных субкода SW1 и SW2. Служебная информация размещаемая в зоне данных субкода может переписываться и дописываться независимо от музыкальной информации записанной в ИКМзоне. Изменить ее не изменяя основных данных невозможно. Это та информация которая записана в символах W1 и W2 кода идентификации ID зоны ИКМданных.
27168. DAT-кассета 167.5 KB
  Так же как у видеокассеты рабочая поверхность ленты защищена шарнирной крышкой и доступ к ней обеспечивается только при откидывании этой крышки. Эта панель к тому же застопоривает шарнирную крышку препятствуя доступу к магнитному слою ленты. Таким образом DATкассета хорошо закрыта со всех сторон и магнитный слой находящейся нее ленты надежно защищен от пыли царапин отпечатков пальцев и других внешних воздействий. В дальнейшем при опускании кассеты внутрь транспортного механизма специальный рычаг ЛПМ поднимает шарнирную крышку...
27169. Система автотрекинга 174.5 KB
  Необходимые условия для вхождения в синхронизм обеспечиваются за счет того что перед началом каждого очередного массива данных на дорожке зона ИКМданных или зона данных субкода размещаются вспомогательные последовательности содержащие большое число переходов IBG вводный и выводной сигналы. Уровни сигналов ошибки при точном а и неточном б следовании головки по дорожке IBG4 F4 F4 F3 F1 F4 F3 IBG3 IBG2 В Нечетный кадр F4 F1 IBG1 1 1 1 2 блока 1 2 2 ATF2 ATF1 Направление движения ленты Направление движения головки 2 2 1 2 1 1 1 IBG1...