Gaussian proposal density using moment matching in SMC methods
|
|
- Whitney Melton
- 5 years ago
- Views:
Transcription
1 Stat Comput (009 19: DOI.07/s Gaussian poposal density using moment matching in SMC methods S. Saha P.K. Mandal Y. Boes H. Diessen A. Bagchi Received: 3 Apil 007 / Accepted: 30 June 008 / Published online: 31 July 008 The Autho(s 008. This aticle is published with open access at Spingelin.com Abstact In this aticle we intoduce a new Gaussian poposal distibution to be used in conjunction with the sequential Monte Calo (SMC method fo solving non-linea filteing poblems. The poposal, in line with the ecent tend, incopoates the cuent obsevation. The intoduced poposal is chaacteized by the exact moments obtained fom the dynamical system. This is in contast with ecent wos whee the moments ae appoximated eithe numeically o by lineaizing the obsevation model. We show futhe that the newly intoduced poposal pefoms bette than othe simila poposal functions which also incopoate both state and obsevations. Keywods Bayesian filteing Nonlinea dynamic system Sequential Monte Calo methods Paticle filteing Impotance sampling Moment matching This wo was suppoted by a eseach gant fom THALES Nedeland BV. S. Saha ( P.K. Mandal A. Bagchi Depatment of Applied Mathematics, Univesity of Twente, 7500 AE, Enschede, The Nethelands s.saha@math.utwente.nl P.K. Mandal p..mandal@math.utwente.nl A. Bagchi a.bagchi@math.utwente.nl Y. Boes H. Diessen THALES Nedeland BV, Haasbegestaat 49, 7554 PA, Hengelo, The Nethelands Y. Boes yvo.boes@nl.thalesgoup.com H. Diessen hans.diessen@nl.thalesgoup.com 1 Intoduction Conside a nonlinea dynamic system given by x f(x 1,w, (1 y h(x,v, 1,,... ( whee (x ae the unobsevable system values (the state with (nown initial pio density p(x 0 p(x 0 x 1 and (y ae the obseved values (the measuements. The pocess noises (w ae assumed to be independent of the measuement noises (v. The poblem is to estimate the unobseved system value x n in some optimal manne fom all the obsevations y 1:n (y 1,y,...,y n,uptotimen, o equivalently, estimate the conditional density (also nown as filteed density p(x n y 1:n. An analytical solution can be found only fo a few special cases such as when both the system and obsevation equations (1 ( ae linea and the noise pocesses ae Gaussian (Kalman filte. Fo geneal models, analytical appoximations such as Extended Kalman filte and Gaussian sum filte (Andeson and Mooe 1979; Jazwinsi 1970; Bagchi 1993 and othe appoximate methods using numeical integation (Kitagawa 1987, the unscented Kalman filte (Julie and Uhlmann 1997; Wan and van de Mewe 000 and the Gaussian quadatue Kalman filte (Ito and Xiong 000 ae poposed in the liteatue. Simulation based sequential Monte Calo (SMC methods, also nown as Paticle Filtes (PF s, povide the taget filte density p(x n y 1:n in the fom of a cloud of paticles (Handschin and Mayne 1969; Aashi and Kumamoto 1975; Godon et al. 1993; West1993; Pitt and Shephad 1999; Doucet et al. 001; Aulampalam et al. 00. The biggest advantage of the SMC method is that it can easily adapt to
2 04 Stat Comput (009 19: nonlineaity in the model and/o non-gaussian noises. The efficiency of the PF algoithm depends on the so-called impotance function π(, also often efeed to as the poposal distibution, used to geneate the paticles. It has been shown in (Doucet et al. 000 that the impotance function of the fom π p(x x 1,y is optimal in a cetain sense. Thee ae, howeve, two majo, in pactice pohibitive, dawbacs fo using this impotance function. Fistly, dawing samples accoding to p(x x 1,y is, in geneal, difficult. Secondly, it is also difficult to get an analytical expession which is needed fo the weight update. In this aticle we popose as the impotance function a Gaussian appoximation of p(x x 1,y by matching the exact moments up to second ode of the distibution of (x,y conditional on x 1. Recently, Doucet et al. (000 and Guo et al. (005 also poposed simila impotance functions. Howeve, Doucet et al. (000 use the lineaised appoximation of the obsevation model to calculate the moments, while Guo et al. (005 appoximate the moments by diffeent numeical methods such as the Gaussian-Hemite quadatue ule o the Julie-Uhlmann quadatue ule. Besides the possibility that the use of exact moments may lead to bette appoximation of the optimal poposal distibution, ou method has one distinct advantage ove the one in Guo et al. (005, namely that, it is computationally less demanding. We also note that, the dynamical systems consideed by Doucet et al. (000 and Guo et al. (005 ae with additive Gaussian noise pocesses, wheeas ou method would wo fo moe geneal models, so long as the obsevation model is a polynomial. Fo compaison we have also consideed unscented paticle filte (van de Mewe et al. 000 which, too, wos fo geneal dynamical systems. Ou expeimental esults with additive Gaussian noise pocesses show that the oveall pefomance of ou poposal function is bette than that of the othe poposals, consideing the tade off between the RMSE and the computational load. The est of the aticle is oganized as follows. In Sect. the geneal SMC method is eviewed vey biefly and the ole of the impotance function is discussed. We descibe ou poposed impotance function using the exact moments in Sect. 3. Constuction of othe impotance functions as poposed by Doucet et al. (000, Guo et al. (005 and van de Mewe et al. (000 ae biefly eviewed in Sect. 4. Implementation issue of ou poposed method is discussed in Sect. 5. Section 6 contains the numeical compaison esults of these methods based on two examples one with polynomial (Sect. 6.1 and the othe with non-polynomial (Sect. 6. obsevation equation. Finally, Sect. 7 concludes the aticle. Geneal SMC method and impotance function Suppose that the system dynamics ae given by (1 (. The sequential Monte Calo method is based on impotance sampling and allows one to estimate ecusively in time the distibution function p(x y 1:. The estimate is given in the fom of a weighted paticle cloud {(x (i, w(i, i 1,...,N}. Given the obsevations up to the cuent time (y 1: and the paticles up to time 1(x (i 0: 1, x(i is dawn accoding to a nomalized impotance function π(x x (i 0: 1,y 1: whose suppot includes that of the tue posteio and subsequently the impotance weights ae updated. Fo a full account we efe the eades to Aulampalam et al. (00. Usually in pactice, the impotance function is taen to be the tansition density, i.e., π(x x (i 0: 1,y 1: p(x x (i 1 because it is easily available fom the model. It is nown that this algoithm suffes fom the degeneacy poblem, that is to say, the vaiance of the impotance weights can only incease ove time. It has been shown by Doucet et al. (000 that an impotance function of the fom π(x x (i 0: 1,y 1: p(x x (i 1,y addesses this issue by minimizing the vaiance of the (unnomalized impotance weight w (i conditional upon x (i 0: 1 and y 1:. In geneal, though, this choice of the impotance function is not pactical as it is difficult to geneate samples fom this distibution. Futhemoe, one needs an analytical expession of the impotance function to be used in the weight update equation, which is also geneally difficult with this choice of the impotance function. 3 Impotance function based on exact moment matching (EMM Suppose the system dynamics ae given by (1 (. We futhe assume the following. Assumption A All the moments of (x,y conditional on x 1 up to second ode, i.e., E(x x 1, E(x x T x 1, E(y x 1, E(y y T x 1, and E(x y T x 1 ae nown. To detemine the impotance function, to be used in conjunction with the paticle filteing algoithm, we poceed as follows. We appoximate the joint distibution of (x,y, conditional on x 1, by a Gaussian distibution with matching moments up to the second ode. Let the coesponding mean μ ( and the covaiance ( be given by μ ( ( ( μ 1 μ ( ( ( and ( 11 ( 1 ( ( 1 T (. (3
3 Stat Comput (009 19: Note that μ ( and ( can be calculated fom the moments which ae assumed to be nown fom Assumption A. Subsequently, we tae the impotance function to be the conditional distibution of x,given(x 1,y, deived fom the appoximated Gaussian distibution above. Thus, we assume π(x x (i 0: 1,y 0: N (m,, whee m μ ( 1 + ( 1 ( 11 ( 1 [ (] 1 ( y μ (, (4 [ ( ] 1 ( ( 1 T. (5 We note hee that a sufficient condition fo Assumption A to hold is: Condition 1 Both E(y x and E(y y T x ae polynomials in x of degee at most m and all conditional moments of x, given x 1, ae nown up to ode m. In the special case when the noise pocesses in (1 (ae additive Gaussian, suppose the system dynamics ae given by x f(x 1 + w, w N (0,Q, (6 y h(x + v, v N (0,R, 1,,... (7 then the conditional moments of all ode of x given x 1 ae nown. If, in addition, h( in (7 is a polynomial, then Condition 1 would be satisfied and hence Assumption A would hold. The pecise fomulas fo the quantities can be found in Appendix A. 4 Othe Gaussian impotance functions Thee ae othe Gaussian impotance functions poposed in the liteatue. We mention thee of them hee. The fist two ae based on the Gaussian appoximation of the optimal impotance function p(x x 1,y which ae ideologically simila to EMM, but the moments ae appoximated in diffeent ways. The thid one, on the othe hand, uses a ban of unscented Kalman filtes to obtain the poposal density. 4.1 Impotance function by lineaization (LIN In an ealie pape Doucet et al. (000 conside a dynamical system with additive Gaussian noise, given by (6 (7. Obseving that the optimal impotance function p(x x 1,y is Gaussian when h( in the obsevation model (7 is linea, the authos lineaize the obsevation equation (7 to obtain y h(f (x 1 + C (x f(x 1 + v (8 whee C h x (f (x 1. Subsequently, they use the coesponding Gaussian distibution as impotance function. This essentially educes to appoximating the conditional distibution of (x,y,givenx 1, by the Gaussian distibution with mean vecto μ and covaiance matix given by μ ( f(x 1 h(f (x 1 and ( Q QC T C Q C QC T + R. 4. Numeically appoximated moment matching In a moe ecent aticle, Guo et al. (005 conside the same dynamical system given by (6 (7. The impotance function poposed by them is also in effect deived fom a Gaussian appoximation of the joint distibution of (x,y conditional on x 1. Guo et al. (005, howeve, appoximate the moments in (3 by vaious numeical techniques, such as the Gauss-Hemite quadatue (GHQ ule and the Julie-Uhlmann quadatue (JUQ ule. We efe the eade to the oiginal aticle fo the details. 4.3 Unscented paticle filte (UPF Unscented paticle filte algoithm of van de Mewe et al. (000 is suitable fo geneal dynamical systems given by (1 (. In this method, fo each paticle a sepaate unscented Kalman filte is popagated to geneate the Gaussian poposal distibution. Once again we efe to the oiginal aticle fo details. 5 Implementation of the EMM Clealy, the EMM as descibed in Sect. 3 can be implemented if the Assumption A holds. A pope classification of models fo which Assumption A holds is not vey easy. Howeve, as mentioned in Sect. 3, if the dynamical system is given by (6 (7 with h( in (7 a polynomial function, then EMM can be implemented. When the exact values of the quantities in (3 cannot be calculated, we popose to appoximate the obsevation equation by one of polynomial fom and implement the EMM to deive the impotance function. Fo instance, conside a eal-valued dynamical system given by (6 (7. We assume futhe that the function h( is n times diffeentiable. We appoximate h( locally by its n-th degee Taylo polynomial aound x f(x 1 to get the following obsevation equation. y a m (x x m + v with a m 1 m! ( m h(x x m x x (9. (
4 06 Stat Comput (009 19: Then the quantities in (3 can be appoximated by the coesponding quantities fo the dynamical system govened by (6 and (. See Appendix B fo details. Note that the appoach poposed above extends the methodology used by Doucet et al. (000 whee the obsevation equation is appoximated by the fist degee Taylo polynomial, wheeas we conside highe degee polynomials. It is also wothwhile to note the diffeence between the appoach followed by Guo et al. (005 and the one poposed above. Guo et al. (005 wo with the given nonlinea model and duing setting up of the Gaussian impotance density, they appoximate the moments. We, on the othe hand, fist appoximate the obsevation equation with a n-th degee polynomial and futhe deive the Gaussian impotance density using the exact moments (based on the appoximated polynomial model. In the following section we pesent two illustative numeical examples. We notice that EMM pefoms bette than the othe methods in the sense that its computational load is consideably less, while the qualities of the filtes (in tems of RMSE ae compaable. 6 Numeical simulation esults In this section we conside two examples one with a polynomial obsevation model and the othe with a nonpolynomial model and compae the filteed estimates obtained by diffeent methods. In both examples we conside additive Gaussian noise pocesses. 6.1 Polynomial obsevation model As in Doucet et al. (000 we conside the system dynamics to be given by (6 (7 with f(x 1 x 1 + 5x x cos(1., (11 h(x x 0. (1 In ou simulations, we set Q, R 1 and geneate a time seies data of length 0 stating with x 0 N (0, 1. Given only the noisy obsevations y, paticle filte algoithm is pefomed with the impotance functions descibed in Sect. 4 (LIN, GHQ, JUQ, UPF and the new one (EMM poposed in Sect. 3. We estimate the state sequence x, 1,,...,0, with all the diffeent methods mentioned above. Note that, in this case, the diffeences in the moments used in EMM and LIN can be clealy seen. The moments used in EMM, as given in (3, ae ( f(x 1 μ ( f and (x 1 ( ( 0 + Q 0 Q f(x 1 Q f(x 1 Q f (x 1 Q 0 + Q 00 + R, while the moments used in LIN, as given in (9, ae ( f(x 1 μ and ( f (x 1 0 Q f(x 1 Q f(x 1 Q f (x 1 Q 0 + R. Fo GHQ, we use the five point quadatue ule and fo JUQ, the thee (n 1 sigma points wee calculated using κ. Fo UPF, the paametes ae taen to be the same as that of van de Mewe et al. (000 with α 1,β 0,κ and P 0 1. Howeve, while esampling we use systematic esampling scheme wheeas van de Mewe et al. (000use esidual esampling scheme. Fo all methods, the initial distibution p(x 0 is taen to be N (0, 5 and esampling was done when the effective sample size became less than one-thid of the oiginal sample size N. Fo each method, we fist calculate the oot mean squaed eo (RMSE ove M 0 uns fo each time point and then the aveage (ove time RMSE, given by ( M 1 Mj1 ( ˆx j xj 1.Heex j is the tue (simulated state fo time in the j-th un and ˆx j is the coesponding (point estimate using a PF method. Each of these methods is implemented with diffeent Monte Calo sample sizes N 0, 50, 500 and 00. In Table 1 the aveage RMSE s ae pesented. Also epoted ae the aveage (ove the 0 uns CPU time, in seconds, to complete a un and the aveage numbe of esampling steps (NRS out of the 0 time steps. Fist of all, we see fom the table that, as expected, the pefomances (as measued by RMSE of all the methods become simila as sample size N inceases. This is in confomity with the fact that fo any poposal distibution the paticle filte conveges to the tue posteio as N. The UPF seems to have a consideably highe computational load than the othe methods. This can be explained by the fact that one needs to un the unscented Kalman filte fo each paticle (at each step to calculate the poposal. Pefomances of GHQ, JUQ and EMM ae moe o less simila (which is bette than LIN, but the time taen to aive at the estimate is less in EMM than that by GHQ and JUQ. It appeas that the numbes of esampling steps ae almost the same fo GHQ, JUQ and EMM, which is slightly bette than LIN. So, the exta computational load fo GHQ and
5 Stat Comput (009 19: Table 1 Compaison of the pefomance of diffeent poposal distibutions with a polynomial obsevation model Table Compaison of the pefomance of diffeent poposal distibutions with a nonpolynomial obsevation model N LIN GHQ JUQ UPF EMM RMSE CPU NRS RMSE CPU NRS RMSE CPU NRS RMSE CPU NRS N LIN GHQ JUQ UPF EMM RMSE CPU NRS RMSE CPU NRS RMSE CPU NRS RMSE CPU NRS JUQ elative to ou EMM method can be constued as a esult of computing the moments numeically. Thus one can conclude that the EMM method is moe efficient compaed to the othe methods as it is computationally less demanding in aiving at the compaable level of efficiency. 6. Non-polynomial obsevation model Let us conside the following model x x 1 + 5x x cos(1. + w, w N (0,, (13 y tan 1 (x + v, v N (0, 1, 1,,... (14 Once again, a time seies data of length 0 was simulated stating with x 0 N (0, 1 and the diffeent paticle filtes wee applied on the obsevation y. Hee the exact moments given in (3 ae unnown. Fo EMM we have consideed a nd degee Taylo polynomial, as descibed in Sect. 5. The othe setup ae as in the pevious example in Sect The pefomances of diffeent methods ae pesented in Table. Again, compaing the RMSE s we obseve that the pefomances of GHQ, JUQ, UPF and EMM ae faily simila, and they ae all bette than LIN. But when CPU times ae compaed, UPF is the wost pefome. A vey close loo eveals that GHQ and JUQ may poduce slightly lowe RMSE compaed to EMM. Howeve, this elative gain is achieved at the expense of high computational load. Thus, consideing the tade of between the RMSE and the computational cost, EMM appeas to povide a pactical and efficient poposal density. 7 Conclusion In this aticle a new impotance function has been poposed which is based on the Gaussian appoximation of the conditional distibution of (x,y,givenx 1, with the fist two moments matched exactly to those of the tue conditional distibution. To use the poposed method one needs to now the moments of the system dynamics up to the second ode. A specific case in which this is satisfied is when the noise pocesses ae additive Gaussian and the obsevation equation is polynomial.when the exact moments ae not nown but the noise pocesses ae additive Gaussian and the obsevation model is smooth, we use a polynomial appoximation of the obsevation model to deive the impotance function. With the help of two examples it has been shown that the poposed EMM method povides a moe pactical and efficient poposal density consideing the tade off between the pefomance (RMSE and the computational load. Open Access This aticle is distibuted unde the tems of the Ceative Commons Attibution Noncommecial License which pemits any noncommecial use, distibution, and epoduction in any medium, povided the oiginal autho(s and souce ae cedited. Appendix A Conside the system dynamics as given in (6 (7. Then, the conditional moments of x given x 1 can be deived as follows. m ( m E(x m x 1 [f(x 1 ] E(w m 0 m ( m [f(x 1 ] μ m, (15 0
6 08 Stat Comput (009 19: whee μ j is the j-th (aw moment of the N (0,Q distibution, given by a m E[(x x m x 1 ] + x μ +1 0 and μ (!! Q, a m μ m+1 + x a m μ m, (1 fo 0, 1,,... (16 When h(x is a polynomial of the fom, h(x n 0 a x, E(y y T x 1 a m a l E[(x x m+l x 1 ] both E(y x and E(y y T x ae polynomials in x, thus l0 satisfying Condition 1. Subsequently, we have a m a l μ m+l, ( E(y x 1 a m E(x m x l0 1 whee μ j s ae given by (16. m ( m a m [f(x 1 ] μ m, 0 (17 Refeences E(x y T x 1 a m E(x m+1 x 1 m+1 ( m + 1 a m 0 [f(x 1 ] μ m+1, (18 E(y y T x 1 a m a l E(x m+l x 1 l0 m+l l0 0 ( m + l a m a l [f(x 1 ] μ m+l. (19 Then (15 (19 ensue the validity of Assumption A. Appendix B Conside the dynamical system descibed by (6 and (. Obseving that the conditional distibution of (x x given x 1 is N(0,Q, the following moments can be calculated as E(y x 1 E(x y T x 1 a m E[(x x m x 1 ] a m μ m, (0 a m E[x (x x m x 1 ] a m E[(x x m+1 x 1 ] Aashi, H., Kumamoto, H.: Constuction of discete-time nonlinea filte by Monte Calo methods with vaiance-educing techniques. Syst. Contol 19(4, 11 1 (1975. (In Japanese Andeson, B., Mooe, J.: Optimal Filteing. Pentice-Hall, Englewood Cliffs (1979 Aulampalam, S., Masell, S., Godon, N., Clapp, T.: A tutoial on paticle filtes fo online nonlinea/non-gaussian Bayesian tacing. IEEE Tans. Signal Pocess. 50(, (00 Bagchi, A.: Optimal Contol of Stochastic Systems. Seies in Systems and Contol Engineeing. Pentice-Hall Intenational, Englewood Cliffs (1993 Doucet, A., de Feitas, J., Godon, N. (eds.: Sequential Monte Calo Methods in Pactice. Spinge, New Yo (001 Doucet, A., Godsill, S., Andieu, C.: On sequential Monte Calo sampling methods fo Bayesian filteing. Stat. Comput., (000 Godon, N., Salmond, D.J., Smith, A.: Novel appoach to nonlinea/ non-gaussian Bayesian state estimation. IEE Poc. F, Rada Signal Pocess. 140(, (1993 Guo, D., Wang, X., Chen, R.: New sequential Monte Calo methods fo nonlinea dynamic systems. Stat. Comput. 15(, (005 Handschin, J.E., Mayne, D.Q.: Monte Calo techniques to estimate the conditional expectation in multistage nonlinea filteing. Int. J. Contol 9(5, (1969 Ito, K., Xiong, K.: Gaussian filtes fo nonlinea filteing poblems. IEEE Tans. Automat. Cont. 45(5, 9 97 (000 Jazwinsi, A.: Stochastic Pocesses and Filteing Theoy. Academic Pess, New Yo (1970 Julie, S., Uhlmann, J.: A new extensions of the Kalman filte to nonlinea systems. In: Poceedings of AeoSense: The 11th Intenational Symposium on Aeospace/Defense Sensing, Simulation and Contols, Multi Senso Fusion, Tacing and Resouce Management, vol. II, pp SPIE, Bellingham (1997 Kitagawa, G.: Non-Gaussian state-space modeling of nonstationay time seies. J. Am. Stat. Assoc. 8(400, 3 63 (1987 Pitt, M.K., Shephad, N.: Filteing via simulation: Auxiliay paticle filtes. J. Am. Stat. Assoc. 94(446, (1999 van de Mewe, R., de Feitas, N., Doucet, A., Wan, E.: The unscented paticle filte. Technical Repot CUED/F-INFENG/TR 380, Cambidge Univesity Engineeing Depatment, Cambidge (000 Wan, E., van de Mewe, R.: The unscented Kalman filte fo nonlinea estimation. In: Poceedings of the IEEE 000 Adaptive Systems fo Signal Pocessing, Communications and Contol Symposium (AS-SPCC, Lae Louise, Albeta, Canada. IEEE Pess, New Yo (000 West, M.: Mixtue models, Monte Calo, Bayesian updating and dynamic models. Comput. Sci. Stat. 4, (1993
Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationAn extended target tracking method with random finite set observations
4th Intenational Confeence on Infomation Fusion Chicago Illinois USA July 5-8 0 An extended taget tacing method with andom finite set obsevations Hongyan Zhu Chongzhao Han Chen Li Dept. of Electonic &
More informationCentral Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution
Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India
More informationGradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function
Intenational Confeence on Infomation echnology and Management Innovation (ICIMI 05) Gadient-based Neual Netwok fo Online Solution of Lyapunov Matix Equation with Li Activation unction Shiheng Wang, Shidong
More informationResearch Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function
Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationOn Polynomials Construction
Intenational Jounal of Mathematical Analysis Vol., 08, no. 6, 5-57 HIKARI Ltd, www.m-hikai.com https://doi.og/0.988/ima.08.843 On Polynomials Constuction E. O. Adeyefa Depatment of Mathematics, Fedeal
More informationMultiple Criteria Secretary Problem: A New Approach
J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and
More informationAPPLICATION OF MAC IN THE FREQUENCY DOMAIN
PPLICION OF MC IN HE FREQUENCY DOMIN D. Fotsch and D. J. Ewins Dynamics Section, Mechanical Engineeing Depatment Impeial College of Science, echnology and Medicine London SW7 2B, United Kingdom BSRC he
More informationHammerstein Model Identification Based On Instrumental Variable and Least Square Methods
Intenational Jounal of Emeging Tends & Technology in Compute Science (IJETTCS) Volume 2, Issue, Januay Febuay 23 ISSN 2278-6856 Hammestein Model Identification Based On Instumental Vaiable and Least Squae
More informationF-IF Logistic Growth Model, Abstract Version
F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth
More informationApplication of Parseval s Theorem on Evaluating Some Definite Integrals
Tukish Jounal of Analysis and Numbe Theoy, 4, Vol., No., -5 Available online at http://pubs.sciepub.com/tjant/// Science and Education Publishing DOI:.69/tjant--- Application of Paseval s Theoem on Evaluating
More informationMagnetometer Calibration Algorithm Based on Analytic Geometry Transform Yongjian Yang, Xiaolong Xiao1,Wu Liao
nd Intenational Foum on Electical Engineeing and Automation (IFEEA 5 Magnetomete Calibation Algoithm Based on Analytic Geomety ansfom Yongjian Yang, Xiaolong Xiao,u Liao College of Compute Science and
More informationAnalytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra
Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba
More informationDetermining solar characteristics using planetary data
Detemining sola chaacteistics using planetay data Intoduction The Sun is a G-type main sequence sta at the cente of the Sola System aound which the planets, including ou Eath, obit. In this investigation
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationLET a random variable x follows the two - parameter
INTERNATIONAL JOURNAL OF MATHEMATICS AND SCIENTIFIC COMPUTING ISSN: 2231-5330, VOL. 5, NO. 1, 2015 19 Shinkage Bayesian Appoach in Item - Failue Gamma Data In Pesence of Pio Point Guess Value Gyan Pakash
More informationA New Method of Estimation of Size-Biased Generalized Logarithmic Series Distribution
The Open Statistics and Pobability Jounal, 9,, - A New Method of Estimation of Size-Bied Genealized Logaithmic Seies Distibution Open Access Khushid Ahmad Mi * Depatment of Statistics, Govt Degee College
More informationModeling and Calculation of Optical Amplification in One Dimensional Case of Laser Medium Using Finite Difference Time Domain Method
Jounal of Physics: Confeence Seies PAPER OPEN ACCESS Modeling and Calculation of Optical Amplification in One Dimensional Case of Lase Medium Using Finite Diffeence Time Domain Method To cite this aticle:
More informationSolving Some Definite Integrals Using Parseval s Theorem
Ameican Jounal of Numeical Analysis 4 Vol. No. 6-64 Available online at http://pubs.sciepub.com/ajna///5 Science and Education Publishing DOI:.69/ajna---5 Solving Some Definite Integals Using Paseval s
More informationA matrix method based on the Fibonacci polynomials to the generalized pantograph equations with functional arguments
A mati method based on the Fibonacci polynomials to the genealized pantogaph equations with functional aguments Ayşe Betül Koç*,a, Musa Çama b, Aydın Kunaz a * Coespondence: aysebetuloc @ selcu.edu.t a
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationIdentification of the degradation of railway ballast under a concrete sleeper
Identification of the degadation of ailway ballast unde a concete sleepe Qin Hu 1) and Heung Fai Lam ) 1), ) Depatment of Civil and Achitectual Engineeing, City Univesity of Hong Kong, Hong Kong SAR, China.
More information6 PROBABILITY GENERATING FUNCTIONS
6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to
More informationState tracking control for Takagi-Sugeno models
State tacing contol fo Taagi-Sugeno models Souad Bezzaoucha, Benoît Max,3,DidieMaquin,3 and José Ragot,3 Abstact This wo addesses the model efeence tacing contol poblem It aims to highlight the encouteed
More informationNumerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method
Numeical Invesion of the Abel Integal Equation using Homotopy Petubation Method Sunil Kuma and Om P Singh Depatment of Applied Mathematics Institute of Technology Banaas Hindu Univesity Vaanasi -15 India
More informationControl Chart Analysis of E k /M/1 Queueing Model
Intenational OPEN ACCESS Jounal Of Moden Engineeing Reseach (IJMER Contol Chat Analysis of E /M/1 Queueing Model T.Poongodi 1, D. (Ms. S. Muthulashmi 1, (Assistant Pofesso, Faculty of Engineeing, Pofesso,
More informationCSCE 478/878 Lecture 4: Experimental Design and Analysis. Stephen Scott. 3 Building a tree on the training set Introduction. Outline.
In Homewok, you ae (supposedly) Choosing a data set 2 Extacting a test set of size > 3 3 Building a tee on the taining set 4 Testing on the test set 5 Repoting the accuacy (Adapted fom Ethem Alpaydin and
More informationFUSE Fusion Utility Sequence Estimator
FUSE Fusion Utility Sequence Estimato Belu V. Dasaathy Dynetics, Inc. P. O. Box 5500 Huntsville, AL 3584-5500 belu.d@dynetics.com Sean D. Townsend Dynetics, Inc. P. O. Box 5500 Huntsville, AL 3584-5500
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationSTATE VARIANCE CONSTRAINED FUZZY CONTROL VIA OBSERVER-BASED FUZZY CONTROLLERS
Jounal of Maine Science and echnology, Vol. 4, No., pp. 49-57 (6) 49 SAE VARIANCE CONSRAINED FUZZY CONROL VIA OBSERVER-BASED FUZZY CONROLLERS Wen-Je Chang*, Yi-Lin Yeh**, and Yu-eh Meng*** Key wods: takagi-sugeno
More informationA Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation
Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas
More informationA Power Method for Computing Square Roots of Complex Matrices
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 13, 39345 1997 ARTICLE NO. AY975517 A Powe Method fo Computing Squae Roots of Complex Matices Mohammed A. Hasan Depatment of Electical Engineeing, Coloado
More informationGoodness-of-fit for composite hypotheses.
Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationAalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández
Aalbog Univesitet Load Estimation fom atual input Modal Analysis Aenlle, Manuel López; Bincke, Rune; Canteli, Alfonso Fenández Published in: Confeence Poceedings Publication date: 005 Document Vesion Publishe's
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationRecent Advances in Chemical Engineering, Biochemistry and Computational Chemistry
Themal Conductivity of Oganic Liquids: a New Equation DI NICOLA GIOVANNI*, CIARROCCHI ELEONORA, PIERANTOZZI ARIANO, STRYJEK ROAN 1 DIIS, Univesità Politecnica delle ache, 60131 Ancona, ITALY *coesponding
More informationA Multivariate Normal Law for Turing s Formulae
A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition
More informationPredictive Filtering for Nonlinear Systems
Pedictive Filteing fo Nonlinea Systems John L. Cassidis F. Landis Makley Abstact In this pape, a eal-time pedictive filte is deived fo nonlinea systems. he majo advantage of this new filte ove conventional
More informationA scaling-up methodology for co-rotating twin-screw extruders
A scaling-up methodology fo co-otating twin-scew extudes A. Gaspa-Cunha, J. A. Covas Institute fo Polymes and Composites/I3N, Univesity of Minho, Guimaães 4800-058, Potugal Abstact. Scaling-up of co-otating
More informationSpherical Solutions due to the Exterior Geometry of a Charged Weyl Black Hole
Spheical Solutions due to the Exteio Geomety of a Chaged Weyl Black Hole Fain Payandeh 1, Mohsen Fathi Novembe 7, 018 axiv:10.415v [g-qc] 10 Oct 01 1 Depatment of Physics, Payame Noo Univesity, PO BOX
More informationBifurcation Analysis for the Delay Logistic Equation with Two Delays
IOSR Jounal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 5 Ve. IV (Sep. - Oct. 05), PP 53-58 www.iosjounals.og Bifucation Analysis fo the Delay Logistic Equation with Two Delays
More informationDuality between Statical and Kinematical Engineering Systems
Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.
More informationNuclear size corrections to the energy levels of single-electron atoms
Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea
More informationTESTING THE VALIDITY OF THE EXPONENTIAL MODEL BASED ON TYPE II CENSORED DATA USING TRANSFORMED SAMPLE DATA
STATISTICA, anno LXXVI, n. 3, 2016 TESTING THE VALIDITY OF THE EXPONENTIAL MODEL BASED ON TYPE II CENSORED DATA USING TRANSFORMED SAMPLE DATA Hadi Alizadeh Noughabi 1 Depatment of Statistics, Univesity
More informationEstimation of the Correlation Coefficient for a Bivariate Normal Distribution with Missing Data
Kasetsat J. (Nat. Sci. 45 : 736-74 ( Estimation of the Coelation Coefficient fo a Bivaiate Nomal Distibution with Missing Data Juthaphon Sinsomboonthong* ABSTRACT This study poposes an estimato of the
More informationA DETAILED DESCRIPTION OF THE DISCREPANCY IN FORMULAS FOR THE STANDARD ERROR OF THE DIFFERENCE BETWEEN A RAW AND PARTIAL CORRELATION: A TYPOGRAPHICAL
Olkin and Finn Discepancy A DETAILED DESCRIPTION OF THE DISCREPANCY IN FORMULAS FOR THE STANDARD ERROR OF THE DIFFERENCE BETWEEN A RAW AND PARTIAL CORRELATION: A TYPOGRAPHICAL ERROR IN OLKIN AND FINN (995
More informationLocalization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matrix
Jounal of Sciences, Islamic Republic of Ian (): - () Univesity of Tehan, ISSN - http://sciencesutaci Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matix H Ahsani
More informationarxiv: v2 [physics.data-an] 15 Jul 2015
Limitation of the Least Squae Method in the Evaluation of Dimension of Factal Bownian Motions BINGQIANG QIAO,, SIMING LIU, OUDUN ZENG, XIANG LI, and BENZONG DAI Depatment of Physics, Yunnan Univesity,
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationInformation Retrieval Advanced IR models. Luca Bondi
Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the
More informationRotor Blade Performance Analysis with Blade Element Momentum Theory
Available online at www.sciencediect.com ScienceDiect Enegy Pocedia 5 (7 ) 3 9 The 8 th Intenational Confeence on Applied Enegy ICAE6 Roto Blade Pefomance Analysis with Blade Element Momentum Theoy Faisal
More informationThis is a very simple sampling mode, and this article propose an algorithm about how to recover x from y in this condition.
3d Intenational Confeence on Multimedia echnology(icm 03) A Simple Compessive Sampling Mode and the Recovey of Natue Images Based on Pixel Value Substitution Wenping Shao, Lin Ni Abstact: Compessive Sampling
More informationA NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM
Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496
More informationLight Time Delay and Apparent Position
Light Time Delay and ppaent Position nalytical Gaphics, Inc. www.agi.com info@agi.com 610.981.8000 800.220.4785 Contents Intoduction... 3 Computing Light Time Delay... 3 Tansmission fom to... 4 Reception
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More informationBayesian Analysis of Topp-Leone Distribution under Different Loss Functions and Different Priors
J. tat. Appl. Po. Lett. 3, No. 3, 9-8 (6) 9 http://dx.doi.og/.8576/jsapl/33 Bayesian Analysis of Topp-Leone Distibution unde Diffeent Loss Functions and Diffeent Pios Hummaa ultan * and. P. Ahmad Depatment
More informationAlternative Tests for the Poisson Distribution
Chiang Mai J Sci 015; 4() : 774-78 http://epgsciencecmuacth/ejounal/ Contibuted Pape Altenative Tests fo the Poisson Distibution Manad Khamkong*[a] and Pachitjianut Siipanich [b] [a] Depatment of Statistics,
More informationMathematisch-Naturwissenschaftliche Fakultät I Humboldt-Universität zu Berlin Institut für Physik Physikalisches Grundpraktikum.
Mathematisch-Natuwissenschaftliche Fakultät I Humboldt-Univesität zu Belin Institut fü Physik Physikalisches Gundpaktikum Vesuchspotokoll Polaisation duch Reflexion (O11) duchgefüht am 10.11.2009 mit Vesuchspatne
More informationMitscherlich s Law: Sum of two exponential Processes; Conclusions 2009, 1 st July
Mitschelich s Law: Sum of two exponential Pocesses; Conclusions 29, st July Hans Schneebege Institute of Statistics, Univesity of Elangen-Nünbeg, Gemany Summay It will be shown, that Mitschelich s fomula,
More informationConjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationA Comparison and Contrast of Some Methods for Sample Quartiles
A Compaison and Contast of Some Methods fo Sample Quatiles Anwa H. Joade and aja M. Latif King Fahd Univesity of Petoleum & Mineals ABSTACT A emainde epesentation of the sample size n = 4m ( =, 1, 2, 3)
More information1 Similarity Analysis
ME43A/538A/538B Axisymmetic Tubulent Jet 9 Novembe 28 Similaity Analysis. Intoduction Conside the sketch of an axisymmetic, tubulent jet in Figue. Assume that measuements of the downsteam aveage axial
More informationImplicit Constraint Enforcement for Rigid Body Dynamic Simulation
Implicit Constaint Enfocement fo Rigid Body Dynamic Simulation Min Hong 1, Samuel Welch, John app, and Min-Hyung Choi 3 1 Division of Compute Science and Engineeing, Soonchunhyang Univesity, 646 Eupnae-i
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationChapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)
Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation
More informationComputers and Mathematics with Applications
Computes and Mathematics with Applications 58 (009) 9 7 Contents lists available at ScienceDiect Computes and Mathematics with Applications jounal homepage: www.elsevie.com/locate/camwa Bi-citeia single
More informationThe Substring Search Problem
The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is
More informationOn the global uniform asymptotic stability of time-varying dynamical systems
Stud. Univ. Babeş-Bolyai Math. 59014), No. 1, 57 67 On the global unifom asymptotic stability of time-vaying dynamical systems Zaineb HajSalem, Mohamed Ali Hammami and Mohamed Mabouk Abstact. The objective
More informationOn the Poisson Approximation to the Negative Hypergeometric Distribution
BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution
More informationHypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events
Intenational Jounal of Contempoay Mathematical Sciences Vol. 12, 2017, no. 5, 243-253 HIKARI Ltd, www.m-hikai.com https://doi.og/10.12988/ijcms.2017.7728 Hypothesis Test and Confidence Inteval fo the Negative
More informationUsing Laplace Transform to Evaluate Improper Integrals Chii-Huei Yu
Available at https://edupediapublicationsog/jounals Volume 3 Issue 4 Febuay 216 Using Laplace Tansfom to Evaluate Impope Integals Chii-Huei Yu Depatment of Infomation Technology, Nan Jeon Univesity of
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationChaos and bifurcation of discontinuous dynamical systems with piecewise constant arguments
Malaya Jounal of Matematik ()(22) 4 8 Chaos and bifucation of discontinuous dynamical systems with piecewise constant aguments A.M.A. El-Sayed, a, and S. M. Salman b a Faculty of Science, Aleandia Univesity,
More informationAppendix B The Relativistic Transformation of Forces
Appendix B The Relativistic Tansfomation of oces B. The ou-foce We intoduced the idea of foces in Chapte 3 whee we saw that the change in the fou-momentum pe unit time is given by the expession d d w x
More informationMathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden
Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty
More informationLiquid gas interface under hydrostatic pressure
Advances in Fluid Mechanics IX 5 Liquid gas inteface unde hydostatic pessue A. Gajewski Bialystok Univesity of Technology, Faculty of Civil Engineeing and Envionmental Engineeing, Depatment of Heat Engineeing,
More informationSafety variations in steel designed using Eurocode 3
JCSS Wokshop on eliability Based Code Calibation Safety vaiations in steel designed using Euocode 3 Mike Byfield Canfield Univesity Swindon, SN6 8LA, UK David Nethecot Impeial College London SW7 2BU, UK
More informationField emission of Electrons from Negatively Charged Cylindrical Particles with Nonlinear Screening in a Dusty Plasma
Reseach & Reviews: Jounal of Pue and Applied Physics Field emission of Electons fom Negatively Chaged Cylindical Paticles with Nonlinea Sceening in a Dusty Plasma Gyan Pakash* Amity School of Engineeing
More informationONE-POINT CODES USING PLACES OF HIGHER DEGREE
ONE-POINT CODES USING PLACES OF HIGHER DEGREE GRETCHEN L. MATTHEWS AND TODD W. MICHEL DEPARTMENT OF MATHEMATICAL SCIENCES CLEMSON UNIVERSITY CLEMSON, SC 29634-0975 U.S.A. E-MAIL: GMATTHE@CLEMSON.EDU, TMICHEL@CLEMSON.EDU
More informationMULTILAYER PERCEPTRONS
Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR
More informationANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE
THE p-adic VALUATION OF STIRLING NUMBERS ANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE Abstact. Let p > 2 be a pime. The p-adic valuation of Stiling numbes of the
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More informationPAPER 39 STOCHASTIC NETWORKS
MATHEMATICAL TRIPOS Pat III Tuesday, 2 June, 2015 1:30 pm to 4:30 pm PAPER 39 STOCHASTIC NETWORKS Attempt no moe than FOUR questions. Thee ae FIVE questions in total. The questions cay equal weight. STATIONERY
More informationModeling of High Temperature Superconducting Tapes, Arrays and AC Cables Using COMSOL
Except fom the Poceedings of the COMSOL Confeence 2010 Pais Modeling of High Tempeatue Supeconducting Tapes, Aays and AC Cables Using COMSOL Oleg Chevtchenko * Technical Univesity of Delft, The Nethelands
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More information2. The Munich chain ladder method
ntoduction ootstapping has become vey popula in stochastic claims eseving because of the simplicity and flexibility of the appoach One of the main easons fo this is the ease with which it can be implemented
More informationMONTE CARLO STUDY OF PARTICLE TRANSPORT PROBLEM IN AIR POLLUTION. R. J. Papancheva, T. V. Gurov, I. T. Dimov
Pliska Stud. Math. Bulga. 14 (23), 17 116 STUDIA MATHEMATICA BULGARICA MOTE CARLO STUDY OF PARTICLE TRASPORT PROBLEM I AIR POLLUTIO R. J. Papancheva, T. V. Guov, I. T. Dimov Abstact. The actual tanspot
More informationarxiv: v2 [astro-ph] 16 May 2008
New Anomalies in Cosmic Micowave Backgound Anisotopy: Violation of the Isotopic Gaussian Hypothesis in Low-l Modes Shi Chun, Su and M.-C., Chu Depatment of Physics and Institute of Theoetical Physics,
More informationarxiv: v1 [math.co] 1 Apr 2011
Weight enumeation of codes fom finite spaces Relinde Juius Octobe 23, 2018 axiv:1104.0172v1 [math.co] 1 Ap 2011 Abstact We study the genealized and extended weight enumeato of the - ay Simplex code and
More informationMEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION
MEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION Walte T. Fedee 337 Waen Hall, Biometics Unit Conell Univesity Ithaca, NY 4853 and Tey P. Speed Division of Mathematics & Statistics,
More informationChem 453/544 Fall /08/03. Exam #1 Solutions
Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law
More informationEncapsulation theory: the transformation equations of absolute information hiding.
1 Encapsulation theoy: the tansfomation equations of absolute infomation hiding. Edmund Kiwan * www.edmundkiwan.com Abstact This pape descibes how the potential coupling of a set vaies as the set is tansfomed,
More informationChapter 6 Balanced Incomplete Block Design (BIBD)
Chapte 6 Balanced Incomplete Bloc Design (BIBD) The designs lie CRD and RBD ae the complete bloc designs We now discuss the balanced incomplete bloc design (BIBD) and the patially balanced incomplete bloc
More informationSolution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere
Applied Physics eseach; Vol. 9, No. 6; 7 ISSN 96-969 E-ISSN 96-9647 Published by Canadian Cente of Science and Education Solution of a Spheically Symmetic Static Poblem of Geneal elativity fo an Elastic
More informationELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS
THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,
More informationPhysics 221 Lecture 41 Nonlinear Absorption and Refraction
Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.
More informationFREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS
Please cite this aticle as: Izabela Zamosa Fee tansvese vibations of non-unifom beams Scientific Reseach of the Institute of Mathematics and Compute Science Volume 9 Issue pages 3-9. The website: http://www.amcm.pcz.pl/
More information