Решение задач оптимального управления с использованием математической системы MATLAB и пакета имитационного моделирования SIMULINK. Сивохин А.В - 110 стр.

UptoLike

110
Таким образом, звено демпфировано, если в его математическом описании
имеется составляющая, пропорциональная первой производной
)(t
λ
.
Составляющая, пропорциональная второй производной, характеризует
инерционные свойства звена. Для уравнения (4.38) эти составляющие
определяются соответственно коэффициентами
1
T и
2
2
T .
4.4 Программная реализация аналитических моделей
M-функция для построения модуля, фазы и годографов частотной функции
)(
ω
iK
на комплексной плоскости и в пространстве имеет вид:
function var= RKiw(A,B,C)
%
%-- ГРАФИКИ ЧАТОТНОЙ ФУНКЦИИ K(iw) :
%
%-- A =
2
2
T ;
%-- B =
1
T ;
%-- C =
δ
;
%-- R =
δ
ρ
2
1
2T
T
=
;
%
%-- 1.Диапазон частот:
%
w=0.01:0.001:10;
%
%-- 2.Коэффициент демпфирования:
%
R=B/(2*sqrt(A)*sqrt(C));
%
%-- 3.Вид колебательного звена:
%
if (R>1)
disp('Сильно-демпфированное звено')
end
if (R==1)
disp('Апериодическое предельно-демпфированное звено')
end
if (1>R && R>sqrt(2)./2)
disp('Нормально-демпфированное звено')
end
if (R==sqrt(2)./2)
disp('Критически-демпфированное звено')
end
    Таким образом, звено демпфировано, если в его математическом описании
имеется составляющая, пропорциональная первой производной             λ (t ) .
Составляющая, пропорциональная второй производной, характеризует
инерционные свойства звена. Для уравнения (4.38) эти составляющие
определяются соответственно коэффициентами T1 и T22 .

                  4.4 Программная реализация аналитических моделей


     M-функция для построения модуля, фазы и годографов частотной функции
K (iω ) на комплексной плоскости и в пространстве имеет вид:

    function var= RKiw(A,B,C)
    %
    %-- ГРАФИКИ ЧАТОТНОЙ ФУНКЦИИ K(iw) :
    %
    %-- A = T22 ;
    %-- B = T1 ;
    %-- C = δ ;
                   T1
    %-- R = ρ =           ;
                  2T2 δ
    %
    %-- 1.Диапазон частот:
    %
       w=0.01:0.001:10;
    %
    %-- 2.Коэффициент демпфирования:
    %
       R=B/(2*sqrt(A)*sqrt(C));
    %
    %-- 3.Вид колебательного звена:
    %
       if (R>1)
          disp('Сильно-демпфированное звено')
       end
       if (R==1)
          disp('Апериодическое предельно-демпфированное звено')
       end
       if (1>R && R>sqrt(2)./2)
          disp('Нормально-демпфированное звено')
       end
       if (R==sqrt(2)./2)
          disp('Критически-демпфированное звено')
       end

                                       110