I need help with the following MATLAB code. I want to pull out $8$ matrices from the nested for loop. I want four matrices for each noise level for guassian and possion distribution defined to a variable. How do I do that?

$\%$ Number of unit vectors

n $=10$;

$\%$ Define noise levels in degrees

noise$\_$levels$\_$deg $=[1/3600,10/3600,1,5]$; $\%$ arcseconds and degrees

$\%$ Define a random rotation matrix $($Direction Cosine Matrix$)$

R $=$ orth$($randn$(3))$; $\%3$x$3$ orthogonal matrix

$\%$ Generate $10$ random unit vectors in the body frame

unit$\_$vectors$\_$body $=$ randn$(3,$ n$)$;

unit$\_$vectors$\_$body $=$ unit$\_$vectors$\_$body $./$ vecnorm$($unit$\_$vectors$\_$body$)$;

$\%$ Define noise distributions $($Gaussian and Poisson$)$

for noise$\_$type $=["$Gaussian$",$ "Poisson"$]$

disp$([\text{'}$Noise Type: $\text{'},$ char$($noise$\_$type$)])$;

for i $=1$:length$($noise$\_$levels$\_$deg$)$

noise$\_$level $=$ noise$\_$levels$\_$deg$($i$)$;

$\%$ Convert noise level to radians

noise$\_$level$\_$rad $=$ deg$2$rad$($noise$\_$level$)$;

disp$([\text{'}$Noise Level: $\text{'},$ num$2$str$($noise$\_$level$),\text{'}$ degrees'$])$;

$\%$ Initialize arrays to hold noisy vectors

noisy$\_$vectors$\_$body $=$ zeros$(3,$ n$)$;

noisy$\_$vectors$\_$inertial $=$ zeros$(3,$ n$)$;

for j $=1$:n

$\%$ Apply noise to the unit vectors

if noise$\_$type $==$ "Gaussian"

noise $=$ noise$\_$level$\_$rad $*$ randn$(3,1)$; $\%$ Gaussian noise

else

noise $=$ noise$\_$level$\_$rad $*$ poissrnd$(1,3,1)$; $\%$ Poisson noise

end

$\%$ Perturb the unit vector's direction

perturbed$\_$vector $=$ unit$\_$vectors$\_$body$($:$,$ j$)+$ noise;

perturbed$\_$vector $=$ perturbed$\_$vector $/$ norm$($perturbed$\_$vector$)$; $\%$ Renormalize

$\%$ Store noisy vector in body frame

noisy$\_$vectors$\_$body$($:$,$ j$)=$ perturbed$\_$vector;

$\%$ Transform to inertial frame using rotation matrix

noisy$\_$vectors$\_$inertial$($:$,$ j$)=$ R $*$ perturbed$\_$vector;

end

$\%$ Display the noisy vectors in the inertial frame for this noise level

disp$(\text{'}$Noisy Vectors in Inertial Frame:$\text{'})$;

disp$($noisy$\_$vectors$\_$inertial$)$;

end

end