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

   w=0.01:0.001:1;
%
%-- 2.Частотная функция:
%
   Kiw=1./((i.*w).*(T.*i.*w+1));
%
%-- 3.Модуль частотной функции:
%
   subplot(2,2,1)
   plot(w,abs(Kiw),'r')
   xlabel('w')
   ylabel('abs(Kiw)')
%
%-- 4.Фаза частотной функции:
%
   subplot(2,2,2)
   plot(w,angle(Kiw),'r')
   xlabel('w')
   ylabel('angle(Kiw)')
%
%-- 5.Годограф на плоскости:
%
   subplot(2,2,3)
   plot(real(Kiw),imag(Kiw),'r')
   xlabel('real(Kiw)')
   ylabel('image(Kiw)')
%
%-- 6.Годограф в трехмерном пространстве (комета):
%
   subplot(2,2,4)
   comet3(real(Kiw),imag(Kiw),w)
   xlabel('real(Kiw)')
   ylabel('image(Kiw)')
   zlabel('w')
%
%-- Конец функции RKiw(T).

function [t1,t2] = Param(T,k,Umax,Xn)
%
%-- ФУНКЦИЯ ДЛЯ РАСЧЕТА ПАРАМЕТРОВ МАТЕМАТИЧЕСКИХ
%--                          МОДЕЛЕЙ:
%
%-- T = 0.62;           %-- постоянная при x'';
%-- k = 0.0023;         %-- постоянная при управляющем воздействии u(t);
%-- Umax = 220;         %-- максимальное значение управляющего
                                    131