Average of multiple matrices to create a new matrix with the same dimension.

Question

Paulo Eduardo Beiral 2019-10-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/484458-average-of-multiple-matrices-to-create-a-new-matrix-with-the-same-dimension

评论： meghannmarie 2019-10-10

Hi all!!

I want to take a total of 2416 variables named T (2112x2688), find the average of element by element of all the 2416 variables to create a single matrix with the avarege values.

I tried the following code, and I was able to get a single matrix with values, but those values doesn't correspond to the true averege values (I proved it myself with a test using only 10 T variables).

I attacched 2 example .mat files if anyone want to check out.

Thank's very much for any help!!!

Regards, Paulo Beiral.

diret = '/data/processamento_beiral/2';
cd (diret);
files = './*.mat';
dir (files);
arrumado = dir(files);
all = [1 355];
%% Average
for i = 1:length(arrumado)
    if str2double(arrumado(i).name(6:8)) >= all(1) && str2double(arrumado(i).name(6:8)) <= all(2) 
        turb(i) = load(arrumado(i).name);
        x = turb(:).T;
        x(isnan(x)) = 0;
        meanMatrix = x/2416;
    end
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

meghannmarie 2019-10-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/484458-average-of-multiple-matrices-to-create-a-new-matrix-with-the-same-dimension#answer_395576

If I am understanding your question right, you are wanting to find average of each element and ignore the nans?

diret = 'New Folder';
cd (diret);
files = './*.mat';
dir(files);
arrumado = dir(files);
all = [1 355];
%% Average
for i = 1:length(arrumado)
    if str2double(arrumado(i).name(6:8)) > all(1) && str2double(arrumado(i).name(6:8)) <= all(2) 
        turb(i) = load(arrumado(i).name);
    end
end
x = {turb.T};
x = cellfun(@(x) shiftdim(x,-1),x,'UniformOutput',false)';
x = cell2mat(x);
meanMatrix = squeeze(mean(x,'omitnan'));

2 个评论
显示无隐藏无

Paulo Eduardo Beiral 2019-10-10

在 MATLAB Online 中打开

Thank's very much meghannmarie!!

That's exactly what i needed, but there's still a doubt in my mind about the following line:

x = cellfun(@(x) shiftdim(x,-1),x,'UniformOutput',false)';

I read the matlab help about cellfun and shiftdim, but could you explain why do they have to be used in this case?...and specifically, why the "@" before the "(x)"?

Again, thank you for helping!

meghannmarie 2019-10-10

So we are shifting the dimension on each matrix so each resulting matrix wil be the size of 1x2112x2688. Now we can turn the cell array into a matrix and the matrixes from cells will be concatonated on that first dimension. This sets us up to get the mean on that dimension.

To apply the shiftdim to all the cells, I use cellfun. Since shiftdim requires 2 inputs (the matrix and -1), I use and anonymous function @(x) shiftdim(x,-1) - (you can look anonymous function up in the help). The anonymous function is applying shiftdim(x,-1) to each cell. If I used @shiftdim, I wouldnt be feeding the second input to the function.

请先登录，再进行评论。