dBEuler - Scilab

Oﬃcial publisher & professional services provider of Scilab, Open Source So+ware for Numerical Computa5on & Simula5on

dBEuler project Statistical pass-by rolling noise measurements in Scilab ECOSURF programme IFSTTAR Guillaume DUTILLEUX

ERA32 Acoustique PCI Acoustique et Vibrations LRPC Strasbourg France Centre d'Études techniques de l'Équipement de l'Est

On January 1st 2014, the 8 CETE facilities, Certu, Cetmef and Setra will merge to form a new structure called Cerema. JTAV AUTUN - 10/10/2012 - Loïc TOUSSAINT Ministère de l'Écologie, du Développement Durable et de la Mer Bruit de roulement au passage : Evolution des outils de mesurage et de capitalisation

1

Introduction ●

●

Burden of disease due to environmental noise –

> 1 M DALYs / year

–

Sleep disturbance

–

Cognitive impairment

–

Cardio-vascular disease

Major role of road traffic noise –

Passengers

–

Goods

●

Mostly rolling noise above 50 km/h

●

Pavements matter !

●

●

Standardized procedures for pavement assessment and classification Software required, why not Scilab ? ScilabTEC - 26/06/2013 - Palaiseau, France - Guillaume Dutilleux - MEDDE/CETE Est/PCI Acoustique et Vibrations dBEuler project - Statistical pass-by rolling noise measurements in Scilab

2

Outline ●

Pass-by measurement procedure

●

(Short) history of SPB systems

●

dBEuler

●

–

General aspects

–

Audio acquisition

–

Post-processing

–

Report generation

–

Customization

–

Help and maintenance

Conclusions / Perspectives

ScilabTEC - 26/06/2013 - Palaiseau, France - Guillaume Dutilleux - MEDDE/CETE Est/PCI Acoustique et Vibrations dBEuler project - Statistical pass-by rolling noise measurements in Scilab

3

Pass-by measurements - 1/2 ●

S(tatistical)PB or C(ontrolled)PB

●

Combined acquistion (speed, acoustic pressure)

●

Statistical analysis


4

Pass-by measurements - 1/2 ●

S(tatistical)PB or C(ontroled)PB

●

Combined acquistion (speed, acoustic pressure)

●



5

Pass-by measurements – 2/2

For each vehicle category ScilabTEC - 26/06/2013 - Palaiseau, France - Guillaume Dutilleux - MEDDE/CETE Est/PCI Acoustique et Vibrations dBEuler project - Statistical pass-by rolling noise measurements in Scilab

6

Systems for SPB in France ●

●

●

1990 dBLRS - Analog recording –

Measurement microphones

–

NAGRA IV

–

Dedicated hardware

–

Proprietary software

2006 dBEuler 1.0 – Digital recording –

Generic hardware

–

Scilab 3.1.1

–

Third party software

2013 dBEuler 2.0 –

Measurement hardware

–

Scilab 5.4.1


7

dBEuler 2.0 features ●

●

Supported standards –

SPB - S 31119:1993

–

SPB - NF EN ISO 11819-1:1997

–

CPB - NF S 31-119-2:2000

Management of a complete SPB/CPB campaign –

Acquisition on max. 4 channels

–

Pass-by post-processing

–


–

Report generation

–

Possible non standard campaigns


8

Design principles ●

Free software

●

Existing software

●

Cross-platform technologies

●

« Low-tech » languages

●

Hardware –

Generic audio devices

–

Acquisition device

●

Campaign = directory tree

●

Working copy and lost+found

●

Standard and human-readable formats


9

Architecture dBEuler (Scilab, tcl, XSLT) Scilab 5.4.1 Shell SoX

Image Magick

Tcl/Tk

Audio dll + exe (C) xsltproc

LaTeX

Portaudio

SinusMess driver

Microsoft Windows XP | Windows 7 32 bits


10

Starting dBEuler


11

Starting dBEuler


12

Recording calibrated audio ●

Built-in audio : client / server connection

Scilab

C dll Client

Record / Stop Audio samples

C Server

●

●

Portaudio

standard audio device

SINUS/Harmonie device

Standalone recorders


13

Acoustic post-processing 1/4

xgetmouse() ScilabTEC - 26/06/2013 - Palaiseau, France - Guillaume Dutilleux - MEDDE/CETE Est/PCI Acoustique et Vibrations dBEuler project - Statistical pass-by rolling noise measurements in Scilab

14



15



16



17

Acoustic post-processing 2/4 ●

Time-weighted sound pressure level

●

Low-pass filtering NF EN 61672-1:2003 function [L,Lmax]=compute_Ltau(sig,fech,tau) //sig : array (n_sig,p_samples), input signal, Pa //fech : scalar, sampling frequency, Hz //tau : scalar, time constant, s. //L : array (n_sig, p_samples), Ltau, dB re 2e-5 Pa //Lmax : array (n_sig,1), max level, dB re 2e-5 Pa. p0=2e-5; fc=(1/tau)/(fech/2);

lptau=iir(1,'lp','butt',[fc 0],[0,0]); lptau=tf2ss(lptau); L=2*dB(1/p0*sqrt(flts(sig.^2,lptau))); Lmax=max(L,'c'); endfunction


18

Acoustic post-processing 3/4 ●

A-weighting NF EN 61672-1:2003

●

Digital filtering –

Bilinear transform

–

Frequency domain


19

Acoustic post-processing 4/4 Class I 1/3 d'octave filter bank NF EN 61260:1995 Template conformance

●

●

Filter type ● Butterworth ● 3rd order Scilab functions used ● iir() ● flts()


20

Statistical analysis ●

Standards: linear least squares → reglin()

●

Standards in practice: « expert » least squares (GUI + reglin())


21

Statistical analysis ●

Standards: linear least squares → reglin()

●

Standards in practice: « expert » least squares (GUI + reglin())

●

dBEuler 2.0: robust regression by m-estimator (α=-0.5)


22

Customization: preferences //warning when insufficient dynamic range checkDynamicRange=checkDynamicRange6dB; //=checkDynamicRange10dB; //=checkDynamicRangeNone; : silent operation //computing spectrum MPspectrumType='L3Fmax'; //=L3Fmax : (complies with NS S 31-119-2) //=L3FmaxISO : (complies with NF EN ISO 11819-1) compute_regression=reglin; //=reglin; : standard Scilab function //=M_estim_SEF_alpha_m0_5 : robust M-estimator

PreferenceCampaign.sce ... //Audio hardware //MPRecorder="Harmonie" //MPRecorder="Server" MPRecorder="None" MPmetrics='ltau'; //MPmetrics='leq' //Extracting pass-bys extractPassBy=extractPassByUser; //=extractPassByUser; //=extraitPassBy4s; //=extraitPassBy10dB; ...

rmOutlier=%T; //%F : no outlier suppression //type of campaign MPtype='SPB'; //Vehicle categories MPcatVehicle=['LV','PL','UT','HGV']; MPcatVref=[90,90,90,80]; //recording mode postProcess='mono'; //n-channel mode //='stereo'; : conventional mode //ID of recording channels : MPchannel=['7m50-1m20',1,'7m50-3m',2,'Comments',3]; ...


23

Report

unix_g() ScilabTEC - 26/06/2013 - Palaiseau, France - Guillaume Dutilleux - MEDDE/CETE Est/PCI Acoustique et Vibrations dBEuler project - Statistical pass-by rolling noise measurements in Scilab

24

Install, Use and Maintain ●

●

Distribution –

Scilab-controled setup

–

Infrastructure

–

dBEuler specific code

–

Sample campaigns

Documentation –

Admin, User, Dev

–

Docbook 5 + DBLaTeX


25

Conclusions ●

●

dBEuler means –

Calibrated audio hardware control

–

Interactive post-processing

–

Outlier-aware statistical analysis

–

Pre-press quality report generation

–

… mostly in Scilab

An experiment in intermediate technology –

From proprietary to free software

–

More responsibility but more flexibility

–

… Scilab offers a suitable feature level / learning curve


26

Perspectives ●

●

Infrastructure –

Scilab 5.x → 6.x ?

–

Fewer 3rd party components

Potential new features –

Interface with radar speedometer

–

Acoustic speedometer

–

Extended filter bank

–

Fully automated post-processing ?


27

Acknowledgements

dBEuler 1.x ScilabTEC - 26/06/2013 - Palaiseau, France - Guillaume Dutilleux - MEDDE/CETE Est/PCI Acoustique et Vibrations dBEuler project - Statistical pass-by rolling noise measurements in Scilab

dBEuler 2.0 28

Thank you for your attention Guillaume DUTILLEUX ERA32 Acoustique PCI Acoustique et Vibrations CETE Est / LRPC de Strasbourg France [email protected] On January 1st 2014, the 8 CETE facilities, Certu, Cetmef and Setra will merge to form a new structure for expertise and engineering in relation to risks, urban and country planning, the environment, and transport which is to be known as Cerema. - CETE: Governmental engineering survey department - Certu: Centre for studies on urban planning, transport and public facilities - Cetmef: Institute for maritim and inland waterways - Sétra : Technical department for transport, roads and bridges

Centre d'Études Techniques de l'Équipement de l'Est JTAV AUTUN - 10/10/2012 - Loïc TOUSSAINT Ministère de l'Écologie, du Développement Durable et de la Mer Bruit de roulement au passage : Evolution des outils de mesurage et de capitalisation

29

dBEuler - Scilab

dBEuler - Scilab

Suggest Documents

SciLab

SCILAB

Scilab

Name - Scilab

Name - Scilab

Scilab manual

Scilab / Xcos

Name - Scilab

Manual Scilab

Description - Scilab

Exercices de seconde - Scilab

Scilab for very beginners

Introduction to Scilab

Scilab for Real Dummies

Programming in Scilab

INTRODUCTION TO SCILAB

LENGUAJE DE PROGRAMACIÓN SCILAB

(GUI)IN SCILAB.

What is Scilab?

Plotting in Scilab

Scilab for very beginners

Scilab for Matlab Info

Scilab Bag Of Tricks

Scilab and modeFRONTIER integration