23028

Задачі ідентифікації динаміки систем з розподіленими параметрами

Лекция

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

Псевдоінверсні методи [2227] обернення алгебраїчних інтегральних та функціональних перетворень дозволяють виконати таку заміну побудувати моделюючі функції в неперервному або дискретному вигляді тільки при відомій функції матриці Гріна в необмеженій просторовочасовій області. Викладена ж в лекції 2 методика побудови функції дозволяє виконати це для систем динаміка яких описана вже диференціальним рівнянням вигляду 1.7 зведеться до знаходження перетворюючої функції функції Гріна в нашому розумінні такої що 15.4 побудови...

Русский

2013-08-04

276.5 KB

2 чел.

6

 Стоян В.А.  

Лекція 15. Задачі ідентифікації динаміки систем з розподіленими

                   параметрами

15.1. Загальна постановка задачі. Розглядувані вище підходи до дослідження динаміки систем з розподіленими параметрами ґрунтувалися на заміні її диференціальної моделі (1.1)(1.4) інтегральною моделлю вигляду (1.9), (1.11)(1.13). Псевдоінверсні методи [22,27] обернення алгебраїчних, інтегральних та функціональних перетворень дозволяють виконати таку заміну (побудувати моделюючі функції в неперервному або дискретному вигляді) тільки при відомій функції (матриці) Гріна  в необмеженій просторово-часовій області. Викладена ж в лекції 2 методика побудови функції  дозволяє виконати це для систем, динаміка яких описана вже диференціальним рівнянням вигляду (1.1).

Для більшості досліджуваних практиками процесів давно побудовані диференціальні моделі і проблем використання наведених нами підходів до моделювання їх динаміки в прямій та оберненій постановці не існує. Технічні ж процеси, які мають місце при науковому вивченні та практичній реалізації сучасних інформаційних технологій, досить часто виходять за межі моделей вигляду (1.1). Нерідкі випадки, коли досліджуваний процес настільки складний і не вивчений, що побудувати диференціальні залежності для нього або принципово неможливо, або вони будуть настільки далекими від реальності, що користуватися ними буде недоцільно. Тут і потрібні не фізико-технологічні підходи до опису динаміки процесу, які ґрунтуються на математичному представленні недостатньо вивчених, або недоступних для формалізації реальних процесів, а принципово нові методики, основані на спостереженнях за станом системи та зовнішньо-динамічними умовами, в яких ця система функціонує. Другими словами: потрібен математичний апарат, який би міг надійно і з бажаним (або наперед відомим) ступенем точності пов’язати спостережуваний стан системи зі спостереженнями за тими фізико-технологічними факторами, які його спричиняють.

Якщо, залишаючись в рамках прийнятих нами термінологій, позначити через  вектор-функцію стану системи в просторово-часовій області , а через вектор-функцію зовнішньо-динамічних впливів на цю систему, то проблема побудови такої бажаної інтегральної моделі вигляду (1.7) зведеться до знаходження перетворюючої функції (функції Гріна в нашому розумінні) такої, що

   (15.1)

де

,           (15.2)

та  – області, в яких виконується спостереження за станом  системи та вектора  зовнішньо-динамічних впливів на неї відповідно.

На підтвердження сказаного в лекції 1 ще раз зауважимо, що в загальному випадку розв’язати задачу побудови множини  неможливо. Тому, як і при побудові та дослідженні моделей вигляду (1.9), (1.11)-(1.13), зупинимося на трьох частинних випадках поставленої задачі.

15.2 Проблеми ідентифікації дискретизованих інтегральних моделей динаміки систем з розподіленими параметрами. Як і при розв’язанні (лекція 3) сформульованої в лекції 1 загальної задачі побудови інтегральної моделі вигляду (1.9), (1.11)-(1.13) зупинимося для початку на випадку, коли вектор-функція стану () та вектор-функція () зовнішньо-динамічних впливів на систему спостерігаються в скінчених множинах точок  та  відповідно.

В цьому випадку шукана перетворююча матрична функція  вироджується в матрицю  блоки-елементи  якої знаходяться з умови, щоб

                                        (15.3)

де - крок дискретизації області . Позначивши через  номер спостереження за системою та вводячи до розгляду вектори

де  для  задачу (15.3) знаходження елементів матриці G замінимо наступною:

                                                            (15.4)

Задача (15.4) побудови матриці G значень перетворюючої матричної функції  значно простіша задачі (15.1), (15.2) побудови множини  матричних функцій . Вона може бути розв’язана з використанням методів лінійної алгебри, покладених в основу побудови (лекція 3) дискретизованого варіанту інтегральної моделі динаміки розглядуваних систем. Дійсно, як неважко бачити, ідентифікаційна задача (15.4) еквівалентна задачі обернення матричної системи

