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

%
%-- 4.Временной интервал:
%

   t = 0.000:0.001:10.00;

%
%-- 5.Функция переходной проводимости:
%
   if R >= 1 %-- Апериодическое звено:
      Lambda = 1./C-A2.*exp(-Beta.*t)+A1.*exp(-Alpha.*t);
   else      %-- Колебательное звено:
      Lambda = -1./2.*exp(-1./2.*(B-(B.^2-4.*C.*A).^(1./2))./...
               A.*t).*(B+(B.^2-4.*C.*A).^(1./2))./(B.^2-...
               4.*C.*A).^(1./2)./C+1./2.*exp(-1./2.*(B+(B.^2-...
               4.*C.*A).^(1./2))./A.*t).*(B-(B.^2-4.*C.*A).^...
               (1./2))./(B.^2-4.*C.*A).^(1./2)./C+1./C;
   end

   subplot(4,2,1)
   plot(t,Lambda,'r')
   xlabel('t')
   ylabel('Lambda')

%
%-- 6.Функция веса:
%

   if R >= 1   %-- Апериодическое звено:
      Weight   = A2.*Beta.*exp(-Beta.*t)-A1.*Alpha.*exp(-Alpha.*t);
   else        %-- Колебательное звено:
      Weight   = 1./4.*(B-(B.^2-4.*C.*A).^(1./2))./A.*exp(-1./2.*...
                 (B-(B.^2-4.*C.*A).^(1./2))./A.*t).*(B+(B.^2-...
                 4.*C.*A).^(1./2))./(B.^2-4.*C.*A).^(1./2)./C-...
                 1./4.*(B+(B.^2-4.*C.*A).^(1./2))./A.*...
                 exp(-1./2.*(B+(B.^2-4.*C.*A).^(1./2))./A.*t).*(B-...
                 (B.^2-4.*C.*A).^(1./2))./(B.^2-4.*C.*A).^(1./2)./C;
   end

   subplot(4,2,2)
   plot(t,Weight,'r')
   xlabel('t')
   ylabel('Weight')


                                        114