In general, if you have a matrix A and want to take the mean across a certain dimension, dim, you should use mean(A,dim).
dim = 1 gives the mean for each column, dim = 2 gives the mean for each row, and so on.
Here's an example for your case:
M = rand(3621,1000);
meanOfEachRow = mean(M,2)