GU=Y,                                                                                           (15.5)

де

згідно критерію

                                                                    (15.6)

дослідженої детально в [].

15.3 Проблеми ідентифікації систем з дискретно спостережуваною функцією стану. Другим частинним випадком ідентифікації задачі (15.1), (15.2) є випадок, коли на відміну від розглядуваного вище, дискретні по  вимірювання значень  вектор-функції стану  виконуються при неперервно спостережуваній вектор-функції  зовнішньо-динамічних збурень системи.

Проблема знаходження вектора  значень  матричної функції , які відповідають спостережуваним в  значенням  вектор-функції  зведеться до обернення системи інтегральних рівнянь:

                                                                     (15.7)

де .

Позначивши через N кількість спостережень за системою, а через  номер такого спостереження, від системи (15.7) перейдемо до наступної:

                                   (15.8)

де матриця Y співпадає з визначеною вище, а

(тут u(n)() n-е вимірювання вектор-функції  зовнішньо-динамічних збурень).

Нижче ми побудуємо та дослідимо множину

            (15.9)

псевдорозвязків ідентифікаційних співвідношень (15.8) узагальнюючи методику розв’язання дискретизованого варіанту задачі, аналогічно тому як в лекції 4 узагальнювалися результати розв’язання задач сформульованих в лекції 3.

15.4 Проблеми ідентифікації систем з дискретно спостережуваною функцією зовнішньо-динамічних впливів. Аналогічно сформульованому вище заслуговує на увагу ще один частинний випадок ідентифікаційної задачі (15.1),(15.2) – випадок, коли вектор-функція  стану системи може бути спостережуваною неперервно по s, а вектор-функція  зовнішньо-динамічних впливів доступна для вимірювань тільки в M точках .

Проблеми знаходження вектора  значень  матричної функції , які спостережувані в точках , зовнішньо-динамічні впливи um=u(sm) перетворює в спостережувану в  вектор-функцію , зведеться до обернення системи функціональних рівнянь

                                        (15.10)

де .

Вимагаючи виконання співвідношення (15.9) при кожному n-му () спостереженні заключаємо, що шукана матрична вектор-функція  задовольнятиме рівнянню

                                        (15.11)

де матриця U співпадає з визначеною вище, а

(тут)n-е вимірювання вектор-функції  стану системи).

Далі, з використанням методики поширення результатів розв’язання дискретного варіанту (лекція 3) задачі моделювання динаміки систем з розподіленими параметрами на дискретно керовані (лекція 5) системирезультати розв’язання дискретизованого варіанту (15.5) ідентифікаційної задачі (15.1),(15.2) будуть поширені і на розглядуваний тут випадок.

Буде побудована та досліджена множина

       (15.12)

псевдорозв'язків ідентифікаційних співвідношень (15.11).

15.5 Ідентифікаційно-псевдоінверсний підхід для побудови матричної функції Гріна динаміки систем з розподіленими параметрами. Розглянуті вище задачі побудови матричних вектор-функцій  та , якщо вони будуть розв’язані, дозволяють досить близько підійти до вирішення проблеми побудови інтегральної моделі (1.7) динаміки систем з розподіленими параметрами. Це можна зробити для систем неперервно спостережуваних як по вхідному зовнішньо динамічному збуренню , так і по вихідному вектору  стану системи.

Дійсно, виконуючи спостереження:

  1.  за вектор-функцією  в точках  при заданій вектор-функції  зовнішньо-динамічних збурень;

б) за вектор-функцією  при вимірюваній в точках  вектор-функції зовнішньо-динамічних збурень можна ставити та розв’язувати задачі (15.8)-(15.9) та (15.11)-(15.12) по побудові матричних вектор-функцій   та   відповідно.

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

15.6. Про загальність псевдоінвертного підходу до ідентифікації систем з розподіленими параметрами. Сформульовані в пунктах 15.3 та 15.4 задачі (15.8)-(15.9) та (15.11)-(15.12) по знаходженню матричних вектор-функцій  та  є допоміжними для розв’язання задачі побудови перерізів двохаргументної матричної функції Гріна  при ідентифікації інтегральної моделі динаміки систем з розподіленими параметрами.

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

  1.  дискретного матричного перетворювача вигляду

,                   (15.13)

де  та  векторні вхід та вихід системи, а  - шукана матриця, якою виконюється перетворення входу x у вихід y;

2) інтегрального перетворювача вигляду

,                  (15.14)

де ,  для , розподілений в S вхід та дискретний вихід системи, а  шукана матрична функція, якою вхід  перетворюється у вихід y;

3) функціонального перетворювача вигляду

(),                  (15.15)

де ,  та  для  дискретний та розподілений в S вхід та вихід системи, а шукана матрична функція, якою вхід x перетворюється у вихід .

