%   .txt-Datei in Matlab einlesen ( Matrixform )

rohdaten1=dlmread('01_white_20kmh_1\23.txt', '\t', 23, 0);

%   Aufteilen in Zeit- u. Wertevektor

time1=rohdaten1(:,1);
soundpressurelevel1=(rohdaten1(:,2))';

%   Plotten der Rohdaten

a = figure ('Name','Rohdaten Fahrt 1 im Vergleich','NumberTitle','off');
    subplot ( 2,1,1 ) ;
        plot ( time1 , soundpressurelevel1 , 'r' );
        xlabel ( 'Zeit [s]' ) ; ylabel ( 'Schalldruck [Pa]' ) ; title( 'Rohsignal' ) ;
        grid on ; hold on ;
        
    subplot(212);
        plot ( time1 , soundpressurelevel1 , 'r' ) ;
        grid on ; hold on ;
        
%   FFT für Fahrt 1 mit Plot

t11=time1;
x11=soundpressurelevel1';

T11=diff(t11(1:2));
fs11=1/T11;
N11=length(x11);

fn=0.5*fs11;
df=fs11/N11;

f11 = 0:df:fn;

H11=fft(x11);
H11_pos0 = abs(H11(1:(N11/2)+1));

%H11_pos01=[H11_pos0(1)/N11 H11_pos0(2:(N11/2))/(N11/2) H11_pos0(end)/N11];
H11_pos01= H11_pos0./(N11/2);
H11_pos01(1)= H11_pos01(1)/(N11);
H11_pos01(end)= H11_pos01(end)/(N11);

b=figure('Name','FFT der 1. Fahrt','NumberTitle','off');
    subplot(3,1,1);
        semilogx(f11,H11_pos01,'r');                  
        xlabel('Frequenz[Hz]'); title('Rohdaten Spektrum');
        grid on; hold on;
        
    subplot(3,1,3);
        semilogx(f11,H11_pos01,'r');
        grid on; hold on;
        

%   A-Bewertung im Frequenzbereich

Ra=( (12200^2.*f11.^4)/((f11.^2+20.6^2) .* (f11.^2+12200^2) .* (sqrt(f11.^2+107.7^2)) .* (sqrt(f11.^2+737.9^2))) );

H11A=H11_pos01 .* Ra';

figure(b);
    subplot(3,1,2);
        semilogx(f11,H11A,'b');
        xlabel('Frequenz [Hz]'); title('Rohdaten Spektrum A-Bewertet');
        grid on; hold on;
        
    subplot(3,1,3);
        semilogx(f11,H11A,'b');
        xlabel('Frequenz [Hz]'); title('Vergleich der Spektren');
        grid on, hold on,
        legend('Spektrum','Spektrum A-Bewertet',0);