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

20. Anhang 201
display('Fehler - Kalibrierfrequenz');
hold;
end %- if-Anweisung
%- Frequenzbereich um die theoretische Kalibrierfrequenz festlegen
von_f_k = frequenz_kali - f_bereich; %- [Hz]
bis_f_k = frequenz_kali + f_bereich; %- [Hz]
%- Frequenzbereich in einen Stützstellenbereich umrechnen
von_index_k = (von_f_k - mod(von_f_k, delta_f)) / delta_f + 1;
bis_index_k = (bis_f_k - mod(bis_f_k, delta_f)) / delta_f + 1;
%- Berechnung des Kalibrierfaktors über den Gesamtpegel
GP_kali = sqrt(1/epsilon_flattop * ↵
sum(aps_lin_mittel(von_index_k:bis_index_k).^2));
kk_bereich = 10^(pegel_kali / 20) * p_null ./ GP_kali;
m_k_bereich (i,1) = kk_bereich;
k_bereich = num2str(kk_bereich,8);
aps(:,1) = F;
aps(:,i+1) = aps_lin_mittel; %- aps(:,2) rechter Kanal
%- aps(:,3) linker Kanal
%---------------------------------------------------------------------
%- Plotten des linearen APS
%--------------------------------------------------------------------figure(1);
set(gcf,'name', datei_string);
subplot(2,2,i);
plot(aps(:,1), aps(:,i+1));
title(['Lineares APS von Kanal: ', kanal, '; Kalibrierfrequenz = ', ↵
num2str(frequenz_kali) ' Hz; k = ', num2str(m_k_bereich(i,1)), '
[Pa/EU]']);
xlabel('f [Hz]');
ylabel(['Amplitude [EU]']);
grid;
%---------------------------------------------------------------------
%- Ablegen des Kalibrierfaktors in einer Datei
%---------------------------------------------------------------------
fid = fopen(textdatei,'a');
fprintf(fid,'%11.8f\t%s \n', kk_bereich, kanal);
fclose(fid);
display(['Datei ', kanal, ' gespeichert']);
kanal = 'links';
maximum (i,1) = max(aps_lin_mittel);
%---------------------------------------------------------------------
%- Anwendung der Kalibrierfaktoren auf das Kalibriersignal
%- Mit dieser Kontrolle kann geprüft werden, ob der Kalibrierfaktor
%- richtig berechnet wurde.
%--------------------------------------------------------------------if
i == 1
amp_rechts = zeitdaten_roh(1:laenge_zeitdaten_roh,2);
kali_fak_re = m_k_bereich(1,1)
amp_rechts_kali = amp_rechts(1:wertezahl) * kali_fak_re;
[B_re_th, F_re_th, T_re] = specgram(amp_rechts_kali, blockgroesse, ↵
abtastrate, window, ueberlappung);
Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz