Чисельні методи рішення систем диференціальних рівнянь на ЦОМ

Точність і погрішність отриманих рішень

При розрахунках на АОМ повстає питання про точність розв'язання. Основна величина погрішності вноситься вирішальними елементами.

Основні види погрішностей вирішальних елементів:

1. Погрішність, внесена кінцевим значенням коефіцієнтів підсилення підсилювача.

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

2. Погрішність опору ланцюгів зворотних зв'язків (ЗЗ) і вхідних ланцюгів.

Вона обумовлюється неточністю номінативних значень елементів, витоками конденсаторів і т.п.

3. Погрішність, обумовлена нестабільністю нульового рівня підсилювачів.

Дрейф нуля ОП обумовлений нестабільністю живлячої напруги, зміною параметрів ОП при нагріванні й т.п. Для зменшення впливу дрейфу нуля перед кожним запуском моделі провадиться установка нулів ОП. Також необхідно обмежувати час роботи електронної моделі на рівні 20 - 30 с.

4. Погрішність нелінійних диференціальних перетворень.

У загальному випадку, погрішність, внесена вирішальними елементами, знаходиться на рівні декількох %, що в більшості технічних рішень є цілком достатнім.

Практично всі завдання, що виникають при моделюванні електроприводів, зводяться до рішення систем нелінійних ДР. Математичний апарат для аналітичного рішення таких завдань або занадто громіздкий, або відсутній. У цих умовах єдиним виходом є застосування чисельних методів інтегрування СДР.

Рішення систем нелінійних диференціальних рівнянь полягає у відшуканні залежностей Y_i(t) при відомих початкових умовах t = t₀ (початковий час) і Y_i0 = Y(t₀) (початкові значення змінних системи).

Всі відомі чисельні методи рішення систем диференціальних рівнянь діляться на два класи - однокрокові й багатокрокові. Реалізація однокрокових методів полягає в розкладанні шуканої функції в ряд Тейлора, у якого втримується те або інше число членів, що визначає порядок методу N. При однокрокових методах рішення кожне нове значення Y_j+1 знаходиться по одному відомому попередньому значенню Y_j. Як приклад реалізації однокрокового методу приведемо найпростіший метод Ейлера (першого порядку)

де h - крок інтегрування.

Сутність даного методу складається в заміні інтегральної кривої на кожному кроці інтегрування відрізком прямій.

Більше складним, але й більше точним є метод Рунге -Кутта 4-го порядку, що задається наступними формулами:

Таким чином, при використанні методу Рунге-Кутта 4-го порядку необхідно робити чотири обчислення правих частин системи диференціальних рівнянь. Тому час обчислень по цьому методі істотно зростає, хоча точність останнього істотно вище.

Однокроковим методам притаманна властивість "самостартування", а також є порівняно проста можливість зміни кроку інтегрування в процесі обчислень.

Реалізація багатокрокових методів полягає в тому, що в їх по ряду попередніх значень Y(t) с точністю екстраполяційного полінома Лагранжа обчислюється певне значення функції. Залежно від числа використовуваних попередніх значень визначається й порядок методу. До даної групи методів відносяться: група методів Адамса, метод Мілна, метод Хеммінга, група методів прогнозу й корекції.

Реалізація методу Адамса другого порядку має вигляд

Реалізація методу Адамса четвертого порядку має вигляд

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

Дата добавления: 2014-01-04