Aufgrund einer Wartung wird GitLab am 28.09. zwischen 10:00 und 11:00 Uhr kurzzeitig nicht zur Verfügung stehen. / Due to maintenance, GitLab will be temporarily unavailable on 28.09. between 10:00 and 11:00 am.

Commit d5405809 authored by Markus Mueller-Trapet's avatar Markus Mueller-Trapet

added low frequency correction term

parent 944914f8
...@@ -15,7 +15,7 @@ sArgs = struct('pos1_data','anything','bandsperoctave',3,'freqVector',[],'create ...@@ -15,7 +15,7 @@ sArgs = struct('pos1_data','anything','bandsperoctave',3,'freqVector',[],'create
%% reference curves %% reference curves
if strcmpi(sArgs.type,'iso') % Reference curve and frequencies according to ISO 717-2 if strcmpi(sArgs.type,'iso') % Reference curve and frequencies according to ISO 717-2
outputStr = 'Ln_w (C_I)'; outputStr = 'Ln_w (C_I, C_{I50-2500})';
roundingFactor = 0.1; roundingFactor = 0.1;
deficiencyLimit = Inf; deficiencyLimit = Inf;
if sArgs.bandsperoctave == 1 if sArgs.bandsperoctave == 1
...@@ -81,7 +81,13 @@ end ...@@ -81,7 +81,13 @@ end
%% adaptation term for ISO %% adaptation term for ISO
if strcmpi(sArgs.type,'iso') if strcmpi(sArgs.type,'iso')
C = round(round((10.*log10(sum(10.^(NISPL./10))))/0.1)*0.1 - 15 - impactInsulationClass); C = round(round((10.*log10(sum(data.freq2value(100,2500).^2))+93.98)/0.1)*0.1 - 15 - impactInsulationClass);
if min(data.freqVector) <= 50 % if we have data
C_50 = round(round((10.*log10(sum(data.freq2value(50,2500).^2))+93.98)/0.1)*0.1 - 15 - impactInsulationClass);
else
C_50 = nan;
end
C = [C C_50];
else else
C = 0; C = 0;
end end
...@@ -100,7 +106,7 @@ if sArgs.createPlot ...@@ -100,7 +106,7 @@ if sArgs.createPlot
bar(gca,deficiencies.freqVector,deficiencies.freq,'hist'); bar(gca,deficiencies.freqVector,deficiencies.freq,'hist');
[maxDef,maxIdx] = max(deficiencies.freq); [maxDef,maxIdx] = max(deficiencies.freq);
if strcmpi(sArgs.type,'iso') if strcmpi(sArgs.type,'iso')
singleNumberString = [outputStr ' = ' num2str(impactInsulationClass) ' (' num2str(C) ') dB']; singleNumberString = [outputStr ' = ' num2str(impactInsulationClass) ' (' num2str(C(1)) ',' num2str(C(2)) ') dB'];
else else
singleNumberString = [outputStr ' = ' num2str(impactInsulationClass) 'dB']; singleNumberString = [outputStr ' = ' num2str(impactInsulationClass) 'dB'];
end end
......
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