Dokument [PDF, 9,1 MB] - FB 4 Allgemein - Fachhochschule ...

20. Anhang 211
%- Es liegt ein lineares APS vor, also muss quadratisch gemittelt
%- werden.
aps_lin_mittel_Pa_gesamt = sqrt(sum(aps_lin_Pa_gesamt.^2, 2)/n);
%---------------------------------------------------------------------
%- Bearbeitung des gesamten Frequenzbereiches
%- Die Ergebnisse zum gesamten Frequenzbereich sind durch
%- "_gesamt" gekennzeichnet
%---------------------------------------------------------------------
%- Berechnung des Gesamtpegels des gesamten APS-Bereiches
gp_aps_lin_mittel_Pa_gesamt = sqrt(1/epsilon * ↵
sum(aps_lin_mittel_Pa_gesamt.^2));
%- Schalldruckpegel berechnen, in dB für den gesamten Bereich
aps_lin_mittel_dB_gesamt = 20. * log10(aps_lin_mittel_Pa_gesamt ./ ↵
p_null);
gp_aps_lin_mittel_dB_gesamt = 20. * log10(gp_aps_lin_mittel_Pa_gesamt ./ ↵
p_null);
%----------------------------------------------------------------------
%- A-Bewertung
aps_lin_mittel_dB_A_gesamt = aps_lin_mittel_dB_gesamt + Lp_A_spl_gesamt;
%- Amplituden aus aps_lin_mittel_dB_A_gesamt für gp_ges in dB(A)
amp_A_Pa_gesamt = 10.^(aps_lin_mittel_dB_A_gesamt ./ 20) * p_null;
gp_lin_A_Pa_gesamt = sqrt(1/epsilon * sum(amp_A_Pa_gesamt.^2));
gp_aps_lin_mittel_dB_A_gesamt = 20. * log10(gp_lin_A_Pa_gesamt ./ p_null);
%----------------------------------------------------------------------
%----------------------------------------------------------------------
%- Zeitabhängiger GP
%- Abstand der abgetasteten Punkte
dt = 1/abtastrate;
%- Anzahl der delta_t im Zeitsignal
n_delta_t = (dauer-mod(dauer, delta_t))/delta_t;
disp(['Berechnung von ', num2str(n_delta_t), ' Gesamtpegeln']);
%- Anzahl der dt in delta_t
%- Über dt wird dann delta_t formuliert damit liegt man direkt im
%- Abtastraster
n_dt_delta_t = (delta_t-mod(delta_t, dt))/dt;
%- Formulierung von delta_t durch dt
delta_t_dt = n_dt_delta_t * dt;
for i = 1 : n_delta_t
i_s = (i-1) * n_dt_delta_t + 1;
i_e = i * n_dt_delta_t;
%- Zeitachse
gp_z(i,1) = i * delta_t_dt;
%- Gesamtpegel
gp_z(i,2) = 20.*log10(sqrt(1/n_dt_delta_t * ↵
sum(zeitdaten_Pa(i_s:i_e).^2))/p_null);
end
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz