%% HW3 Problem 4 Part E. Eric Aliotta. 01.25.2017

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