Навигация
Алгоритми розрахунку періодичного режиму в нелінійній схемі
18104
знака
0
таблиц
1
изображение
АЛГОРИТМИ РОЗРАХУНКУ ПЕРІОДИЧНОГО РЕЖИМУ В НЕЛІНІЙНІЙ СХЕМІ
1. Коротка характеристика методів розрахунку
Вузли, проектування яких нас цікавить, в більшості випадків використовуються в режимах, де нелінійність виявляється досить сильно. Тому не розглядають методи, які знаходять практичне застосування лише для аналізу слабких нелінійних схем. До їх числа відноситься, наприклад, метод рядів Вольтера.
Методи, які дозволяють розраховувати довільні нелінійні кола, поділяють на дві великі групи – часові та спектральні.
Характерна властивість методів першої групи – це інтегрування нелінійних диференційних рівнянь до отримання усталеного періодичного режиму. Головний недолік їх в тому, що цікавлячись усталеним режимом, ми повинні розрахувати ще і перехідні процеси. Із цим недоліком можна миритись, якщо зусилля і затрати на незаплановані розрахунки невеликі, тобто якщо, наприклад, перехідний режим продовжується недовго. Звісно, що перехідний процес в схемі тим коротший, чим менше її вибірність.
Необхідність визначення періодичного режиму у вибірних пристроях створила ряд прийомів, скоротивши розрахункову процедуру. Пояснимо їх на прикладі. Нехай у нелінійній схемі період усталеного режиму відомий, задана величина і, крім цього, можна говорити, що перехідний процес практично завершується за L періодів. Таким чином, щоб знайти періодичний режим, треба інтегрувати диференційне рівняння схеми протягом L періодів. Перший прийом складається в зменшені часу інтегрування на кожному періоді. Це змушує використовувати такі чисельні методи, які, зберігаючи потрібну точність, дозволяють вести інтегрування з максимальним кроком. Ідея другого прийому – виконувати інтегрування не на кожному періоді, а із пропусками. Для її реалізації формується функція незв’язності, котра характеризує ступінь досягнення усталеного режиму. За допомогою цієї функції, із початковими умовами на якомусь періоді, визначаються початкові умови для інтегрування на наступному періоді.
Очевидно, якщо перехідний процес закінчиться, то початкові умови, використані на попередньому періоді, співпадуть з початковими умовами, обчисленими для наступного. Виявилось, що за початковими умовами для k–го періоду можна приблизно знайти початкові умови для (k+m)-го періоду, де m- ціле число, більше одиниці. В результаті, число періодів, протягом яких треба інтегрувати рівняння, скоротиться в m разів.
У спектральних методах розрахунку визначається періодичне рішення нелінійних диференційних рівнянь, записаних у формі ряду Фур’є. Відносно спектральних компонент цього ряду утворюється система нелінійних (трансцендентних) рівнянь, котра вирішується за допомогою ітерацій. Різновиди методів цієї групи визначається тим, як побудовано ітераційний процес.
Для схемотехнічного проектування розрахунок періодичного режиму потрібен як у випадку, коли період процесу відомий, так і коли період повинен бути знайдений. Перша ситуація характерна для підсилювачів потужності, помножувачів та дільників частоти, тобто для схем, в яких є зовнішня дія. В таких схемах в якості робочого використовують періодичний режим із періодом, рівним періоду зовнішнього сигналу або в ціле число разів більшим за цей період. Подібні схеми звуть неавтономними. Другий клас схем – автономні, наприклад, автогенератори. В них період коливань визначається внутрішніми параметрами і знаходиться разом із амплітудою коливань.
У зв’язку з тим, що розрахунок процесу з відомим періодом простіше, спочатку розглядається цей випадок, а потім вказується, які зміни потрібно ввести в розрахунок в іншій ситуації.
2. Про чисельні методи інтегрування звичайних диференційних рівнянь
Маючи на увазі задачу розрахунку періодичного режиму часовим методом, звернемось до чисельного інтегрування диференційних рівнянь, які описують процеси в електронних схемах. Це цікавить нас з точки зору скорочення часу розрахунків.
Відповідними методами інтегрування можна вважати такі, котрі потребують малих витрат на обчислення при дотриманні заданої точності. Під витратами розумітимемо загальний об’єм розрахунків (число арифметичних операцій), необхідний об’єм пам’яті ЕОМ, машинний час і т.ін. Цілком зрозуміло суперечність потреб високої точності і малих витрат. Дійсно, якщо не вдаватися в подробиці, то здається, що точність чисельного інтегрування тим вище, чим менший крок. З іншого боку, із зменшенням кроку зростає час розрахунку. Справа ускладнюється ще й тим, що в деяких ситуаціях зменшення кроку не вирішує проблему точності. Це зустрічається при інтегруванні “жорстких” диференційних рівнянь, в яких коефіцієнти значно різняться через великий розкид постійних часу кола. В результаті, в перехідному процесі є швидкі та повільні складові, із правильним розрахунком яких впорається далеко не всякий метод інтегрування. Таким чином, в обчислювальній математиці з’явилась необхідність вияву властивостей чисельних методів інтегрування, які впливають на витрати та точність розрахунків. В даний час прийнято характеризувати ці методи точністю і стійкістю.
Точність визначається помилками, які виникають під час розрахунків. Для порівняння методів вводиться поняття “локальна помилка”, яка стосується одного кроку інтегрування. При цьому якщо метод спирається на результати попередніх кроків, припускається, що перші дані точні. Локальна помилка складається з методичної та помилки округлення. Перша залежить від методу, друга від арифметичної точності ЕОМ. В практично використаних методах локальна методична помилка повинна допускати оцінку. Крім того, має місце локальна перехідна помилка, яка зобов’язана своїм існуванням похибкам на попередньому кроці.
Для користувача важлива загальна, глобальна помилка, відповідальна за весь інтервал інтегрування. Зв’язок локальної помилки із загальною складний. Тому при аналізі методу глобальна помилка не розглядається, а уяву про неї одержують за допомогою тестових прикладів.
Друга характеристика методів чисельного інтегрування – їх стійкість. Практичний прояв її виглядає так. Припустимо, інтегруємо тестове диференційне рівняння, для якого відомо точне рішення. Обчислення проводимо кілька разів з постійним кроком, збільшуючи крок від одного розрахунку до іншого. До деякої величини кроку, яку звуть критичною, похибка інтегрування постійно зростає. Але як тільки крок перейде критичне значення, похибка різко зростає (виникає чисельна нестійкість) і отримані результати значно відрізняються від справжніх. Із сказаного випливає, що нестійкий метод не придатний для задач, де виникає потреба інтегрування диференційних рівнянь протягом тривалого часу: щоб уникнути чисельної нестійкості, необхідно змінювати крок, а це збільшує кількість кроків і призводить до зростання загальної похибки через накопичення усіх видів помилок – перехідної, округленої, методичної.
Методи чисельного інтегрування розрізняються в точності, стійкості і ряді інших властивостей. Наведемо прийняту в літературі класифікацію і вкажемо властивості методів, належних до окремих класів.
Методи поділяються на дві великі групи. Відмінна особливість обох груп – спосіб апроксимації заданої функції. В першій групі використовується розклад в ряд Тейлора, в другій – апроксимація функції поліномом з тейлорівським розкладанням. Широке розповсюдження отримали методи другої групи, котрі у вітчизняній літературі звуться кінцево-різницевими.
спектральний інтегрування нелінійний періодичний
Формула довільного кінцево-різницевого методу, відносно до рішення диференційного рівняння 
при початкових умовах 
, записується так
, (1)
де 
,
- крок інтегрування.
Число p задає кількість попередніх кроків, які визначають значення 
шуканої функції. При р=0 метод зветься однокроковим. Якщо 
, то метод явний, при 
- неявний. В останньому випадку шукане значенні входить до правої частини (1) як аргумент нелінійної функції. Співмножники 
, 
(їх число дорівнює 2р+3) шукається методом невизначених коефіцієнтів при поліміальній апроксимації невідомої функції 
. Число m називають порядком методу. За допомогою невизначених коефіцієнтів складають m+1 рівняння відносно , . Коли 2p+3 > m+1, тоді частиною співмножників задаються.
Відомі формули Ейлера – явна
та неявна
відносяться до однокрокового методу першого порядку. Його локальна методи- чна похибка оцінюється величиною 
. Якщо зберегти порядок методу і зробити його багатокроковим, то підбором співмножників в (1) можна зменшити методичну похибку. Таким чином, точність кінцево –різницевих методів зростає із збільшенням їх порядку і числа попередніх кроків, які враховуються. Однак треба мати на увазі, що підвищення порядку супроводжується зменшенням області стійкості.
На перший погляд уявляється, що явні методи мають перевагу над неявними тому, що в останніх значення виходить із рішення нелінійного рівняння, а в явних розраховується за аналітичним виразом. Але як свідчить аналіз, неявні методи більш стійкі. Отож, вони допускають при заданій точності більший крок.
В даний час в алгоритмах чисельного інтегрування проблемно-орієнтованих програм використовується кінцево-різницеві методи, які мають бажану стійкість та дозволяють оцінювати локальну методичну похибку на кожному кроці. За допомогою цієї оцінки підтримується максимальний розмір кроку і вибирається мінімальний порядок методу. Для зменшення об’єму розрахунків в неявних методах розраховується спочатку за відповідною явною формулою (прогноз), а потім уточнюється за допомогою неявної (корекція). Після вибору методу чисельного інтегрування програміст основні зусилля направляє на створення ефективного алгоритму, який визначає розмір кроку.
Відносно методів інтегрування, спираючись на розклад невідомої функції в ряд Тейлора, наприклад методом Рунге-Кутта різних порядків, можна зазначити, що вони знаходять обмежене використання. Пов’язано це з двома обставинами: по-перше, ускладнюється оцінка локальної методичної похибки на кожному кроці інтегрування; по-друге, для визначення треба m разів обчислити значення першої частини диференційних рівнянь (m – порядок методу), причому ці значення неможливо використовувати на інших кроках. Друга властивість знижує ефективність розрахунків порівняно з кінцево-різницевими формулами. Методом Рунге-Кутта зручно починати чисельне інтегрування, якщо воно ведеться за багатокроковими різницевими формулами, для отримання необхідних початкових значень. Справа в тому, що методи Рунге-Кутта виявляються явними та однокроковими. Тому використання їх на початковій стадії обчислення не дуже позначиться на загальних часових втратах, а необхідна точність забезпечується правильним вибором порядку.
Похожие работы
, коли розглядається стійкість положення рівноваги, і виявляються періодичними функціями часу для збурень періодичного процесу. Очевидно, що для виділених станів схеми аналіз стійкості у великому найбільш складний, оскільки він зв’язаний з рішенням диференційного рівняння. Наступним за складністю буде вивчення локальної стійкості періодичного режиму, а самим простим – аналіз локальної стійкості – ...
... і паразитних спектральних складових. По-друге, ці коливання можуть бути інтенсивними та приводити до виходу з ладу елементів схеми (пробій та вигорання транзистора, згорання опору і т. і.). Отже, при проектуванні нелінійних вузлів треба приймати заходи по запобіганню паразитного збудження. 2 Що означає «опрацювати» схему нелінійного вузла? Серед розробників РЕА є термін «опрацьована» схема. ...
... 4 виконавців, тривалість виконання робіт складає 153 днів, ціна теми складає 7. Забезпечення життєдіяльності і питання екології 7.1 Аналіз шкідливих і небезпечних чинників В даній дипломній роботі і досліджений сервопривод урахуванням нелінійності. Основним місцем роботи була науково-виробнича дослідницька лабораторія (НВДЛ) №402, кафедри №301. При проведенні робіт в лабораторіях поді ...
... Висновки по розділу 3 У даному розділі диплома була розроблена автоматизована інформаційна система розрахунку прибутку на гірничо-збагачувальному підприємстві. Дана система була розроблена для підвищення ефективності роботи підприємства. В основу алгоритмів обробки даних покладені методи математичної статистики й оптимізаційні моделі. Для проектування і реалізації автоматизованої інформаційної ...
0 комментариев