23030

Проблеми моделювання динаміки систем з розподіленими параметрами

Лекция

Экономическая теория и математическое моделирование

4 і модель ця адекватно описує динаміку фізикотехнічного обєкту процесу то можна ставити і розвязувати: Прямі задачі динаміки визначення векторфункції стану ys при заданих зовнішньодинамічних факторах ; Обернені задачі динаміки визначення векторфункцій які б згідно певного критерію дозволяли отримувати задану картину змін векторфункції ys або наближатися до неї.4 побудовані апробовані практикою а відповідні математичні теорії дозволяють розвязувати як прямі так і обернені задачі динаміки таких систем....

Русский

2013-08-04

1.64 MB

3 чел.

8

Курс лекцій по моделюванню динаміки систем з розподіленими параметрами

Лекція 1. Проблеми моделювання динаміки систем з

                 розподіленими параметрами

1.1. Диференціальна модель динаміки. Будемо розглядати систему (процес) стан якої (якого) в просторово-часовій області , обмеженій контуром Г, визначається вектор-функцією Для більшості фізико-технічних систем (процесів) після вивчення і формалізації їх суті вдається побудувати залежність вектор-функції y(s) від вектор-функції  зовнішньо-динамічних факторів, які в області  діють на систему. В загальному випадку залежність ця записується системою диференціальних рівнянь, яку в подальшому будемо записувати у вигляді:

, (1.1)
де — вектор частинних похідних по просторових координатах  та похідна по часу
t, а

  (1.2)
задані матричні диференціальні оператори.

Опис динаміки системи співвідношенням (1.1) буде не повним, якщо не врахувати початковий стан системи та вплив оточуючого її середовища (довкілля), який проявляється через граничний стан контуру Γ області . Вплив початково-крайових збурень на стан системи в подальшому будемо описувати наступними співвідношеннями:

                                                (1.3)

  (1.4)
де  – задані матричні оператори, а  та  – задані вектори та векторні функції, розмірність яких, як правило, узгоджується з розмірністю та порядком матричних диференціальних операторів  та
.

Якщо модель системи описана співвідношеннями (1.1) - (1.4), і модель ця адекватно описує динаміку фізико-технічного об’єкту (процесу), то можна ставити і розв’язувати:

  1.  Прямі задачі динаміки – визначення вектор-функції стану y(s) при заданих зовнішньо-динамічних факторах ,  ,  ;
  2.  Обернені задачі динаміки - визначення вектор-функцій  ,  ,  , які б згідно певного критерію дозволяли отримувати задану картину змін вектор-функції y(s) (або наближатися до неї).

Для більшості класичних систем моделі вигляду (1.1) - (1.4) побудовані, апробовані практикою, а відповідні математичні теорії дозволяють розв’язувати як прямі, так і обернені задачі динаміки таких систем.

Проблеми виникають там, де в силу складності системи (процесу), модель (1.1) - (1.4) не повна (нема узгодженості між розмірністю і порядком системи (1.1) та розмірністю і кількістю співвідношень (1.3), (1.4)), не адекватно описує динаміку розглядуваної системи (процесу), або не зовсім вивчене і формалізоване у вигляді (1.3), (1.4) середовище, в якому функціонує система (протікає процес). Розв’язок проблеми – в настройці моделі (1.1) та умов (1.3), (1.4) на реальну динамічну картину, яка спостерігається для системи (процесу). Останнє виконується двояко:

  1.  Шляхом ідентифікації параметрів співвідношень (1.1), (1.3), (1.4) на експериментальних даних, якщо система (процес) формалізовані моделлю (1.1) - (1.4);
  2.  Моделюванням зовнішньо-динамічної обстановки, якщо співвідношення (1.3), (1.4) важко формалізуються, або характеристики , , зовнішньо-динамічного оточення системи не доступні для вимірювання.

Іншими словами кажучи: практика роботи з моделями (1.1) - (1.4) вимагає уміння ставити та розв’язувати наступні задачі:

;                                                 (1.5)

(1.6)
де  –спостережуваний стан системи,  вектор параметрів моделі, від яких залежать матричні оператори
, , , а  та  – фіктивні зовнішньо-динамічні фактори, якими моделюється ефект впливу на вектор-стану y(s) системи початково-крайових умов (1.3), (1.4).

Методи постановки та розв’язання задач (1.5), (1.6) для різних режимів і умов функціонування системи (1.1) - (1.4) і будуть вивчатися нижче.

