Частина тексту файла (без зображень, графіків і формул):
Міністерство освіти та науки України
Національний університет «Львівська політехніка»Кафедра «Телекомунікацій»
/
Практична-лабораторна робота 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=4; omega=2*pi*4;
t=[0:0.001:1.5];
u=1-exp(-t/0.001)+0.03*(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');
Графік перехідного процесу для коливного розв’язку:
δ/2<ω. При δ=4 та ω=2*pi*4
Графік перехідного процесу для аперіодичного розв’язку:
δ/2>ω. При δ=40 та ω=2*pi*3
Графік перехідного процесу для граничного розв’язку:
δ/2=ω. При δ=12 та ω=2*pi*6:
Шуми стають проявлятись на графіки при збільшенні коефіцієнта амплітуди шуму до значення 0,03
Висновок: виконавши лабораторну роботу я показав зміну форми вихідного сигналу системи лінійних диференціальних рівнянь ІІ-го порядку. Також дослідив вплив шумової добавки вхідного сигналу, підібравши мінімальне значення його множника, при якому стають помітним шуми на графіку вихідного сигналу.
Ви не можете залишити коментар. Для цього, будь ласка, увійдіть
або зареєструйтесь.
Ділись своїми роботами та отримуй миттєві бонуси!
Маєш корисні навчальні матеріали, які припадають пилом на твоєму комп'ютері? Розрахункові, лабораторні, практичні чи контрольні роботи — завантажуй їх прямо зараз і одразу отримуй бали на свій рахунок! Заархівуй всі файли в один .zip (до 100 МБ) або завантажуй кожен файл окремо. Внесок у спільноту – це легкий спосіб допомогти іншим та отримати додаткові можливості на сайті. Твої старі роботи можуть приносити тобі нові нагороди!