Приложение 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