load -mat SRS003.MAT
load -mat SRS005.MAT
load -mat SRS006.MAT
load -mat SRS007.MAT
load -mat SRS008.MAT

h = 6.626E-34;
mu = 3E8/1.604E-6;
deltaT = 1.0;

freq = srs005(:,1);
Gain = 0.75;
Respons = 0.7;
Po = 1.11/Gain/Respons

Gain = 0.75E3;
Po = 1.11/Gain/Respons;

sigmaE = sqrt(h*mu*Po*deltaT);

sigmaI = sigmaE*Respons;

sigmaV = sigmaI*Gain;

subplot(2,2,1)
shotnoise = ones(801,1)*sigmaV;
loglog(srs005(:,1),srs005(:,2),srs005(:,1),shotnoise,srs003(:,1),srs003(:,2))
hold on;

subplot(2,2,2)
Po = 2.12/Gain/Respons

sigmaE = sqrt(h*mu*Po*deltaT);

sigmaI = sigmaE*Respons;

sigmaV = sigmaI*Gain;
shotnoise = ones(801,1)*sigmaV;
loglog(srs006(:,1),srs006(:,2),srs006(:,1),shotnoise,srs003(:,1),srs003(:,2))
hold on;

subplot(2,2,3)
Po = 3.47/Gain/Respons

sigmaE = sqrt(h*mu*Po*deltaT);

sigmaI = sigmaE*Respons;

sigmaV = sigmaI*Gain;
shotnoise = ones(801,1)*sigmaV;
loglog(srs007(:,1),srs007(:,2),srs007(:,1),shotnoise,srs003(:,1),srs003(:,2))
hold on;

subplot(2,2,4)
Po = 4.86/Gain/Respons

sigmaE = sqrt(h*mu*Po*deltaT);

sigmaI = sigmaE*Respons;

sigmaV = sigmaI*Gain;
shotnoise = ones(801,1)*sigmaV;
loglog(srs008(:,1),srs008(:,2),srs008(:,1),shotnoise,srs003(:,1),srs003(:,2))
hold on;


grid on
legend('Total Noise','Shot Noise','Electronics Noise')
xlabel('Frequency (Hz)')
ylabel('Vrms/Hz^{0.5}')
title('Intensity Noise at 4.86 V')
axis([1 1000 10E-10 10E-4])


subplot(2,2,3)
grid on
legend('Total Noise', 'Shot Noise','Electronics Noise')
xlabel('Frequency (Hz)')
ylabel('Vrms/Hz^{0.5}')
title('Intensity Noise at 3.47 V')
axis([1 1000 10E-10 10E-4])


subplot(2,2,2)
grid on
legend('Total Noise','Shot Noise','Electronics Noise')
xlabel('Frequency (Hz)')
ylabel('Vrms/Hz^{0.5}')
title('Intensity Noise at 2.12 V')
axis([1 1000 10E-10 10E-4])


subplot(2,2,1)
grid on
legend('Total Noise','Shot Noise','Electronics Noise')
xlabel('Frequency (Hz)')
ylabel('Vrms/Hz^{0.5}')
title('Intensity Noise at 1.11 V')
axis([1 1000 10E-10 10E-4])