Title: Robust statistics from multiple pings improves noise tolerance in sonar.

Size: px
Start display at page:

Download "Title: Robust statistics from multiple pings improves noise tolerance in sonar."

Transcription

1 Title: Robust sttistics fro ultiple pings iproves noise tolernce in sonr. Authors: Nicol Neretti*, Nthn Intrtor, Leon N Cooper This work ws supported in prt by ARO (DAAD ), nd ONR (N--C- 96). N. Neretti is with the Institute for Brin nd Neurl Systes nd Physics Deprtent, Brown University, Providence RI 9. N. Intrtor is with the Institute for Brin nd Neurl Systes nd Physics Deprtent, Brown University, Providence RI 9, nd is on leve fro Tel-Aviv University. L. N Cooper is with the Physics Deprtent, Institute for Brin nd Neurl Systes nd Neuroscience Deprtent, Brown University, Providence RI 9. Edics: -DETE - Detection, Estition, Clssifiction, nd Trcking; -NPAR - Pretric nd Nonpretric Sttisticl Signl Processing. Corresponding Author: Nicol Neretti, Box 843, Brown University, Providence, RI 9 Tel. (4) , Fx (4) , E-il: Nicol_Neretti@brown.edu Abstrct The bility of sonr to detect objects is strongly influenced by the operting signl-to-noise rtio (SNR). As sound plitude decys very fst in wter this sensitivity reduces the effective sonr rnge. It is well known tht the rnge ccurcy decys for incresing levels of noise until brekpoint is reched fter which ccurcy deteriortes by severl orders of gnitude. In this pper we present robust fusion of tie-dely estites fro ultiple pings tht reduces the SNR corresponding to the ccurcy brekpoint. We show tht siple verge of the tiedely estites does not shift the brekpoint to lower SNR. The ethod cn iprove the resilience to noise of sonr syste, hence incresing its potentil rnge of opertion. I. INTRODUCTION The theory of optil receivers studies the design of pulses nd receivers to obtin optil detection in the presence of noise. Considerble work on the theoreticl ccurcy of rnge esureents hs been done in the pst strting with the Woodwrd eqution, which hs been derived using different ethods. A coprehensive description of erlier work cn be found in []. The vlidity of the Woodwrd eqution depends on vrious ssuptions in prticulr the ssuptions of very low signl-to-noise rtios (SNR), therefore it ust be reexined for the cse of low SNR's. The theory of optil receivers shows tht the tched filter receiver xiizes the output pek-signl-toen-noise (power) rtio [, 3], nd is the optiu ethod for the detection of signls in noise. Infortion bout the distnce of the trget is extrcted by coputing the tie t which the crosscorreltion between the echo nd replic of the pulse is xiu. This dely is converted into distnce by ens of the sound velocity in the prticulr ediu in considertion (e.g. wter or ir). This type of receiver is generlly referred to s coherent receiver.

2 The clssicl theory of optil receivers describes the rnge ccurcy of sonr syste vi the wellknown Woodwrd eqution, which cn be derived by using vriety of ethods [, 4-7]. For sll SNR's, one of the preters in the clssicl eqution i.e. the bndwidth hs to be odified, nd the receiver is then clled seicoherent. In [8] it ws shown tht the trnsition between the two types of behviors occurs t different SNR's depending on chrcteristics of the pulses such s bndwidth nd center frequency. With this observtion, novel syste bsed on n dptive choice of the pulse ws proposed [8]tht cn iprove ccurcy in the cse of reltively low SNR, when biguity in the choice of the correct pek of the cross-correltion function cnnot be voided. Due to the nonliner nture of the tie-dely estition proble, when the SNR drops below certin criticl vlues threshold effects tke plce. Threshold effects cn be chrcterized by shrp deteriortion of the tie-dely estitor vrince. Severl sttisticl bounds hve been used in the pst to describe the ccurcy of tched filter receiver perfornce in between the boveentioned SNR criticl vlues. These include the Crer-Ro lower bound [9-], the Brnkin bound [], nd the Ziv-Zki bound [3-6]. Such bounds hve been pplied both to the proble of tie-dely estition [5, 7-3] nd of frequency estition [4-6]. In prticulr, the Brnkin bound hs been used to define the SNR brekpoints corresponding to the chnge in behvior of the optil receiver s the SNR decreses for the cse of single ping nd single echo [7, 9, ], single ping nd ultiple echoes [, 3], nd ultiple pings nd single echo [8, ]. In this pper we nlyze the signl-to-noise brekpoint, studying the probbility of choosing the correct pek fro the noisy cross-correltion function with ethod siilr to the one in [5], where the threshold effect ws relted to the existence of highly probble outliers fr fro the true tie-dely vlue. This will enble us to extend the result to the cse of ultiple pings without priori knowledge on the tie-dely itself. This pproch is different fro tht used in previous work on the ultiple pings nd single echo cse [8, ], where the ultiple echoes for single object re obtined rtificilly vi ultiple receivers nd unique ping. In fct, the bounds found in [8, ] re vlid only if the noise t the different receivers is totlly uncorrelted or if the distnce between trnsducer nd receiver is constnt, both conditions difficult to relize in prctice. Section II introduces two working odels for the tie dely estition tht re used to copute the probbility of error for single observtion. In section III we show tht while verging of ultiple observtions iproves ccurcy, it does not increse noise tolernce. In section IV we derive theoreticl bound for the probbility of correct tie-dely estition within certin ccurcy tolernce bsed on the ode of the observtions. Section V describes set of Monte Crlo siultions confiring the theoreticl predictions. II. SINGLE PING BREAKPOINT A. A siplified odel for the utocorreltion function: δ-function. We first present siple odel where the noiseless cross-correltion function of the sonr ping is δ-like function tht is zero everywhere except fro t = where it is equl to the pulse energy (Fig. ). This odel requires pulse with infinite bndwidth; therefore it is not relizble

