%   .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 ( 4,1,1 ) ;
        plot ( time1 , soundpressurelevel1 , 'r' );
        xlabel ( 'Zeit [s]' ) ; ylabel ( 'Schalldruck [Pa]' ) ; title( 'Rohsignal' ) ;
        grid on ; hold on ;
        
    subplot(414);
        plot ( time1 , soundpressurelevel1 , 'r' ) ;
        grid on ; hold on ;
        
%   FFT für Fahrt 1 mit Plot

t11=time1;
x11=soundpressurelevel1';

T11=diff(t11(2:1));
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((N11/2)+1)/N11];

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

%   A-Bewertung im Frequenzbereich

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

H11A=H11_pos01.*Ra;

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