// Leistungsspektrum1 Mittelwert
mode(0);
// Hier werden WAV Dateien analysiert
// Die Abtastrate muss 48000 sein
// Es wird nur 1 Kanal analysiert
// Überdeckung der Analyse 87,5 //
//  
//
mclose('all');
clear;
xdel();
// WAV Datei auswählen   
[dateiname,pfad]=uigetfile('*.wav','D:\Mdaten');    //WAV Datei wählen
if isequal(dateiname,0);    //Ist eine Date gewählt?
    //
else
    //
    fe=pfad+'\'+dateiname;  //Vollständigen Pfad erstellen 
    info = wavread(fe,'info');  //WAV Info`s laden
    fd=fileinfo(fe);            //Datei Info`s laden
    disp('////////////////////////////////////////////////////////////')
    disp('WAV-Datei: '+fe); //Ausgabe vollständiger Pfad
    p='Kanäle: '+string(info(2));   //Ausgabe Anzahl der Kanäle
    disp(p);                        //Ausgabe Anzahl der Kanäle
    p='Samplerate: '+string(info(3));   //Ausgabe Samplerate
    disp(p);                            //Ausgabe Samplerate
    p='Dateilänge : '+string(fd(1));    //Ausgabe Dateilänge
    disp(p);                            //Ausgabe Dateilänge
    p='Laufzeit : '+string(info(8)/info(3));    //Laufzeit in Sek.
    disp(p);                                            //Laufzeit in Sek.
        disp('////////////////////////////////////////////////////////////')
    if info(3)==48000 then  //Stimmt die Samplerate der Audiodatei?
        samplmax=info(8);   //Dateilänge als Obergrenze erfassen
        /////////////////////////////////////////////
        start1=1;       //Erstes Byte des Rechenzyklus
        start2=48000;   //Letztes Byte des Rechenzyklus************************
        ueberdeckung=6000;//Überdeckungsabstand
        teiler=2400;//Amplitudengröße
        //===============================================================
        kanal=1;    //Nummer des Kanals der analysiert wird. 
        frequ1=1;   //Die FFT wird von dieser Frequenz an berechnet**********
        frequ2=500; //Die FFT wird bis zu dieser Frequenz berechnet************
        //===============================================================
        yF=wavread(fe,[start1 start2]); //Bytes für nächsten Rechnzyklus
        //
        // window funktion berechenen hanning
        yW=window('hn',start2);
        yFW=yW.*yF(kanal,:);
        //
        yFFTL=fft(yFW); //Die FFT wird berechnet
        // powerspektrum
        w=1;    //Zähler für das Ergebnis Array
        Py = (yFFTL).* conj(yFFTL);   //Powerspektrum der FFT berechnen
        Py=Py/teiler;
        //
        alles(w,frequ1:frequ2)=Py(frequ1:frequ2);  //Ergebnisse ins Array 
        start1=start1 + ueberdeckung;   //Der Startpunkt für den nächsten Zyklus
        start2=start2 + ueberdeckung;   //Der Endpunkt für den nächsten Zyklus
        // alle folgenden berechnen
        while start2<samplmax   //Prüfe ob Ende der Audiodaten erreicht
            //
            w=w+1;  //Zähler für das Ergebnis Array
            yF=wavread(fe,[start1 start2]); //Audiodaten werden gelesen
            //
            yFW=yW.*yF(kanal,:);
            //
            yFFTL=fft(yFW); //Die FFT wird berechnet
            start1=start1+ueberdeckung; //Der Startpunkt für den nächsten Zyklus
            start2=start2+ueberdeckung; //Der Endpunkt für den nächsten Zyklus
            // powerspektrum
            Py = (yFFTL).* conj(yFFTL);//Powerspektrum der FFT berechnen
            Py=Py/teiler;
            alles(w,frequ1:frequ2)=Py(frequ1:frequ2);  //Ergebnisse ins Array 
        end
    else   
    end
//
    alles2=mean(alles,'r'); //Mittelwert berechnen
    plot(alles2(1:100)); //abgebildete Frequenzen
////////////////////////////////////////////////////////
//Handle of type with properties:
//========================================
/////////////////////////////////////////////////////////
a=gca();
a.visible = "on";
a.axes_visible = ["on","on"];
a.axes_reverse = ["on","off"];
a.auto_ticks = ["off","off"];

a.grid = [1,1];
a.grid_thickness = [1,1];
a.grid_style = [1,1];
a.grid_position = "background";

a.title.font_size=2;
a.title.text="FFT - Mittelwert";
a.x_label.text="Frequenz";
a.x_label.visible="on";
a.y_label.font_size=2;
a.y_label.text="Amplitude";
a.y_label.visible="on";

a.auto_scale="off";

a.children.children(1).thickness=3;
a.children.children(1).foreground=2;
//a.x_ticks=tlist(['locations','labels'],[0,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,280,290,300],["0","10","20","30","40","50","60","70","80","90","100","110","120","130","140","150","160","170","180","190","200","210","220","230","240","250","260","270","280","290","300"]);
a.x_ticks=tlist(['locations','labels'],[0,5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100],["0","5","10","15","20","25","30","35","40","45","50","55","60","65","70","75","80","85","90","95","100"]);

a.y_ticks=tlist(['locations','labels'],[0,50,100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000],["0","50","100","150","200","250","300","350","400","450","500","550","600","650","700","750","800","850","900","950","1000"]);

a.sub_tics=[4,4];

end
mclose('all')
f=gcf();           
f.figure_position = [30,30];
f.figure_size = [1400,900];

////////////////////////////////////////////////////