PHYSICS EXPERIMENTS BASED ON DATA ACQUISITION THE LONGITUDINAL ACOUSTIC DOPPLER EFFECT

Similar documents
PHYSICS EXPERIMENTS BASED ON DATA ACQUISITION THE LONGITUDINAL ACOUSTIC DOPPLER EFFECT

Doppler effect (Item No.: P )

STUDIES ON WAVES AND OSCILLATIONS WITH DATA ACQUISITION SYSTEMS *

AVAL AUDIO-VISUAL ACTIVE LOCATOR. Faculty Sponsor: Professor Kathleen E. Wage Kelly Byrnes Rony Alaghbar Jacob Cohen

Inquiring activities on the acoustic phenomena at the classroom using sound card in personal computer

CONSIDERATIONS ON USING THE VIRTUAL INSTRUMENTS FOR THE ACQUISITION AND ANALYSIS OF EXPERIMENTAL DATA FROM DYNAMIC SYSTEMS

Lab VIEW Programming for Vibration Analysis

LOW LEVEL CURRENTS AND VOLTAGES MEASUREMENTS

Enhancing Physics Teaching with Technology.

ni.com Sensor Measurement Fundamentals Series

HEAD. Advanced Filters Module (Code 5019) Overview. Features. Module with various filter tools for sound design

Microphone Array Measurements for High-speed Train

Acoustic Doppler Effect

Chapter 16 Sound. Copyright 2009 Pearson Education, Inc.

AN AUTOMATED ALGORITHM FOR SIMULTANEOUSLY DETERMINING ULTRASONIC VELOCITY AND ATTENUATION

INTERACTIVE TECHNOLOGY LABORATORIES

Individually configurable system. Microphone Arrays.

3D Intermodulation Distortion Measurement AN 8

Development of a sonic boom measurement system at JAXA

THE PATTERNS OF THE SOUND INTENSITY DISTRIBUTION OF MIDRANGE LOUDSPEAKER

INTERACTIVE SKETCHING OF THE URBAN-ARCHITECTURAL SPATIAL DRAFT Peter Kardoš Slovak University of Technology in Bratislava

On the accuracy reciprocal and direct vibro-acoustic transfer-function measurements on vehicles for lower and medium frequencies

Investigation of Noise Spectrum Characteristics for an Evaluation of Railway Noise Barriers

Applying Virtual Oscilloscope to Signal Measurements in Scintillation Detectors

Non-Data Aided Doppler Shift Estimation for Underwater Acoustic Communication

Instrumentation (ch. 4 in Lecture notes)

Waves-Wave Behaviors

SigCal32 User s Guide Version 3.0

Directivity Controllable Parametric Loudspeaker using Array Control System with High Speed 1-bit Signal Processing

Experiments with wave, using low-cost amplitude modulated ultrasonic techniques

Beam Dynamics + Laser Micro Vibrometry 1

Speech Intelligibility Enhancement using Microphone Array via Intra-Vehicular Beamforming

Supplementary User Manual for BSWA Impedance Tube Measurement Systems

Biomechanical Instrumentation Considerations in Data Acquisition ÉCOLE DES SCIENCES DE L ACTIVITÉ PHYSIQUE SCHOOL OF HUMAN KINETICS

Acoustic Projector Using Directivity Controllable Parametric Loudspeaker Array

Summary. Page SEG SEG Denver 2014 Annual Meeting

NON CONTACT VIBRATION MEASUREMENTS ON PARABOLIC SURFACE ANTENNA. Dorin Simoiu 1, Liviu Bereteu 1

Lab 12 Laboratory 12 Data Acquisition Required Special Equipment: 12.1 Objectives 12.2 Introduction 12.3 A/D basics

Virtual Instrument for FPGA based Spectrum Analyzer

Investigating Electromagnetic and Acoustic Properties of Loudspeakers Using Phase Sensitive Equipment

Measurement & Control of energy systems. Teppo Myllys National Instruments

Clarification of the Effect of High-Speed Train Induced Vibrations on a Railway Steel Box Girder Bridge Using Laser Doppler Vibrometer

Multi-spectral acoustical imaging

Resonance Mode Acoustic Displacement Transducer

Considerations for Analog Input and Output

SENSING OF METAL-TRANSFER MODE FOR PROCESS CONTROL OF GMAW

Speech Intelligibility Enhancement using Microphone Array via Intra-Vehicular Beamforming

Influence of Vibration Amplitude Oscillations on the Conical Sieve Suspended Dorel STOICA, Gheorghe VOICU, Carmen RUSĂNESCU

Laboratory 1: Motion in One Dimension

Transfer Function (TRF)

LAB Week 7: Data Acquisition

A Road Traffic Noise Evaluation System Considering A Stereoscopic Sound Field UsingVirtual Reality Technology

WELCOME TO PHYC 493L Contemporary Physics Lab

Experiments with wave, using low-cost amplitude modulated ultrasonic techniques

CONTRIBUTION REGARDING NOISE MEASUREMENT ACOUSTIC PROCEDURES ON BOARD

Residual Phase Noise Measurement Extracts DUT Noise from External Noise Sources By David Brandon and John Cavey

Panca Mudji Rahardjo, ST.MT. Electrical Engineering - UB

EE 300W 001 Lab 2: Optical Theremin. Cole Fenton Matthew Toporcer Michael Wilson

Enhanced Resonant Inspection Using Component Weight Compensation. Richard W. Bono and Gail R. Stultz The Modal Shop, Inc. Cincinnati, OH 45241

Lab M6: The Doppler Effect

MEMS. Platform. Solutions for Microsystems. Characterization

[4] (b) Fig. 6.1 shows a loudspeaker fixed near the end of a tube of length 0.6 m. tube m 0.4 m 0.6 m. Fig. 6.

Acoustic Signature of an Unmanned Air Vehicle - Exploitation for Aircraft Localisation and Parameter Estimation

Waves-Wave Behaviors

Measurement System for Acoustic Absorption Using the Cepstrum Technique. Abstract. 1. Introduction

GENERATION OF SIGNALS USING LABVIEW FOR MAGNETIC COILS WITH POWER AMPLIFIERS

f n = n f 1 n = 0, 1, 2.., (1)

Measuring Voltage and Time Quantities of a Signal Through a Virtual Oscilloscope

Acoustic Yagi Uda Antenna Using Resonance Tubes

PerSec. Pervasive Computing and Security Lab. Enabling Transportation Safety Services Using Mobile Devices

A mechanical wave is a disturbance which propagates through a medium with little or no net displacement of the particles of the medium.

Audio Engineering Society. Convention Paper. Presented at the 115th Convention 2003 October New York, New York

Intelligent traffic control system using PLC

Computer based experiments for off-campus teaching and learning of AC electricity

The Development of Laser Ultrasonic Visualization Equipment and its Application in Nondestructive Inspection

Measuring the Speed of Sound in Air Using a Smartphone and a Cardboard Tube

Transient high sensitive soot measurement. AVL Micro Soot Sensor. Manfred Linke

INTERFEROMETRIC VIBRATION DISPLACEMENT MEASUREMENT

EENG473 Mobile Communications Module 3 : Week # (12) Mobile Radio Propagation: Small-Scale Path Loss

Control and Signal Processing in a Structural Laboratory

Fumiaki UEHAN, Dr.. Eng. Senior Researcher, Structural Mechanics Laboratory, Railway Dynamics Div.

EE 233 Circuit Theory Lab 3: First-Order Filters

NEW YORK STATE TEACHER CERTIFICATION EXAMINATIONS

BSWA Impedance Tube Solutions

Development of 4/16-Channel Data Acquisition System Using Lab VIEW

Using Signal Express to Automate Analog Electronics Experiments

Proceedings of Meetings on Acoustics

Sound Waves Speed Intensity Loudness Frequency Pitch Resonance Sound Waves

Localization of underwater moving sound source based on time delay estimation using hydrophone array

JOHANN CATTY CETIM, 52 Avenue Félix Louat, Senlis Cedex, France. What is the effect of operating conditions on the result of the testing?

Harmonic Distortions Analyzer for Power Rectifiers

Ascent Ground and Satellite Demonstration

Measurement of Amplitude Modulation AN 6

CHAPTER 7 HARDWARE IMPLEMENTATION

SUMMARY. ) f s Shock wave Sonic boom UNIT. Waves transmit energy. Sound is a longitudinal mechanical wave. KEY CONCEPTS CHAPTER SUMMARY