1.2. Інтегральна модель динаміки. Практика роботи з системами, що функціонують в певній просторово-часовій області, показала доцільність вибору моделі динаміки у формі (1.1) - (1.4) і довела дієвість математичних методів розв’язання прямих та обернених задач динаміки таких систем. Задачі ж (1.5) та (1.6) для систем, динаміка яких описується моделлю (1.1) - (1.4), розв’язуються важко.

Зважаючи на це, для розв’язання задач (1.5), (1.6) виберемо тут іншу — інтегральну — модель, яку запишемо у вигляді:

, (1.7)
де  – відома матрична функція.

Явна аналітична залежність y(s) від u(s) спрощує роботу з моделлю, а отже і алгоритми розв’язання задач (1.5), (1.6).

Зауважимо відразу, що співвідношення (1.7) на відміну від співвідношень (1.1) - (1.4) з фізико-технічних міркувань побудувати не просто. Не просто модель (1.7) отримати і із моделі (1.1) - (1.4) – для цього до вигляду (1.7) необхідно привести аналітичний розв’язок початково-крайової задачі (1.1) - (1.4). Як показує практика, це можливо тільки для деяких нескладних моделей і для досить простих умов їх функціонування.

Розглядувані ж нижче методи дозволять в якійсь мірі розв’язати і цю проблему. Для початку будуть запропоновані методи побудови моделі (1.7) для систем, які функціонують в необмеженій просторово-часовій області.

. (1.8)
Побудована для цього випадку залежність

(1.9)
вектор-функції стану системи від зовнішньо-динамічних збурень
u(s) буде не чим іншим, як розв’язком рівняння (1.1) без врахування початково-крайових умов(1.3), (1.4). Матрична функція  в цьому розв’язку буде функцією Гріна такою, що

,                          (1.10)

де , а (s-s’) так звана –функція, або функція одиничного джерела.

Враховуючи, що

бачимо, що дійсно  формі (1.9) задовольнятиме рівнянню (1.1) в області .

Функцію стану y(s) системи (1.1) - (1.4), що функціонує в замкненій просторово-часовій області , або (як ми назвали тут) інтегральну модель динаміки системи побудуємо у вигляді

                                          (1.11)

де  та  – складові вектор-функції y(s) системи, якими визначається внесок початково-крайових умов (1.3), (1.4) в загальну картину динаміки системи. При цьому вектор-функції  та  виберемо у вигляді

                                     (1.12)

                                  (1.13)

де , а  та  – фіктивні зовнішньо-динамічні збурення, які (як видно з (1.12) - (1.13)) діють за межами часової та просторової областей відповідно і визначаються з умови, щоб з врахуванням (1.11) виконувалися співвідношення (1.3), (1.4) при заданих . Співвідношення ж (1.1) з представленням вектор-функції стану y(s) у вигляді (1.10) будуть виконуватися. Останнє легко перевіряється з врахуванням властивостей (1.10) функції, -функції  та особливостей визначення функцій  та  в (1.12) та (1.13).

1.3. Проблеми переходу від диференціальної форми моделі динаміки системи до інтегральної. Заміна класичної диференціальної форми (1.1) моделі динаміки систем з розподіленими параметрами інтегральною формою (1.9), (1.11) - (1.13), яка може бути зручнішою при розв’язанні задач (1.5), (1.6), викликає чималі труднощі. Заключаються вони в:

  1.  Побудові матричної функції Гріна  рівняння (1.1) для необмеженої просторово-часової області ;
  2.  Побудові вектор-функцій  та , якими б моделювалися початково-крайові умови (1.3), (1.4).

Математично перша проблема зводиться до розв’язання рівняння (1.10). Розв’язок другої визначатиметься наступною системою інтегральних рівнянь:

       (1.14)

                       (1.15)

При цьому мається на увазі, що функція зовнішньо-динамічних збурень u(s) відома.

Варіанти розв’язання цих двох проблем ми розглянемо в подальшому, однак зауважимо, що в деяких випадках функціонування системи (1.1) інтегральні співвідношення (1.14), (1.15) спрощуються. Це випадки, коли при розгляді динаміки системи (1.1) можна нехтувати початковим, або крайовим збуренням. Для цих двох випадків співвідношення (1.14), (1.15) спрощуються до наступних:

а) для динаміки системи без врахування крайових умов (в необмеженій просторовій області):

      (1.16)

б) для динаміки системи без врахування початкового стану(в необмеженій часовій області):

                         (1.17)

В подальшому ми розглянемо три варіанти обернення співвідношень (1.14) - (1.17):

