function [A,B]=ex(C,D)
LN=length(C);
SN=length(D);
A=zeros(1,ceil((LN-2*SN+2)/2));
B=zeros(1,floor((LN-2*SN+2)/2));
for i=1:2:LN-2*SN+2
    for j=1:SN;
        A((i+1)/2)=A((i+1)/2)+C(i-2+2*j)*D(j);
    end;
end;
for i=2:2:LN-2*SN+2
    for j=1:SN;
        B(i/2)=B(i/2)+C(i-2+2*j)*D(j);
    end;
end;