Приложение 3. Содержимое файла data_k.m:

function [R, Q, area] = DATA(a)

% a(a1, a2)  - вектор варьируемых параметров системы

% R, Q  - числитель и знаменатель передаточной функции

% Задание области вариации параметров a1,a2. Формат: [минимум, шаг, максимум]

area = [ 0 0.01 1; 0 0.001 0.2 ];

% Варьируемые параметры можно переобозначить

k=a(1);

T=a(2);

% Задание коэффициентов числителя и знаменателя п.ф., зависящих от параметров

% Записываются по убыванию степеней

R = [100*k*T, 100*k];                         % Числитель передаточной функции

Q = [4, 9, 9+100*k*T, 1+100*k];               % Характеристический полином

Приложение 4. Matlab, Листинг. Программа построения фазовых траекторий и передаточных характеристик в среде.

% function_ode.m

function dXdt = fun_ode(t, X)

    global T01;

    global T02;

    global K0;

    x = X(1);

    y = X(2);

    dXdt(1) = y;

    dXdt(2) = ( K0*Saturation( x ) - x - (T01+T02)*y )/(T01*T02);   

    dXdt = dXdt(:);   

end

% main.m

global T01;

    global T02;

    global K0;

    T01 = 2; T02 = 6; K0 = 10; Kos = 0.5;

    t0 = 0; tf = 15;  

figure(1)   

    X0 = [1; 1];

    [t, X] = ode45('fun_ode', [t0 tf] , X0);   

    plot(X(:, 1), X(:, 2));

    hold on

    X0 = [-1; -1];

    [t, X] = ode45('fun_ode', [t0 tf] , X0);   

    plot(X(:, 1), X(:, 2));

figure(2)

    plot(t, X(:, 1));

% saturation.m