looping through 2 .mat files

1 次查看（过去 30 天）

C.G. 2021-10-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1573043-looping-through-2-mat-files

评论： Mathieu NOE 2021-11-19

I have 2 matlab files with data in. Im trying to get my code to loop through both files and put the data into the 6 pre-assigned matrix.

I am getting the error:

Unable to perform assignment because the left and right sides have a different number of elements.
timestep = height(vel);
time_sec = timestep/100;
inputrate = 1;
%% Files  
files = dir('*.mat');
num_files = length(files);
%% Data output
resGT = zeros(100000,2);
GTflucc = zeros(100000,2);
mean_GT = zeros(numel(files),1);
std_GT = zeros(numel(files),1);
max_GT = zeros(numel(files),1);
min_GT = zeros(numel(files),1);
for a = 1:100000
load(files(a).name)
v = vel;
resGT(a) = sqrt(v(:,1).^2 +v(:,2).^2 + v(:,3).^2);
GT2 = resGT.^2;
GTflucc(a) = sqrt(cumsum(GT2)./time_sec);
%stats
mean_GT(a) = mean(resGT);
std_GT(a) = std(resGT);
max_GT(a) = max(resGT);
min_GT(a) = min(resGT);
end 

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

请先登录，再进行评论。

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

回答（1 个）

Mathieu NOE 2021-10-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1573043-looping-through-2-mat-files#answer_817963

在 MATLAB Online 中打开

hello

had to modify a bit your code to make it work

try this :

clc
clearvars
% timestep = height(vel); % are not yet defined (after load only)
% time_sec = timestep/100; % idem
inputrate = 1;
%% Files  
files = dir('*.mat');
num_files = length(files);
%% Data output
resGT = zeros(100000,2);
GTflucc = zeros(100000,2);
mean_GT = zeros(numel(files),1);
std_GT = zeros(numel(files),1);
max_GT = zeros(numel(files),1);
min_GT = zeros(numel(files),1);
for a = 1:num_files
    load(files(a).name)
    v = vel;
    timestep = height(vel);
    time_sec = timestep/100;
    resGT = sqrt(v(:,1).^2 +v(:,2).^2 + v(:,3).^2);
    GT2 = resGT.^2;
    GTflucc = sqrt(cumsum(GT2)./time_sec);
    %stats
    mean_GT(a) = mean(resGT);
    std_GT(a) = std(resGT);
    max_GT(a) = max(resGT);
    min_GT(a) = min(resGT);
end 