clear; close all;
Fs = 1600;% Częstotliwości próbkowania
Ts = 1/Fs;% Czas próbkowania
L = 128;% Liczba próbek
t_d = (0:L-1)*Ts;% Oś czasu - wartości dyskretne
t_c = linspace(0,t_d(end),100);    % Oś czasu - wartości ciągłe
A = [230, 115, 75, 35, 15];% Parametry sygnały - Amplituda
f = [50, 100, 250, 400, 800];% Parametry sygnały - Częstotliwości
Phi = [10, 20, 30, 40, 50];% Parametry sygnały - Przesunięcie
a1 = 2.5;
a0 = 3;
b1 = 2.5;
b0 = 3;

y_d = zeros(1:length(L));
x = zeros(1,L);
y_d(1) = 0;

y2_d=0;

for i = 1:5
    y_d = A(i)*sin(2*pi*f(i)*t_d+Phi(i));

    y2_d = y2_d + y_d;
end

for n = 2:length(t_d)
    y_d(n) = (a1/a0) * y_d(n-1) + (b1/b0) * y2_d(n);
end

% Wykres
figure
set(gcf,'DefaultLineLineWidth',2);
set(gcf,'DefaultAxesFontSize',14);
% stem(t_d,y_c,'s')
hold on;
stem(t_d,y_d,'s')
hold off
grid on
xlabel('Czas[s]')
ylabel('Amplituda')
legend('- d -')