Второй сем элтеха / Novikov_CM9_71_eltex1_matlab

clear all; %#ok<CLALL>
close all;

K = 255000; lK = log10(K);
f = 2e6; lf = log10(f); % 2 МГц
R1 = 18e3; % 18 кОМ
R2 = 1.4e3; % 1.4 кОМ
R3 = 163e3; % 163 кОМ
C1 = 5e-6; % 5 мкФ
C2 = 26e-6; % 26 мкФ
C3 = 710e-12; % 710 пФ

T_0 = R1 * C1;
f_0 = 1 / (2 * pi * T_0); lf_0 = log10(f_0);

T_1 = C2 * (R2 * R3) / (R2 + R3);
T_2 = C3 * (R2 + R3);
T_3 = C2 * (R2 + R3);
T_4 = C3 * (R2 * R3) / (R2 + R3);
f_1 = 1 / (2 * pi * T_1); lf_1 = log10(f_1);
f_2 = 1 / (2 * pi * T_2); lf_2 = log10(f_2);
f_3 = 1 / (2 * pi * T_3); lf_3 = log10(f_3);
f_4 = 1 / (2 * pi * T_4); lf_4 = log10(f_4);

lf_c = lf - lK;
f_pos = [-3, lf_3, lf_0, lf_1, lf_c, lf_2, lf_4, lf, 7];
K_OY = [20 * lK, 20 * lK, 20 * lK, 20 * lK, 20 * lK, 20 * lK - 20 * (lf_2-lf_c), 20 * lK - 20 * (lf_4-lf_c), 0, 20 * lK - 20 * (7 - lf_c)];
K_OC = [0, 0, -20 * (lf_0-lf_3), -20 * (lf_1-lf_3), -20 * (lf_1-lf_3), -20 * (lf_1-lf_3), 0, 0, 0];
K_OC_mirror = [0, 0, 20 * (lf_0-lf_3), 20 * (lf_1-lf_3), 20 * (lf_1-lf_3), 20 * (lf_1-lf_3), 0, 0, 0];
K_RC = [20*(-lf_0-3), 20*(-lf_0+lf_3), 0, 0, 0, 0, 0, 0, 0];
fi = [90, 180, 90, 0, 0, -90, 0, -90, -90];

K_res = [K_RC(1) , K_RC(2), K_OC_mirror(3), K_OC_mirror(4), K_OC_mirror(5), ...
K_OC_mirror(6), K_OC_mirror(7), K_OC_mirror(8), K_OY(9)];

f1 = figure('Name', 'ЛАЧХ ОУ', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([0, 1.05*max(K_OY)]); % Границы на горизонтальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([0, 1.05*max(f_pos)]); % Границы осей на вертикальной оси.
plot(f_pos, K_OY);
plot([0, 1], [20, 0], "--"); % LineStyle - линия - пример из методического пособия

f2 = figure('Name', 'ЛАЧХ входной цепи', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([-8, 6]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([0, 1]); % % Границы на горизонтальной оси.
plot(f_pos, K_RC);

f3 = figure('Name', 'ЛАЧХ обратной связи', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([-50, 10]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([-3, 7]); % Границы на горизонтальной оси.
plot(f_pos, K_OC);

f4 = figure('Name', 'Зеркальная ЛАЧХ обратной связи', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([-10, 50]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([-3, 7]); % Границы на горизонтальной оси.
plot(f_pos, K_OC_mirror);

f5 = figure('Name', 'ЛАЧХ усилителя с учетом ООС', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([-20, 115]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([-2, 7]); % Границы на горизонтальной оси.
plot(f_pos, K_OC_mirror);
plot(f_pos, K_OY);
legend('ЛАЧХ ООС','ЛАЧХ ОУ');

f6 = figure('Name', 'ЛАЧХ активного фильтра', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([-70, 115]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([-3, 7]); % Границы на горизонтальной оси.
plot(f_pos, K_OC_mirror);
plot(f_pos, K_OY);
plot(f_pos, K_RC);
legend('ЛАЧХ ООС','ЛАЧХ ОУ', 'ЛАЧХ RC');

f7 = figure('Name', 'Суммарная ЛАЧХ активного фильтра', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Коэффициент усиления K, дБ');
ylim([1.05*min(K_res), 1.05*max(K_res)]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([-3, 7]); % Границы на горизонтальной оси.
plot(f_pos, K_res);

f8 = figure('Name', 'Суммарная ЛФЧХ активного фильтра', 'NumberTitle', 'off');
hold on;
grid on;
ylabel('Фаза \phi, градусы.');
ylim([-190, 190]); % Границы на вертикальной оси.
xlabel('Частота единичного усиления f_{ед}, дек.'); % TeX markup
xlim([-3, 7]); % Границы на горизонтальной оси.
stairs(f_pos, fi);