3 3 in prctice. However, its study will give soe insight on the echnis leding to the shrp decy in rnge ccurcy experienced by tch-filtered receivers. When white Gussin noise is dded to the returning sonr ping, then the cross correltion vector hs Gussin distribution with ultidiensionl centers t zero for ll vlues tht re outside of the bin t=, nd center t =E > for the vlue of the cross-correltion t t=, where E is the energy of the noise-free echo (see Appendix B for discussion on the vlidity of this pproxition). Let x i be the vlue t ech bin of the cross-correlted signl nd x the vlue t t=. We then hve: r r x v σ N e r T p( x, x, K, xn ) =, v = [,,, K,] () ( σ π ) Using the Gussin distribution, we cn clculte the probbility for the event x > xi, i, which corresponds to the correct echo dely estition, for given noise level σ. Denoting this probbility by α, we hve: (, ) x x x P X > X i = α = dx dx dx L dxn N e ( σ π ) i 3 + ( x ) x x σ σ = dx e e dx ( σ π ) N ( x ) N σ = e + erf( x) dx N + + π [ ] r r x v σ N () Note tht this probbility depends only on / σ, which is proportionl to the signl to noise rtion (SNR). β = α denotes the probbility of error, i.e., the probbility tht the plitude of t lest one of the points outside of the correct bin is lrger thn the plitude correct pek. B. A ore relistic odel for the utocorreltion function For ping with finite bndwidth it is ore relistic to odify the odel in eqution () to include the finite extension in tie of the utocorreltion function (Fig. ). In this odel the utocorreltion function is pproxited by piecewise constnt function, with n plitude equl to within the centrl intervl I of length, nd zero elsewhere. In this cse, we need to consider the width of the priori window of the cross-correltion. The width of the window corresponds to the sonr rnge. While the potentil error in dely estition is reduced when the width is reduced, so is the sonr rnge. If the priori window hs length of L nd the spling frequency is f s, then there will be N = N + N = L fspoints, N = f s of which will be within the centrl bin ( correct bin ), nd outside the centrl bin but within the priori window. Without loss of generlity, we consider the cse where nd re integers. N N N

4 4 We define rndo vector such tht the first N rndo vribles correspond to the plitudes of the points within the correct bin, while the lst N correspond to the plitudes of the points outside. As ws derived bove, for white Gussin noise, the joint probbility density function for the vector of n rndo vribles is given by: r r x = r σ p( x ), v =, x, K, xn N e,,,,,,, ( σ π ) 443 K 443 K N N T (3) The desired probbility of tie dely estition within the correct bin of the center of the crosscorreltion function is given by N P rg x { X j} I = P Xi > k, j N i= = N + ( X k i) x x x i i i dx k L dx dx dx N e N ( σ π ) i k k i= r r x v σ N ( x + ) x ( y ) x y σ σ σ = N dx N e ( ) e dy e dy σ π + N ( x σ ) = N e + erf x N π + erf ( x σ ) N ( ) [ ] N dx (4) Thus, for given noise level, the probbility tht the correct bin is selected is: + N ( x ) N N σ ( σ ) x N e erf [ erf x ] dx π ( ) (5) α = + + C. Accurcy brekpoint Let T be the R.V. whos probbility distribution is given by Eq. 4. For SNR vlue tht is bove SNR, we hve seen (Eq. ) tht the probbility for flling in the correct bin α is bove α. Let σ be the STD of the distribution of the echo loction in the centrl (correct) bin, nd let σ be the STD of the distribution which is outside the centrl bin. Now, suppose we sple fro the originl distribution T whose cuultive function is given by eqution (5). For n observtions, frction of α of the n observtions flls in the correct bin on verge, while βn fll outside. The stndrd devition of the distribution will then be give by: std ( T) = α std ( T ) + β std ( T ) = ασ + βσ (6)

5 5 We define the brekpoint s the level of noise for which the contribution of T to the totl error becoes doinnt. Thus, the RMSE will be significntly lrger thn the one given by the unifor distribution on I lone when σ α < σ + σ (7) We then define the probbility brekpoint to be: ( σ σ ) ( σ σ ) α = (8) + it is possible to find the SNR brekpoint s the SNR vlue for which eqution (5) equls the vlue in (8). III. WHY DOES THE MEAN FAIL? The clssicl wy to fuse infortion fro ultiple observtions is to verge. As the observtions re independent nd identiclly distributed, the centrl liit theore (CLT) iplies tht the stndrd devition (error) of the verged R.V. should be n ties sller thn the error de by ech of the n observtions seprtely. This is indeed the cse before the brekpoint. However, this process does not iprove the sitution fter the brekpoint nd, in prticulr, does not shift the brekpoint to lower SNR s. Thus, while verging iproves ccurcy, it does not increse noise tolernce. Below we provide theticl nlysis which explins why the brekpoint does not chnge. The esureent process described bove is equivlent to spling for unifor distribution on the intervl I with probbility α nd fro unifor distribution on the intervl I, with gp corresponding to I, with probbility β. Suppose we sple n ties to obtin T, T, K, Tn nd use the sple en s our estite for the dely: T n i n i = = T (9) δ δ δ On verge, α n vlues will be in the correct bin, T, T, K, Tα n, while β n will be spled fro the unifor distribution, T u, T u, K, u. Then the sple en cn be decoposed into two prts: T T β n = + n αn β n Ti Ti i= i= ()

