Міністерство освіти та науки України Національний університет «Львівська політехніка» Кафедра КСА Курсова робота з дисципліни «КОМП’ЮТЕРНІ МЕТОДИ ДОСЛІДЖЕННЯ СИСТЕМ КЕРУВАННЯ» АНАЛІЗ ДИНАМІЧНИХ ПРОЦЕСІВ НЕЛІНІЙНИХ КОРЕГУЮЧИХ ЛАНОК СИСТЕМ КЕРУВАННЯ Варіант 3 Зміст І - Теоретичні відомості. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3 ІІ - Завдання. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 ІІІ - Виведення системи диференціальних рівнянь. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6 IV - Блок-схема алгоритму. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 V - Код програми. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .9 VI - Графічні результати перехідних процесів. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .13 VII - Висновок. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .16 VIII - Список використаної літератури. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .16 I. Теоретичні відомості Метод Рунге-Кутта Методи Рунге-Кутта — важлива група чисельних методів розв’язування (систем) звичайних диференціальних рівнянь. Названі на честь німецьких математиків Карла Рунге і Мартіна Кутта,які відкрили ці методи. Цей метод айбільше часто вживається при чисельному відшуканні розв’язку задачі Коші
(1), при умові
і дозволяє одержати наближення високої точності. Геометрично цей метод для задачі Коші також полягає в тому, що на малому відрізку [х; х+h] інтегральна крива у=у(х) рівняння (1) заміняється відрізком прямої, що проходить через точку (х; у(х)). Однак в основу методу покладений більше тонкий, чим у методах Ейлера, підхід до визначення напрямку цього відрізка прямій. Нехай відрізок
розділений на п рівних частин точками
,
і визначені наближені значення
розв’язку диференціального рівняння відповідно в точках
. Переходимо до відрізка
й відшукання
(рис. 1). / Рис. 1 Визначаємо
− напрямок дотичної до інтегральної кривої в точці
, і точку перетину прямих
і
, тобто точку
. Знаходимо напрямок дотичної в точці
:
і із точки
проводимо пряму з кутовим коефіцієнтом
:
до перетину із прямою
. Одержуємо точку
. Знаходимо напрямок дотичної в точці
:
і із точки
проводимо пряму з кутовим коефіцієнтом
:
до перетину із прямою
. Одержуємо точку
. Далі визначаємо напрямок дотичної в точці
:
. Остаточний напрямок відрізка ламаної, що представляє наближений розв’язок задачі, буде рівним
і проводимо із точки
пряму
, до перетинання із прямої
в точці
, де
вважаємо наближеним значенням розв’язку в точці
(див. рис. 1). Алгоритм методу Рунге-Кутта для диференціального рівняння першого порядку: Передбачаються заданими рівняння
, початкова умова
і відрізок
. 1. Задаємо число п точок поділу відрізка
й обчислюємо крок
. Вважаємо відомими
й переходимо до дії 2. 2. Нехай знайдені
. Визначаємо
,
,
,
,
,
,
. Якщо
(k+1=n), то процес закінчений. Числа
представляють наближені значення шуканого розв’язку в точках
. Якщо ж
(k+1<n), то повторюємо дію 2, вважаючи вихідним
. Всі розрахунки по алгоритму зручно оформляти у вигляді таблиці. Обчислення по методу Рунге-Кутта значно ускладнені в порівнянні з методом Ейлера, але за рахунок цього він дає меншу похибку при заміні точного розв’язку
наближеним
. З теорії наближених методів відомо, що при кроці інтегрування h має місце оцінка
, так що похибка одного кроку обчислень (визначення
по
) має порядок
(або
). Сумарна похибка за п кроків, тобто похибка приблизного наближеного розв’язку
в точці
буде порядку
(або
). Звідси, якщо збільшити п у два рази, похибка приблизно зменшиться в 16 разів. Тому для оцінки наближеного розв’язку
, отриманого із кроком h, повторюють обчислення із кроком 2h і за абсолютну похибку приймають число
, де
− наближений розв’язок із кроком 2h. Наведена оцінка є оцінкою методу й не враховує похибку, отриману при округленні. Поліноміальна апроксимація нелінійних характеристик елементів Метод вибраних точок Як і кубічний сплайн цей метод використовують, коли функція задана у вигляді таблиці значень аргументу і функції. Спочатку необхідно вибрати вираз для апроксимації. Наприклад, можна взяти поліном третього порядку
. (1) Тут необхідно визначити чотири коефіцієнти
. Для цього необхідно мати систему чотирьох алгебричних рівнянь. Візьмемо з таблиці два значення функції
(2) і доповнимо їх ще двома умовами, а саме значеннями похідної в точках
і

,
. (3) Похідна кубічного поліному (1) буде мати вигляд
. (4) Для значень (2), (3), (4) складемо систему лінійних алгебричних рівнянь, підставивши значення вузлів апроксимації (2) в (1), а значення похідних у точках
і
(3) у (4)
(5) Розв’язавши систему рівнянь (5), визначимо коефіцієнти апроксимації
на відрізку
ІІ. Завдання / Метод: Рунге-Кутта(формули в):
Дані параметрів схеми:

Апроксимація залежності
виконується з вибором розрахункової формули
Напруга живлення заданана рисунку, де
. ІІІ. Виведення системи диференціальних рівнянь
Рівняння струмів і напруг Вузол І: i1=i2+i3 Контур І: