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

20. Anhang 204
tic
%---------------------------------------------------------------------
%- Wahl des Kanals der WAV-Datei, der analysiert werden soll;
%- diese Wahl ist notwendig, weil das Programm nur einen Kanal
%- analysiert.
%--------------------------------------------------------------------kanal
= input('Bitte den zu analysierenden Kanal angeben ("rechts" oder ↵
"links"):', 's');
%---------------------------------------------------------------------
%---------------------------------------------------------------------
%- Programmeinstellungen
%- Standardeinstellung für ein Signal, das mit 44100 Hz abgetastet wurde:
%- Blockgröße von 32768 Werten
%- Divisor von 2.56
%---------------------------------------------------------------------
format long g
divisor = 2.56; %- Es wird eine kleine Frequenzspanne als
%- 44100/2 Linien betrachtet.
fenster_typ = 'hanning'; %- oder 'flattop'; für das Flattop-Fenster
%- muss die Flattop-Funktion [Haukap]
%- existieren (vgl.
%- Kalibrierfaktor_Kanal_bereich.m.
ueberlappung = 0;
p_null = 2e-5; % [Pa]
blockgroesse = 32768;
delta_t = 0.3; %- Zeitspanne für die GP(t)-Berechnug
f_mn_a = 16; %- Erste beachtete Terz (Nenn-Frequenz),
f_mn_e = 12500; %- letzte beachtete Terz (Nenn-Frequenz);
%- Dadurch das die Frequenzspanne durch
%- den Divisor dividiert wird, ist am Ende der f-Spanne darauf zu achten,
%- welche Terz noch gefüllt werden kann. Bei einem Divisor von 2,56 beträgt
%- die letzte Frequenz 17226,5625 Hz, das reicht nicht für das Terzband mit
%- f_m = 16000 Hz, wo f_2 = 17783 Hz beträgt. Die hier eingestellte
%- Frequenzspanne ist auch bei der Schmalbandanalyse zu beachten.
q = 10^(1/10); %- Stufensprung für die Normreihe R10
%- Oktavmittenfrequenzen für die Achsbeschriftung
oktav_m = {'16'; ''; ''; '31,5'; ''; ''; '63'; ''; ''; '125'; ''; ''; ↵
'250'; ''; ''; '500'; ''; ''; '1000'; ''; ''; '2000'; ''; ''; '4000'; ''; ↵
''; '8000'; ''; ''; '16000'};
terz_string = ['Terzbereich [', num2str(f_mn_a), '; ', num2str(f_mn_e), '] ↵
Hz'];
freq_string_1 = ['Frequenzspanne [', num2str(f_mn_a), '; ', ↵
num2str(f_mn_e), '] Hz'];
%---------------------------------------------------------------------
%- Vorbereiten der Terzbänder in dem angegebenen Frequenzbereich.
%-
%- Es wird gemäß DIN EN 61260 mit den exakten Terz-Bandmittenfre-
%- quenzen gerechnet. Aus diesen werden dann die oberen und unteren
%- Bandeckfrequenzen errechnet, nach denen die Frequenzlinien den
%- einzelnen Terzbändern zugeordnet werden.
%- Die Terzmittenfrequenzen sind nach der Normreihe R10 formuliert,
Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz