UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA

Size: px
Start display at page:

Download "UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA"

Transcription

1 UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA Naval Uiversity of Gdyia Gdyia, Śmidowicza 69, Polad Gdańsk Uiversity of Techology 8-95 Gdańsk, Narutowicza 11/1, Polad The paper deal with the basic problem coected with the free vibratio of fiite cylidrical steel shell immersed i water. At the first there is described the theoretical model of this problem by meas of the Helmholtz equatio with the vo Neuma boudary coditio o the radiated surface. Next the results of a umerical ivestigatio are preseted. The fial results are the acoustic characteristics obtaied i the sea. These oes were compared to the umerical results.. INTRODUCTION Uderwater object detectio, idetificatio ad classificatio immersed or buried by meas of a acoustic method has bee of great iterest for a few decades by acousticias maries ad archaeologists. There are ofte used the objects e.g. mies, amphoras i a fiite cylidrical shape or similar. The method of active locatios (echolocatios) is very popular i detectio of the immersed or buried objects. I the shape of reflected soudig pulse oe ca observe the rigded part of soudig pulse ad so-called tail. I this tail we ca fid a substatial feature of the object, that scattered the soud pulse. This paper described the basic method that allows to determie the shape, volume of the cylidrical steel shell immersed or buried i the sea. 1. PROBLEM FORMULATION For describig the pressure distributio i a far field we used the Helmholtz equatio i the followig form: r r r Φ R) + k Φ( R) = f ( ) (1) (

2 1 1 = ( q ) + q q q q q = z + y, R = ϕ q + x + x, The solutio of this equatio is as follow: r Φ( r ) Φ( R) = r G( R, r N S G( R, r is a radiated surface ad ) is a Gree fuctio. S ) ds () O the source surface for r = a should be satisfied the boudary coditio so-called vo Neuma coditio i this form: Φ( r ) r r = w& ( ) (3) N r = a w& ( r ) is ormal velocity of the shell vibratio. The pressure distributio is described by the followig formulae: acρ iπ / ( + 1) i( ωt kr) p = W& m me e [ J ( kasi Θ) + ie ( kasi Θ)] ψ (4) 4R si( kl / cos Θ) ( kl cos Θ) ψ = [1 + ] cos Θ (mπ ) ( kl cos Θ) c - speed of soud; ρ - medium desity; J ( ), E ( ) - Bessel ad Weber-Lommel fuctios order respectively. Geometric cofiguratio of radiatig cylider is preseted o figure 1.

3 Fig. 1. Geometric cofiguratio of radiatig cylider We assume that cylidrical surface is covered by cotiuous distributio of sigle sources give by: r Wm mπ w& m ( ) = & [1 cos (x + L)] cos ϕ (5) L m, are modal ad odal umber respectively, is amplitude of ormal velocity. W& m Aother relatio is obtaied whe we take ito accout the pressure wave radiatio by a part of a ifiite cylider. It is give by these formulae: ixξ = iρωm iϕ A ( ξ, ϕ ) H ( τr) e p e π d τh τa ξ m (6) ( ) ξ is a wave umber i x directio. After simplificatio we obtai: cρε ( ) iπ / i( ω t kr) p = W& ψ m m me e 4πkR (si Θ) H ( ka si Θ) (7) ε ( ) = 1 for = ε ( ) = for 1 Distributio of soud pressure ca be fid by meas of the relatio (4) or (7) for the cases whe we kow ω ad coected with them W &. O the figure the chart of p m ad p m as fuctio of kasi Θ was preseted. m m 1

4 Fig.. The chart of p ad m 1 p as fuctio of kasi Θ m To ormalized these fuctio we divided it by factor: acρ W& m ϕ. (8) 4R. EXPERIMANTAL INVESTIGATION As a first step we are carried out examiatios of the ormal velocity distributio o the cylidrical steel shell. The cylidrical shell was immersed i water i a aechoic basi. O its surface was mouted velocity trasducer as is show i figure 3. Fig. 3. Positio of mouted velocity trasducer at cylidrical shell These sigals after a pulse excitatio were recorded ad aalysed. O the base of this examiatio we ca estimate fudametal frequecies coected with free oscillatios frequecies. Figure below show a example of velocity sigal observed at poit VIII.

5 Fig. 4. Sigals observed at poit VIII O the figure 5 are show the form of fuctio for differet free vibratios frequeces. The velocity distributio alog the x - axis ad ϕ - agle. Fig. 5. The forms of fuctio for differet free vibratio frequeces

6 3. MEASUREMENTS CARRIED OUT IN THE SEA The ext step of these ivestigatios were measuremets of scattered soud pressure by a fiite cylidrical steel shell excitig by soud pulse geerated by the impulsive source. This kid of source allows to excite low free oscillatios of the cylidrical shell for e.g. less tha f = Hz. I figure 6 are show two curves as a results of compariso of two differet relatios betwee the values obtaied by meas of theoretical relatio ad measuremet oe. The better approximatio gives the relatio (7). Fig. 6. Results of compariso of two differet relatios 4. CONCLUSIONS Fidig a free oscillatio of a immersed fiite cylidrical steel shell we ca estimate its volume ad shape ad some times also iside cotet. These ivestigatios are a first step to fid a method for iterpretatios some relatio betwee geometrical relatio of the target ad scattered soud field. The ext step of the ivestigatio will be devoted to the proper sigal processig allowig to extract the iterestig iformatio from scattered sigal. REFERENCES [1] Proceedigs of The applicatio of recet advaces i uderwater.detectio ad survey techiques to uderwater archeology, Bodrum, May 4. [] A. Kołodziejski, E. Kozaczka, Broadbad source of uderwater disturbaces, Archives of Acoustics 1,,, 197-1, 1986

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,

More information

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Damped Vibration of a Non-prismatic Beam with a Rotational Spring Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a No-prismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS. ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics

More information

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + 62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of

More information

Modeling the Radiation Field of Normal-Beam Ultrasonic Transducers by Multi-Gaussian Beam Technique

Modeling the Radiation Field of Normal-Beam Ultrasonic Transducers by Multi-Gaussian Beam Technique Modelig the Radiatio Field of Normal-Beam Ultrasoic Trasducers by Multi-Gaussia Beam Techique Ami Yaghootia, Farhag Hoarvar, ad Mehdi Zeighami NDE Lab., Faculty of Mechaical Egieerig, K. N. Toosi Uiversity

More information

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray

More information

Types of Waves Transverse Shear. Waves. The Wave Equation

Types of Waves Transverse Shear. Waves. The Wave Equation Waves Waves trasfer eergy from oe poit to aother. For mechaical waves the disturbace propagates without ay of the particles of the medium beig displaced permaetly. There is o associated mass trasport.

More information

577. Estimation of surface roughness using high frequency vibrations

577. Estimation of surface roughness using high frequency vibrations 577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia

More information

Lecture 3. Digital Signal Processing. Chapter 3. z-transforms. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. 2016

Lecture 3. Digital Signal Processing. Chapter 3. z-transforms. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. 2016 Lecture 3 Digital Sigal Processig Chapter 3 z-trasforms Mikael Swartlig Nedelko Grbic Begt Madersso rev. 06 Departmet of Electrical ad Iformatio Techology Lud Uiversity z-trasforms We defie the z-trasform

More information

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity

More information

FINALTERM EXAMINATION Fall 9 Calculus & Aalytical Geometry-I Questio No: ( Mars: ) - Please choose oe Let f ( x) is a fuctio such that as x approaches a real umber a, either from left or right-had-side,

More information

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare Spectral Aalysis This week i lab Your ext experimet Homework is to prepare Next classes: 3/26 ad 3/28 Aero Testig, Fracture Toughess Testig Read the Experimets 5 ad 7 sectios of the course maual Spectral

More information

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series. .3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(

More information

Motor Acoustics Overview

Motor Acoustics Overview Motor Acoustics Overview Ceter for Acoustics & Vibratio May 5, 9 Dr. Marti L. Pollack APS Backgroud Research, Developmet, & Egieerig Small Busiess Approximately 57 Egieers & Scietists ~ ½ hold PhDs i Egieerig,

More information

MAT 271 Project: Partial Fractions for certain rational functions

MAT 271 Project: Partial Fractions for certain rational functions MAT 7 Project: Partial Fractios for certai ratioal fuctios Prerequisite kowledge: partial fractios from MAT 7, a very good commad of factorig ad complex umbers from Precalculus. To complete this project,

More information

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Molecular ad Quatum Acoustics vol. 7, (6) 79 ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Jerzy FILIPIAK 1, Lech SOLARZ, Korad ZUBKO 1 Istitute of Electroic ad Cotrol Systems, Techical Uiversity of Czestochowa,

More information

Acoustic Field inside a Rigid Cylinder with a Point Source

Acoustic Field inside a Rigid Cylinder with a Point Source Acoustic Field iside a Rigid Cylider with a Poit Source 1 Itroductio The ai objectives of this Deo Model are to Deostrate the ability of Coustyx to odel a rigid cylider with a poit source usig Coustyx

More information

Analysis Methods for Slab Waveguides

Analysis Methods for Slab Waveguides Aalsis Methods for Slab Waveguides Maxwell s Equatios ad Wave Equatios Aaltical Methods for Waveguide Aalsis: Marcatilis Method Simple Effective Idex Method Numerical Methods for Waveguide Aalsis: Fiite-Elemet

More information

Chapter 6 Sampling Distributions

Chapter 6 Sampling Distributions Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to

More information

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,

More information

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the

More information

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series

More information

Math 113 Exam 3 Practice

Math 113 Exam 3 Practice Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you

More information

Question 1: The magnetic case

Question 1: The magnetic case September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to

More information

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions .03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio

More information

Chapter 1. Complex Numbers. Dr. Pulak Sahoo

Chapter 1. Complex Numbers. Dr. Pulak Sahoo Chapter 1 Complex Numbers BY Dr. Pulak Sahoo Assistat Professor Departmet of Mathematics Uiversity Of Kalyai West Begal, Idia E-mail : sahoopulak1@gmail.com 1 Module-2: Stereographic Projectio 1 Euler

More information

Solution of EECS 315 Final Examination F09

Solution of EECS 315 Final Examination F09 Solutio of EECS 315 Fial Examiatio F9 1. Fid the umerical value of δ ( t + 4ramp( tdt. δ ( t + 4ramp( tdt. Fid the umerical sigal eergy of x E x = x[ ] = δ 3 = 11 = ( = ramp( ( 4 = ramp( 8 = 8 [ ] = (

More information

Probability, Expectation Value and Uncertainty

Probability, Expectation Value and Uncertainty Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such

More information

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c 4th Iteratioal Coferece o Advaced Materials ad Iformatio Techology Processig (AMITP 06) Directio of Arrival Estimatio Method i Uderdetermied Coditio Zhag Youzhi a, Li eibo b, ag Hali c Naval Aeroautical

More information

Signal Processing in Mechatronics

Signal Processing in Mechatronics Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal

More information

Vector Quantization: a Limiting Case of EM

Vector Quantization: a Limiting Case of EM . Itroductio & defiitios Assume that you are give a data set X = { x j }, j { 2,,, }, of d -dimesioal vectors. The vector quatizatio (VQ) problem requires that we fid a set of prototype vectors Z = { z

More information

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

Signals & Systems Chapter3

Signals & Systems Chapter3 Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of

More information

n 3 ln n n ln n is convergent by p-series for p = 2 > 1. n2 Therefore we can apply Limit Comparison Test to determine lutely convergent.

n 3 ln n n ln n is convergent by p-series for p = 2 > 1. n2 Therefore we can apply Limit Comparison Test to determine lutely convergent. 06 微甲 0-04 06-0 班期中考解答和評分標準. ( poits) Determie whether the series is absolutely coverget, coditioally coverget, or diverget. Please state the tests which you use. (a) ( poits) (b) ( poits) (c) ( poits)

More information

Name: Math 10550, Final Exam: December 15, 2007

Name: Math 10550, Final Exam: December 15, 2007 Math 55, Fial Exam: December 5, 7 Name: Be sure that you have all pages of the test. No calculators are to be used. The exam lasts for two hours. Whe told to begi, remove this aswer sheet ad keep it uder

More information

Chapter 2 Feedback Control Theory Continued

Chapter 2 Feedback Control Theory Continued Chapter Feedback Cotrol Theor Cotiued. Itroductio I the previous chapter, the respose characteristic of simple first ad secod order trasfer fuctios were studied. It was show that first order trasfer fuctio,

More information

11 Correlation and Regression

11 Correlation and Regression 11 Correlatio Regressio 11.1 Multivariate Data Ofte we look at data where several variables are recorded for the same idividuals or samplig uits. For example, at a coastal weather statio, we might record

More information

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION Hisao Izuchi, Pricipal Egieerig Cosultat, Egieerig Solutio Uit, ChAS Project Operatios Masato Nishiguchi, Egieerig Solutio Uit, ChAS

More information

Appendix: The Laplace Transform

Appendix: The Laplace Transform Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio

More information

[ ] sin ( ) ( ) = 2 2 ( ) ( ) ( ) ˆ Mechanical Spectroscopy II

[ ] sin ( ) ( ) = 2 2 ( ) ( ) ( ) ˆ Mechanical Spectroscopy II Solid State Pheomea Vol. 89 (003) pp 343-348 (003) Tras Tech Publicatios, Switzerlad doi:0.408/www.scietific.et/ssp.89.343 A New Impulse Mechaical Spectrometer to Study the Dyamic Mechaical Properties

More information

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu

More information

Phys 6303 Final Exam Solutions December 19, 2012

Phys 6303 Final Exam Solutions December 19, 2012 Phys 633 Fial Exam s December 19, 212 You may NOT use ay book or otes other tha supplied with this test. You will have 3 hours to fiish. DO YOUR OWN WORK. Express your aswers clearly ad cocisely so that

More information

Chapter 10: Power Series

Chapter 10: Power Series Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

CHAPTER 2. Mean This is the usual arithmetic mean or average and is equal to the sum of the measurements divided by number of measurements.

CHAPTER 2. Mean This is the usual arithmetic mean or average and is equal to the sum of the measurements divided by number of measurements. CHAPTER 2 umerical Measures Graphical method may ot always be sufficiet for describig data. You ca use the data to calculate a set of umbers that will covey a good metal picture of the frequecy distributio.

More information

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Discrete Orthogonal Moment Features Using Chebyshev Polynomials Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical

More information

Chapter 7 z-transform

Chapter 7 z-transform Chapter 7 -Trasform Itroductio Trasform Uilateral Trasform Properties Uilateral Trasform Iversio of Uilateral Trasform Determiig the Frequecy Respose from Poles ad Zeros Itroductio Role i Discrete-Time

More information

INTRODUCTION J. Acoust. Soc. Am. 101 (6), June /97/101(6)/3313/9/$ Acoustical Society of America 3313

INTRODUCTION J. Acoust. Soc. Am. 101 (6), June /97/101(6)/3313/9/$ Acoustical Society of America 3313 A combied itegro-modal approach for predictig acoustic properties of irregular-shaped cavities J. Missaoui ad L. Cheg Departmet of Mechaical Egieerig, Laval Uiversity, Québec, Québec G1K 7P4, Caada Received

More information

c. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area?

c. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area? 1. Research Methodology a. What is meat by the supply chai (SC) coordiatio problem ad does it apply to all types of SC s? Does the Bullwhip effect relate to all types of SC s? Also does it relate to SC

More information

(# x) 2 n. (" x) 2 = 30 2 = 900. = sum. " x 2 = =174. " x. Chapter 12. Quick math overview. #(x " x ) 2 = # x 2 "

(# x) 2 n. ( x) 2 = 30 2 = 900. = sum.  x 2 = =174.  x. Chapter 12. Quick math overview. #(x  x ) 2 = # x 2 Chapter 12 Describig Distributios with Numbers Chapter 12 1 Quick math overview = sum These expressios are algebraically equivalet #(x " x ) 2 = # x 2 " (# x) 2 Examples x :{ 2,3,5,6,6,8 } " x = 2 + 3+

More information

Wind Energy Explained, 2 nd Edition Errata

Wind Energy Explained, 2 nd Edition Errata Wid Eergy Explaied, d Editio Errata This summarizes the ko errata i Wid Eergy Explaied, d Editio. The errata ere origially compiled o July 6, 0. Where possible, the chage or locatio of the chage is oted

More information

The coalescent coalescence theory

The coalescent coalescence theory The coalescet coalescece theory Peter Beerli September 1, 009 Historical ote Up to 198 most developmet i populatio geetics was prospective ad developed expectatios based o situatios of today. Most work

More information

10-701/ Machine Learning Mid-term Exam Solution

10-701/ Machine Learning Mid-term Exam Solution 0-70/5-78 Machie Learig Mid-term Exam Solutio Your Name: Your Adrew ID: True or False (Give oe setece explaatio) (20%). (F) For a cotiuous radom variable x ad its probability distributio fuctio p(x), it

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departmet of Electrical Egieerig ad Computer Sciece 6.34 Discrete Time Sigal Processig Fall 24 BACKGROUND EXAM September 3, 24. Full Name: Note: This exam is closed

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Formation of A Supergain Array and Its Application in Radar

Formation of A Supergain Array and Its Application in Radar Formatio of A Supergai Array ad ts Applicatio i Radar Tra Cao Quye, Do Trug Kie ad Bach Gia Duog. Research Ceter for Electroic ad Telecommuicatios, College of Techology (Coltech, Vietam atioal Uiversity,

More information

Analysis of composites with multiple rigid-line reinforcements by the BEM

Analysis of composites with multiple rigid-line reinforcements by the BEM Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100

More information

Math 113 Exam 3 Practice

Math 113 Exam 3 Practice Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This

More information

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING Mechaical Vibratios FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING A commo dampig mechaism occurrig i machies is caused by slidig frictio or dry frictio ad is called Coulomb dampig. Coulomb dampig

More information

Random Variables, Sampling and Estimation

Random Variables, Sampling and Estimation Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig

More information

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations Geeraliig the DTFT The Trasform M. J. Roberts - All Rights Reserved. Edited by Dr. Robert Akl 1 The forward DTFT is defied by X e jω = x e jω i which = Ω is discrete-time radia frequecy, a real variable.

More information

A second look at separation of variables

A second look at separation of variables A secod look at separatio of variables Assumed that u(x,t)x(x)t(t) Pluggig this ito the wave, diffusio, or Laplace equatio gave two ODEs The solutios to those ODEs subject to boudary ad iitial coditios

More information

Introduction to Signals and Systems, Part V: Lecture Summary

Introduction to Signals and Systems, Part V: Lecture Summary EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

More information

Infinite Sequences and Series

Infinite Sequences and Series Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet

More information

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty

More information

Monte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem

Monte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem Australia Joural of Basic Applied Scieces, 5(): 097-05, 0 ISSN 99-878 Mote Carlo Optimizatio to Solve a Two-Dimesioal Iverse Heat Coductio Problem M Ebrahimi Departmet of Mathematics, Karaj Brach, Islamic

More information

Eton Education Centre JC 1 (2010) Consolidation quiz on Normal distribution By Wee WS (wenshih.wordpress.com) [ For SAJC group of students ]

Eton Education Centre JC 1 (2010) Consolidation quiz on Normal distribution By Wee WS (wenshih.wordpress.com) [ For SAJC group of students ] JC (00) Cosolidatio quiz o Normal distributio By Wee WS (weshih.wordpress.com) [ For SAJC group of studets ] Sped miutes o this questio. Q [ TJC 0/JC ] Mr Fruiti is the ower of a fruit stall sellig a variety

More information

Math 257: Finite difference methods

Math 257: Finite difference methods Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...

More information

BIOS 4110: Introduction to Biostatistics. Breheny. Lab #9

BIOS 4110: Introduction to Biostatistics. Breheny. Lab #9 BIOS 4110: Itroductio to Biostatistics Brehey Lab #9 The Cetral Limit Theorem is very importat i the realm of statistics, ad today's lab will explore the applicatio of it i both categorical ad cotiuous

More information

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio

More information

Streamfunction-Vorticity Formulation

Streamfunction-Vorticity Formulation Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first

More information

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? Harold G. Loomis Hoolulu, HI ABSTRACT Most coastal locatios have few if ay records of tsuami wave heights obtaied over various time periods. Still

More information

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to: 2.003 Egieerig Dyamics Problem Set 9--Solutio Problem 1 Fid the equatio of motio for the system show with respect to: a) Zero sprig force positio. Draw the appropriate free body diagram. b) Static equilibrium

More information

Chapter 4 : Laplace Transform

Chapter 4 : Laplace Transform 4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic

More information

THE KALMAN FILTER RAUL ROJAS

THE KALMAN FILTER RAUL ROJAS THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider

More information

MATH 10550, EXAM 3 SOLUTIONS

MATH 10550, EXAM 3 SOLUTIONS MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios

More information

Written exam Digital Signal Processing for BMT (8E070). Tuesday November 1, 2011, 09:00 12:00.

Written exam Digital Signal Processing for BMT (8E070). Tuesday November 1, 2011, 09:00 12:00. Techische Uiversiteit Eidhove Fac. Biomedical Egieerig Writte exam Digital Sigal Processig for BMT (8E070). Tuesday November, 0, 09:00 :00. (oe page) ( problems) Problem. s Cosider a aalog filter with

More information

Part I: Covers Sequence through Series Comparison Tests

Part I: Covers Sequence through Series Comparison Tests Part I: Covers Sequece through Series Compariso Tests. Give a example of each of the followig: (a) A geometric sequece: (b) A alteratig sequece: (c) A sequece that is bouded, but ot coverget: (d) A sequece

More information

Chapter 10 Partial Differential Equations and Fourier Series

Chapter 10 Partial Differential Equations and Fourier Series Math-33 Chapter Partial Differetial Equatios November 6, 7 Chapter Partial Differetial Equatios ad Fourier Series Math-33 Chapter Partial Differetial Equatios November 6, 7. Boudary Value Problems for

More information

OPTIMAL PIECEWISE UNIFORM VECTOR QUANTIZATION OF THE MEMORYLESS LAPLACIAN SOURCE

OPTIMAL PIECEWISE UNIFORM VECTOR QUANTIZATION OF THE MEMORYLESS LAPLACIAN SOURCE Joural of ELECTRICAL EGIEERIG, VOL. 56, O. 7-8, 2005, 200 204 OPTIMAL PIECEWISE UIFORM VECTOR QUATIZATIO OF THE MEMORYLESS LAPLACIA SOURCE Zora H. Perić Veljo Lj. Staović Alesadra Z. Jovaović Srdja M.

More information

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL

More information

Orthogonal Gaussian Filters for Signal Processing

Orthogonal Gaussian Filters for Signal Processing Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios

More information

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle Similarity betwee quatum mechaics ad thermodyamics: Etropy, temperature, ad Carot cycle Sumiyoshi Abe 1,,3 ad Shiji Okuyama 1 1 Departmet of Physical Egieerig, Mie Uiversity, Mie 514-8507, Japa Istitut

More information

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B. Amme 3500 : System Dyamics ad Cotrol Root Locus Course Outlie Week Date Cotet Assigmet Notes Mar Itroductio 8 Mar Frequecy Domai Modellig 3 5 Mar Trasiet Performace ad the s-plae 4 Mar Block Diagrams Assig

More information

Goodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen)

Goodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen) Goodess-of-Fit Tests ad Categorical Data Aalysis (Devore Chapter Fourtee) MATH-252-01: Probability ad Statistics II Sprig 2019 Cotets 1 Chi-Squared Tests with Kow Probabilities 1 1.1 Chi-Squared Testig................

More information

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation Temperature Distributio i a Cylidrical Core ad Coil Assembly with Heat Geeratio P. S. Yeh, Ph.D. Abstract A residetial distributio trasformer cosists of a cylidrical steel tak ad a core-ad-coil assembly

More information

Avoidance of numerical singularities in free vibration analysis of Euler-Bernoulli beams using Green functions

Avoidance of numerical singularities in free vibration analysis of Euler-Bernoulli beams using Green functions WSEAS TRASACTIOS o APPLIED ad THEORETICAL MECHAICS Goraka Štimac Ročević, Braimir Ročević, Ate Skoblar, Saji Braut Avoidace of umerical sigularities i free vibratio aalysis of Euler-Beroulli beams usig

More information

EXPERIMENT OF SIMPLE VIBRATION

EXPERIMENT OF SIMPLE VIBRATION EXPERIMENT OF SIMPLE VIBRATION. PURPOSE The purpose of the experimet is to show free vibratio ad damped vibratio o a system havig oe degree of freedom ad to ivestigate the relatioship betwee the basic

More information

2. Neutronic calculations at uranium powered cylindrical reactor by using Bessel differential equation

2. Neutronic calculations at uranium powered cylindrical reactor by using Bessel differential equation Trasworld Research Network 37/661 (), Fort P.O. Trivadrum-695 03 Kerala, Idia Nuclear Sciece ad Techology, 01: 15-4 ISBN: 978-81-7895-546-9 Editor: Turgay Korkut. Neutroic calculatios at uraium powered

More information

FIR Filter Design: Part II

FIR Filter Design: Part II EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak

More information

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c. 5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =

More information

Lesson 03 Heat Equation with Different BCs

Lesson 03 Heat Equation with Different BCs PDE & Complex Variables P3- esso 3 Heat Equatio with Differet BCs ( ) Physical meaig (SJF ) et u(x, represet the temperature of a thi rod govered by the (coductio) heat equatio: u t =α u xx (3.) where

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Indian Institute of Information Technology, Allahabad. End Semester Examination - Tentative Marking Scheme

Indian Institute of Information Technology, Allahabad. End Semester Examination - Tentative Marking Scheme Idia Istitute of Iformatio Techology, Allahabad Ed Semester Examiatio - Tetative Markig Scheme Course Name: Mathematics-I Course Code: SMAT3C MM: 75 Program: B.Tech st year (IT+ECE) ate of Exam:..7 ( st

More information

Progress In Electromagnetics Research, PIER 51, , 2005

Progress In Electromagnetics Research, PIER 51, , 2005 Progress I Electromagetics Research, PIER 51, 187 195, 2005 COMPLEX GUIDED WAVE SOLUTIONS OF GROUNDED DIELECTRIC SLAB MADE OF METAMATERIALS C. Li, Q. Sui, ad F. Li Istitute of Electroics Chiese Academy

More information

6.003 Homework #12 Solutions

6.003 Homework #12 Solutions 6.003 Homework # Solutios Problems. Which are rue? For each of the D sigals x [] through x 4 [] (below), determie whether the coditios listed i the followig table are satisfied, ad aswer for true or F

More information