%Z wykladu
function y=splot(x1,x2,p1,p2)

%oblicz minimalna dlugosc wektora czasu dyskretnego
len=length(x1)+length(x2); %wektor dlugosci sygna?ów ograniczonych
minlen=len-1;
mv=p1+p2
dt=(p1+p2):minlen-1+(p1+p2);               % skala dyskretnego czasu
% aby ?atwiej rysowa? i liczy? dope?nij wektory zerami
x11=zeros(1,minlen-length(x1));
x1=[x1 x11];
x22=zeros(1,minlen-length(x2));
x2=[x2 x22];
y=zeros(1,minlen);
for n=1:minlen          %petla zewnetrzna
    for k=1:n           % pętla wewn?trzna
        y(n)=y(n)+x1(k)*x2(n-k+1);
    end
end
%stem(dt,y)
subplot(4,1,1); stem(dt-mv+p1,x1)
title('Sygnal pierwszy')
subplot(4,1,2); stem(dt-mv+p2,x2)
title('Sygnal drugi')
subplot(4,1,3:4); stem(dt,y); grid on
title('Splot sygnalów')