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

20. Anhang 209
display('Fehler');
hold
end %- switch
%---------------------------------------------------------------------
switch fenster_typ
case 'hanning'
fenster_typ
window = hanning(blockgroesse);
fenster = 2;
epsilon = 1.5;
case 'flattop'
fenster_typ
windos = flattop(blockgroesse);
fenster = 4.688600988;
epsilon = 3.83;
end %- if-Anweisung
%- Aus den Zeitdaten werden nun soviele Werte übernommen, wie ganze Blöcke
%- hineinpassen. Es passen n_block mal ganze Blöcke in die Zeitdaten-
%- Matrix.
laenge_zeitdaten_Pa = length(zeitdaten_Pa);
dauer = laenge_zeitdaten_Pa / abtastrate; % [s]
dauer_block = blockgroesse / abtastrate; % [s]
delta_f = 1 / dauer_block; % [Hz]
mittelungsanzahl = n_block;
%- Frequenzspanne ermitteln
f_min = 0;
%- Theoretische Maximalfrequenz
f_max_th = abtastrate / divisor;
%- Maximalfrequenz, die mit delta_f erfasst wird
f_max = f_max_th - mod(f_max_th, delta_f);
%- delta_f passt n-mal in f_max, es gibt aber n+1 Stützstellen
n_f = f_max / delta_f + 1;
%- Matrix mit den abgedeckten Frequenzen
frequenzen = linspace(f_min, f_max, n_f);
%- Spline für die A-Bewertung erzeugen
%- Es wird aufgrund der dbA-Stützstellen und des Frequenzbereiches
%- "frequenzen" eine Ausgleichskurve erzeugt und für jeden Frequenz-
%- wert ein dBA-Korrekturwert abgelegt
Lp_A_spl_gesamt = (spline(f_mn, dba_stuetz, frequenzen))';
freq_string_2 = ['Frequenzspanne [', num2str(f_min), '; ',num2str(f_max), ↵
'] Hz'];
%---------------------------------------------------------------------
%---------------------------------------------------------------------
%- Berechnung des linearen APS über die Funktion
%- Im Nachfolgenden wird zwischen dem Frequenzbereich (_bereich) und
%- dem gesamten Frequenzbereich (_gesamt) unterschieden. Der vorgege-
%- bene Frequenzbereich (_bereich) wird durch den angegebenen Terz-
Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz