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

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

mv.fh.duesseldorf.de
von mv.fh.duesseldorf.de Mehr von diesem Publisher
22.07.2013 Aufrufe

20. Anhang 210 %- bereich bestimmt. %--------------------------------------------------------------------- [B_th, F_th, T] = specgram(zeitdaten_Pa, blockgroesse, abtastrate,window, ↵ ueberlappung); %--------------------------------------------------------------------- %- Indizes für die Bereichsbeschreibung index_f_start = find(F_th >= f_start); index_f_ende = find(F_th

20. Anhang 211 %- Es liegt ein lineares APS vor, also muss quadratisch gemittelt %- werden. aps_lin_mittel_Pa_gesamt = sqrt(sum(aps_lin_Pa_gesamt.^2, 2)/n); %--------------------------------------------------------------------- %- Bearbeitung des gesamten Frequenzbereiches %- Die Ergebnisse zum gesamten Frequenzbereich sind durch %- "_gesamt" gekennzeichnet %--------------------------------------------------------------------- %- Berechnung des Gesamtpegels des gesamten APS-Bereiches gp_aps_lin_mittel_Pa_gesamt = sqrt(1/epsilon * ↵ sum(aps_lin_mittel_Pa_gesamt.^2)); %- Schalldruckpegel berechnen, in dB für den gesamten Bereich aps_lin_mittel_dB_gesamt = 20. * log10(aps_lin_mittel_Pa_gesamt ./ ↵ p_null); gp_aps_lin_mittel_dB_gesamt = 20. * log10(gp_aps_lin_mittel_Pa_gesamt ./ ↵ p_null); %---------------------------------------------------------------------- %- A-Bewertung aps_lin_mittel_dB_A_gesamt = aps_lin_mittel_dB_gesamt + Lp_A_spl_gesamt; %- Amplituden aus aps_lin_mittel_dB_A_gesamt für gp_ges in dB(A) amp_A_Pa_gesamt = 10.^(aps_lin_mittel_dB_A_gesamt ./ 20) * p_null; gp_lin_A_Pa_gesamt = sqrt(1/epsilon * sum(amp_A_Pa_gesamt.^2)); gp_aps_lin_mittel_dB_A_gesamt = 20. * log10(gp_lin_A_Pa_gesamt ./ p_null); %---------------------------------------------------------------------- %---------------------------------------------------------------------- %- Zeitabhängiger GP %- Abstand der abgetasteten Punkte dt = 1/abtastrate; %- Anzahl der delta_t im Zeitsignal n_delta_t = (dauer-mod(dauer, delta_t))/delta_t; disp(['Berechnung von ', num2str(n_delta_t), ' Gesamtpegeln']); %- Anzahl der dt in delta_t %- Über dt wird dann delta_t formuliert damit liegt man direkt im %- Abtastraster n_dt_delta_t = (delta_t-mod(delta_t, dt))/dt; %- Formulierung von delta_t durch dt delta_t_dt = n_dt_delta_t * dt; for i = 1 : n_delta_t i_s = (i-1) * n_dt_delta_t + 1; i_e = i * n_dt_delta_t; %- Zeitachse gp_z(i,1) = i * delta_t_dt; %- Gesamtpegel gp_z(i,2) = 20.*log10(sqrt(1/n_dt_delta_t * ↵ sum(zeitdaten_Pa(i_s:i_e).^2))/p_null); end %----------------------------------------------------------------------- %----------------------------------------------------------------------- Fachhochschule Düsseldorf Diplomarbeit 2002/03, Terence Klitz

20. Anhang 210<br />

%- bereich bestimmt.<br />

%---------------------------------------------------------------------<br />

[B_th, F_th, T] = specgram(zeitdaten_Pa, blockgroesse, abtastrate,window, ↵<br />

ueberlappung);<br />

%---------------------------------------------------------------------<br />

%- Indizes für die Bereichsbeschreibung<br />

index_f_start = find(F_th >= f_start);<br />

index_f_ende = find(F_th

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!