SW Series Impedance Tube Solutions

Rocking Drones with Intentional Sound Noise on Gyroscopic Sensors

CHAPTER 12 SOUND ass/sound/soundtoc. html. Characteristics of Sound

HDTV Mobile Reception in Automobiles

CORRECTION NOTICE SOUTH AFRICAN CIVIL AVIATION AUTHORITY CIVIL AVIATION ACT, 2009 (ACT NO. 13 OF 2009)

Transcription:

Romanian Reports in Physics, Vol. 68, No. 3, P. 1312 1325, 2016 PHYSICS EXPERIMENTS BASED ON DATA ACQUISITION THE LONGITUDINAL ACOUSTIC DOPPLER EFFECT M. OPREA 1,2, CRISTINA MIRON 1* 1 University of Bucharest, Faculty of Physics, Bucharest, Romania 2 Mihai Viteazul School, Călăraşi, Romania * Corresponding author: cmiron_2001@yahoo.com Received August 18, 2014 Abstract. The presence of a data acquisition device in a modern physics laboratory sets the ground for an advanced scientific investigation performed by the students, due to the fact that the studied physical phenomena can be qualitatively and quantitatively explored in an exhaustive manner. A physics experiment includes the acquisition, processing, storing and interpreting of the experimental data taken with the help of sensors associated with the studied physical quantities (temperature, pressure, force etc.). The analog signals received by the sensors connected to a data acquisition device are turned into digital signals. These are later on processed with the help of a software application in order to monitor the variation in the physical quantities correlated with the studied phenomena. In this paper we focus on the study on the longitudinal acoustic Doppler Effect by performing the spectral analysis of the sound picked up by a microphone, using the NIDAQ acquisition device and the graphical programming environment LabView. Key words: Doppler effect, Doppler shift, virtual experiment, data acquisition, signal frequency. 1. INTRODUCTION Recent studies in the didactics of Physics highlight the students attitude towards the scientific content of this discipline and the computer-based strategies that a teacher may employ in order to capture the attention and interest of his students [1 3]. On quite a significant number of occasions they show a certain reluctance in approaching scientifical matters and a considerable lack of interest. In such circumstances, the teacher must raise the qualitative standards of the experiment [4]. Regarded as a bidirectional communication interface teacherstudent, the experiment can bring substantial qualitative influences upon the actors involved in the didactic process. A successful physics experiment must necessarily gain the attention of the students: it has to be dynamic, spectacular and efficient. Only this way the students involvement, enthusiasm and creativity will become visibly enhanced.

2 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1313 With minimal hardware (data acquisition device, computer) and software res (LabView) used in an intelligent and didactically relevant way this goal can be achieved, as some authors have pointed out [5 9]. Based on these considerations, we set out to design an experiment for the illustration of the acoustic Doppler effect, which we extensively elaborate on in this paper. 2. THEORETICAL BACKGROUND To begin with, we are going to highlight some theoretical aspects concerning the acoustic Doppler effect [10]. In a given physical environment (atmospheric air), we have an acoustic S and a receiver R which travel one towards another at the speed of vs and v R. The S generates a sound with a frequency f 0 which travels with the speed c (Fig. 1). Fig. 1 Source and receiver moving towards one another. The colored versions can be accessed at http://www.infim.ro/rrp/. The apparent frequency perceived by the receiver R is c v f f 0 c v f0 is the frequency of the signal emitted by the acoustic. We will analyse two particular cases of the longitudinal Doppler effect: R S, where a) Stationary receiver, moving ( v R 0 ). In this c case f f 0, when the moves towards the receiver, and c v S c f f 0, when the moves away from the receiver c ( vs vs v S ). If we know the values of f and f 0, we can determine the

1314 M. Oprea, Cristina Miron 3 f f0 speed of the based on the relation: v S c, where + and f designate the cases where the approaches the receiver and moves away from it. b) Moving receiver, stationary ( v 0 ). The observer perceives a c vr frequency f f 0, when the approaches him, and c c vr f f 0, when the moves away from him ( vr vr ). If we c know the values of f and f 0, we can determine the speed of the receiver f f0 based on the relation: v R c, where + and designate the cases f where the receiver approaches the and moves away from it. The present paper illustrates the experimental aspects of the two situations a) and b) described above for the determination of the Doppler shift in case of approach and recession, as well as of the speeds v S and vr in these cases. As we can see in Fig. 2 and Fig. 3, the acoustic signal generated by an audio will be captured by a microphone with a parabolic reflector which is connected to the sound card of a laptop. From here, the signal is amplified and applied on an analog input of a data acquisition device NIDAQ6008, connected to the laptop through a USB port. The data flow is processed with a software application made in LabView, which offers information on the spectral content of the signal detected by the microphone. This way one can observe the Doppler shift of an analysed spectral component. S Fig. 2 Experimental setup. The colored versions can be accessed

4 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1315 Fig. 3 Experimental schema. Before discussing the experimental results, we will briefly illustrate the application we designed for the analysis of the experimental data. 3. LABVIEW APPLICATION The Front Panel of the application contains two diagrams in which we can observe the amplitudes of the signals acquired by the NIDAQ6008 and the acoustic frequency spectrum recorded by it (Fig. 4). Fig. 4 Front Panel of LabView application. The colored versions can be accessed Two groups of indicators display the maximal and minimal values of the signal s amplitude and of the frequencies associated with these values. A numerical indicator lists the values of the signal s amplitude recorded in the experiment. Studying the diagram of the application, one can observe that the dynamic data flow is converted into a numeric data flow and graphically displayed on the

1316 M. Oprea, Cristina Miron 5 Front Panel after being processed by the Power Spectrum element. The array of numeric values can be saved to a text-based measurement file (.lvm), from where it can be exported for analysis in MS Excel (Fig. 5). Fig. 5 Block diagram of LabView application. The colored versions can be accessed at http://www.infim.ro/rrp/. 4. EXPERIMENTAL RESULTS As we previously stated, the experimental part of this paper includes two situations: 1. Moving, stationary receiver; 2. Stationary, moving receiver. As an acoustic we used a speaker connected to a minilaptop on which we started a software application for generating audio signals (Fig. 6). We then proceeded to select the frequency of 1 khz. Based on the LabView application we designed for the signal acquisition, we verified the amplitude and frequency of the audio signal generated by the speaker (Fig. 7). Fig. 6 Acoustic. The colored versions can be accessed at http://www.infim.ro/rrp/. Fig. 7 signal f = 1kHz. The colored versions can be accessed at http://www.infim.ro/rrp/.

6 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1317 A. Moving, stationary receiver. The audio was placed in a car which drove on the highway in a linear motion over a distance of 100m, at the end of which we installed a receiver composed of a parabolic microphone connected to the sound card of a laptop. From the output of the sound card, the signal amplified through dedicated software (Audacity) was applied on an analog input of the data acquisition device NIDAQ6008 (Fig. 8). Fig. 8 Experimental setup (speaker, parabolic mirophone, DAQ device). The colored versions can be accessed The acoustic signal was started while the car was moving towards the receiver with a constant speed v S. For a better resolution of the measurements, the recordings of the signal detected by the microphone were made when the car was in its vicinity, after previously choosing two reference markers simetrically placed relative to the receiver (10 m) (Fig. 9). Fig. 9 Source moving towards receiver with v S between the reference markers R1 and R2 on a distance d R1R2 = 10 m. The colored versions can be accessed The experimental results indicate a sonic frequency of 1050 Hz when the vehicle approaches the receiver and a frequency of 950 Hz when the vehicle moves

1318 M. Oprea, Cristina Miron 7 away from the receiver (Fig. 10, Fig. 11). Given the average air temperature of 29 C, the measured Doppler shift was 50 Hz when the approached the receiver and when it moved away from the receiver (Table 1). Fig. 10 Case A. Moving, stationary receiver; Doppler shift measurement approach (Labview left, Ms Excel right). The colored versions can be accessed Fig. 11 Case A. Moving, stationary receiver; Doppler shift measurement recession (Labview left, Ms Excel right). The colored versions can be accessed type Table 1 measurements Moving stationary receiver, f = 1 khz, approaching v S, t = 29 C receding Experimental values 1050 950 Doppler shift 50 50 The recording of the air s temperature in order to establish a precise value of the sound speed was performed with a thermistor NTC, which performed the role of a sensor and was connected to the data acquisition device (Fig. 12). During the time interval 15:00 15:20 of the experiment there weren t any notable changes of atmospheric pressure and temperature and the thermic variation interval was situated between the values of 28.57 C and 28.91 C, which made us use the rounded value t = 29 C.

