Proportionalspalten in Matrix finden

Ich habe eine große Matrix (1.000 Zeilen und 50.000 Spalten). Ich weiß, dass einige Spalten korreliert sind (der Rang ist nur 100) und ich vermute, dass einige Spalten sogar proportional sind. Wie finde ich solche proportionalen Spalten? (Ein Weg wäre eine Schleifecorr(M(:,j),M(:,k))), aber gibt es etwas effizienteres?