- дискретизованих по s та ;

- дискретизованих по s;

- дискретизованих по .

Варіанти ці зводяться до побудови розв’язку наступних алгебраїчних та функціональних рівнянь відповідно:

                                                                          (1.18)

                                                          (1.19)

.                                                                 (1.20)

Випишемо вирази матриці С, векторів  та, а також матричних функцій та  для задачі (1.14), (1.15).

1.У рівнянні (1.18):

                             (1.21)

а

                     

       (1.22)

Тут  задані значення кількості точок дискретизації , ,  

2. У рівнянні (1.19):

                                           (1.23)

де  та  моделюючі вектор-функції, введені в (1.12), (1.13); та  визначені в (1.22), а

                  (1.24)

.

Тут  при i=1 та  при i=2.

3. У рівнянні (1.20)

                                           (1.25)
де значення вектор-функції стану для
; та визначені в (1.22), а

            (1.26)

Вирази (1.21), (1.23) та (1.25) для матриці С векторів , матричних функцій  спростяться для задач (1.16) та (1.17) і з врахуванням позначень (1.22), (1.24), (1.26) запишуться у вигляді:

1.Для задачі (1.16):

                                             (1.27)

;

2.Для задачі (1.17):

                                            (1.28)

;

Отже проблема переходу від диференціального запису моделі динаміки у формі (1.1)-(1.4) до інтегрального представлення у формі (1.9), (1.11)-(1.13) буде розв’язана, якщо:

  1.  буде побудована матрична функція Гріна  для необмеженої просторово-часової області ;
  2.  будуть методи побудови розв’язків (або наближень до них) рівнянь (1.18)-(1.20).

Ці два питання будуть детально вивчені нами в подальшому.

1.4. Проблеми ідентифікації параметрів моделі динаміки систем з розподіленими параметрами. Як відзначалося вище, не завжди вдається адекватно описати динаміку системи, особливо це стосується складних і важких для формалізації процесів. В таких випадках диференціальні оператори рівняння (1.1), а інколи і співвідношень (1.3) та (1.4), можуть залежати від невідомих параметрів. А це значить, що в загальному випадку вигляд цих операторів може бути таким:


Тут  – вектори параметрів моделі, які знаходяться з експерименту.

Враховуючи, що від  залежатиме функція Гріна , за умови, що вдалося обернути співвідношення (1.18)-(1.20) так, що

(1.29)
та позначаючи через  набір експериментальних даних, проблему ідентифікації параметрів  зведемо до розв’язання наступної задачі:

, (1.30)
де

, (1.31)
або

                                            (1.32)

Зрозуміло, що це задача непроста, тим паче, що і співвідношення (1.29) є дуже спрощеним представленням розв’язку рівнянь (1.18)-(1.20). Однак алгоритм розв’язання задачі після побудови обернень типу (1.29) ми запропонуємо.

1.5. Проблеми моделювання зовнішньо-динамычного оточення динамыки систем з розподіленими параметрами. Часто виникає необхідність, знаючи структуру моделі, відновити значення , , зовнішньо-динамічного збурення, початкових та крайових збурюючих факторів – усіх разом, чи деяких з них. Виходячи із згаданих вище експериментальних даних  та враховуючи, що вплив початково-крайових збурень на стан системи згідно (1.12), (1.13) моделюється функціями  та , проблему цю зводимо до розв’язання наступної задачі

.                                         (1.33)

Тут y(s) – вектор-функція стану системи, яка співвідношеннями (1.11)-(1.13) пов’язана з .

З використанням методики обернення співвідношень (1.18)-(1.20) і ця задача буде розв’язана нижче.

1.6. Задача оптимального розміщення спостерігачів та керувань. В процесі переходу від систем (1.14)-(1.17) до рівнянь (1.18)-(1.20) виконувалася дискретизація початково-крайових умов і керуючої вектор-функції u(s) точками  та  відповідно. Точками  дискретизувався і експеримент.

Вибір цих точок через розв’язки (1.29) рівнянь (1.18)-(1.20) впливатиме і на величину

,
тому має смисл постановка задачі оптимізації процедури вибору точок дискретизації шляхом мінімізації величини . Зрозуміло, що задача ця буде розв’язана після побудови алгоритму розв’язання задач (1.30) та (1.33).

Це і є коло тих проблем і задач, розв’язання яких буде розглянуте нами в подальшому.

7

                                   Стоян В.А.


 

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

