Satellites in Orbit

%% Satellites in Orbit
% Assessment:  08
% Calculates the angular momentum, rotational kinetic energy,
%   and the angle between those two values for problem 1. Calculates the
%       specific angular momentum for problem 2.

%% Problem 1

I = [5000 0 0               %[kg-m^2] Moment of inertia
    0 5000 0
    0 0 9000];
omega = [0.2
    0.4
    2.3];                       %[rad/s] Angular velocity
omegaT = omega';
H = I * omega;                  %[kg-m^2/s]Angular momentum
Trot = .5 * omegaT * H;             %[J] Rotational kinetic energy
theta = acos(dot(H, omega)/(norm(H)*norm(omega)));        % [rad] Angle between the Angular velocity and momentum vectors
fprintf('The angular momentum is: %i kg*m^2 \n', H)
fprintf('The rotational kinetic energy is: %e Joules \n', Trot)
fprintf('The angle between the angular velocity and momentum is: %i rad \n', theta)

%% Problem 2

r = [1150               %[km] Position vector of the satellite
    9700
    -2400];
v = [-6.49              %[m/s] Velocity vector of the satellite
    -1.28
    1.65];
h = cross(r,v);          %[km^2/s] Specific angular momentum vector
fprintf('The specific angular momentums are: %i1 km^2/s, %i2 km^2/s, and %i3 km^2/s', h)