#39;k_1', 'k_2', 'k_3', 'k_4', 'k_5');
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
for k = 1: len_K
A_ (:,:, k) = A + B * K (k, :);
end
size (A_);
% -------------------------------------------------- ----------------------%
% Знаходження фазових координат
X (:, 1) = X_0;
h = 0.01;
time_X (1) = 0;
for k = 1: len_K
X (:, k +1) = X (:, k) + h * A_ (:,:, k) * X (:, k);
time_X (k +1) = time_X (k) + h;
end
X (:, k +1) = [];
time_X (k +1) = [];
% -------------------------------------------------- ----------------------%
% Знаходження управління
for k = 1: len_K
u (k) = K_pr (k, :) * X (:, k);
end
% -------------------------------------------------- ----------------------%
% Побудова u (t) і X (t)
figure (3);
plot (time_X, u, 'r-', 'LineWidth', 2)
title ('u (t)');
xlabel ('t')
grid on
В
figure (4);
plot (time_X, X (1, ​​:), 'LineWidth', 2)
hold on
title ('x_1 (t)');
xlabel ('t')
grid on
В
figure (5);
plot (time_X, X (2, :), 'LineWidth', 2)
title ('x_2 (t)');
xlabel ('t')
grid on
В
figure (6);
plot (time_X, X (3, :), 'LineWidth', 2)
title ('x_3 (t)');
xlabel ('t')
grid on
В
figure (7);
plot (time_X, X (4, :), 'LineWidth', 2)
title ('x_4 (t)');
xlabel ('t')
grid on
В
figure (8);
plot (time_X, X (5, :), 'LineWidth', 2)
title ('x_5 (t)');
xlabel ('t')
grid on
В
save Sravnenie_stabilizacii_2 time_X X u
Sravnenie_stabilizacii.m
close all
В
load Sravnenie_stabilizacii_1 X21 X22 X23 X24 X25 uu1
load Sravnenie_stabilizacii_2 time_X X u
В
figure (31);
plot (time_X, u, time_X, uu1, 'LineWidth', 2)
title ('u (t)');
xlabel ('t')
hl = legend ('u (t) - Управління з змін. коеф. ',' u (t) - управління з пост. коеф. ');
set (hl, 'FontName', 'Courier');
grid on
В
figure (41);
plot (time_X, X (1, ​​:), time_X, X21, 'LineWidth', 2)
hold on
title ('x_1 (t)');
xlabel ('t')
hl = legend ('x_1 (t) - З змін. коеф. ',' x_1 (t) - з пост. коеф. ');
set (hl, 'FontName', 'Courier');
grid on
В...