6 6 If σδ is the stndrd devition of the δ-like distribution, nd σ u is the stndrd devition of the unifor distribution, then, pplying the centrl liit theore to the two sus in eqution () we obtin: αn std Ti = αnσ i= β n std T = βnσ i= i () The root-en-squre error will be significntly lrger thn the one given by the δ-like distribution lone when βnσ > αnσ, i.e. when σ α < σ + σ () It is observed tht this bound does not iprove with the nuber of pings nd is equl to the bound found for single ping in eqution (7). This explins why verging the echo dely estites fro ultiple pings hve not been found useful, nd probbly led to the wrong conclusion tht no iproveent cn be chieved fro ultiple pings dt. It should be noted tht it is possible to shift the brekpoint by estiting the tie dely for the verged cross-correltion functions of ll the observtions, process tht would require n extreely good lignent of the cross-correltion functions. This is bsiclly equivlent to reducing the noise level by verging the echoes [4] (see lso discussion in Appendix A), which is not relistic in sitution where the sonr is not copletely still with respect to the trget, or where it is not prcticl to store the entire echo wvefors for off line processing. However, successful ipleenttion of this verging cn be chieved by using ultiple receivers [8, ]. IV. USING THE MODE Suppose we divide the priori window into intervls of length equl to the size of the correct bin, to obtin = L/ intervls B, B, K, B, with B = I representing the correct bin (Fig. 3). If p, p, K, p re the probbilities for n estite to fll in ech of the intervls nd Y, Y, K, Y re rndo vribles representing the nuber of estites flling in ech intervl, then

7 7 j= j= Y p j j = n = (3) The joint probbility distribution for the nuber of estites in ech bin is given by the ultinoil distribution n! k PY k Y k Y k p p L p (4) k k ( =, =, K, = ) = k! k! Lk! The probbility of choosing the correct bin using the ode is the probbility tht the nuber of estites flling in the correct bin, is greter thn the nuber of estites flling in ny other bin k, i : i k n! k k P = P( Y > Yj, j ) = p p L p correct k, k, K, k k! k! Lk! k > kj, j k+ k+ L+ k = n k (5) The su in eqution (5) cn be decoposed into two prts: ) the probbility tht ore tht hlf of the n points flls into the correct bin; ) the probbility P <5% tht even if less thn hlf of the n points fll in the correct bin, the nuber of points in it is greter thn tht of ny other bin. P >5% ( correct bin) = 5% +P 5% P P > < (6) The first ter in (6) cn be written s: P k k > 5% k, k, K, k k! k! Lk! k > n/ k + k + L+ k = n k > n/ = n ( ) p k ( p p 3 L p k ) n! = p p n k L p k (7) The probbility of n estite to fll outside the correct bin is unifor over the priori window with probbilityβ, so tht the probbility for it to fll in ny intervl of size is β ( ) :

8 8 p p = ( α) ( ), j j = α (8) Substituting (8) into eqution (7) we obtin P > 5% n k ( ) ( ) n k k α α, (9) k> n/ = whereα is function of the SNR through Eq. (5). The coputtion of is ore coplicted. However, it is possible to derive n upper bound ( SNR >5% ) on the SNR brekpoint for the tie-dely ccurcy coputed by using the ode of n estites, s the SNR for which P> 5% = α, where α is given in Eq. (8): n k ( ) ( SNR> ) ( SNR> ) P <5% n k k α α 5% 5% = α () k> n/ The tighter bound corresponding to the totl probbility of choosing the correct bin given by Eq. (6) will lwys be lower thn the one derived fro Eq. (): SNRBP SNR > 5%, () thus, the brekpoint which the ode cn chieve will be for lower SNR thn the one clculted bove, which is lredy significntly better thn the brekpoint chieved by either single ping or by verging of the echo dely estites of ultiple pings (see Figure 4). V. EXPERIMENTAL RESULTS To test the theticl results presented in the previous sections, we developed set of Monte Crlo siultions using cosine pcket. We first nlyzed the histogrs of the errors in the dely estite of the idel receiver for different SNR's (figure 5). For high SNR ( db) ll the errors re sll nd follow the Woodwrd eqution tht corresponds to vlues within the centrl bin in figure 5. As the level of the noise increses, lrge errors in the estites pper. The errors re uniforly distributed over the entire -priori window, nd the reltive rtio between the correct estites (centrl bin) nd the level of the unifor distribution decreses with SNR (figures 5b, 5c, nd 5d). However, even for high levels of noise the centrl pek is significntly lrger tht the rest of the distribution. Figure 6 shows the perfornce of idel receiver for single ping. For high SNR the ccurcy follows the Woodwrd eqution corresponding to coherent idel s expected fro the theory of optil receivers. The perfornce breks for low SNR round 7 db. Figures 6b, 6c nd 6d

9 9 show the nlysis of the ccurcy brekpoint for different nuber of pings, 5 nd respectively. The blue line describes the optil ccurcy tht cn be chieved using crosscorreltion fro ultiple pings. Its breking point represents the optil breking point tht could hve been chieved using sttionry sonr nd trget, nd tht could be predicted by using the Brnkin bound s in [8, ] (see lso discussion in Appendix A). This breking point however is not ttinble, s it relies on creful registrtion of returns fro different pings. Such creful registrtion cn only be done if the distnce between object to trget is kept constnt, or if it is known for ech ping in dvnce. It cn be seen tht robust fusion of ultiple pings bsed on the ode (light blue, nd gent lines) iproves noise resiliency while retining close to optil chievble ccurcy under ultiple pings. In generl there is no significnt iproveent in the resiliency to noise when siple en of the observtions is used due the strong contintion of the distribution fro outliers (red lines). This confirs the theticl result presented in section III. Figure 7, shows sury of the results for the different ethods. The brekpoint for the verged cross-correltion function (blue squres) follows the idel curve obtined by reducing the level of the noise x explined in Appendix B (solid blue line). The brekpoint of the estite obtined fro the en does not substntilly chnge s the nuber of pings in incresed. A ore robust sttistics such s the edin iproves the resiliency to noise s the nuber of pings increses (green tringles). The best results re obtined by using the ode of the estites fro the ultiple pings (gent dionds). APPENDIX A In this section we show tht estiting the tie-dely fro the verged cross-correltion functions shifts the SNR brekpoint s if the noise level ws reduced by verging the echoes. We odel ech echo s su of two coponents: y r = u r + r η () r where u is n ttenuted replic of the pulse, nd r η is white noise with stndrd devition equl to σ η. We define the signl to noise rtion to be: E d = (3) N where E is the energy of the noise-free coponent u r esured in Ws, nd N is the spectrl density of the noise esured in W/Hz. If f = / dt is the spling frequency, then: s

10 r η dt σ N = = η fs fs r r u E = u dt = fs (4) nd the SNR cn be expressed s: d = r u σ η (5) The SNR cn be expressed in db s: SNR( db) = log d = log r u σ η (6) When verging n echoes, the stndrd devition of the gussin process generting the noise for ech echo is reduced by n due to the centrl liit theore. Thus, the signl to noise rtio becoes: r u SNRn( db) = log dn = log = ση n r u = log log n= SNR( db) log n ση (7) If BP is the brekpoint corresponding to single ping, then the brekpoint for the verge of n echoes will be: BPn = BP logn (8) APPENDIX B If x( t) is the pulse nd yt ( ) = ut ( ) + η( t) is the echo, then their cross-correltion φ xy cn be expressed s:

11 φ ( τ) = xt ( ) yt ( + τ) dt= xy = (9) = x( tut ) ( + τ) dt+ xt ( ) η( t+ τ) dt = φ ( τ) + φ ( τ) xu xη Moreover if the noise is white then the utocorreltion of the noise: R ( t, t ) = σδ( t t ) (3) ηη η If we tret x( t ) s filter cting on the noise η( t) to generte filtered version of the noise φ() t = φx η, then the utocorreltion R φφ is: η η ( ) ( ) ( ) R ( t, t ) = R t α, t β x α x β dαdβ = φφ ηη + + η xx ( ) ( ) ( ) = σ δ t β t + α x α x β dαdβ = ( ) ( ) ( t) = σ x β + t x β dβ = = σφ (3) Thus, the noise on top of the cross-correltion function between the pulse nd the noise-free coponent of the echo is not white, since its utocorreltion function is not δ-function. However, since the support of φ xx is sll copred to the length of the -priori window, we pproxite R φφ with σδ( t), nd consider the noise to be white. η REFERENCES [] M. I. Skolnik, Introduction to Rdr Systes, st ed: McGrw-Hill Book Copny, 96. [] M. I. Skolnik, Introduction to Rdr Systes, 3 rd ed: McGrw-Hill,. [3] D. O. North, "An Anlysis of the Fctors which Deterine Signl/Noise Discriintion in Pulse-crrier Systes," RCA, Tech. Rept. PTR-6C, June [4] A. J. Mllinckrodt nd T. E. Sollenberger, "Optiu-pulse-tie Deterintion," IRE Trns., vol. PGIT-3, pp. 5-59, 954. [5] D. Slepin, "Estition of Signl Preters in the Presence of Noise," IRE Trns., vol. PGIT-3, pp , 954.

