РОЗДІЛ 3 ХАРАКТЕРИСТИКИ САК

Часові характеристики

Диференціальні рівняння незалежно від форми подання є самою загальною формою опису САК і не дають наочного зображення її властивостей. Більш наочно характеризують ці властивості функції y(t) , що є рішеннями диференціальних рівнянь.

Відомо, що те саме диференціальне рівняння має безліч рішень, конкретний вигляд яких залежить від початкових умов і від характеру функції x(t). Тому в ТАК властивості систем і їхніх елементів характеризують рішеннями, що відповідають нульовим початковим умовам і одному з типових впливів на вході, що називаються часовими характеристиками.

Найбільш широке використання при описі динамічних властивостей одержала перехідна функція h(t) . Перехідною функцією називають функцію, що описує зміну вихідної величини, що виникає після подачі на вхід одиничного східчастого впливу l(t) при нульових початкових умовах. Графік перехідної функції називається перехідною характеристикою.

Другою часовою характеристикою є імпульсна перехідна функція w(t) . Під цією функцією мають на увазі функцію, що описує зміну вихідної величини, яка виникає після подачі на вхід дельта-функції при нульових початкових умовах. Графік w(t) називають імпульсною перехідною характеристикою.

З попереднього викладу виходить, що лінійні САК описуються диференціальними рівняннями вигляду

 

де x(t) і y(t) – відповідно вхідна і вихідна величини; ai , bj – коефіцієнти;
n – порядок рівняння.

З курсу вищої математики відомо, що інтегрування рівняння (3.1) зводиться до знаходження суми загального рішення однорідного рівняння без правої частини yс(t) і якого-небудь часткового рішення неоднорідного рівняння yв (t) , тобто

Зміна вихідної величини, обумовлена складовою yс(t) , називається вільним рухом, тому що залежить тільки від вигляду лівої частини рівняння (3.1), тобто від внутрішніх властивостей самого об'єкту. Складова yв (t) , навпаки, залежить від характеру вхідного впливу, тому відповідна зміна називається змушеним рухом.

Складову yс(t) шукаємо у вигляді

де p – деяке раціональне число.

Підставивши (3.3) у рівняння (3.1) при нульовій правій частині, одержимо:

або

Останнє рівняння називається характеристичним. Таким чином, вираз (3.3) є рішенням вихідного рівняння за умови, що p є коренем рівняння (3.4). Оскільки це рівняння має n коренів, маємо і n лінійно незалежних рішень
yi (t) . Скористаємося відомою теоремою математики, що зтверджує, що коли n лінійно незалежних функцій yi (t) є рішеннями однорідного рівняння, то загальне рішення цього рівняння визначається виразом

де Ci – довільні постійні інтегрування.

Відмітимо, що вираз (3.5) справедливий тільки у випадку, якщо всі корені pi є простими. Якщо ж який-небудь корінь pj має кратність r , то в (3.5) замість r доданків вигляду (3.3) треба включити складову вигляду

Часткове рішення yв(t) звичайне шукається в тому ж вигляді, в якому задана права частина, тобто залежно від вигляду функції x(t) .

Розглянемо приклад. Приклад 3.1. САК описується диференціальним рівнянням першого по-рядку

Визначити часові характеристики h(t) і w(t) . Вирішення. Спочатку знайдемо h(t) . Маємо характеристичне рівняння:

Його єдиний корінь p = -1T. Отже

Змушену складову hв(t) шукатимемо у вигляді hв(t) = C2 . Підставивши це рішення у вихідне рівняння, одержимо C2 = – k . Тоді використаємо початкову умову h(0) = 0. Для цього запишемо рівняння:

Звідки С1 = -k

Остаточно одержимо:

 

 

Графік отриманого рішення представлений на рис. 3.1

Рис.3.1 – Характеристика h(t)

 

Для визначення w(t) вихідне рівняння перетворимо до вигляду


 

і проінтегруємо отриманий вираз

або

 

Введемо позначення

Останнє рівняння ідентичне вихідному за умови, що x(t) = 1(t) . Отже

Остаточно одержуємо:

Таким чином

w(t) = h(t)

 

Графік отриманого рішення поданий на рис. 3.2.

Рис. 3.2 – Характеристика w(t)

На рис. 3.3,а для умов прикладу 3.1 зображена реакція системи при подачі на вхід лінійного сигналу x(t) = 2t , а на рис. 3.3,б - гармонійного сигналу

Рис. 3.3 – Реакція системи при подачі на вхід інших типових сигналів

Застосування перетворення Лапласа значно спрощує визначення тимчасових характеристик.

Хід вирішення при цьому наступний:

1. Перетворимо вихідне рівняння (3.1) за Лапласом при нульових початкових умовах:

 

 

2. Вирішимо алгебраїчне рівняння (3.7) відносно Y(s) при заданому X (s) :

3. Визначимо оригінал вирішення y(t) .

У загальному випадку для знаходження y(t) використовують зворотне перетворення Лапласа (L-1 - перетворення), обумовлене формулою
Рімана-Мелліна:

де a=Re s >c0 може бути будь-яким постійним числом >c0 .

Більш простим методом є використання довідкових таблиць, в яких
наводяться зображення F (s) і відповідні їм оригінали y(t).

Більш простим методом є використання довідкових таблиць, в яких
наводяться зображення F (s) і відповідні їм оригінали y(t) .


У разі, якщо зображення є дрібно-раціональною функцією, тобто причому l<r , а коефіцієнти ci , d j - дійсні числа, застосовується формула розкладання Хевісайда:

