function [SNRF, SNRUV] = SNG(noOfRuns)
NoOfPhotons = 1000000;
Q = 0.6;
b = 1;
E = 10000;
c = 0.00001;
stD = 1000;
Mean = 1000000;
for i = 1:noOfRuns
extPhotons = NoOfPhotons + randn(noOfRuns,1) * stD;
F = 0.01*extPhotons*2.303*E*b*c*Q;
stD2 = sqrt(F);
detPhotonsF = F + randn(noOfRuns,1)* stD2;
mean_detPhotonsF = mean(detPhotonsF);
stdev = std(detPhotonsF);
SNGF (i) = mean_detPhotonsF/stdev
end
Best Answer