// Leistungsspektrum2 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(); // Daten für die eigene Colormap cm2 = [0,0,0.5625; 0,0.3333,0.7083 0,0.6667,0.8542; 0,1,1; 0.0625,1,0.9375; 0.125,1,0.875; 0.1875,1,0.8125; 0.25,1,0.75; 0.3125,1,0.6875; 0.375,1,0.625; 0.4375,1,0.5625; 0.5,1,0.5; 0.5625,1,0.4375; 0.625,1,0.375; 0.6875,1,0.3125; 0.75,1,0.25; 0.8125,1,0.1875; 0.875,1,0.125; 0.9375,1,0.0625; 1,1,0; 1,0.9565,0; 1,0.913,0; 1,0.8696,0; 1,0.8261,0; 1,0.7826,0; 1,0.7391,0; 1,0.6957,0; 1,0.6522,0; 1,0.6087,0; 1,0.5652,0; 1,0.5217,0; 1,0.4783,0; 1,0.4348,0; 1,0.3913,0; 1,0.3478,0; 1,0.3043,0; 1,0.2609,0; 1,0.2174,0; 1,0.1739,0; 1,0.1304,0; 1,0.087,0; 1,0.0435,0; 1,0,0; 0.9762,0,0; 0.9524,0,0; 0.9286,0,0; 0.9048,0,0; 0.881,0,0; 0.8571,0,0; 0.8333,0,0; 0.8095,0,0; 0.7857,0,0; 0.7619,0,0; 0.7381,0,0; 0.7143,0,0; 0.6905,0,0; 0.6667,0,0; 0.6429,0,0; 0.619,0,0; 0.5952,0,0; 0.5714,0,0; 0.5476,0,0; 0.5238,0,0; 0.5,0,0]; // 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=1;//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************ //=============================================================== w=0; //Zähler für das Ergebnis Array // yW=window('hn',start2); // fenster fsetlegen // while start2