clear; close all;
Fs = 50;% Częstotliwości próbkowania
Ts = 1/Fs;% Czas próbkowania
L = 36;% 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 = 10;% Parametry sygnały - Amplituda
f = 50;% Parametry sygnały - Częstotliwości
Phi = 0;% Parametry sygnały - Przesunięcie fazowe
a1 = 2.5;
a0 = 3;
b1 = 2.5;
b0 = 3;

y_d=zeros(1:length(L));
x=zeros(1,L);
y_d(1)=b1/b0;

for n = 2:length(t_d)
    y_d(n) = (a1/a0) * y_d(n-1) + (b1/b0)*x(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 -')