function [ rms ] = RMS( s )
% Calculate the root mean square 
% Input
%   s = digital signal
% Output
%   rms = root mean square of s
rms = sqrt( sum( s .* s ) / length(s) );
end

