%
% Beispiel zur Nutzung von DSPntelOktave
%
%
%   Dipl.-Ing. Armin Rohnen LbA
%   Labor für Maschinendynamik
%   Hochschule für angewandte Wissenschaften München
%   University of applied sciences Munich
%
%   31. Juli 2016
%
%
%%
% Parameter die als Übergabeparameter erforderlich sind

clear;
tic


%filename = '15ccw_01.wav';
%[data, fs] = audioread(filename);
load('spieluhr.mat');
ntel = 12;                      % ntelOktav
filterord = 8;                  % Steilheit der ntelOktavfilter 
tdata = OData(:,1);             % Zeitsignal für die ntelOktav-Analyse
zeitk = 125;                    % Zeitkonstante in ms für die RMS-Pegelberechnung
overlap = 50;                   % Overlap für die RMS-Pegelberechnung 0 = kein Overlap

[dataFilt, dataOut, timeOut, frequenzen] = DSPntelOktave(ntel, filterord, fs, tdata, zeitk, overlap);


toc

%%
%
%

mesh(frequenzen, timeOut, dataOut);
set(gca,'XScale','log');
axis([100 10000 0 25 0.1 0.6])
colorbar;
view(60,30);