84684. Виды аналитических и информационных документов, направляемых Посольством в Центр 14.29 KB
  В ней освещается широкоформатная тема например внешняя торговля страны пребывания за год или полгода вопросы состояния отношений страны пребывания с третьими странами участие страны пребывания в работе международных организаций и т. Оно касается крупных внутриполитических или внешнеполитических проблем страны пребывания. Грамотно составленный отчет является базой для внесения корректив во внешнеполитическую линию России в отношении данной страны. Составляются на видных государственных и общественных деятелей страны пребывания.
84685. Структура Посольства, распределение обязанностей, штатное расписание 13.41 KB
  Структура Посольства распределение обязанностей штатное расписание. Посол представляет РФ непосредственно руководит работой Посольства несёт персональную ответственность за выполнение возложенных на Посольство задач и осуществление им функций определяет в соответствии с нормативными актами МИДа России структуру Посольства распределяет должностные обязанности между его сотрудниками. Сотрудников посольства. Обеспечение безопасности посольства его сотрудников и членов их семей организация защиты государственной и иной охраняемой законом...
84686. Роль иностранного языка в работе дипломата 12.06 KB
  Роль иностранного языка в работе дипломата. Роль иностранного языка в работе чрезвычайно важна. Профессиональное владение иностранными языками может пригодиться дипломату во многих сферах его дипломатической деятельности. Знание в совершенстве языка страны пребывания значительно расширяет количество источников информации печатные издания на иностранном языке интернетресурсы телевидение радио и т.
84687. МГИМО: история, структура, кадровая база МИД России 15.87 KB
  МГИМО: история структура кадровая база МИД России. МГИМО является одним из старейших университетских центров страны по подготовке специалистов международного профиля. Первый набор в МГИМО составил 200 студентов. С 1946 года на учебу в МГИМО стали направляться студенты из зарубежных стран.
84688. Статус профессиональной дипслужбы во внешнеполитической системе, дипломаты и политическое руководство 13.63 KB
  Служба в аппарате внешнеполитического ведомства это особая разновидность федеральной государственной службы Российской Федерации. Это профессиональное осуществление целей и функций внешней политики России посредством исполнения государственных должностей федеральной государственной службы утвержденных в: центральном аппарате МИД РФ диппредставительствах и консульских учреждениях за рубежом представительствах при международных организациях представительствах МИД на территории РФ на отдельных государственных должностях госслужбы в...
84689. Порядок приема в МИД и требования, предъявляемые к выпускникам вузов, впервые поступающим на государственную службу 14.85 KB
  К дипломатической службе не относятся технические сотрудники. В России долгое время не было специального закона о дипломатической службе. Закон 79ФЗ о гражданской службе от 2004 года распространялся и на дипломатическую службы. В этих законах формулируются принципы государственной службы: равных доступ граждан к службе профессионализм и компетентность стабильность.
84690. Дипломатические ранги и порядок их присвоения 15.73 KB
  Присвоение дип рангов дип работникам производится в соответствии с установленными ФЗ Об основах гос службы РФ и иными нормативными правовыми актами РФ квалификационными требованиями к профессиональному образованию стажу и опыту работы по специальности знанию Конституции РФ Федеральных Законовнов применительно к исполнению своих должностных обязанностей а также с учетом срока пребывания в дип ранге результатов служебной деятельности и при наличии сертификата о соответствующем уровне владения иностранным языком. Присваиваются следующие...
84691. Порядок отбора кандидатов на замещение дипломатических должностей в загранучреждениях 17.24 KB
  Форма прохождения дипломатической службы в каждой стране своя; она определяется национальными особенностями и традициями этой страны а также материальными и кадровыми возможностями. Однако под влиянием глобализации происходит стандартизация дипломатической службы во всем мире. Это крайне ценимые качества на дипломатической службе. Главные требования которые предъявляются к лицам претендующим на должность дипломатической службы: сильная политологическая правовая и общегуманитарная подготовка классического образца воспитанность и хорошие...
84692. Прохождение службы в центральном аппарате МИДа и за рубежом. Ротация дипломатических кадров 14.15 KB
  Основные элементы системы прохождения дипслужбы можно представить следующим образом: А отбор на службу оценка претендентов с точки зрения их профессиональнокфалификационной психологической и духовнонравственной пригодности к работе в структурах МИД и загранучреждений РФ; Б процедура утверждения в должности замещение госдолжности дипслужбы; В профессиональная адаптация путём испытательного срока стажировок повышения квалификации и других мероприятий; Г оценка служебной деятельности формирование здоровых ценностных установок на...