function y = running_rms(x, fs, T)
%RUNNING_RMS Running rms approximation

Krms = T*fs;
alpha = exp(-1/Krms);
y = sqrt(filter(1/Krms, [1, -alpha], abs(x).^2));

end
