clc,close all,clear,beep off;
% Czestotliwosci próbkowania
Fs = 50;
% Czas próbkowania
Ts = 1/Fs;
% Liczba próbek
L = 36;
% Os czasu - wartosci dyskretne
t_d = (0:L-1)*Ts;

% Sygna³ - wymuszenie
x1 = zeros(1,length(t_d));
x1(1) = 1;

% Sygna³ - wartosci osi Y
a1 = 2.5;
a0 = 3;
b1 = 2.5;
b0 = 3;

y_d1i = zeros(1,length(t_d));
y_d1i(1) = b1/b0*x1(1);

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

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

% Czestotliwosci próbkowania
Fs = 1;
% Czas próbkowania
Ts = 1/Fs;
% Liczba próbek
L = 20;
% Os czasu - wartosci dyskretne
t_d1 = (0:L-1)*Ts;

% Sygna³ - wymuszenie
x1 = t_d1 >= 10;

% Sygnal - wartosci osi Y
a1 = 2.5;
a0 = 3;
b1 = 2.5;
b0 = 3;

y_du1 = zeros(1,length(t_d1));
y_du1(1) = b1/b0*x1(1);

for n = 2:length(t_d1)
    y_du1(n) = a1/a0*y_du1(n-1) + b1/b0*x1(n);
end

% Wykres
figure
set(gcf,'DefaultLineLineWidth',2);
set(gcf,'DefaultAxesFontSize',14);
stem(t_d1,y_du1,'s')
grid on
xlabel('Czas[s]')
ylabel('Amplituda')
legend('- d -')