Моделирование САУ. Работа с пакетом Control System Toolbox. Навыки моделирования САУ с использованием М-файлов и в пакете Simulink, страница 3

grid;                               %..сетка

figure(11);                         %..создаем_графическое_окно

clf();                              %..очистка_графического_окна

bode(w4);                                                         %..График_диограммы_Боде_колебательного_звена

title('График диограммы Боде Колебательного звена');              %..служит_для_вывода_название_графика

grid;                               %..сетка

figure(12);                         %..создаем_графическое_окно

clf();                              %..очистка_графического_окна

nyquist(w4);                                                      %..График_диограммы_Найквиста_колебательного_звена

title('График диограммы Найквиста Колебательного звена');         %..служит_для_вывода_название_графика

grid;                               %..сетка

%..3. Для колебательного звена определить по АЧХ значение максимальной амплитуды и частоты, при которой она достигается.

[A_CH_X,CH_F_X,A_F_CH_X]=bode(w4);

[A,index_A]=max(A_CH_X(1,1,:));

disp(' ____________Результат_вычислений_________ ');          %..вывод_в_командное_окно

disp('|Максимальная амплитуда|      Частота     |');          %..вывод_в_командное_окно

fprintf('| %20.13g | %16.8g |\n',A,A_F_CH_X(index_A)); %..вывод_в_командное_окно_значений

disp('|______________________|__________________|');          %..вывод_в_командное_окно

%..4. С помощью обозревателя ltiview  построить временную и частотные характеристики апериодического и колебательного звена в одном окне.

ltiview(w3,w4);

%..5. Система описывается передаточной функцией вида: Необходимо

%подобрать такие значения A, B, C, чтобы система была а) устойчива, б) неустойчива. Построить график переходной характеристики и карту полюсов, подтверждающие правильность выбора коэффициентов.Сделать вывод о том, как влияет параметр K на вид переходной характеристики системы.

A1=7;                     %..Параметр_системы

B1=-15;                   %..Параметр_системы

C1=70;                    %..Параметр_системы

A2=-7;                    %..Параметр_системы

B2=-15;                   %..Параметр_системы

C2=-70;                   %..Параметр_системы

w_p1=tf(k,[A1 B1 C1]); %..Придаточная_функция_неустойчивой_системы

w_p2=tf(k,[A2 B2 C2]); %..Придаточная_функция_устойчивой_системы

figure(13);                         %..создаем_графическое_окно

clf();                              %..очистка_графического_окна

subplot(1,2,1);

step(w_p1);                         %..Графие_переходной_характеристики_придаточной_функции_w_p1

title('Графие переходной характеристики');         %..служит_для_вывода_название_графика

grid;                               %..сетка

subplot(1,2,2);

pzmap(w_p1);                                             %..Карта_полюсов_придаточной_функции_w_p1

title('Карта полюсов');             %..служит_для_вывода_название_графика

figure(14);                         %..создаем_графическое_окно

clf();                              %..очистка_графического_окна

subplot(1,2,1);

step(w_p2);                         %..Графие_переходной_характеристики_придаточной_функции_w_p2

title('Графие переходной характеристики');         %..служит_для_вывода_название_графика

grid;                               %..сетка

subplot(1,2,2);