Розглядувані нижче підходи до розв’язання задач (15.5)-(15.6), (15.8)-(15.9) та (15.11)-(15.12) дозволяють розв’язати задачі ідентифікації матриці А та матричних функцій, – перетворювачів вигляду (15.13)-(15.15) згідно критеріїв:

,                  (15.16)

,                  (15.17)

,                  (15.18)

де N – кількість спостережень за відповідним перетворювачем, а n – номер такого спостереження.

Враховуючи сказане, а також те, що задачі (15.13)-(15.18), в силу своєї абстрактності простіші для сприйняття, ніж сформульовані вижче задачі (15.5)-(15.6), (15.8)-(15.9) та (15.11)-(15.12), при побудові ідентифікаційних методів дослідження динаміки розглядуваних систем будемо виходити з останніх.

Для випадків, коли точність, з якою знайдені згідно (15.16)-(15.18) матриця А та матричні функції, , задовольняють перетворенням (15.13)-(15.15), недостатня - розглянемо варіанти нелінійного (поліноміального) перетворення вхідного сигналу у вихідний. Другими словами: крім матриці  А та матричних функцій, ,  побудуємо ще і поліноміальне перетворення  , , де , , таке, щоб

,

,

,

де

,

.

Зауважимо, що всі ці задачі стосовно алгебраїчного перетворення вигляду (15.5) розв’язані в [ ]. Запропонована там методика розв’язання сформульованих вище задач і буде викладена нами нижче.

5

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


 

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

51541. СЦЕНАРІЙ ДО ДНЯ ВЧИТЕЛЯ 54.5 KB
  Ведуча: Вклонiться люди вчителю своєму Володаревi слова i пера Крiзь бiль i радiсть роздуми i втому любов його безсмертна проросла. Ведуча: Дорогi вчителi До Вас сьогоднi слова щирої шани сердечної вдячностi. Ведуча: Щоранку поспiшає до школи наш директор.I сьогоднi дорогий наш директор даруємо Вам...
51542. 1 вересня 91.5 KB
  І ліцей відкриває гостинно всі двері назустріч новій дітворі. Сміх навкруги розсипає Дружна ліцейна сімя Сьогодні знов у наше братство проводжає Дітей своїх Україна моя Ідуть вони в ліцей знання щоб здобувати Щоб стати найкращими в світі І щоб Україну свою будувати Й прославить у цьому столітті. Знов гомінливо у ліцейному дворі Це значить що стартує рік навчальний Ми зустрічаємо із відкритим серцем першачків Ліцей співає пісню величальну Хор співає Зелене жито Слова палкі мелодія врочиста Державний гімн і прапор майорить...
51543. Планирование на предприятии. Е.О. Ушакова 574.5 KB
  В методических указаниях изложены основные практические вопросы, входящие в стандартный курс дисциплины «Планирование на предприятии». В данном издании приведена программа курса, практические задания, темы курсовых работ, методические указания по выполнению курсовой работы, темы семинарских занятий, контрольные вопросы, тестовые задания для проверки знаний студентов и список рекомендуемой литературы.
51544. Сценарій свята День знань для першокласників 11.86 MB
  Сценарій свята Перший урок у першому класі Дійові особи: Коза Дереза Телесик вчителька першокласники учні 5ого класу. Вчителька. Шановна пані Може ви знаєте що тут за свято І хто ви така Вчителька. Здрастуйте добродійко Я вчителька.
51545. Знову до школи 54 KB
  How mny units re there in this textbook Wht re the nmes of these units Wht re we going to study in them Wht sections re there t the end of the book Wht informtion cn you find in these sections Учні ланцюжком відповідають на запитання вчителя. Look t the picture nd tell wht we re going to tlk bout during our lesson tody. Т: sk your friends wht they did during the summer holidys. HO1: Wht did your friends do during June: ________________________________________________________________________________ July:...
51548. Клімат . 8 клас 337 KB
  Красоткіна Е 6 Відповіді: А 6 циклон Б В Г Д Е Завдання №2. Красоткіна А 1 Все зацвілопора натхненна Така чарівна у природі мить Убрались в зелень явори і клени І в світллі сонця річечка блищить. Красоткіна Б 2 Пора летить бентежна ніжна У квітах росах чарівна. Красоткіна...
51549. Знайомство з країною Фантазією 40.5 KB
  Обладнання: альбом плакат із зображенням замку малюнок Королеви кольорові олівці фломастери клубок ниток для гри фотографії дітей яблуко гриб шматочки льоду. Як ви думаєте де живе Королева Звичайно у замку. На дошці плакат із зображенням замку. Яку гарну країну ми створили який чудовий замок збудували А про господиню Королеву замку забули Тож давайте познайомимось із нею ближче.