Частина тексту файла (без зображень, графіків і формул):
МІНІСТЕРСТВО ОСВІТИ І НАУКИ, МОЛОДІ ТА СПОРТУ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ „ЛЬВІВСЬКА ПОЛІТЕХНІКА”
ІКТА
Кафедра Захист інформації
З В І Т
До лабораторної роботи №4
з курсу:
„ Засоби прийому та обробки інформації в системах технічного захисту ”
на тему:
„ ВИЗНАЧЕННЯ ХАРАКТЕРИСТИК ЛІНІЙНИХ СТАЦІОНАРНИХ СИСТЕМ ЗА ДОПОМОГОЮ ІНСТРУМЕНТАЛЬНИХ ЗАСОБІВ MatLab ”
Мета роботи : вивчити засоби МАТЛАБ для визначення характеристик лінійних стаціонарних систем та отримати навики проведення аналізу динамічних систем.
Завдання:
Номер
Параметри динамічної системи
варіанту
Поліном чисельника
Поліном знаменника
2
1.0
(((1/5)s)^2+2*1/5*0.1s+1)(3s+1)
Лістинг програми:
clc
clf
S=tf('s');
H=((1.0)/((((1/5)*S)^2+2*1/5*0.1*S+1)*(3*S+1)))
[Num,Den]=tfdata(H,'V');
sys_tf=tf(Num,Den) % форма tf
[z,p,k]=tf2zp(Num,Den);
sys_zpk=zpk(z,p,k) % форма zpk
figure(1)
step(sys_tf) % перехідна функція
damp(sys_tf) % полюси системи
figure(2)
impulse(sys_tf) % імпульсна функція
[X,w]=freqs(Num,Den); % комплексний коефіцієнт передачі
figure(3)
subplot(2,1,1)
plot(w,abs(X)),grid on % АЧХ
subplot(2,1,2)
plot(w,unwrap(angle(X))*180/pi),grid on % ФЧХ
figure(4)
bode(sys_tf),grid on % ЛАЧХ
nyquist(sys_tf) % АФЧХ
Transfer function:
1
--------------------------------
0.12 s^3 + 0.16 s^2 + 3.04 s + 1
Transfer function:
1
--------------------------------
0.12 s^3 + 0.16 s^2 + 3.04 s + 1
Zero/pole/gain:
8.3333
-------------------------
(s+0.3333) (s^2 + s + 25)
Eigenvalue Damping Freq. (rad/s)
-3.33e-001 1.00e+000 3.33e-001
-5.00e-001 + 4.97e+000i 1.00e-001 5.00e+000
-5.00e-001 - 4.97e+000i 1.00e-001 5.00e+000
/
Перехідна функція
/
Імпульсна функція
/
Амплітудо-частотна характеристика
/
Висновок:
На цій лабораторній роботі я ознайомився із засобами Matlab для визначення характеристик лінійних стаціонарних систем та отримати навики проведення аналізу динамічних систем.