12 [6] P. M. Woodwrd, Probbility nd Infortion Theory, with Applictions to Rdr. New York: McGrw-Hill Book Copny, Inc., 953. [7] M. I. Skolnik, "Theoreticl Accurcy of Rdr Mesureents," IRE Trns., vol. ANE-7, pp. 3-9, 96. [8] N. Neretti, N. Intrtor, nd L. N. Cooper, "Adptive pulse optiiztion for iproved sonr rnge ccurcy," IEEE Signl Processing Letters, 3. [9] H. Crer, Mtheticl ethods of sttistics. Princeton, NJ: Princeton Univ. Press, 946. [] C. R. Ro, "Infortion nd ccurcy ttinble in the estition of sttisticl preters," Bull. Clcutt Mth. Soc., vol. 37, pp. 8-9, 945. [] R. A. Fisher, "On the theticl foundtions of theoreticl sttistics," Phil. Trns. Roy. Soc., vol., pp. 39, 9. [] E. W. Brnkin, "Loclly best unbised estites," Ann. Mth. Stt., vol., pp , 946. [3] S. Bellini nd G. Trtr, "Bounds on error in signl preter estition," IEEE Trns. Coun., vol. COM-, pp , 974. [4] D. Chzn, M. Zki, nd J. Ziv, "Iproved lower bound on signl preter estition," IEEE Trns. Infortion Theory, vol. IT-, pp. 9-93, 975. [5] L. P. Seidn, "Perfornce liittions nd error clcultions for preter estition," Proc. IEEE, vol. 58, pp , 97. [6] J. Ziv nd M. Zki, "Soe lower bounds on signl preter estition," IEEE Trns. Infortion Theory, vol. IT-5, pp , 969. [7] I. Reuven nd H. Messer, "A Brnkin-type lower bound on the estition error of hybrid preter vector," IEEE Trns. Infortion Theory, vol. 43, pp , 997. [8] S.-K. Chow nd P. M. Schultheiss, "Dely estition using nrrow-bnd processes," IEEE Trns. ASSP, vol. ASSP-9, pp , 98. [9] R. J. McAuly nd E. M. Hofstetter, "Brnkin bounds on preter estition," IEEE Trns. Infortion Theory, vol. IT-7, pp , 97. [] R. J. McAuly nd L. P. Seidn, "A useful for of the Brnkin lower bound nd its ppliction to PPM threshold nlysis," IEEE Trns. Infortion Theory, vol. IT-5, pp , 969. [] J. Tbrikin nd J. L. Krolik, "Brnkin bounds for source locliztion in n uncertin ocen environent," IEEE Trns. Signl Processing, vol. 47, pp , 999. [] A. Zeir nd P. M. Schultheiss, "Relizble lower bounds for tie dely estition: Prt - Threshold phenoen," IEEE Trns. Signl Processing, vol. 4, pp. -7, 994. [3] A. Zeir nd P. M. Schultheiss, "Relizble lower bounds for tie dely estition," IEEE Trns. Signl Processing, vol. 4, pp. 3-33, 993. [4] L. Knockert, "The Brnkin bound nd threshold behvior in frequency estition," IEEE Trns. Signl Processing, vol. 47, pp , 997. [5] D. C. Rife nd R. R. Boorstyn, "Single-tone preter estition fro discrete-tie observtions," IEEE Trns. Infortion Theory, vol. IT-, pp , 974. [6] B. Jes, B. D. O. Anderson, nd R. C. Willison, "Chrcteriztion of threshold for single tone xiu likelihood frequency estition," IEEE Trns. Signl Processing, vol. 43, pp. 87-8, 995.

13 3 Figure Cptions Fig. - Siple odel where the noiseless cross-correltion function of the sonr ping is δ-like function tht is zero everywhere (within the priori window of length L) except fro t = where it is equl to the pulse energy (top). The botto figure shows the probbility of selecting given tie loction in the cross-correltion function predicted by this odel. Fig. - A ore relistic odel for the noiseless cross-correltion function includes the finite extension in tie of the utocorreltion function. In this odel the noiseless cross-correltion function is pproxited by piecewise constnt function with n plitude equl to within the centrl intervl I of length nd zero elsewhere (within the priori window of length L). The botto figure shows the probbility of selecting given tie loction in the crosscorreltion function predicted by this odel. Fig. 3 Model used in the coputtion of the probbility of correct bin choice in the cse of ultiple pings. We divide the priori window into intervls of length equl to the size of the correct bin, to obtin = L/ intervls B, B, K, B, with B = I representing the correct bin. Fig. 4 Probbility of king the correct choice s function of SNR for different nubers of pings. The rrows indicte upper bounds on the SNR brekpoints. Fig. 5 Histogrs of the errors in the dely estite in Monte Crlo siultion for different SNR's. For high SNR ( db) ll the errors re sll nd follow the Woodwrd eqution tht corresponds to vlues within the centrl bin in figure (). As the level of the noise increses, lrge errors in the estites pper. The errors re uniforly distributed over the entire -priori window, nd the reltive rtio between the correct estites (centrl bin) nd the level of the unifor distribution decreses with SNR, see figures (b), (c), nd (d). However, even for high levels of noise the centrl pek is significntly lrger tht the rest of the distribution. Fig. 6 RMSE s function of SNR nd nuber of pings,, nd respectively for the Cosine Pcket. Notice how the SNR brekpoint for the verge of ultiple pings (red line) does not decrese with the nuber of pings. Fig. 7 Brekpoint in db s function of nuber of pings for different ethods. The solid line corresponds to the noise reduction due to verging s discussed in Appendix C.

14 4 Figures L α α N L Fig. L α α N L Fig.

15 5 B B 3 B 4 B 5 B Fig. 3.8 Ping (BLACK) Pings (RED) 5 Pings (BLUE) 5 Pings (GREEN) PROB. CORRECT SNR (db) Fig. 4

16 6 () (b) (c) (d) Percent (%) Error in Pek Loction (s) Fig. 5

17 7 Ping Pings Pings Pings -4-4 RMSE (s) -6-8 Ave. XC Men Medin Mode SNR(dB) Fig. 6

18 8 5 Men Brekpoint (db) 5 Ave. XC Medin Mode Idel Nuber of pings Fig. 7

PULSE-TRAIN BASED TIME-DELAY ESTIMATION IMPROVES RESILIENCY TO NOISE

PULSE-TRAIN BASED TIME-DELAY ESTIMATION IMPROVES RESILIENCY TO NOISE PULSE-TRAIN BASED TIME-DELAY ESTIMATION IMPROVES RESILIENCY TO NOISE 1 Nicola Neretti, 1 Nathan Intrator and 1,2 Leon N Cooper 1 Institute for Brain and Neural Systes, Brown University, Providence RI 02912.

More information

ECONOMETRIC THEORY. MODULE IV Lecture - 16 Predictions in Linear Regression Model

ECONOMETRIC THEORY. MODULE IV Lecture - 16 Predictions in Linear Regression Model ECONOMETRIC THEORY MODULE IV Lecture - 16 Predictions in Liner Regression Model Dr. Shlbh Deprtent of Mthetics nd Sttistics Indin Institute of Technology Knpur Prediction of vlues of study vrible An iportnt

