Это зависит от размеров ваших кватернионных данных. Если один кватернион равен (4x1), то:
orientation = zeros(4,n)
for ii=1:size(acc,1)
% quaternion data
qahrs = ifilt(acc(ii,:), gyro(ii,:), mag(ii,:));
% supposed to be an array of quaternions
orientation(:,ii) = qahrs;
end
и если это (1x4)
orientation = zeros(n,4)
for ii=1:size(acc,1)
% quaternion data
qahrs = ifilt(acc(ii,:), gyro(ii,:), mag(ii,:));
% supposed to be an array of quaternions
orientation(ii,:) = qahrs;
end
Я предположил, что n - это число кватернионов.