clear all

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

fileName = 'DeploymentNov24-09\BTNoiseData\gndedJ1BT\R00000\D000000';
segLength = 2^16;
overLapPercent = 0;
fileSize = 48000 * 60 * 2;

hSpectrum = spectrum.welch("WindowName', '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(ADCData));

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