polar (fi_w, A_w, 'k');
grid on
xlabel ('Re (W (jw))')
ylabel ('Im (W (jw))')
title ('Function AFCHX (fi_w, A_w)')
% ---------------------------------------------- ---------------------------%
В
% ---------------------- Побудова ЛАЧХ ---------------------- --------------%
figure ('Name', '[0,100]');
w = -100: 0.01: 100;
LA_w = 20 * log (sqrt ((b0 ^ 2 + b1 ^ 2. * w. ^ 2)./((-a1 * w. ^ 2 + a3 * w. ^ 4). ^ 2 + (a0 * w-a2 * w. ^ 3 + a4 * w. ^ 5). ^ 2)));
plot (w, LA_w, 'k', 'LineWidth', 2);
grid on
xlabel ('w')
ylabel ('L (w)')
title ('Function L (w)')
% ---------------------------------------------- ---------------------------%
В
% ---------------------- Побудова ФАЧХ ------------------------------------%
% ---------------------------------------------- ---------------------------%
В
% ---------------------- Побудова h (t) ------------------------------------%
figure ('Name', '[0,50]');
t = 0: 0.01: 50;
h_t = 0.0024 * exp (-13.5832. * t) - 0.2175 * exp (-0.6848. * t) ...
+ 0.1452 * exp (-0.5850. * t). * cos (2.7677. * t) ...
- 0.2217 * exp (-0.5850. * t). * sin (2.7677. * t) ...
+ 0.5825. * t + 0.0699;
plot (t, h_t, 'k', 'LineWidth', 2);
grid on
xlabel ('t')
ylabel ('h (t)')
title ('Function h (t)')
% ---------------------------------------------- ---------------------------%
В
% ---------------------- Побудова k (t) ------------------- -----------------%
figure ('Name', '[0,50]');
t = 0: 0.01: 50;
k_t = - 0.0329 * exp (-13.5832. * t) + 0.1489 * exp (-0.6848. * t) ...
- 0.6986 * exp (-0.5850. * t). * cos (2.7677. * t) ...
- 0.2721 * exp (-0.5850. * t). * sin (2.7677. * t) ...
+ 0.5826;
plot (t, k_t, 'k', 'LineWidth', 2);
grid on
xlabel ('t')
ylabel ('k (t)')
title ('Function k (t)')
% ---------------------------------------------- ---------------------------%
В
x1 = tf ([b1 b0], [a4 a3 a2 a1 a0 0]);
ltiview (x1)
ProstranstvoSostoyanii.m
clc
clear all
В
% format rational
В
b1 = 9;
b0 = 5;
В
a5 = 0.1153;
a4 = 1.78;
a3 = 3.92;
a2 = 14.42;
a1 = 8.583;
a0 = 0;
В
% 1. Матриця Фробеніуса
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%
A = [0 1 0 0 0;
0 0 1 0 0;
0 0 0 1 0;
0 0 0 0 1;
0 -a1/a5 -a2/a5 -a3/a5 -a4/a5]
В
B = [0, 0, 0, 0, 1/a5]
В ...
