Conv2 explanation for specific example

4 vues (au cours des 30 derniers jours)
Leo
Leo le 26 Juil 2011
Hey Guys, i need some help for the conv2 implementation for given problem:
K>> conv2([1 2 3], [1 0 0], 'full')
ans =
1 2 3 0 0
Why are the zeros on the right side? Does varies with the border depending on the kernel-center?
For this example its clear:
K>> conv2([1 2 3], [1 0 1], 'full')
ans =
1 2 4 2 3
Thanks for your help, Leo

Réponses (4)

Jan
Jan le 26 Juil 2011
Because you have vectors, your example is equivalent to:
conv([1, 2, 3], [1, 0, 0])
The algorithm of CONV is described exhaustively in: doc conv
E.g. the last element of the result is "w(2*n-1) = u(n)*v(n)", which is 0 in your example. Perhaps your are searchng for the 'same' or 'valid' shapes?
  2 commentaires
the cyclist
the cyclist le 26 Juil 2011
I noticed with a bit of amusement that conv() calls conv2().
Jan
Jan le 26 Juil 2011
Indeed?
In Matlab 6.5 it called FILTER after reordering the inputs: CONV(A, B) equals CONV(B, A), but FILTER(B, 1, A) is faster, if length(A) < length(B).
I think I should read through all of the toolbox functions again...

Connectez-vous pour commenter.


the cyclist
the cyclist le 26 Juil 2011
Isn't the convolution calculation essentially doing what I have written below? Does that make it clearer where the zeros come from?
a = [1 2 3]
b = [1 0 0]
conv_a_b = b(1)*[a 0 0] + b(2)*[0 a 0] + b(3) * [0 0 a]

Leo
Leo le 26 Juil 2011
thank you for your really quick responses! @Jan: yeah, i read conv2, but there is no information about how MATLAB set the paddings. In my work I have to apply a shift_kernel on a matrix to center the kernel, and the shif-kernel could be a matrix too. @cyclist: not really.. should 'a' replaced by [1 2 3] -> b(1) * [1 2 3 0 0] + ... ?
best, leo
  1 commentaire
the cyclist
the cyclist le 26 Juil 2011
In response to what you said to Jan: There is no "padding". MATLAB is treating the trailing zeros in "b" exactly how it would treat them if they were non-zero, and giving you the result.
In response of what you said to me: Yes. What I typed is functional code, using the "a" that you defined.

Connectez-vous pour commenter.


Leo
Leo le 2 Août 2011
Hey guys,
thanks for your answers. I had an error in reasoning. I forgot that the kernel should be mirrored.
Best!

Catégories

En savoir plus sur Performance and Memory dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by