Localizando colunas proporcionais na matriz

Eu tenho uma grande matriz (1.000 linhas e 50.000 colunas). Sei que algumas colunas estão correlacionadas (a classificação é de apenas 100) e suspeito que algumas colunas sejam até proporcionais. Como posso encontrar essas colunas proporcionais? (uma maneira seria loopcorr(M(:,j),M(:,k))), mas há algo mais eficiente?