Working with a set of matrices

Question

Lorcan O'Connor 2021-4-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/791349-working-with-a-set-of-matrices

评论： Lorcan O'Connor 2021-4-3

I want to generate a finite matrix group by starting with the S ={identity} and multiplying by some generators until S no longer grows in size.

Is there an efficient way to deal with a "set" like this, without manualling checking for duplicates? For example, it seems awkward if you use cell arrays to store the elements of S.

2 个评论
显示无隐藏无

Matt J 2021-4-3

For example, it seems awkward if you use cell arrays to store the elements of S.

Why?

Lorcan O'Connor 2021-4-3

Well say I multiply everything in the cell array by a matrix and only want to take the new elements, how do I efficiently take a “set Union” (or even check if a matrix is in the array), without just looping and checking everything?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2021-4-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/791349-working-with-a-set-of-matrices#answer_665899

编辑：Matt J 2021-4-3

在 MATLAB Online 中打开

I'm not sure I will fully understand the task without seeing an example. However, if S is organized as an MxNxP array,

[M,N]=deal(3,3);
S=cat(3,eye(M,N), randi(5,M,N), eye(M,N) )
S = 
S(:,:,1) =

     1     0     0
     0     1     0
     0     0     1


S(:,:,2) =

     4     5     4
     2     1     3
     1     5     2


S(:,:,3) =

     1     0     0
     0     1     0
     0     0     1

then you can get rid of duplicate slices with

s=reshape(S,[],size(S,3));
S=reshape( unique(s.','rows').', M,N,[])
S = 
S(:,:,1) =

     1     0     0
     0     1     0
     0     0     1


S(:,:,2) =

     4     5     4
     2     1     3
     1     5     2