Commit 6b501e0d authored by Mueller-Trapet's avatar Mueller-Trapet

include SNR in analytic result

parent 19f3b348
......@@ -129,7 +129,7 @@ if sArgs.analytic
% RT for analytic result
RT = ita_roomacoustics(ir,'T30','freqRange',[125 8000],'bandsperoctave',1);
RT = RT.T30.freq;
m_analytic = 1./sqrt(1 + (2*pi.*bsxfun(@times,fm,RT)./13.8).^2);
m_analytic = 1./(bsxfun(@times,sqrt(1 + (2*pi.*bsxfun(@times,fm,RT)./13.8).^2),(1+10.^(-SNR.'./10))));
m_analytic = min(bsxfun(@times,m_analytic,correctionTerm),1);
SNR_eff_analytic = min(max(real(10.*log10(m_analytic./(1-m_analytic))),-15),15);
TI_analytic = (SNR_eff_analytic + 15)./30;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment