Correlations of matrix
matrix1 = rand(10, 100);
matrix2 = rand(10, 100);
for ind1=1:size(matrix1, 2)
    for ind2=1:size(matrix2, 2)
        [xc{ind1, ind2},lags{ind1, ind2}] = xcorr(matrix1(:,1),matrix2(:,2));
        % Here you can add any other operation you want
    end
end