8 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1319 Fig. 12 Experimental setup (left) and LabView temperature measurement (right). The colored versions can be accessed Five consecutive experimental measurements were performed during the 20 minutes dedicated to the data acquisition process (Table 2). Table 2 measurements Moving stationary receiver, f = 1 khz, v S, t = 29 C Time of the recording Freq. appr. Freq. reced. 15:00 15:05 15:10 15:15 15:20 Average values Approx. values 1050 1050 1050 1049 1049 1049.6 1050 949 949 950 950 950 949.6 950 As we stated previously, knowing the Doppler shift enables us to calculate the speed of the, which we proceeded on doing with the help of a LabView application that we designed (Fig. 13). By averaging the values obtained when the approached and moved away from the receiver, we obtained a speed of v S = 59.72 km/h. Fig. 13 Case A. Moving, stationary receiver; Speed determination based on Doppler shift ( approach left, recession right). The colored versions can be accessed

1320 M. Oprea, Cristina Miron 9 B. Stationary, moving receiver. In this case, the parabolic microphone was placed in the car and the was put sideways from the road. The acoustic signal was started as the car moved towards the at a constant speed v. Just like in case A (moving stationary S receiver), the recordings of the signal detected by the microphone were performed when the car drove near the two reference markers symmetrically placed relative to the. During the experiment, which was performed in the time interval 15:30 15:50, the temperature was situated around the average value of 29 C. The experimental values which were obtained are present in Table 3. These results are identical to the ones reported in the previously studied case A. Table 3 measurements Stationary (1000 Hz), moving receiver ( v R ), t = 29 C type approaching receding receiver receiver Experimental values 1050 950 Doppler Shift 50 50 We calculated the value of the speed v R = 59.58 km/h (Fig. 14). Fig. 14 Case B. Stationary, moving receiver; Speed determination based on Doppler shift (receiver approach left, receiver recession right). The colored versions can be accessed Additional determinations. Changing the frequency of the acoustic to 2000 Hz and 3000 Hz, we performed additional experiments in order to determine the Doppler shift and the speed of the and the receiver (Table 4 and Table 5).

10 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1321 Table 4 Determination of frequency and speed Moving, stationary receiver Initial experimental conditions: v S t = 29 C Source frequency= 2000 Hz Source frequency= 3000 Hz Experimental values approaching receding approaching Doppler Shift receding Calculated v S (km/h) 2100 1909 100 91 56.77 3151 2864 151 136 56.84 Table 5 Determination of frequency and speed Stationary, moving receiver Initial experimental conditions: v R t = 29 C Source frequency = 2000 Hz Source frequency = 3000 Hz Experimental values approaching receiver receding receiver approaching receiver Doppler Shift receding receiver Calculat v ed R (km/h) 2096 1905 96 95 56.89 3143 2857 143 143 56.79 Maintaining the frequency of the constant, we determined the Doppler shift for different speeds. This time, the air temperature was situated around the value of 26 C (Table 6 and Table 7). Calculated speed type Table 6 measurements Moving (1000 Hz), stationary receiver, t = 26 C Freq. app. 57.63 km/h 66.53 km/h 75.95 km/h Freq. reced. Freq. app. Freq. reced. Source Freq. app. Freq. reced. Experimental values 1051 954 1059 947 1068 940 Doppler shift 51 46 59 53 68 60

1322 M. Oprea, Cristina Miron 11 Calculated speed type Experiment al values Doppler shift Table 7 measurements Stationary (1000 Hz), moving receiver, t = 26 C Freq. app. receiver 57.19 km/h 66.72 km/h 76.26 km/h Freq. reced. receiver Freq. app. receiver Freq. reced. receiver Freq. app. receiver Freq. reced. receiver 1051 954 1059 947 1068 940 51 46 59 53 68 60 At the end, the dependency characteristics were drawn: detected frequency frequency (at constant speed) (Table 8, Figs. 15, 16) and detected frequency speed (at constant frequency) (Table 9, Figs. 17, 18). Table 8 measurements (t = 29 C) Moving v S = 59.72 km/h; Stationary receiver Stationary ; Moving receiver v R = 59.58km/h Source frequency Approach frequency Receding frequency 1000 1050 950 100 2000 2100 1909 191 3000 3151 2864 287 Approach freq. receding freq. 1000 1050 950 100 2000 2096 1905 191 3000 3143 2857 286 Fig. 15 Moving, stationary receiver. The colored versions can be accessed

