Частина тексту файла (без зображень, графіків і формул):
Міністерство освіти та науки України
Національний університет «Львівська політехніка»Кафедра «Телекомунікацій»
/
Практична робота 1
з дисципліни „Основи теорії систем”.
Дослідження розв’язків лінійної системи другого порядку.
Роботу
Львів – 2017
Теоретичний вступДосліджуємо лінійну систему другого порядку з диференціальним рівнянням
d2y(t)/dt2+δ∙dy(t)/dt+ω2∙y(t)=ω2∙u(t) (1)
та початковими умовами y(0)=y0; dy(0)/dt=y10.
Відповідне однорідне диференціальне рівняння:
d2y(t)/dt2+δ∙dy(t)/dt+ω2∙y(t)=0; y(0)=y0; dy(0)/dt=y10. (2)
Якщо коефіцієнти δ, ω – константи, то загальний розв’язок рівняння (1) складається із загального розв’язку однорідного рівняння (2) та частинного розв’язку неоднорідного рівняння (1):
y(t)=a1exp(p1t)+a2exp(p2t)+u(t). (3)
Коефіцієнти a1 і a2 є сумісними розв’язками рівнянь початкових умов однорідного диференціального рівняння a1+a2=y0; a1p1+a2p2=y10:
a1=(y10–y0 p2)/(p1–p2); a2=y0(1+p2/(p1–p2))–y10/(p1–p2).
Характеристичне рівняння утворюється формальними замінами у рівнянні (2) d/dt → p; d2/dt2 → p2; та скорочення y(t):
p2+δp+ω2=0. (4)
Показники експонент у (3) p1 і p2 є розв’язками характеристичного рівняння (4):
. (5)
Є два суттєво відмінні види розв’язків лінійної системи другого порядку: коливний та аперіодичний. Коливний розв’язок є тоді, коли p1 і p2 комплексні спряжені, тобто підкореневий вираз у формулі (5) від’ємний: δ/2<ω. Аперіодичні розв’язки є при δ/2>ω. Граничному випадку відповідає рівність: δ/2=ω.
Розв’язки залежать також від вхідного сигналу u(t) та від початкових умов.
Лінійні системи зручно записувати за допомогою передавальних функцій (transfer functions), які утворюються формальною заміною оператора диференціювання у рівнянні системи (1) символом p (d/dt → p, d2/dt2 → p2):
(p2+δp+ω2) y(t)=ω2 u(t).
Звідси дробово-раціональна передавальна функція формально є відношенням вихідного сигналу до вхідного
W(p)=y(t)/u(t)=ω2/(p2+δp+ω2). (6)
Текст програми дослідження розв’язків системи другого порядку на мові MATLAB наведено далі.
clear;
delta=8; omega=2*pi*5;
t=[0:0.001:1.5]; u=1-exp(-t/0.001)+0.0*(rand(size(t))-0.5);
num=[0 0 omega^2]; den=[1 delta omega^2];
y=lsim(num, den, u, t);
plot(t, u, 'r', t, y, 'b');
Графіки числових експериментів
/
Графік коливного розв’язку
/
Графік аперіодичного розв’язку
/
Графік граничного розв’язку
/
Графік вхідного сигналу з шумовою добавкою 0.01
Висновок:
Під час виконання практичної роботи я досліджував розв’язки лінійного диференціального рівняння 2-го порядку, побудував графіки для трьох випадків, визначив вплив шумової добавки вхідного сигналу, а також знайшов значення шумового коефіцієнту при якому шум помітний на графіку.
Ви не можете залишити коментар. Для цього, будь ласка, увійдіть
або зареєструйтесь.
Ділись своїми роботами та отримуй миттєві бонуси!
Маєш корисні навчальні матеріали, які припадають пилом на твоєму комп'ютері? Розрахункові, лабораторні, практичні чи контрольні роботи — завантажуй їх прямо зараз і одразу отримуй бали на свій рахунок! Заархівуй всі файли в один .zip (до 100 МБ) або завантажуй кожен файл окремо. Внесок у спільноту – це легкий спосіб допомогти іншим та отримати додаткові можливості на сайті. Твої старі роботи можуть приносити тобі нові нагороди!