Most optimization functions in Octave such as fminunc expect vectors the parameter to be optimized. However, sometimes (e.g. in the field of neural networks) the vector representation is somewhat unnatural and clumsy to work with. Converting between vector and matrix representation is referred to as unrolling (matrix -> vector) and reshaping (vector -> matrix).
Unrolling
Unrolling is quite simple. Given a matrix A, you have two different possibilities of addressing the elements of A. Either you address rows and colums separately, i.e., A(i,j) or you use one single index or index range A(k) which indexes A column by column (Remember that the colon stands for the range from 1 to end):
|
1 2 3 4 5 6 7 8 9 10 11 12 |
>> A = [1 2 3; 4 5 6] A = 1 2 3 4 5 6 >> v = A(:) v = 1 4 2 5 3 6 |
As you can see, A(:) is exactly what we want: A vector representation of the matrix.
Reshaping
The reverse way, reshaping, is performed by a command called reshape. It takes as arguments the vector to be reshaped and the row and column count of the target matrix. Continuing the example from above, we get the following:
|
1 2 3 4 5 |
>> Ares = reshape(v, 2, 3) Ares = 1 2 3 4 5 6 |
Unrolling and reshaping with multiple matrices
Say you have multiple matrices, A,B, and C. Then the described procedure works analogously:
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 |
>> A = [1 1; 1 1]; B = [2 2]; C = [3;3]; >> v = [A(:); B(:); C(:)] v = 1 1 1 1 2 2 3 3 >> Ares = reshape(v(1:4), 2, 2) Ares = 1 1 1 1 >> Bres = reshape(v(5:6), 1, 2) Bres = 2 2 >> Cres = reshape(v(7:8), 2, 1) Cres = 3 3 |