C = convn(A,B)
C = convn(A,B,'shape')
Returns the full N-dimensional convolution (default).
Returns the central part of the result that is the same size as A.
Returns only those parts of the convolution that can be computed without assuming that the array A is zero-padded. The size of the result is
max(size(A)-size(B) + 1, 0)