How to write a 'horizontal direct product' (*~) from Gauss in Matlab?

Hi,

I am trying to write the 'horizontal direct product' (*~) from Gauss in Matlab. Anyone can help me?

E.g. Horizontal direct product z = x*~y;

x = [1 2; 3 4]; y = [5 6; 7 8];

z = [5 6 10 12; 21 24 28 32]

NB. The input matrices x and y must have the same number of rows. The result will have cols( x) * cols( y) columns.

Thanks!

Answer 1

z = repmat(y,1,size(x,2)).*kron(x,ones(1,size(y,2)))

Answer 2

That is a subset of the output of kron(x,y)

Thanks, Walter. I noted that, but do you know how to extract this subset?

Row indexing, but I do not know what the general pattern is (I would need a few additional examples of different lengths to go on.)

Answer 3

Thanks, Matt. That's perfect.