More information

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY YIFEI PAN, MEI WANG, AND YU YAN ABSTRACT We estblish soe uniqueness results ner 0 for ordinry differentil equtions of the

More information

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY The Atwood Mchine OBJECTIVE To derive the ening of Newton's second lw of otion s it pplies to the Atwood chine. To explin how ss iblnce cn led to the ccelertion of the syste. To deterine the ccelertion

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

Exponents and Powers

Exponents and Powers EXPONENTS AND POWERS 9 Exponents nd Powers CHAPTER. Introduction Do you know? Mss of erth is 5,970,000,000,000, 000, 000, 000, 000 kg. We hve lredy lernt in erlier clss how to write such lrge nubers ore

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

OXFORD H i g h e r E d u c a t i o n Oxford University Press, All rights reserved.

OXFORD H i g h e r E d u c a t i o n Oxford University Press, All rights reserved. Renshw: Mths for Econoics nswers to dditionl exercises Exercise.. Given: nd B 5 Find: () + B + B 7 8 (b) (c) (d) (e) B B B + B T B (where 8 B 6 B 6 8 B + B T denotes the trnspose of ) T 8 B 5 (f) (g) B

More information

Monte Carlo method in solving numerical integration and differential equation

Monte Carlo method in solving numerical integration and differential equation Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The

More information

APPROXIMATE INTEGRATION

APPROXIMATE INTEGRATION APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY YIFEI PAN, MEI WANG, AND YU YAN ABSTRACT We study ordinry differentil equtions of the type u n t = fut with initil conditions

More information

PHYS 601 HW3 Solution

PHYS 601 HW3 Solution 3.1 Norl force using Lgrnge ultiplier Using the center of the hoop s origin, we will describe the position of the prticle with conventionl polr coordintes. The Lgrngin is therefore L = 1 2 ṙ2 + 1 2 r2

More information

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.) MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give

More information

r = cos θ + 1. dt ) dt. (1)

r = cos θ + 1. dt ) dt. (1) MTHE 7 Proble Set 5 Solutions (A Crdioid). Let C be the closed curve in R whose polr coordintes (r, θ) stisfy () Sketch the curve C. r = cos θ +. (b) Find pretriztion t (r(t), θ(t)), t [, b], of C in polr

More information

15. Quantisation Noise and Nonuniform Quantisation

15. Quantisation Noise and Nonuniform Quantisation 5. Quntistion Noise nd Nonuniform Quntistion In PCM, n nlogue signl is smpled, quntised, nd coded into sequence of digits. Once we hve quntised the smpled signls, the exct vlues of the smpled signls cn

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph

More information

Acceptance Sampling by Attributes

Acceptance Sampling by Attributes Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

PHY 5246: Theoretical Dynamics, Fall Assignment # 5, Solutions. θ = l mr 2 = l

PHY 5246: Theoretical Dynamics, Fall Assignment # 5, Solutions. θ = l mr 2 = l PHY 546: Theoreticl Dynics, Fll 15 Assignent # 5, Solutions 1 Grded Probles Proble 1 (1.) Using the eqution of the orbit or force lw d ( 1 dθ r)+ 1 r = r F(r), (1) l with r(θ) = ke αθ one finds fro which

More information

Continuous Random Variables

Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Lecture 3 Gaussian Probability Distribution

Lecture 3 Gaussian Probability Distribution Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Contact Analysis on Large Negative Clearance Four-point Contact Ball Bearing

Contact Analysis on Large Negative Clearance Four-point Contact Ball Bearing Avilble online t www.sciencedirect.co rocedi ngineering 7 0 74 78 The Second SR Conference on ngineering Modelling nd Siultion CMS 0 Contct Anlysis on Lrge Negtive Clernce Four-point Contct Bll Bering

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve. Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F

More information

Ideal Gas behaviour: summary

Ideal Gas behaviour: summary Lecture 4 Rel Gses Idel Gs ehviour: sury We recll the conditions under which the idel gs eqution of stte Pn is vlid: olue of individul gs olecules is neglected No interctions (either ttrctive or repulsive)

More information

Scientific notation is a way of expressing really big numbers or really small numbers.

Scientific notation is a way of expressing really big numbers or really small numbers. Scientific Nottion (Stndrd form) Scientific nottion is wy of expressing relly big numbers or relly smll numbers. It is most often used in scientific clcultions where the nlysis must be very precise. Scientific

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

Tests for the Ratio of Two Poisson Rates

Tests for the Ratio of Two Poisson Rates Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

The steps of the hypothesis test

The steps of the hypothesis test ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

different methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).

different methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s). Mth 1A with Professor Stnkov Worksheet, Discussion #41; Wednesdy, 12/6/217 GSI nme: Roy Zho Problems 1. Write the integrl 3 dx s limit of Riemnn sums. Write it using 2 intervls using the 1 x different

More information

An Introduction to Trigonometry

An Introduction to Trigonometry n Introduction to Trigonoetry First of ll, let s check out the right ngled tringle below. The LETTERS, B & C indicte the ngles nd the letters, b & c indicte the sides. c b It is iportnt to note tht side

More information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth 3 Exm Prctice Februry 8, 03 Exm will cover 7.4, 7.5, 7.7, 7.8, 8.-3 nd 8.5. Plese note tht integrtion skills lerned in erlier sections will still be needed for the mteril in 7.5, 7.8 nd chpter 8. This

