calibrate.m

来自「计算声质中的Aures Tonality」· M 代码 · 共 12 行

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