Моделирование_4
.docxНикитиной С.
ИВТ-44
Лабораторная работа №5
Вариант 11
Код Matlab состоит из функций ltiview, tf и функций рисования подобно ЛР№2.
Упражнение 1
Переходная характеристика звена:
Из графика видно, что
Прямые показатели:
Время регулирования
Пусть , тогда . Исходя из графика
Перерегулирование
По графику .
Частота колебаний
Число колебаний за время
Время достижения первого максимума
Время нарастания переходного процесса
Декремент затухания
Упражнение 2
Косвенные показатели:
Показатель колебательности
Резонансная (собственная) частота
Полоса пропускания
Частота среза
Запас устойчивости по модулю и по фазе
Упражнение 3
Из графика видно, что
Прямые показатели:
Время регулирования
Так как , то , тогда . Исходя из графика
Перерегулирование
По графику , но , тогда
Частота колебаний
Колебания непериодические, так как система является сложной составной и период непостоянен.
Число колебаний за время
Время достижения первого максимума
Время нарастания переходного процесса
Декремент затухания
Script:
Упражнение 1:
clc, clear
hold on
grid on
delta = 5;
h_u = 1;
t=0:0.01:3000;
h=1-(exp(-0.002.*t).*cos((0.099).*t)+exp(-0.002.*t).*sin((0.099).*t).*0.02);
h1=h_u+delta/100;
h2=h_u-delta/100;
line([0 3000],[1 1],'Color','k')
line([0 3000],[h1 h1],'Color','r','LineStyle','--')
line([0 3000],[h2 h2],'Color','r','LineStyle','--')
plot(t,h)
xlabel('t')
ylabel('h(t)')
Упражнение 2:
clear,clc
f=tf([10],[1 0.4 11])
bode(f)
grid on
Упражнение 3:
clear,clc,close all
k = 0.1;
T = 1.4;
ksi = 0.2;
T2 = 0.35;
delta = 0.01;
w1 = tf([k 0],[T 1])
w2 = tf([1],[T2^2 2*T2*ksi 1])
w = w1*w2
ltiview({'step'},w);
>> h_u = 0;
>> h1=h_u+delta;
>> h2=h_u-delta;
>> figure(1) (график ltiview с Print to figure)
>> hold on
>> line([0 12],[1 1],'Color','k')
>> line([0 12],[h1 h1],'Color','r','LineStyle','--') >> line([0 12],[h2 h2],'Color','r','LineStyle','--')