де sj - корінь рівняння D(s) = 0 ; N - число різних корінь;
kj - кратність j -го кореня.

Диференціальні рівняння реальних САК звичайно мають прості корені s j і отже для них kj =1. Тоді вираз (3.10) з урахуванням співвідношення

 

матиме більш простий вигляд :

(3.11)

Якщо поліном D(s) має q1 кратних і q2 простих коренів, то (3.11) записується у вигляді

Оскільки визначення часових характеристик САК проводиться при
типових впливах, наведемо зображення цих впливів

 

Найменування впливу Оригінал Зображення
Східчаста функція a ×1(t) а s
Функція-дельта-функція d(t)
Розглянемо приклади.    

Приклад 3.2. Визначити часові характеристики h(t) і w(t) для САК з прикладу 3.1 операторним методом. Вирішення.

Визначимо h(t) . Для цього перетворимо за Лапласом вихідне рівняння з урахуванням того, що x(t) =1(t) :

Звідки

Отриманий вираз є дрібно-раціональною функцією, до якої можна

 

застосувати формулу розкладання Хевісайда. Тоді:

Рівняння D(s) = Ts2 + s = 0 має два корені: s1 = 0 і s2 = -1T. Скориставшись формулою (3.11), остаточно одержимо:

Аналогічно визначимо w(t) з огляду на те, що x(t) =d(t) . Маємо:

Скориставшись формулою (3.11), остаточно одержимо:

Приклад 3.3. Рівняння САК має вигляд

Визначимо часову характеристику h(t) при T = 0,3 з; x= 0,5 ; k =10 . Вирішення.

Перетворимо вихідне рівняння за Лапласом при нульових початкових умовах:

Звідки

Використаємо формулу розкладання Хевісайда. Маємо: C(s) =10

 


Рівняння

має три корені: s1=0

 

 

 

Скориставшись формулою (3.11), остаточно одержимо:

Введемо позначення: 10 = Asinj0 ; 5,774 = Acosj0.

Вирішивши отримані рівняння, одержимо:

Тоді

 

Графік характеристики h(t) наведений на рис. 3.4

Рис. 3.4 – Характеристика h(t)

 

 

Частотні характеристики

Частотні характеристики описують передаточні властивості САК в ре-жимі сталих гармонійних коливань, викликаних зовнішнім гармонійним впливом. Ці характеристики широко використовують в ТАК, тому що реальні зовнішні впливи можуть бути представлені у вигляді суми гармонійних сигналів. Вони визначаються змушеною складовою рішення диференціального рівняння при подачі на вхід впливу:

Представимо вплив (3.13) за допомогою формули Ейлера у вигляді суми двох експонентних впливів:

де

і

Вирішимо (3.1), підставивши в праву частину вираз (3.14). При цьому будемо шукати тільки змушену складову рішення yв(t) .

Використовуючи принцип суперпозиції, рішення yв(t) можна подати у вигляді двох складових

де у1(t) –рішення при


 


 


 

Будемо шукати y1(t) у вигляді:

 

З останнього виразу маємо:

W( jw) називають частотною передаточною функцією . Зрівнявши (3.18) з виразом для передаточної функції W(s) , можна зробити висновок
про те, що W( jw) є частковим випадком W(s) при s = jw.

Скориставшись прямим перетворенням Фур,є

можна зробити наступне визначення: частотною передатною функцією називається відношення вихідної величини до вхідної, перетворених за Фур'є при нульових початкових умовах.

W( jw), як і будь-яка функція комплексної змінної, може бути представлена в алгебраїчній і показовій формах.

Алгебраїчна форма:

де P(w) і Q(w) - речовинна і мнима частини відповідно.
Показова форма:

де - модуль, а - аргумент

Підставивши (3.20) в (3.17), одержимо:

(3.21)

Аналогічно одержимо складову у2(t)

Склавши (3.21) і (3.22) остаточно маємо

 

Таким чином, при гармонійному впливі на вході вихідна величина після закінчення перехідного процесу (yc(t) = 0 ) також змінюється за гармонійним законом, але з іншою амплітудою і фазою. При цьому відношення амплітуд вихідної і вхідної величин дорівнює модулю, а зміщення фаз – аргументу
W( jw). Крива, що описує кінець вектора частотної передатної функції на комплексній площині при зміні частоти від 0 до ¥, називається
амплітудно-фазовою частотною характеристикою (АФЧХ).

Крім АФЧХ, що є самою загальною частотною характеристикою, розрізняють наступні види частотних характеристик:

­ амплітудна частотна характеристика (АЧХ) – графік функції

A(w) = W( jw) ;

 

­ фазова частотна характеристика (ФЧХ) – графік функції j(w) = Arg W( jw) ;

­ речовинна частотна характеристика – графік функції P(w) = Re W( jw) ;

­ мнима частотна характеристика
– графік функції

Q(w) = Im W( jw) .

З порівняння (3.23) і (3.13) випливає важлива властивість частотних характеристик - можливість їхнього експериментального визначення на реальному об'єкті.

Приклад 3.4. Визначити частотні характеристики для умов прикладу 3.3.

Вирішення.

Перетворимо вихідне рівняння за Лапласом при нульових початкових

 

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

Зробивши заміну s = jw,маємо:

Одержимо алгебраїчну форму подання W(jw):

 

Звідси

Відповідні графіки подані на рис. 3.5.

Рис. 3.5 – Частотні характеристики