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

20. Anhang 200
B = B_th(1:n_f,:);
F = F_th(1:n_f);
[m, n] = size(B);
%---------------------------------------------------------------------
%- Lineares APS
%---------------------------------------------------------------------
%- Das APS kann aus der Multiplikation der FFT mit ihrer konjugiert
%- komplexen berechnet werden. Daraus geht dann das quadratische APS hervor
%- [V^2]. Das lineare APS erhält man, indem man aus dem quadratischen APS
%- die Wurzel zieht, oder das APS über die Betragsbestimmung der komplexen
%- Elemente in der FFT-Matrix B ermittelt.
%- Der zweite Weg wird hier angewendet.
%---------------------------------------------------------------------
aps_lin = (2.0 * fenster_flattop * abs(B) / blockgroesse) / sqrt(2.0);
%- Der Faktor 2.0 ersetzt die theoretisch durchzuführende Addition
%- der Amplituden des negativen und positiven Frequenzbereiches,
%- was bei reellen Signalen zulässig ist.
%- Der Fensterfaktor (fenster_flattop) macht die Normierung des
%- angewendeten Wichtungsfensters rückgänig.
%- Mit der Division durch sqrt(2.) werden RMS-Amplituden berechnet.
%- Gesonderte Beachtung des Gleichanteils
%- Der Gleichanteil kommt nur einmal vor, darf also nicht mit
%- dem Faktor 2. vervielfacht werden.
aps_lin(1,1:n) = aps_lin(1,1:n) ./ 2.0;
%- Mitteln des APS (quadratisch)
aps_lin_mittel = sqrt(sum(aps_lin.^2, 2)/n);
%- Berechnung des Gesamtschalldruckpegels
GP_aps_lin_mittel = sqrt(1/epsilon_flattop * sum(aps_lin_mittel.^2));
%---------------------------------------------------------------------
%- Berechnung des Kalibrierfaktors
%---------------------------------------------------------------------
%- Ermittlung der Kalibrierfrequenz und des Kalibrierpegels
%---------------------------------------------------------------------
%- Vorgabe, bei welchem Index der Kalibrierpegel zu finden ist
index_1000 = round(1000 / delta_f + 1); %- ungefährer Index des Pegels
index_250 = round(250 / delta_f + 1); %- ungefährer Index des Pegels
delta_index = round(20 / delta_f); %- Bereich um den Pegel
%- Maximum im vorgegebenen Frequenzbereich berechnen
max_1000 = max(aps_lin_mittel((index_1000 - delta_index):(index_1000 + ↵
delta_index)));
max_250 = max(aps_lin_mittel((index_250 - delta_index):(index_250 + ↵
delta_index)));
%- Das größere Maximum kennzeichnet den Kalibrierpegel
if max_1000 > max_250
frequenz_kali = 1000
pegel_kali = 94
elseif max_1000 < max_250
frequenz_kali = 250
pegel_kali = 124
else
Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz