title('S\{x1[n] + x2[n]\}');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 32 0 4]);

% b)
figure;
subplot(2,2,1);
plot(x1);
title('x1[n]');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 33 -1 4]);
subplot(2,2,2);
plot(x2);
title('x2[n]');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 33 -1 4]);

bSx1=x1+1;
bSx2=x2+1;
S3=bSx1+bSx2; %S{x1[n]} +S{x2[n]} suma odpowiedzi sygnalow
subplot(2,2,3);
plot(S3);
title('S\{x1[n]\} +S\{x2[n]\}');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 32 -1 4]);


S4=(x1+x2)+1; % S{x1[n] + x2[n]} odpowiedz sumy sygnalow
subplot(2,2,4);
plot(S4);
title('S\{x1[n] + x2[n]\}');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 32 0 4]);

% c)
figure;
subplot(2,2,1);
plot(x1);
title('x1[n]');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 33 -1 4]);
subplot(2,2,2);
plot(x2);
title('x2[n]');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 33 -1 4]);
cSx1=zeros(1,33);
cSx2=zeros(1,33);
S5=zeros(1,33);

for i=1:31
    cSx1(i)=x1(i+1)-x1(i);
    cSx2(i)=x2(i+1)-x2(i);
    S5(i) = cSx1(i)+cSx2(i);
end
%S5(32)=0.2;
 %S{x1[n]} +S{x2[n]} suma odpowiedzi sygnalow
subplot(2,2,3);
plot(S5);
title('S\{x1[n]\} +S\{x2[n]\}');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 30 -0.3 0.3]);
sygnal=zeros(1,32);
for i=1:31
    sygnal(i)=(x1(i+1)+x2(i+1))-(x1(i)+x2(i));
end
S6=(x1+1)-x1+(x2+1)-x2; % S{x1[n] + x2[n]} odpowiedz sumy sygnalow
subplot(2,2,4);
plot(sygnal);
title('S\{x1[n] + x2[n]\}');
xlabel('Numer próbki');
ylabel('Amplituda');
axis([0 30 -0.3 0.3]);