calibrate.m

计算声质中的Aures Tonality

function[ycal,err]=calibrate(y)
err=1;

Pref=20e-6;
RMS=sqrt(mean(y.^2));
SPLcalc=20*log10(RMS/Pref);
disp([10,'The RMS SPL, calculated as SPLmax=20*log10(RMS(y)/Pref) is', num2str(SPLcalc)])
SPLwant=input('Please enter the RMS SPL as determined during the measurement');
c=10^((SPLwant-SPLcalc)/20);ycal=c*y;
err=0;
return