Міністерство освіти та науки України
Національний університет «Львівська політехніка»
Кафедра КСА
Курсова робота
з дисципліни
«КОМП’ЮТЕРНІ МЕТОДИ
ДОСЛІДЖЕННЯ СИСТЕМ КЕРУВАННЯ»
АНАЛІЗ ДИНАМІЧНИХ ПРОЦЕСІВ НЕЛІНІЙНИХ КОРЕГУЮЧИХ ЛАНОК СИСТЕМ КЕРУВАННЯ
Варіант 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
Контур І: