outputconv

s = [3 2 5 7 9 8 17 2];
h = [2 0 1 4];

% zero out the output value
y = zeros(1, length(s)+length(h)-1);

for i=1:length(y)
    for j=1:length(h)
        if (i-j+1 < 1)
            continue;
        end
        if (i-j+1 > length(s))
            continue;
        end
        y(i) = y(i) + h(j) * s(i-j+1);
    end
end
disp(y);