clear all

Fs = 48000;
dT = 1/Fs;
overlapPercent = 50;

fileName = 'H:\DeploymentNov24-09\BTNoiseData\gndedJ1\R00000\';
hSpectrum = spectrum.welch('Hamming', segLength, overlapPercent);
ADCData = zeros( fileSize, 1, 'int16');

fid = fopen(fileName);
fprintf(1, 'Reading %s \n', fileName);
ADCData = fread(fid, 'int16=>int16');
fclose(fid);

Vref = 1.5;
LSB = Vref / (2^15 - 1);
data = double(LSB * double(rawData));

rms = sqrt(mean(single(data).^2))
    
pxx = psd(hSpectrum, single(data), 'Fs', Fs);
Pxx(pxxIndex,:) = 10 * log10(pxx.Data);