1
Download Diffuse Tomography of Absorbing and
Transcript
function [sys_val] = calc_repeatability(data_fn,num_reps,source,pixel);
% Used to calculate baselines to subtract from recorded spectra.
% s davis 2007
for i = 1:num_reps
% load basis data for one rep
data_temp = load([data_fn,sprintf('_s%g_rep%g.raw',source,i)]);
data_temp = data_temp(pixel,:);
data(i,:) = data_temp;
end
data_avg = median(data,1)';
data_dev = std(data,1,1)';
perc_error = 100*data_dev./data_avg;
SNR = 20*log10(data_avg./data_dev);
sys_val = [data_avg, data_dev, perc_error, SNR];
% find(data_avg >=
figure
semilogy(1:length(data_avg),data_avg)
title('averaged data')
figure
plot(1:length(data_dev),perc_error)
title('Percent error based on stdev')
figure
plot(1:length(SNR),SNR)
title('SNR')
382
