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

20. Anhang 199
textdatei = strcat(pfad, dateiname(1:(laenge_dateiname- ↵
12)),'_Kalibrierfaktor.txt');
fid = fopen(textdatei,'w');
fprintf(textdatei,' \n');
end
%---------------------------------------------------------------------
laenge_zeitdaten_roh = length(zeitdaten_roh);
%- Kanalwahl, rechts oder links
switch lower(kanal)
case {'rechts'}
zeitdaten_kanal = zeitdaten_roh(1:laenge_zeitdaten_roh,2);
case {'links'}
zeitdaten_kanal = zeitdaten_roh(1:laenge_zeitdaten_roh,1);
otherwise
display('Fehler');
hold
end
window = flattop(blockgroesse);
fenster_flattop = 4.688600988;
epsilon_flattop = 3.83;
ueberlappung = 0;
p_null = 2e-5; %- [Pa]
%- 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_roh = length(zeitdaten_roh);
n_block = (laenge_zeitdaten_roh - mod(laenge_zeitdaten_roh, ↵
blockgroesse)) / blockgroesse;
wertezahl = n_block * blockgroesse;
zeitdaten = zeitdaten_kanal(1:wertezahl,1);
laenge_zeitdaten = length(zeitdaten);
dauer = laenge_zeitdaten / abtastrate; %- [s]
dauer_block = blockgroesse / abtastrate; %- [s]
delta_f = 1 / dauer_block; %- [Hz]
mittelungsanzahl = n_block;
f_min = 0;
% theoretische Maximalfrequenz
f_max_th = abtastrate / divisor;
% Maximalfrequenz, die mit delta_f abgedeckt werden kann.
f_max = f_max_th - mod(f_max_th, delta_f);
% delta_f passt n-mal in f_max rein, es gibt aber n+1 Stützpunkte
n_f = f_max / delta_f + 1;
frequenzen = (linspace(f_min, f_max, n_f))';
%---------------------------------------------------------------------
%- Berechnung des linearen APS über die Funktion specgram
%---------------------------------------------------------------------
[B_th, F_th, T] = specgram(zeitdaten, blockgroesse, abtastrate, window,
ueberlappung);
%- Übernahme des berücksichtigten Frequenzbereiches.
Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz