% Example of power spectrum, including normalisation of data by root Hz.  

[intenf,intenmag] = HPpower('tmpsdf\new\inten.dat'); % read data

% Calculate normalisation matrix, imaginatively called 'normalisation'
for n = 1:(length(intenf) - 1)
    normalisation(n) = sqrt(intenf(n+1) - intenf(n))
end
    normalisation(length(intenf)) = normalisation(length(intenf)-1)

% Divide signal by normalisation matrix; 0.672 was the maximum value in this 
% case and is used to get relative intensity values. 
intenvolts = intenmag ./ normalisation' ./ 0.672;
figure(16);
loglog(intenf, intenvolts);
grid on;
title('Relative intensity fluctuations in MC transmitted beam, normalised by root Hz')
xlabel('Frequency / Hz')
ylabel('Normalised intensity signal / Relative intensity (rms) per root Hz')