More information

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1) 3e. Introduction Lecture 3e Rectngulr wveguide So fr in rectngulr coordintes we hve delt with plne wves propgting in simple nd inhomogeneous medi. The power density of plne wve extends over ll spce. Therefore

More information

Chapter 5 : Continuous Random Variables

Chapter 5 : Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll

More information

AN ORIGINAL METHOD FOR URBAN TRAFFIC NOISE PREDICTION

AN ORIGINAL METHOD FOR URBAN TRAFFIC NOISE PREDICTION AN ORIGINA METHOD FOR URBAN TRAFFIC NOISE REDICTION Federico Rossi 1, Uberto Di Mtteo, Sofi Sioni 3 1 University of erugi Industril Engineering Deprtent oc. enti Bss 1, 05100, Terni, Itly frossi@unipg.it

More information

Discussion Introduction P212, Week 1 The Scientist s Sixth Sense. Knowing what the answer will look like before you start.

Discussion Introduction P212, Week 1 The Scientist s Sixth Sense. Knowing what the answer will look like before you start. Discussion Introduction P1, Week 1 The Scientist s Sith Sense As scientist or engineer, uch of your job will be perforing clcultions, nd using clcultions perfored by others. You ll be doing plenty of tht

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

7 - Continuous random variables

7 - Continuous random variables 7-1 Continuous rndom vribles S. Lll, Stnford 2011.01.25.01 7 - Continuous rndom vribles Continuous rndom vribles The cumultive distribution function The uniform rndom vrible Gussin rndom vribles The Gussin

More information

7-1: Zero and Negative Exponents

7-1: Zero and Negative Exponents 7-: Zero nd Negtive Exponents Objective: To siplify expressions involving zero nd negtive exponents Wr Up:.. ( ).. 7.. Investigting Zero nd Negtive Exponents: Coplete the tble. Write non-integers s frctions

More information

Chapter 1: Fundamentals

Chapter 1: Fundamentals Chpter 1: Fundmentls 1.1 Rel Numbers Types of Rel Numbers: Nturl Numbers: {1, 2, 3,...}; These re the counting numbers. Integers: {... 3, 2, 1, 0, 1, 2, 3,...}; These re ll the nturl numbers, their negtives,

More information

fractions Let s Learn to

fractions Let s Learn to 5 simple lgebric frctions corne lens pupil retin Norml vision light focused on the retin concve lens Shortsightedness (myopi) light focused in front of the retin Corrected myopi light focused on the retin

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Deprtment 8.044 Sttisticl Physics I Spring Term 03 Problem : Doping Semiconductor Solutions to Problem Set # ) Mentlly integrte the function p(x) given in

More information

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014 SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 014 Mrk Scheme: Ech prt of Question 1 is worth four mrks which re wrded solely for the correct nswer.

More information

Applications of Bernoulli s theorem. Lecture - 7

Applications of Bernoulli s theorem. Lecture - 7 Applictions of Bernoulli s theorem Lecture - 7 Prcticl Applictions of Bernoulli s Theorem The Bernoulli eqution cn be pplied to gret mny situtions not just the pipe flow we hve been considering up to now.

More information

1 Probability Density Functions

1 Probability Density Functions Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

Discussion Question 1A P212, Week 1 P211 Review: 2-D Motion with Uniform Force

Discussion Question 1A P212, Week 1 P211 Review: 2-D Motion with Uniform Force Discussion Question 1A P1, Week 1 P11 Review: -D otion with Unifor Force The thetics nd phsics of the proble below re siilr to probles ou will encounter in P1, where the force is due to the ction of n

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the

More information

Operations with Polynomials

Operations with Polynomials 38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils

More information

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

Influence of Mean Stress

Influence of Mean Stress Influence of Men tress Discussion hs been liited to copletely reversible stress thus fr. Mening = 0 However, there re ny instnces of dynic loding when en stress is nonzero. Men tresses Incresing en stress

More information

The Spring. Consider a spring, which we apply a force F A to either stretch it or compress it

The Spring. Consider a spring, which we apply a force F A to either stretch it or compress it The Spring Consider spring, which we pply force F A to either stretch it or copress it F A - unstretched -F A 0 F A k k is the spring constnt, units of N/, different for different terils, nuber of coils

More information

INVESTIGATION OF THERMAL PROPERTIES OF SOIL BY IMPULSE METHOD Juraj Veselský

INVESTIGATION OF THERMAL PROPERTIES OF SOIL BY IMPULSE METHOD Juraj Veselský THERMOPHYSICS 6 October 6 5 INVESTIGATION OF THERMAL PROPERTIES OF SOIL BY IMPULSE METHOD Jurj Veselský Fculty of Civil Engineering STU Brtislv, Rdlinského, 8 68 Brtislv Abstrct The therl properties of

More information

V E L O C I T Y a n d V E L O C I T Y P R E S S U R E I n A I R S Y S T E M S

V E L O C I T Y a n d V E L O C I T Y P R E S S U R E I n A I R S Y S T E M S V E L O C I T Y n d V E L O C I T Y R E S S U R E I n A I R S Y S T E M S A nlysis of fluid systes using ir re usully done voluetric bsis so the pressure version of the Bernoulli eqution is used. This

More information

Math 42 Chapter 7 Practice Problems Set B

Math 42 Chapter 7 Practice Problems Set B Mth 42 Chpter 7 Prctice Problems Set B 1. Which of the following functions is solution of the differentil eqution dy dx = 4xy? () y = e 4x (c) y = e 2x2 (e) y = e 2x (g) y = 4e2x2 (b) y = 4x (d) y = 4x

