figure;
count = 0;
for Nrot = [586, 292] % measured precession cycles
count = count +1;
noise_mean = 0;
noise_sd = 50;
% repeat experiment 10,000 times
Nrun = 1e4;
measured_sig = zeros(Nrun,1);
for j=1:Nrun
Mag = 100*ones(Nrot,1);
Noise = noise_sd*randn(Nrot,1) + noise_mean;
Sig = Mag + Noise;
measured_sig(j) = mean(Sig);
end
subplot(1,2,count);
plot(measured_sig);
ylim([80 120]); title([ int2str(Nrot) ' cycles per measurement']);
ylabel('Signal'); xlabel('Simulation Runs');
fprintf('%d Cycles \n',Nrot);
fprintf('Measured signal: %2.2f+/-%2.2f \n',mean(measured_sig),std(measured_sig));
fprintf('SNR : %2.2f\n\n',mean(measured_sig)/std(measured_sig));
end