toe kursa4 4var+3var / toe_ebanina

% Данные ряда Фурье входного сигнала
amplitudes_input = [1.274, 1.273, 0.424, 0, 0.255, 0.424, 0.182, 0, 0.142];
phases_input = [0, -1.572, 3.14, 1.569, 0, -1.572, 3.14, 1.569, 0];
harmonic_numbers = [1, 2, 3, 4, 5, 6, 7, 8, 9];

% Базовая частота
T = 4;
omega0 = 2 * pi / T;

% Передаточная функция H(jw)
H_jw = @(omega) (1.75 * (1j * omega).^2) ./ ((1j * omega).^2 + 1.8 * 1j * omega + 2.07);

% Вычисляем |H(jkω)| и arg(H(jkω)) для каждой гармоники
H_magnitudes = zeros(1, length(harmonic_numbers));
H_phases = zeros(1, length(harmonic_numbers));

for k = 1:length(harmonic_numbers)
H_val = H_jw(harmonic_numbers(k) * omega0);
H_magnitudes(k) = abs(H_val);
H_phases(k) = angle(H_val);
end

% Новые амплитуды и фазы выходного сигнала
amplitudes_output = amplitudes_input .* H_magnitudes;
phases_output = phases_input + H_phases;

% Время для построения графика
t = linspace(0, 4, 4000);

% Восстановление выходного сигнала
y_t = zeros(size(t));

for k = 1:length(harmonic_numbers)
y_t = y_t + amplitudes_output(k) * cos(harmonic_numbers(k) * omega0 * t + phases_output(k));
end

% Построение графика выходного сигнала
figure;
plot(t, y_t, 'g', 'LineWidth', 1.5);
title('Выходной сигнал y(t) через разложение Фурье');
xlabel('Время, с');
ylabel('Амплитуда');
grid on;
legend('Выходной сигнал через ряд Фурье');