23028

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

Лекция

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

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

Русский

2013-08-04

276.5 KB

1 чел.

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

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


 

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

47758. ОСНОВИ ФІНАНСІВ ПІДПРИЄМСТВ 1.28 MB
  Субєктами таких відносин можуть бути підприємства та організації банківські установи та страхові компанії позабюджетні фонди інвестиційні фонди аудиторські організації інші суб'єкти господарювання які є юридичними особами. Формування фінансових ресурсів на підприємствах відбувається під час формування статутного фонду а також у процесі розподілу грошових надходжень у результаті повернення авансованих коштів у основні та оборотні фонди використання доходів на формування резервного фонду фонду споживання і фонду накопичення. Формування...
47759. Криминалистика. Конспект лекций 549 KB
  Понятие принципы элементы и значение планирования расследования преступлений. Техника планирования расследования. Информационно-компьютерное обеспечение раскрытия и расследования преступлений. Роль автоматизированных информационно-поисковых систем АИПС в информационном обеспечении раскрытия и расследования преступлений.
47760. Криминалистика. Курс лекций 3.67 MB
  Криминалистика юридическая наука прикладного характера дает в руки следователей оперативно-розыскных работников и экспертов-криминалистов научно обоснованные и проверенные практикой средства приемы и методы раскрытия расследования и предупреждения любых самых сложных преступлений. Современная криминалистика продолжает активно развиваться впитывая в себя все новые знания из различных областей других наук для разработки на этой основе более совершенных средств и методов расследования. вышло в свет первое в России пособие...
47761. Кримінальний процес. Конспект лекцій 1.57 MB
  Порядок застосування інших заходів процесуального примусулекція №6 порушення кримінальної справи1. Поняття і завдання стадії порушення кримінальної справи2. Приводи і підстави для порушення кримінальної справи. Обставини що виключають порушенння кримінальної справи.
47762. Становлення культурології як самостійної галузі гуманітарного знання 133.03 KB
  Без світу культури важко собі уявити світ особистості. До культури в цілому відноситься широкий діапазон людських почуттів і думок від пошуку смислу життя до естетики. Той чи інший досягнутий рівень культури людства визначає кожен раз заново окультурення кожної народженої людини в результаті чого врештірешт відбувається окультурення людської природи. Потрібні були століття розвитку культури щоб рука сучасної людини змогла навчитися виготовляти складну техніку створювати високодосконалі твори скульптури твори живопису чи...
47763. Культура мови. Опорний конспект лекцій 2.73 MB
  Кожна людина повинна знати мову своєї держави, душею відчувати причетність до культури нації, а через свою культуру - причетність до вселюдської культури
47764. Управление персоналом и трудовые отношения 1.68 MB
  Оценка выполняемой работы. Введение схемы оценки выполняемой работы. Кроме того что данный курс слишком краток Вы просто и не должны быть таким специалистом если не посвятили себя целиком вопросам работы с кадрами. Каждая из вышеперечисленных функций включает в себя: планирование: постановка целей и стандартов разработка правил и последовательности действий разработка планов и прогнозирование некоторых возможностей в будущем; организация: постановка определенных задач перед каждым подчиненным разделение на отделы...
47765. ОРГАНИЗАЦИЯ И ТЕХНОЛОГИЯ РЕКЛАМНОЙ ДЕЯТЕЛЬНОСТИ 1.9 MB
  Данное учебное пособие по дисциплине «Организация и технология рекламной деятельности» предназначено для обучения студентов СПО третьего курса повышенного уровня обучения
47766. Курс лекцій. Історія української культури 3.28 MB
  Вивчення цього предмету повязане з необхідністю гуманізації університетської освіти поєднаних глибоких професійних звань з опануванням багатої історії культури України та людства в цілому Історія української культури†належить до обов’язкових до вивчення в вищому навчальному закладі дисциплін. Зокрема акцентується увага на формуванні української народності та нації на етнічному складі населення сучасної України. Етнічний склад населення України