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

                    %-- воздействия u(t);
%-- Xn = 2.09;      %-- конечное значение управляемой величина x(t);
%
%-- 1.Расчет параметров математических моделей:
%
   z = exp(Xn/(k*Umax*T))%-- промежуточный расчетный параметр;
   y = z + sqrt(z^2-z)      %-- промежуточный расчетный параметр;
   t1 = T*log(y)            %-- время переключения управления;
   t2 = 2*t1-Xn/(k*Umax) %-- время перехода системы в конечную
                            %-- точку;
%
%-- 2.Конец функции Param(T,k,Umax,Xn).

function [Lambda, Weight,XSin]= LWSin(T,a0,w)
%
%-- ГРАФИКИ ПЕРЕХОДНЫХ ФУНКЦИЙ Lambda(t), Weight(t) и
%--                            Sin(A, B, C):
%
%-- 1.Временной интервал:
%
    t = 0:1:10;
%
%-- 2.Функция переходной проводимости:
%
    Lambda = -T+t+T.*exp(-t./T);
   subplot(4,2,1)
    plot(t,Lambda,'r')
    xlabel('t')
    ylabel('Lambda')
%
%-- 3.Функция веса:
%
    Weight = 1-exp(-t./T);
    subplot(4,2,2)
    plot(t,Weight,'r')
    xlabel('t')
    ylabel('Weight')
%
%-- 4.Реакция звеньев на синусоидальные воздействия:
%
      XSin = (a0.*w./T).*(T./w.^2-(T.^3./...
                (1+w.^2.*T.^2)).*exp(-t./T)-...
                (1./(w.^2.*(w.^2+1./T.^2))).*...
                (w.*sin(w.*t)+1./T.*cos(w.*t)));
    subplot(4,2,3)
                                  132