MAXIMUM A POSTERIORI ESTIMATION OF ROOM IMPULSE RESPONSES. Dinei Florencio and Zhengyou Zhang. Microsoft Research, Redmond, WA, USA
|
|
- Merilyn Laurel Harris
- 6 years ago
- Views:
Transcription
1 MAXIMUM A POSTERIORI ESTIMATION OF ROOM IMPULSE RESPONSES Dinei Florencio and Zengyou Zang Microsoft Researc, Redmond, WA, USA ABSTRACT Estimating room impulse responses RIRs) as a number of applications, including personalized audio, analyzing and improving acoustic beavior of concert alls, listening room compensation, sound source localization, and many oters. RIRs ave been estimated in essentially te same fasion for te last 50 years: Compute te cross correlation between a signal played at point A, and te signal received at point B. Best results are obtained wen te signal played is wite noise, or a maximum lengt sequence. No prior knowledge is exploited in computing te RIR, wic is simply assumed to be te cross correlation between played and received signals. In contrast, researc in adaptive RIR estimation a.k.a. adaptive Acoustic Eco Cancellation) as made uge progress by among oter tings) incorporating models for te RIR. In tis paper we propose a new RIR estimation tecnique, based on a maximum a posteriori formulation. More specifically, we estimate te room reverberation time, as well as te room noise level, and use tose as priors for te RIR estimation. Comparison wit ground trut sows an average improvement of 12 db compared to traditional metods. Index Terms RIR, room acoustics, room impulse response, AEC, eco cancellation. 1. INTRODUCTION Wen a sound is played witin an enclosure e.g., a room), it propagates in very complex ways. It reflects off walls and oter surfaces, it diffracts, it is absorbed by air and by tese surfaces, and does all tat in a frequency dependent manner [1, 2]. In general, tis is furter complicated by te directionality of te sound source as well as te directionality of te capturing device e.g., a micropone or te uman ear) [3]. Tus, typically, even if te exact geometry of a room is known, modelling te transfer function between two points in a room is extremely ard, except for crude approximations. Knowledge of tat transfer function is, owever, extremely useful in a number of scenarios. Common application include acoustic eco cancellation AEC) [4, 5, 6, 7], personalized audio [8, 9], analysing and improving acoustic beaviour of concert alls [10, 11], listening room compensation [12, 13, 14], sound source localization [15, 16, 17], and many oters [18, 19, 20, 19, 21]. Due to te aforementioned difficulty in modeling acoustic beavior, room impulse responses RIR) are typically experimentally measured directly in te real environment [22]. More specifically, a sound source and a micropone are placed at eac end-point, a sound is played and recorded, and tat recorded signal is used to compute an estimate of) te RIR between tose two points 1. Furtermore, 1 In applications like AEC, te pat of interest in is between te loudif we need estimates of te RIR in points oter tan te measured ones, RIR interpolation may be needed, wic is, in itself, anoter callenging problem[23]. In tis paper we introduce a tecnique tat makes use of prior knowledge of acoustic beavior to improve te quality of te RIR estimates. Our maximum a posteriori MAP) estimate incorporates a simple model for ambient noise and an exponential decaying model for reverberation. Yet, results are quite appealing: an average 12dB improvement over traditional RIR estimation. A closely related problem is tat of adaptive RIR estimation typically referred as adaptive AEC). In adaptive AEC, one assumes te RIR is canging. Tat is muc more callenging problem, and one as to balance te adaptation parameter according to te dynamics of te system, as well as te ambient noise. State of te art adaptive AEC use models of te RIR to control step size and related parameters [24, 25, 26, 27, 28]. Wile tey greatly improve adaptation, for a fixed system as in our problem), tey underperform a simple matced filter cross correlation), wic as us) uses te wole available signal at once. In a way, our paper brings te RIR estimation in line wit te last 50 years of progress in adaptive AEC. Te remainder of tis paper is organized as follows: Section 2 introduces our notation and explains te traditional RIR estimation metod. Section 3 presents our MAP formulation, and Section 4 presents and discusses te experimental results. Finally, Section 5 presents our main conclusions and future directions. 2. TRADITIONAL RIR ESTIMATION To estimate te room impulse response RIR) between a loudspeaker and a micropone, a signal is played troug te loudspeaker and recorded at te micropone. Witout loss of generality, let us assume te first K samples of an RIR are a reasonable approximation of te RIR for te intended application. Let us denote te training signal as te N 1 vector x, and te corresponding signal received at te micropone as y. To simplify notation and witout loss of generality), we assume x as unit variance. Neglecting boundary effects to simplify te notation, we model te signal y as: y = X + a 1) were is te K 1 vector corresponding to te room impulse response, X is te N K matrix composed of K time-sifted versions and x, and a is te ambient noise, wic we assume to be independent from x. speaker and a micropone, and bot are under te control of te system. In oter scenarios, loudspeakers and/or micropones may ave to be placed at one of bot te desired end points.
2 Te traditional way of estimating te RIR is by playing a wite training sequence x, and correlating it wit te received signal y [29]. More specifically, by computing: ĥ = X T y. 2) Te reasoning for tis approximation can be observed by inserting te expression for y from 1) into 2), i.e.: ĥ = X T X + a) 3) = X T X) + X T a. 4) Mismatc db) db SNR 40dB SNR By taking te expected value of 4), we get: E{ĥ} = Rxx + rxa 5) were R xx = E{X T X} i.e., te K-lag autocorrelation matrix of x), and r xa = E{X T a} i.e., te K-lag cross-correlation vector between x and a). Note tat, if te training sequence is wite, tan R xx = I, and, since te noise is assumed independent of te signal r xa = 0. Tus, X T y is an unbiased estimate of te RIR. Furtermore, as te lengt of te training sequence increases, te estimate becomes increasingly close to te true, converging to te correct value as N. For finite training sequences, owever, tese are just approximations. Indeed, we can rewrite 4) as: ) ĥ = + X T X I + X T a 6) Tus, we see tat tere are two noise terms affecting te traditional way of estimating ĥ, bot related to te fact tat x is finite. Te first one reflects te non-diagonal nature of a finite signal covariance matrix, wile te second term originate from te random correlations between two independent) finite signals. Te first can be alleviated by making te training sequence longer, and te second can be alleviated by making te loudspeaker signal stronger, or te sequence longer. Note tat te non-diagonal elements of te covariance matrix can also be reduced by using a special kind of wite noise, called maximum lengt sequences MLS) [30, 31, 32]. Tese special pseudo-random sequences, do improve te quality of te estimates. However, tey sound like wite noise, and cannot be modified muc. As suc, tey are not appropriate for most applications were users will be present during te measurements including our target application of initial AEC filter estimation). Figure 1 sows te ĥ estimation error as a function of training sequence duration, for ambient noise levels of 0 db and 40 db SNR plotted as 10 log 10 i ĥ[i] [i])2 / i [i])2 )). Te quality of estimates vary widely wit room caracteristics, distance, etc. Te data in Figure 1 was obtained using te parameters tat are closest to our target application in AEC, and are described in more detail in Section 4. As expected, te quality of te estimate increases by rougly 3dB for every doubling of te training sequence lengt. However, to get to reasonable estimates, reasonably long sequences ave to be used. For example, at 40dB SNR level, a sequence wit 262k samples, yields only a 23dB approximation of te true RIR. In oter words, we need 16 seconds of a training sequence to get a just reasonable 23 db approximation. In many cases, owever, we cannot afford too long sequences. Wile for acoustic analysis of a concert alls, it migt be reasonable Lengt of training sequence x10 5 samples) Fig. 1. Mismatc between true and conventional RIR estimation as a functions of training sequence lengt. Results correspond to a room wit T 60 = 300ms, a micropone at 2.2m from te loudspeaker, and F s = 16KHz. Estimated filter is also 300ms long. to play several minute-long training sequences, tis is certainly not te case wen te RIR is being used as initialization for a acoustic eco cancellation in a communication system. 3. MAP ESTIMATION We now propose a maximum a posteriori MAP) formulation to estimate. Under te MAP formulation, we ave: MAP = arg max {f y)} 7) = arg max {f)fy )/fy)} 8) = arg max {f)fy )} 9) were f ) denotes te probability density function. Note tat te probability fy ) can be computed by estimating te te probability of te noise a being y X. We assume tat te ambient noise a, and te RIR are multivariate Gaussian random variables, i.e., a N 0, Σ a), and N 0, Σ ), were Σ a and Σ are te covariance matrices for a and, respectively. Note tat, wile a Gaussian is probably not te best model for, it does capture its key statistical caracteristics, and it makes te matematical treatment muc simpler tan if we assume a more complex model. Tus, based on te Gaussian assumption, we ave: f) = 1 2π)k detσ ) exp 1 2 T Σ 1 a ), 10) and 1 fy ) = 2π)k detσ exp 1 ) a) 2 y X)T Σ 1 y X). 11) Inserting 10) and 11) in 9), taking te log, and disregarding te constants and te log itself since we only care about te at te maximum, but not te maximum value per se), we get:
3 { MAP = arg max T Σ 1 } y X) T Σ 1 y X) 12) Taking te derivative of te function being maximized in te rigt side of 12) and making it equal to zero, we get: Σ 1 MAP X T Σ 1 a y X MAP ) = 0 13) Wic can be rewritten as: ) 1 MAP = X T Σ 1 a X + Σ 1 X T Σ 1 a y 14) Tus, to estimate MAP, besides te received signal y, we need estimates of te covariance matrices of te ambient noise and of te taps of te RIR. Tese are, in effect, our priors on te te ambient noise, and on te RIR itself. We discuss our metod for obtaining tem next Estimating te priors Σ a and Σ Te covariance matrix Σ a captures te caracteristics of te ambient noise. Assuming suc noise is stationary, a reasonable estimate can be easily computed during silence periods, e.g., before playing te training sequence x. Tus, in effect it is truly a priori knowledge. As a proof of concept, and a demonstration of te potential of te tecnique, we ave cosen to use a reasonably simple model for Σ. We already simplified te model by assuming te taps of are jointly Gaussian. Wile te result in 14) would be valid for a generic Σ, we furter simplify our metod by assuming te taps are independent, and tat tey decay at an exponential rate wit a time constant proportional to te room reverberation time T 60. Te independence assumption makes Σ a diagonal matrix. Ten, eac element of te diagonal corresponds to te estimate of te variance of tat tap [33, 34]: σ 2 [k] = A.10 6k/FsT 60, 15) were A is a scaling constant to account for te starting sound level. Note tat tis simplified model disregards te fact tat te direct pat does not obey te same exponential decay as te rest of te reverberation, and te sparseness of early reverberation. Te constants A and T 60 are estimated by computing te cross correlation between x and y i.e., te RIR using te traditional metod), and fitting a function in te form of 15) to te later part of te correlation to avoid te direct pat and early reflections). More specifically, we compute ĥ using Eq. 2), and estimate A and T 60 as: {A, T 60} = arg min A,T 60 K k=k s ĥ[k])2 A.10 6k/FsT 60) 2, 16) were k s is te starting point, cosen to elp ignore te direct pat and te early reflections. Our estimation of te covariance matrix of te taps of te RIR, Σ, is based on te actual received signal y. Tus, rigorously speaking, it is not prior knowledge. Te prior is effectively te exponential decay assumption. In oter words, tere are two estimation steps. In te first, we use te received signal to estimate Σ. Ten, on a second step, we estimate MAP using our estimate of Σ as prior. a Because te first step is simply an estimation of two real numbers A and T 60), we assume it is reasonably robust, and don t consider tis as a significant source of error, and focus mostly on te second step. Altoug subtle, we make tat distinction to clarify our use of te MAP term in naming our estimation algoritm. 4. EXPERIMENTAL RESULTS As we mentioned earlier, obtaining true RIRs is callenging, so te standard way of evaluating RIR estimation metods is based on syntetic data. We compare te results of our MAP metod wit te traditional correlation metod, bot using a random training sequence as well as a maximum lengt sequence. We denote tese by TRAD and MLS, respectively, wile refer to our metod as MAP. For syntetic data, we can directly compare te estimated RIR wit te actual one, and precisely measure te estimation error. Table 1 sows te results for ambient SNRs varying from 0 to 30dB, and sequences lengts of 4095, 8191, and Note tat MLS sequences are restricted to certain lengts e.g., 2 N 1), so we evaluate results only at tose lengts, even toug te MAP and TRAD metods do not ave suc restriction. Te simulated room is a rectangular box measuring meters, wit a T 60 of 300ms, and a micropone at 2.2 m from te loudspeaker. Sampling frequency is 16KHz, and te estimated filter lengt is 300ms tus K = 4800). Te true RIR was obtained using te image metod [35]. For TRAD and MAP, a different training sequence is generated eac trial, as a pseudo-random gaussian noise wit te desired lengt. For te MLS te training sequence is computed as in [31]. For all cases, te received signal y was ten computed by Eq. 4), and used as input for te estimation metod. Results are averaged over 100 trials wit different environment wite) noise. For te MAP metod, we compute A and T 60 by using 16) wit k s = 800 and K = 4800, and use a quadratic error parametric fitting based on te minsearc function from Matlab 2013Ra. As it can be seen in Table 1, te MAP estimation is better tan TRAD and MLS in all cases. Te gains over TRAD vary from 5.5 to 19.4 db, wit an average improvement of 11.8 db. Compared to MLS, te gains are smaller, varying from 1.3 to 6.8 db, wit an average improvement of 2.8 db. Figure 2 sows an example of te true and estimated RIRs by te different metods. As expected, te TRAD metod is muc noisier, and tat estimation noise does not decrease towards te tail of te RIR. Te estimation based on MLS as significantly less noise tan TRAD. However, altoug ard to see in te picture, tat noise does not decrease significantly towards te tail eiter. Finally, te MAP estimation is clearly closer to te ground trut, and te estimation noise is decreasing wit k. Tis is important, particularly in cases were te lengt of te true impulse response is not known a priori RIR estimates and AEC Wile tis paper focuses mostly on te RIR estimation itself, in tis section we quickly comment on our underlying application, as it as some implications on our metod. AEC is typically done by adaptive systems. Tis is partially due to te need to adapt to environment canges, and partially due to te need for a reasonable start up time. In systems were te direct pat between loudspeaker and micropone is known a priori e.g., speakerpones, etc) tat pat is
4 Lengt = 4095 Lengt = 8191 Lengt = SNR TRAD MLS MAP TRAD MLS MAP TRAD MLS MAP Table 1. Estimation mismatc, in db, between te true RIR and te corresponding estimates using traditional correlation TRAD), using correlation based on a MLS, and using te proposed metod MAP). MAP MLS TRAD TRUE Fig. 2. RIR and corresponding estimations for a room wit T 60 = 300ms, and 20dB SNR. From bottom to top: RIR ground trut, traditional RIR estimate, MLS RIR estimate, and MAP RIR estimate. typically measured in te factory or during design), and used as te starting point for te adaptive filter. In systems were te location of loudspeaker is not controlled e.g., PCs, game consoles, etc), tis initial estimate as to be computed after deployment. Indeed, modern communication systems were te placement of te micropone in relation to te loudspeaker is not constrained, typically play a calibration tone to obtain an initial estimate of te AEC filter. Suc calibration tone may be disguised e.g., te sinusoidal sweep at te startup sound in Skype), or not e.g., te Mozart snippet played during te initial Kinect calibration). Regardless of being disguised or not, wite noise is not a desirable sound to be played at loud volumes. For tat reason, as muc as MLS provide improvements over pseudo-random sequences, tey are very specific and cannot be modified. Wile te results reported for TRAD and MAP are also based on wite noise sequences, tese can be easily modified to accommodate virtually any sequence, by witening te received and reference) signals before computing te RIR estimates. Wile some degradation can be expected, te overall trend is similar. In oter words, wile we include results for using MLS in Table 1, we do not consider tem as a valid alternative to many applications. Furtermore, we point out tat, for cases were MLS are valid, te MAP results can be furter improved by using tem as te training sequence as well. 5. CONCLUSIONS Estimating Room impulse responses find application in a number of areas. Estimation time, owever, is a concern, as it can be too long for many applications. Tis time can be reduced by te use of MLS sequences, but tese imply additional constraints on te sequences. In tis paper we presented a new metod for estimation of room impulse responses. Te metod is based on a MAP formulation, using te observed ambient noise and an exponential reverberation decay as priors. Wile priors ave been used in te past in adaptive AEC, tey are used to control te adaptation, and do not extend naturally to fixed or initial) RIR estimation. Using our metod, improvements of 11.8 db over traditional metods, and of 2.8 db over MLS sequences were acieved. In contrast to MLS, te proposed tecnique can be easily modified to support non-wite training sequences, making it particularly appropriate for applications were users will be present, as for example in initial AEC filter estimation. We are currently using tis MAP estimation to speed up product development[36]. Additionally, we are also investigating extensions to multicannel adaptive AEC[7], implications to micropone arrays [37]and evaluating subjective quality of te results using Crowd- MOS [38, 39]
5 6. REFERENCES [1] M Junger and D Feit, Sound, structures, and teir interaction, vol. 240, MIT press Cambridge, MA, [2] M Scroeder, Computer models for concert all acoustics, American Journal of Pysics, vol. 41, pp. 461, [3] L Beranek and T Mellow, Acoustics: Sound Fields and Transducers, Academic Press, [4] A Yellepeddi and D Florencio, Sparse array-based room transfer function estimation for eco cancellation, Signal Proc. Letters, IEEE, vol. 21, no. 2, pp , [5] S Kuo and D Morgan, Active noise control: a tutorial review, Proceedings of te IEEE, vol. 87/6, [6] J Wung, T Wada, B-H Juang, B Lee, M Trott, and R Scafer, A system approac to acoustic eco cancellation in robust ands-free teleconferencing, in Proc. of WASPAA, [7] Z Zang, Q Cai, and J Stokes, Multicannel acoustic eco cancelation in multiparty spatial audio conferencing wit constrained kalman filtering, in Proc. of IWAENC, [8] M Song, C Zang, D Florencio, and H Kang, Personal 3D audio system wit loudspeakers, in Proc. of ICME, [9] Y Huang, J Cen, and J. Benesty, Immersive audio scemes, Signal Processing Magazine, IEEE, vol. 28, no. 1, pp , Jan [10] J Pätynen, S Tervo, and T Lokki, Analysis of concert all acoustics via visualizations of time-frequency and spatiotemporal responses, Te Journal of te Acoustical Society of America, vol. 133, pp. 842, [11] H Morgenstern and B Rafaely, Analysis of acoustic mimo systems in enclosed sound fields, in Proc. of ICASSP, [12] S Goetze, E Albertin, M Kallinger, A Mertins, and K Kammeyer, Quality assessment for listening-room compensation algoritms, in Proc. of ICASSP, [13] J Antons, S Arndt, R Scleicer, S Moller, D O Saugnessy, T Falk, et al., Cognitive, affective, and experience correlates of speec quality perception in complex listening conditions, in Proc. of ICASSP, [14] F Xiong, J Appell, and S Goetze, System identification for listening-room compensation by means of acoustic eco cancellation and acoustic eco suppression filters, in Proc. of ICASSP, [15] Y Rui, D Florencio, W Lam, and J Su, Sound source localization for circular arrays of directional micropones, in Proc. of ICASSP, [16] R Parisi, R Russo, M Scarpiniti, and A Uncini, Localization of audio sources by multiple binaural sensors, in Digital Signal Processing DSP), 18t IEEE Int. Conf, on, [17] S Tervo, J Patynen, and T Lokki, Acoustic reflection localization from room impulse responses, Acta Acustica, vol. 98, no. 3, pp , [18] S Tervo and T Tossavainen, 3D room geometry estimation from measured impulse responses, in Proc. of ICASSP, [19] F Ribeiro, D Florencio, D Ba, and C Zang, Geometrically constrained room modeling wit compact micropone arrays, Audio, Speec, and Language Processing, IEEE Trans. on, vol. 20, no. 5, pp , [20] D Ba, F Ribeiro, C Zang, and D Florencio, L1 regularized room modeling wit compact micropone arrays, in Proc. of ICASSP, [21] F Ribeiro, D Florencio, P Cou, and Z Zang, Auditory augmented reality: Object sonification for te visually impaired, in Proc. of MMSP, [22] J Klein, M Pollow, P Dietric, and M Vorländer, Room impulse response measurements wit arbitrary source directivity, in 40t Italian AIA) Annual Conference on Acoustics, [23] R. Mignot, L. Daudet, and F. Ollivier, Room reverberation reconstruction: Interpolation of te early part using compressed sensing, Audio, Speec, and Language Processing, IEEE Trans. on, vol. PP, no. 99, pp. 1 1, [24] C Contan, M Zeller, W Kellermann, and M Topa, Excitationdependent stepsize control of adaptive volterra filters for acoustic eco cancellation, in Proc. of EUSIPCO, [25] E Habets, S Gannot, and I Coen, Robust early eco cancellation and late eco suppression in te stft domain, in Proc. of IWAENC, [26] K Helwani, H Bucner, and S Spors, Multicannel adaptive filtering wit sparseness constraints, in Acoustic Signal Enancement; Int. Worksop on, [27] S Makino, Y Kaneda, and N Koizumi, Exponentially weigted stepsize nlms adaptive filter based on te statistics of a room impulse response, Speec and Audio Processing, IEEE Trans. on, vol. 1, no. 1, pp , [28] M Rupp, Convergence properties of adaptive equalizer algoritms, Signal Processing, IEEE Trans. on, vol. 59, no. 6, pp , [29] G Turin, An introduction to matced filters, Information Teory, IRE Trans. on, vol. 6, no. 3, pp , [30] J Vanderkooy, Aspects of mls measuring systems, Journal of te Audio Engineering Society, vol. 42, no. 4, pp , [31] M Vorländer and M Kob, Practical aspects of mls measurements in building acoustics, Applied Acoustics, vol. 52, no. 3, pp , [32] I Mateljan, Signal selection for te room acoustics measurement, in Proc. of WASPAA, [33] E Habets, Multi-cannel speec dereverberation based on a statistical model of late reverberation, in Proc. of ICASSP, [34] K Lebart, J Boucer, and P Denbig, A new metod based on spectral subtraction for speec dereverberation, Acta Acustica, vol. 87, no. 3, pp , [35] S Scimmel, M Muller, and N Dillier, A fast and accurate soebox room acoustics simulator, in Proc. of ICASSP, [36] P Moquin, K Venalainen, and D Florencio, Determination of room impulse response for syntetic data acquisition and ASR testing, Te Journal of te Acoustical Society of America, vol. 136, no. 4, pp. 2265, [37] D A Florencio and H Malvar, Multicannel filtering for optimum noise reduction in micropone arrays, in Proc. of ICASSP, [38] F Ribeiro, D Florencio, C Zang, and M Seltzer, CrowdMOS: An approac for crowdsourcing mean opinion score studies, in Proc. of ICASSP, [39] F Ribeiro, D Florencio, and V Nascimento, Crowdsourcing subjective image quality evaluation, in Proc. of ICIP, 2011.
Combining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More information3D NUMERICAL MODELING OF PARAMETRIC SPEAKER USING FINITE-DIFFERENCE TIME-DOMAIN
3D NUMERICAL MODELING OF PARAMETRIC SPEAKER USING FINITE-DIFFERENCE TIME-DOMAIN Lijun Zu Center for Energy and Geo Processing Georgia Institute of Tecnology Atlanta, GA 333 Dinei Florencio Microsoft Researc
More informationNONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura
NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL Georgeta Budura Politenica University of Timisoara, Faculty of Electronics and Telecommunications, Comm. Dep., georgeta.budura@etc.utt.ro Abstract:
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationBayesian ML Sequence Detection for ISI Channels
Bayesian ML Sequence Detection for ISI Cannels Jill K. Nelson Department of Electrical and Computer Engineering George Mason University Fairfax, VA 030 Email: jnelson@gmu.edu Andrew C. Singer Department
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationFinancial Econometrics Prof. Massimo Guidolin
CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis
More informationTHE hidden Markov model (HMM)-based parametric
JOURNAL OF L A TEX CLASS FILES, VOL. 6, NO. 1, JANUARY 2007 1 Modeling Spectral Envelopes Using Restricted Boltzmann Macines and Deep Belief Networks for Statistical Parametric Speec Syntesis Zen-Hua Ling,
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationLearning based super-resolution land cover mapping
earning based super-resolution land cover mapping Feng ing, Yiang Zang, Giles M. Foody IEEE Fellow, Xiaodong Xiuua Zang, Siming Fang, Wenbo Yun Du is work was supported in part by te National Basic Researc
More informationLAPLACIAN MATRIX LEARNING FOR SMOOTH GRAPH SIGNAL REPRESENTATION
LAPLACIAN MATRIX LEARNING FOR SMOOTH GRAPH SIGNAL REPRESENTATION Xiaowen Dong, Dorina Tanou, Pascal Frossard and Pierre Vandergeynst Media Lab, MIT, USA xdong@mit.edu Signal Processing Laboratories, EPFL,
More informationRegularized Regression
Regularized Regression David M. Blei Columbia University December 5, 205 Modern regression problems are ig dimensional, wic means tat te number of covariates p is large. In practice statisticians regularize
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationMultichannel Speech Separation and Enhancement Using the Convolutive Transfer Function
1 Multicannel Speec Separation and Enancement Using te Convolutive Transfer Function Xiaofei Li, Laurent Girin, Saron Gannot and Radu Horaud arxiv:1711.07911v2 [cs.sd] 26 Feb 2018 Abstract Tis paper addresses
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationMIMO decorrelation for visible light communication based on angle optimization
MIMO decorrelation for visible ligt communication based on angle optimization Haiyong Zang, and Yijun Zu Citation: AIP Conference Proceedings 80, 09005 (07); View online: ttps://doi.org/0.03/.4977399 View
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationLines, Conics, Tangents, Limits and the Derivative
Lines, Conics, Tangents, Limits and te Derivative Te Straigt Line An two points on te (,) plane wen joined form a line segment. If te line segment is etended beond te two points ten it is called a straigt
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More informationExperimental investigation of multiple-input multiple-output systems for sound-field analysis
Acoustic Array Systems: Paper ICA2016-158 Experimental investigation of multiple-input multiple-output systems for sound-field analysis Hai Morgenstern (a), Johannes Klein (b), Boaz Rafaely (a), Markus
More information[db]
Blind Source Separation based on Second-Order Statistics wit Asymptotically Optimal Weigting Arie Yeredor Department of EE-Systems, el-aviv University P.O.Box 3900, el-aviv 69978, Israel Abstract Blind
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationA Reconsideration of Matter Waves
A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,
More informationCORRELATION TEST OF RESIDUAL ERRORS IN FREQUENCY DOMAIN SYSTEM IDENTIFICATION
IAC Symposium on System Identification, SYSID 006 Marc 9-3 006, Newcastle, Australia CORRELATION TEST O RESIDUAL ERRORS IN REQUENCY DOMAIN SYSTEM IDENTIICATION István Kollár *, Ri Pintelon **, Joan Scouens
More informationCubic Functions: Local Analysis
Cubic function cubing coefficient Capter 13 Cubic Functions: Local Analysis Input-Output Pairs, 378 Normalized Input-Output Rule, 380 Local I-O Rule Near, 382 Local Grap Near, 384 Types of Local Graps
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationGRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES
GRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES Cristoper J. Roy Sandia National Laboratories* P. O. Box 5800, MS 085 Albuquerque, NM 8785-085 AIAA Paper 00-606 Abstract New developments
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationRobust Parametric Modeling of Speech in Additive White Gaussian Noise
Journal of Signal and Information Processing, 05, 6, 99-08 Publised Online May 05 in SciRes. ttp://www.scirp.org/journal/jsip ttp://d.doi.org/0.436/jsip.05.600 Robust Parametric Modeling of Speec in Additive
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationarxiv: v1 [physics.flu-dyn] 3 Jun 2015
A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows arxiv:156.132v1 [pysics.flu-dyn] 3 Jun 215 S. Dong Center for Computational & Applied Matematics Department
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationPredicting the Intelligibility of Noisy and Non-Linearly Processed Binaural Speech
1 Predicting te Intelligibility of Noisy and Non-Linearly Processed Binaural Speec Asger Heidemann Andersen, Jan Mark de Haan, Zeng-Hua Tan, and Jesper Jensen Abstract Objective speec intelligibility measures
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationAdaptive Neural Filters with Fixed Weights
Adaptive Neural Filters wit Fixed Weigts James T. Lo and Justin Nave Department of Matematics and Statistics University of Maryland Baltimore County Baltimore, MD 150, U.S.A. e-mail: jameslo@umbc.edu Abstract
More informationImpact of Lightning Strikes on National Airspace System (NAS) Outages
Impact of Ligtning Strikes on National Airspace System (NAS) Outages A Statistical Approac Aurélien Vidal University of California at Berkeley NEXTOR Berkeley, CA, USA aurelien.vidal@berkeley.edu Jasenka
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationDerivatives. By: OpenStaxCollege
By: OpenStaxCollege Te average teen in te United States opens a refrigerator door an estimated 25 times per day. Supposedly, tis average is up from 10 years ago wen te average teenager opened a refrigerator
More informationNew Distribution Theory for the Estimation of Structural Break Point in Mean
New Distribution Teory for te Estimation of Structural Break Point in Mean Liang Jiang Singapore Management University Xiaou Wang Te Cinese University of Hong Kong Jun Yu Singapore Management University
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationAn Adaptive Model Switching and Discretization Algorithm for Gas Flow on Networks
Procedia Computer Science 1 (21) (212) 1 1 1331 134 Procedia Computer Science www.elsevier.com/locate/procedia International Conference on Computational Science, ICCS 21 An Adaptive Model Switcing and
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationQuantization of electrical conductance
1 Introduction Quantization of electrical conductance Te resistance of a wire in te classical Drude model of metal conduction is given by RR = ρρρρ AA, were ρρ, AA and ll are te conductivity of te material,
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationThe Complexity of Computing the MCD-Estimator
Te Complexity of Computing te MCD-Estimator Torsten Bernolt Lerstul Informatik 2 Universität Dortmund, Germany torstenbernolt@uni-dortmundde Paul Fiscer IMM, Danisc Tecnical University Kongens Lyngby,
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,
More informationAn Empirical Bayesian interpretation and generalization of NL-means
Computer Science Tecnical Report TR2010-934, October 2010 Courant Institute of Matematical Sciences, New York University ttp://cs.nyu.edu/web/researc/tecreports/reports.tml An Empirical Bayesian interpretation
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationEstimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution
International Journal of Clinical Medicine Researc 2016; 3(5): 76-80 ttp://www.aascit.org/journal/ijcmr ISSN: 2375-3838 Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationOn the computation of wavenumber integrals in phase-shift migration of common-offset sections
Computation of offset-wavenumber integrals On te computation of wavenumber integrals in pase-sift migration of common-offset sections Xiniang Li and Gary F. Margrave ABSTRACT Te evaluation of wavenumber
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More informationDIRECTION ESTIMATION BASED ON SOUND INTENSITY VECTORS. Sakari Tervo
7th European Signal Processing Conference (EUSIPCO 9) Glasgow, Scotland, August 4-8, 9 DIRECTION ESTIMATION BASED ON SOUND INTENSITY VECTORS Sakari Tervo Helsinki University of Technology Department of
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationPart C : Quantum Physics
Part C : Quantum Pysics 1 Particle-wave duality 1.1 Te Bor model for te atom We begin our discussion of quantum pysics by discussing an early idea for atomic structure, te Bor model. Wile tis relies on
More informationUse of fin analysis for determination of thermal conductivity of material
RESEARCH ARTICLE OPEN ACCESS Use of fin analysis for determination of termal conductivity of material Nea Sanjay Babar 1, Saloni Suas Desmuk 2,Sarayu Dattatray Gogare 3, Snea Barat Bansude 4,Pradyumna
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationAn optimum design of robotic food handling by using Burger model
DOI 10.1007/s11370-008-0032-5 ORIGINAL RESEARCH PAPER An optimum design of robotic food andling by using Burger model Naoki Sakamoto Mitsuru Higasimori Tosio Tsuji Makoto Kaneko Received: 28 February 2008
More informationLearning to Reject Sequential Importance Steps for Continuous-Time Bayesian Networks
Learning to Reject Sequential Importance Steps for Continuous-Time Bayesian Networks Jeremy C. Weiss University of Wisconsin-Madison Madison, WI, US jcweiss@cs.wisc.edu Sriraam Natarajan Indiana University
More informationOptimization of Distribution Parameters for Estimating Probability of Crack Detection
Optimization of Distribution Parameters for Estimating Probability of Crack Detection Alexandra Coppe 1, Rapael T. Haftka 2 and Nam-Ho Kim 3 University of Florida, Gainesville, FL, 32611 and Palani Ramu
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationCharacterization of Relay Channels Using the Bhattacharyya Parameter
Caracterization of Relay Cannels Using te Battacaryya Parameter Josepine P. K. Cu, Andrew W. Ecford, and Raviraj S. Adve Dept. of Electrical and Computer Engineering, University of Toronto, Toronto, Ontario,
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationImproved Algorithms for Largest Cardinality 2-Interval Pattern Problem
Journal of Combinatorial Optimization manuscript No. (will be inserted by te editor) Improved Algoritms for Largest Cardinality 2-Interval Pattern Problem Erdong Cen, Linji Yang, Hao Yuan Department of
More informationNotes on Neural Networks
Artificial neurons otes on eural etwors Paulo Eduardo Rauber 205 Consider te data set D {(x i y i ) i { n} x i R m y i R d } Te tas of supervised learning consists on finding a function f : R m R d tat
More informationExcursions in Computing Science: Week v Milli-micro-nano-..math Part II
Excursions in Computing Science: Week v Milli-micro-nano-..mat Part II T. H. Merrett McGill University, Montreal, Canada June, 5 I. Prefatory Notes. Cube root of 8. Almost every calculator as a square-root
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More information