Algorithm for a commutation matrix

Hi everybody, need help to design a Matlab code function that creates a commutation matrix K of size (mn,mn) such that, for a (m,n) matrix A:

vec(A')=K vec(A)

where the vec operator is such that, for example, if:

Thanks for your help.

best,

Answer 1

[Luciano Campos]

I've worked this out. No further help needed.

Thanks anyway.

best,

Answer 2

[Walter Roberson]

You have not defined Commutation Matrix correctly. You have defined it in terms of Conjugate Transpose, but it is instead defined in terms of Transpose.

https://en.wikipedia.org/wiki/Commutation_matrix

[Luciano Campos]

Hi Walter,

Thanks for your answer, but I don't really see the difference between the way I defined the commutation matrix and how it's defined in Wikipedia. In any case, I'd really appreciate if you can help me coding it.

best,