12 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1323 Fig. 16 Stationary, moving receiver. The colored versions can be accessed Table 9 measurements (f = 1000 Hz, t = 26 C) Moving ; Stationary receiver Stationary ; Moving receiver Speed (km/h) Approach frequency Receding frequency Approach freq. receding freq. 57.63 1051 954 97 66.53 1059 947 112 75.95 1068 940 128 57.19 1048 952 96 66.72 1056 944 112 76.26 1064 936 128 Fig. 17 Moving stationary receiver. The colored versions can be accessed

1324 M. Oprea, Cristina Miron 13 Fig. 18 Stationary moving receiver. The colored versions can be accessed Analysing Figs. 15 and 16 one can notice an increase in the offset f as the frequency of the sound increases, while its speed approach f recession remains constant. From Figs. 17 and 18 one can observe an increase in the offset f as the speed of the increases, while its frequency remains approach f recession constant. Furthermore, we designed additional experiments for low speeds in order to determine the Doppler shift and speed, when: a) a person is running and b) a person is going on a bicycle. We inserted below an extract from the experimental measurements for the standard frequency of 2000 Hz. Table 10 measurements, moving stationary receiver f = 2000Hz, t = 30 C Person running Person on bicycle type Approching frequency Receding frequency Approching frequency Experimental values 2020 1980 2025 1975 Doppler Shift 20 20 25 25 Calculated speed (km/h) 11.91 14.91 Receding frequency The LabView application we designed also enables us to go beyond the usual experimental situations which can be studied in class. Thus, we can simulate what happens when the speed of the is close to the speed of sound (sonic boom) and when the frequency of the acoustic reaches very high values (ultrasound spectrum).

14 Experiments based on data acquisition the longitudinal acoustic Doppler effect 1325 5. CONCLUSIONS The didactic efficiency of an experimental approach such as the one illustrated in this paper is guaranteed. The students have shown an authentic desire for research and investigation and have managed to develop a qualitative understanding of the phenomenon. The experiment enhaced their ability to communicate scientifically relevant content, their teamwork spirit, their practical and digital abilities. These observations confirm other researchers findings [3, 4] regarding the importance of having a combined didactic approach: complementing the real experiment with the virtual one. REFERENCES 1. S. Moraru, I. Stoica, F.F. Popescu, Rom. Rep. Phys. 63, 2, 577 586 (2011). 2. C. Kuncser, A. Kuncser, G. Maftei, S. Antohe, Rom. Rep. Phys. 64, 4, 1119 1130 (2012). 3. L. Dinescu, M. Dinica, C. Miron, E.S. Barna, Rom. Rep. Phys. 65, 2, 578 590 (2013). 4. M. Dems, Experimental Methods in Science, Course support, 2 nd semester, Technical University of Łódz Science and Technology. 5. A.A.J. Glean, J.A. Judge, J.F. Vignola, P.F. O Malley and T.J. Woods, Journal of the Acoustical Society of America 129, 4, 2648 (2011). 6. I.K. Lau, Design of Measurement Techniques of Acoustic Glass Shattering System, Final Year Project, UTAR, 2013. 7. P. Dhara, A. Roy, P. Maity, P. Singhai, P.S. Roy, Design of the Data Acquisition System for The Nuclear Physics Experiments at Vecc, 9 th International Workshop on Personal Computers and Particle Accelerator Controls, held at VECC, 2012. 8. D. Jiang, J. Xiao, H. Li, Q. Dai, Eur. J. Phys. 28, 977 982 (2007). 9. D. Amrani, P. Paradis, Lat. Am. J. Phys. Educ. 4, 3, 511 514, (2010). 10. A. Hristev, Mechanics and Acoustics (in Romanian), Edit. Didactică şi Pedagogică, Bucharest, 1982.