% mv240_splMeter.mlx
% Signalanalyse von Luftschalldaten mit  splMeter
%
% Dipl.-Ing. Armin Rohnen LbA
% Hochschule München - Munich University of Applied Sciences
% Faktultät 03 - Labor für Schwingungstechnik und Maschinendanmik
% Dachauer Straße 98b, 80335 München
% eMail: rohnen@hm.edu oder armin@rohnen.net
%
% München, 26. Februar 2022
%
% Lizenzbedingungen:
% Dieses Programm ist freie Software. Sie können es unter den Bedingungen der GNU General Public License, wie von der Free Software Foundation veröffentlicht, weitergeben und/oder modifizieren, entweder gemäß Version 3 der Lizenz oder (nach Ihrer Option) jeder späteren Version.
% Die Veröffentlichung dieses Programms erfolgt in der Hoffnung, daß es Ihnen von Nutzen sein wird, aber OHNE IRGENDEINE GARANTIE, sogar ohne die implizite Garantie der MARKTREIFE oder der VERWENDBARKEIT FÜR EINEN BESTIMMTEN ZWECK. Details finden Sie in der GNU General Public License.
% GNU General Public License http://www.gnu.org/licenses/    http://www.gnu.de/documents/gpl.de.html
%

clear;
load('20220225_Casablanca.mat');

plot(time, data, 'LineWidth', 2)
set(gca, 'FontSize', 16)
title('Zeitsignal der Luftschallmessung')
ylabel('Schalldruck in Pa')
xlabel('Zeit in s')
SPL = splMeter('FrequencyRange', [22 22050], ...
    'FrequencyWeighting', 'A-weighting', ...
    'TimeWeighting', 'Fast', ...
    'PressureReference', 20e-6, ...
    'CalibrationFactor', 1, ...
    'TimeInterval', 1, ...
    'SampleRate', fs);

[Lt,Leq] = SPL(data);

release(SPL);

plot(time, Lt, 'LineWidth', 2);
hold on
plot(time,Leq, 'LineWidth', 2);
set(gca, 'FontSize', 16);
title('Pegelverlauf');
ylabel('Pegel in dB(A)');
xlabel('Zeit in s');
legend('L_t : Momentanpegel', 'L_{eq} : zeitlicher Mittelwert')
axis([0 20 20 100])