More information

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8 Mth 3 Fll 0 The scope of the finl exm will include: Finl Exm Review. Integrls Chpter 5 including sections 5. 5.7, 5.0. Applictions of Integrtion Chpter 6 including sections 6. 6.5 nd section 6.8 3. Infinite

More information

221A Lecture Notes WKB Method

221A Lecture Notes WKB Method A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e

More information

Terminal Velocity and Raindrop Growth

Terminal Velocity and Raindrop Growth Terminl Velocity nd Rindrop Growth Terminl velocity for rindrop represents blnce in which weight mss times grvity is equl to drg force. F 3 π3 ρ L g in which is drop rdius, g is grvittionl ccelertion,

More information

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System SPIE Aerosense 001 Conference on Signl Processing, Sensor Fusion, nd Trget Recognition X, April 16-0, Orlndo FL. (Minor errors in published version corrected.) A Signl-Level Fusion Model for Imge-Bsed

More information

SEA based underwater noise radiation prediction method for ship s conception design

SEA based underwater noise radiation prediction method for ship s conception design INTER-NOISE 016 SEA bsed underwter noise rdition prediction ethod for ship s conception design Xueren WANG 1 ; Chung WU ; Xuhong IAO 3 ; Kifu YE 4 1 Nvl Acdey of Arent, Chin Hrbin Engineering University,

More information

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018 Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.

More information

Math 31S. Rumbos Fall Solutions to Assignment #16

Math 31S. Rumbos Fall Solutions to Assignment #16 Mth 31S. Rumbos Fll 2016 1 Solutions to Assignment #16 1. Logistic Growth 1. Suppose tht the growth of certin niml popultion is governed by the differentil eqution 1000 dn N dt = 100 N, (1) where N(t)

More information

4.1. Probability Density Functions

4.1. Probability Density Functions STT 1 4.1-4. 4.1. Proility Density Functions Ojectives. Continuous rndom vrile - vers - discrete rndom vrile. Proility density function. Uniform distriution nd its properties. Expected vlue nd vrince of

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

4.4 Areas, Integrals and Antiderivatives

4.4 Areas, Integrals and Antiderivatives . res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order

More information

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30 Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function

More information

Sufficient condition on noise correlations for scalable quantum computing

Sufficient condition on noise correlations for scalable quantum computing Sufficient condition on noise correltions for sclble quntum computing John Presill, 2 Februry 202 Is quntum computing sclble? The ccurcy threshold theorem for quntum computtion estblishes tht sclbility

More information

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d Interntionl Industril Informtics nd Computer Engineering Conference (IIICEC 15) Driving Cycle Construction of City Rod for Hybrid Bus Bsed on Mrkov Process Deng Pn1,, Fengchun Sun1,b*, Hongwen He1, c,

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

EFFECTIVE BUCKLING LENGTH OF COLUMNS IN SWAY FRAMEWORKS: COMPARISONS

EFFECTIVE BUCKLING LENGTH OF COLUMNS IN SWAY FRAMEWORKS: COMPARISONS IV EFFETIVE BUING ENGTH OF OUMN IN WAY FRAMEWOR: OMARION Ojectives In the present context, two different pproches re eployed to deterine the vlue the effective uckling length eff n c of colun n c out the

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

Math 113 Exam 1-Review

Math 113 Exam 1-Review Mth 113 Exm 1-Review September 26, 2016 Exm 1 covers 6.1-7.3 in the textbook. It is dvisble to lso review the mteril from 5.3 nd 5.5 s this will be helpful in solving some of the problems. 6.1 Are Between

More information

The graphs of Rational Functions

The graphs of Rational Functions Lecture 4 5A: The its of Rtionl Functions s x nd s x + The grphs of Rtionl Functions The grphs of rtionl functions hve severl differences compred to power functions. One of the differences is the behvior

More information

Conducting Ellipsoid and Circular Disk

Conducting Ellipsoid and Circular Disk 1 Problem Conducting Ellipsoid nd Circulr Disk Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 (September 1, 00) Show tht the surfce chrge density σ on conducting ellipsoid,

More information

Second degree generalized gauss-seidel iteration method for solving linear system of equations. ABSTRACT

Second degree generalized gauss-seidel iteration method for solving linear system of equations. ABSTRACT Ethiop. J. Sci. & Technol. 7( 5-, 0 5 Second degree generlized guss-seidel itertion ethod for solving liner syste of equtions Tesfye Keede Bhir Dr University, College of Science, Deprtent of Mthetics tk_ke@yhoo.co

More information

Discrete Least-squares Approximations

Discrete Least-squares Approximations Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve

More information

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:

More information

Section 14.3 Arc Length and Curvature

Section 14.3 Arc Length and Curvature Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

Name Solutions to Test 3 November 8, 2017

Name Solutions to Test 3 November 8, 2017 Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier

More information

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3 2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is

More information

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction Lesson : Logrithmic Functions s Inverses Prerequisite Skills This lesson requires the use of the following skills: determining the dependent nd independent vribles in n exponentil function bsed on dt from

More information

Section 6: Area, Volume, and Average Value

Section 6: Area, Volume, and Average Value Chpter The Integrl Applied Clculus Section 6: Are, Volume, nd Averge Vlue Are We hve lredy used integrls to find the re etween the grph of function nd the horizontl xis. Integrls cn lso e used to find

More information

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of

More information

MATH 144: Business Calculus Final Review

MATH 144: Business Calculus Final Review MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information