Adjustment of Sampling Locations in Rail-Geometry Datasets: Using Dynamic Programming and Nonlinear Filtering
|
|
- Asher Long
- 5 years ago
- Views:
Transcription
1 Systems and Computers in Japan, Vol. 37, No. 1, 2006 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J87-D-II, No. 6, June 2004, pp Adjustment of Samplin Locations in Rail-Geometry Datasets: Usin Dynamic Prorammin and Nonlinear Filterin Masako Kamiyama 1,2 and Tomoyuki Hiuchi 3 1 Track Technoloy Division, Railway Technical Research Institute, Kokubunji, Japan 2 Department of Statistical Science, The Graduate University for Advanced Studies (Sokendai), Tokyo, Japan 3 Department of Statistical Modelin, The Institute of Statistical Mathematics, Research Oranization of Information and Systems, Tokyo, Japan SUMMARY A track inspection car, which measures the shape of railway tracks (hereafter, rail eometry) while it is runnin on rails, discretizes the measurement results at nearly fixed spatial intervals. However, the distance between the discretized locations (spatial samplin intervals) may shorten or lenthen locally due to slippin or slidin of the car wheel, and this prevents the samplin locations from alinin with those of a dataset obtained with another measurin run. The authors developed an alorithm for approximately alinin the samplin locations of the measurement datasets obtained with different runs. First, they considered this problem as the selection of the series of data correspondin to each supervised data from a trainin dataset, which was constructed by interpolation in order to minimize the evaluation function of a number sequence representin data points. Next, they used the maximum likelihood method to identify the unknown parameters contained in the evaluation function. This problem uses two features of the evaluation function. The first is that the evaluation function is minimized by dynamic prorammin, and the obtained optimum sequence is equivalent to a maximum a posteriori (MAP) estimate in the Bayesian framework. The second is that by convertin the evaluation function to a eneral state space representation, the lo likelihood of the model that includes the parameters is obtained by a nonlinear filterin method. Also, to simplify the search for the identification, they devised a parameter search procedure for the parameters in the autoreressive (AR) model Wiley Periodicals, Inc. Syst Comp Jpn, 37(1): 61 70, 2006; Published online in Wiley InterScience ( wiley.com). DOI /scj Key words: dynamic prorammin; state space representation; nonlinear filter; maximum likelihood method. 1. Introduction The special train called a track inspection car measures the shape of railway tracks (rail eometry) reularly since this eometry varies slihtly under the load of daily passin trains. The measurement results are discretized at equidistant intervals and stored for processin by computers. To obtain temporal variations of the rail eometry, the measurement results ideally should always be discretized at the same locations on the track. However, since the pulse sinals for specifyin the samplin locations are linked with the rotation of the wheel of the track inspection car as shown in Fi. 1, the pulse eneration locations cannot be reproduced. Therefore, the samplin locations on the track Wiley Periodicals, Inc.
2 Fi. 1. Scheme for spatial discretization of measured railway track eometry accordin to wheel rotation pulses. vary for each measurin run. To alin two sets of samplin locations obtained with difference runs, we must estimate the location aps from the discretized datasets, but this is difficult. If the spatial samplin intervals between all of the data were uniform, these aps between two sets of the samplin locations could be easily estimated by calculatin the cross correlation coefficient. However, because the samplin intervals locally shorten or lenthen as shown in Fi. 2 due to the slippin or slidin of the car wheel, this method cannot be applied. (In Fi. 2, the measurement results are discretized once per rotation of the wheel to simplify the explanation.) Fiure 3 shows an example of the measurement results. The datasets (A) and (B) for the upper two raphs show the values of the track aue (spacin between the left and riht rails) for the same railway section measured 1 day Fi. 2. Scheme for variations in spatial samplin intervals accompanyin wheel rotation states. Fi. 3. Example of datasets obtained from two measurements for the same railway section and results of adjustin samplin locations. and 8 days later (0 means that the aue is equal to the reference value). Both datasets were discretized at equal intervals (approximately 30 cm; 8 measurements per wheel rotation). In addition, for both (A) and (B), the location where the first data was obtained is somewhere within one samplin interval in the car movement direction of a certain device that is fixed at a specific location on the track. Similarly, the data for the 3201st point of (A) and 3197th point of (B) were also measured within one samplin interval in the car movement direction from another location detection device. It is known empirically that the aue hardly varies over an 8-day period or so, and the waveforms for (A) and (B), which are shown in the first and second raphs from the top of Fi. 3, are actually similar. Also, the datasets close to the 1st data in (A) and (B) are obtained at similar locations. This is apparent from the difference between the two waveforms, which is shown in the third raph from the top of Fi. 3. Similarly, the fourth raph from the top of Fi. 3 shows that the datasets close to the end data in (A) and (B) are also obtained at similar locations. In other words, despite the fact that (A) and (B) are measurement values of the aue for the same railway section, the number of data points constitutin (B) is smaller than the number constitutin (A) by 4 points (correspondin to approximately 1 m). This is because the samplin interval shortened or lenthened locally due to slippin or slidin of the car wheel durin one or the other measurin runs. 62
3 We developed an alorithm that uses dynamic prorammin (DP [1]) and nonlinear filterin for a eneral state space representation to approximately alin the samplin locations for datasets of this kind. The fifth raph in Fi. 3 shows the results when the samplin locations for dataset (B) were alined with those of dataset (A) by usin the proposed alorithm. This paper describes the location adjustment alorithm in Sections 2 and 3, presents a discussion in Section 4, and presents conclusions in Section Adjustment of Samplin Locations by Usin Dynamic Prorammin and Nonlinear Filterin An overview of the adjustment alorithm that we developed is presented below: (1) Model a mechanism to be empirically considered to yield the nonuniform samplin. (2) Formulate this model in a nonlinear optimization problem that can be solved by dynamic prorammin. However, since this form contains unknown parameters, the solution cannot be uniquely obtained. (3) Represent the above nonlinear form with a eneral state space representation and identify the included unknown parameters throuh a nonlinear filterin alorithm based on the maximum likelihood method. (4) Substitute the estimated parameters to solve the optimization problem of step (2). (End) These steps are explained in detail below Modelin the problem We decided to consider this location adjustment problem as a problem for selectin the most suitable data points correspondin to individual supervised data {Y 1, Y 2,..., Y T 1, Y T } (where T is the number of the supervised data) from the trainin dataset that was interpolated so that the number of points was a multiple of α as shown in Fi. 4. Since the aue dataset is smoothed by an analo low-pass filter before the spatial discretization, the oriinal trainin dataset is discretized so that the oriinal spatial frequency component does not chane. The spatial frequency characteristic of this filter varies with time, since this characteristic chanes with the speed of the inspection car. The interpolated dataset hereafter is called the trainin dataset {X 1, X 2,..., X N 1, X N }. Note that N ( αt) is the number of trainin data after interpolation. To distinuish between the dataset as a set and the individual data, we hereafter denote {Y 1, Y 2,..., Y T 1, Y T } collectively as Y 1:T. The other variables are denoted in a similar manner. Fi. 4. Scheme for the basic idea for adjustin the discretized location aps. Let n t be a data point index of the trainin set selected to match the t-th supervised data Y t ; the value of the correspondin data is X nt. This problem becomes a problem of findin the optimal number sequence (this is denoted as n 1:T ) from amon the feasible number sequences n 1:T. Also, let Y t X nt be denoted by e t (e t ~ N(0, τ 2 )). Next, we discuss the characteristics of n 1:T. Althouh the accurate distribution of the samplin interval is unknown, the runnin distance due to slippin or slidin is empirically determined to be approximately 0.1% of the total distance or less. Therefore, we assume that n def t = n t n t 1 of n 1:T, which corresponds to the samplin interval, is α for almost every t. Also, for the rane of values of n t, this article employs 1 n t 2α 1. This is based on the empirical knowlede that was obtained from conventional runnin experiments. Add to this assumption, we decide that the probability distribution of n t would be constant for all values other than α. We also assume that the second-order difference 2 n def t = n t 2n t 1 + n t 2 of n 1:T is 0 for almost every t. This is because 2 n t = 0 is naturally zero when the wheel rotates reularly; n t = n t 1 = α. Similarly, even when the wheel slips or slides, the wheel in the course of slippin or slidin seems to continue rotatin at the same rate ( 2 n t = 0), or to recover normal rotation ( n t = α) Adjustin the location by usin dynamic prorammin From the discussion above, we model the characteristics of the number sequence n 1:T indicatin the optimal sequence of the data points as follows: n t = def n t n t 1 is α for almost every t. 2 n t = def n t 2n t 1 + n t 2 is 0 for almost every t. 63
4 e t ~ N(0, τ 2 ) (where e t = Y t X nt ). At this time, n 1:T is obtained by optimizin the followin evaluation function F(n t ): (1) Calculate (x t ) = min[(x t 1 ) + f (x t, x t 1 ) x t 1 S t 1 ] for each element x t S t (where min[ ] represents the minimum value of the elements within [ ]). (2) Increment t by 1 and return to step (1). where S def t = [S t S t 1 ] T. When (x t ) is rearded as a state that transitions as t increases, the value of (x t ) depends only on (x t 1 ), x t, and x t 1 alone. (End) Since n 1:T varies with the values of parameters µ 1 and µ 2, we should carefully determine the parameter values. where (1) 2.3. Dynamic prorammin and MAP estimation in the Bayesian framework To determine the values of µ 1 and µ 2, n 1:T, which was obtained in the previous section, should be interpreted statistically. Assume that e t = Y t X nt is a normal random variable with a zero mean and a standard deviation τ. By multiplyin (1) by 1/(2τ 2 ) and exponentiatin it, we obtain (2) (3) (4) The factors µ 1 0 and µ 2 0 are penalties that are added to the evaluation function when n t α and 2 n t 0 respectively. In addition, as mentioned above, the data points at both ends of the supervised and trainin datasets are always measured within a fixed distance from specific locations on thetrack.therefore,theboundaryconditionsforn 1:T are as follows: Assume that both the penalties n 1 and n T have uniform distributions. For obtainin the n 1:T that satisfies the conditions shown in (1) to (4), dynamic prorammin (DP) can be applied [1], since this optimization problem observes the principle of optimality. Substitute x def t = [n t n t 1 ] T (where T represents transpose) in Eq. (1) and let then (x t ) def t = Σ j=1 F(n j ) can be replaced by the followin recursive formulation: The former term of the riht side (first { }) can be interpreted as a certain multiple of the conditional density p(y 1:T n 1:T ) of the supervised dataset Y 1:T assumin normality when all n 1:T are iven. Similarly, the latter term (second { }) can also be interpreted as a certain multiple of the prior distribution p(n 1:T ) of n 1:T in the Bayesian framework. Therefore, Eq. (1) can be interpreted as a search for the n 1:T that maximizes the posterior distribution p(n 1:T Y 1:T ) p(y 1:T n 1:T )p(n 1:T ) [this is called the maximum a posteriori (MAP) estimate]. Therefore, µ 1, µ 2, and τ 2 are hyperparameters (in the Bayesian framework) that ive the prior distribution [2]. Optimal solutions that are obtained by dynamic prorammin, not just for this problem, can be interpreted as MAP estimates [3]. Hence, the hyperparameter values can be evaluated usin the lo likelihood LL(µ 1, µ 2, τ 2 ) def = lo p(y 1:T µ 1, µ 2, τ 2 ) of Y 1:T. In other words, the hyperparameters for maximizin LL (denoted by µ ~ 1, µ ~ 2, and τ ~2 ) are interpreted as the optimal hyperparameters within the maximum likelihood method. The method for computin LL is described in the next section. (5) 64
5 2.4. Estimatin hyperparameters by the maximum likelihood method Lo likelihood LL(µ 1, µ 2, τ 2 ) is obtained as the followin sum of conditional probability distributions: Therefore, to obtain the LL, we should obtain p(y t Y 1:t 1, µ 1, µ 2, τ 2 ) for all t {1, 2,..., T}. Next, to obtain p(y t Y 1:t 1, µ 1, µ 2, τ 2 ), we transform the statistical model (5) into a eneralized state space representation [4]. By lettin x t = def [n t n t 1 ] and y t = def Y t, we can obtain the statistical model as follows. [Observation model] where e t ~ N(0, τ 2 ). [System model] where since q(v t ) = 1. Fiure 5 shows sample distributions for q( ). If we assume that e t and v t are white noises, we can define the prediction and filterin operations as follows. The p(y t y 1:t 1 ) values (where t = 1, 2,..., T 1, T) are obtained as the by-product of the computation for p(n 1:T y 1:T ), and the lo likelihood LL is obtained as the sum of their loarithms [4]. Since Y t / y t = 1, p(y t y 1:t 1 ) is equal to p(y t Y 1:t 1 ). [Prediction] where v t ~ q( x t 1, µ 1, µ 2, τ 2 ) and the distribution of q( x t 1, µ 1, µ 2, τ 2 ) is as follows. When n t 1 n t 2 = α, then where since q(v t ) = 1. When n t 1 n t 2 α, then Fi. 5. Typical distribution q( ) for system noise v t when α = 5, µ 1 = 1, µ 2 = 2, and τ 2 = 1. 65
6 [Filterin] p(x 1 y 0:0 ) is a uniform distribution from 1 to (2α 1). Some filterin operations in a eneralized state space representation involve complicated numerical interation. However, since the system vector x t consists of discrete values, while y t are continuous real values in this model, p(y t y 1:t 1 ) can be computed arithmetically without interation. Note that we used a rid search to maximize LL(µ 1, µ 2, τ 2 ). We now substitute (µ ~ 1, µ ~ 2, and τ ~2 ) in Eq. (1) and denote the n 1:T obtained by usin dynamic prorammin to optimize that equation by n ~ 1:T hereafter. Fi. 6. Fixed-interval smoothin distributions p(n t y 1:T ) (closeup view of t 520) (characterizin the differences e ~ 1:T as a Gaussian white noise sequence) Results and discussion Table 1 shows the parameters that were estimated by the alorithm described above when the actual measurement values for (A) shown in Fi. 3 were taken as the supervised dataset and the values for (B) were taken as the trainin dataset (before adjustment). (Note that α = 5 is used in this paper.) Also, n ~ 1:T that was estimated by dynamic prorammin is shown as (Y 1:T X n ~ 1:T ) in the fifth raph and as n ~ t in the sixth raph from the top of Fi. 3. From the sixth raph, it is apparent that if it is assumed that the track inspection car did not slip or slide durin the measurement run for dataset (A), then it is assumed to have slid for a distance correspondin to four points (approximately 1 m) near t = 1900 durin the measurement run for dataset (B). Also, the track inspection car s runnin speed when measurin dataset (B), which is shown in the bottom raph in Fi. 3 (this is calculated back from the number of distance pulses enerated within a unit time and normally cannot be referred), reveals that the train s speed dropped unnaturally in the vicinity of where the train was assumed to have slid and suests that it is hihly likely that the wheel actually slid and the eneration of distance pulses was temporarily reduced. Next, Fi. 6 shows the results when p(x t y 1:T ) is obtained by fixed-interval smoothin as follows: and then the followin is obtained as the marinal smoothin distribution: It is apparent that the marinal smoothin distribution does not convere and the reliability of n ~ 1:T is low. The major cause is that e 1:T is assumed to be a Gaussian white noise sequence even thouh an analo low-pass filter has smoothed the actual measurement data. As shown in Fi. 7, the autocorrelation of e ~ 1:T = def Y t X n ~ t is hih. Therefore, to improve the model, we apply an autoreressive (AR) model to e 1:T in the next section. Table 1. Estimated parameters obtained usin the actual measurement values in Fi. 3 (characterizin the differences e ~ 1:T as a Gaussian white noise sequence) Fi. 7. Autocorrelation coefficients of a difference sequence e ~ 1:T (characterizin e ~ 1:T as a Gaussian white noise sequence). 66
7 3. Remodelin the Difference Sequence 3.1. Introduction of an autoreressive (AR) model Assume that the difference component e t = Y t X nt is colored noise that can be described as follows usin an autoreressive model: Fi. 8. Autocorrelation coefficients of a residual sequence w ~ 1:T (characterizin e ~ 1:T as a first-order autoreressive process). (w t ~ N(0, σ 2 ) and k is the order of the autoreressive model, which is appropriately determined). If we substitute the followin for the T t=1 (Y t X nt) 2 in Eq. (1), and then optimize it usin dynamic prorammin, n 1:T is obtained where x def t = [n t n t 1 n t k ] T. This model can be transformed into a eneralized state space representation by lettin y def t = Y t k i=1 a i Y t i and X nt 0 = 0, and the maximum likelihood of the parameters can be estimated in a similar manner as in the previous section. LL to be maximized to estimate the parameters is LL(µ 1, µ 2, a 1:k, σ 2 ). Hereafter, the estimated parameters are denoted by (µ ~ 1, µ ~ 2, a ~ 1:k, σ ~2 ). Table 2 shows the estimated parameters when the above alorithm was applied to datasets (A) and (B), which were shown in Fi. 3 when k = 1 was assumed. It is apparent that by introducin the autoreressive model, the lo likelihood LL increased sinificantly from 608 to 3294 and that the reliability of the estimates improved. Next, Fi. 8 shows the autocorrelation coefficients of the residual sequence w ~ def t = X n ~ t a ~ 1X n ~ t 1. This is closer to a white noise autocorrelation than in Fi. 7. However, by introducin the autoreressive model, the number of points at which n ~ t α increased from 4 to 42 (fiure is omitted). In other words, althouh introducin the autoreressive model increased the likelihood, an implausible n ~ 1:T based on conventional knowlede was obtained. Also, by introducin the autoreressive model, the time required for the (µ ~ 1, µ ~ 2, a ~ 1, σ ~2 ) search increased. If Table 2. Estimated parameters obtained usin the actual measurement values in Fi. 3 (characterizin the differences e ~ 1:T as a first-order autoreressive process) (G 1, G 2, G a, G σ ) denote the rid counts that were set for (µ 1, µ 2, a 1, σ 2 ), respectively, the required time is a G 1 G 2 G a G σ multiple of the time for calculatin the LL value once. In our environment (Pentium MHz), since it took approximately 100 seconds per LL calculation (when T = 3201), when all rid point counts were set to 5, introducin the autoreressive model increased the calculation time from approximately 3.5 hours to 17.4 hours. However, aue eometry datasets are massive dependin on the track lenth (hundreds of kilometers), and fast location adjustments are desirable Improvement of the estimation alorithm As we have seen in the previous section, it is apparent that the alorithm must be improved so that an n ~ 1:T that conforms to conventional knowlede is obtained and parameters havin hih likelihood can be estimated. Therefore, we used the fact that (a ~ 1:k,σ ~2 ) are autoreressive parameters representin the difference sequence e ~ t to develop the followin alorithm. (0) (1) Set the initial value a 1:k of a 1:k to 0. (2) Maximize the lo likelihood LL(µ 1, µ 2, σ 2 a 1:k = a (0) 1:k ) accordin to a rid search to obtain provisional values of µ ~ 1 and µ ~ (0) 2 (which are described as µ 1 and (0) µ 2, respectively). In addition, use (0) (µ1, µ (0) 2, a (0) 1:k ) to estimate n ~ 1:T accordin to dynamic prorammin. (3) Use the Yule Walker method to estimate a ~ 1:k and σ ~2 from the estimated difference sequence e ~ 1:k (these are (1) described as a 1:k and σ 2(1) ) [5]. (4) Substitute (a (1) 1:k, σ 2(1) ) and maximize the lo likelihood LL aain to obtain (µ (1) 1, µ (1) 2 ), and estimate n ~ 1:T accordin to dynamic prorammin. (5) Compare the newly obtained n ~ 1:T with the one obtained previously and end the calculation if all entries are equal. If any entries differ, return to step (3). (End of the alorithm) 67
8 Table 3. Transitions of parameter estimates usin the actual measurements shown in Fi. 3 (characterizin differences e ~ 1:T as a first-order autoreressive process and applyin the alternate estimate alorithm) Maximum likelihood estimates of µ 1 (0) and µ2 (0) Estimates of σ (1) and a 1 (1) by dynamic prorammin Maximum likelihood estimates of µ 1 (1) and µ 2 (1) Estimates of σ (2) and a 1 (2) by dynamic prorammin µ ~ µ ~ LL a ~ 1 0 (fixed) σ ~ Remark Estimation ended Hereafter, this alorithm is termed the alternate estimate alorithm. Note that the (µ ~ 1, µ ~ 2, a ~ 1:k, σ ~2 ) that are obtained by the alternate estimate alorithm are not maximum likelihood estimates. (µ ~ 1, µ ~ 2) are maximum likelihood estimates when (a ~ 1:k, σ ~2 ) are already known. As a result of applyin the alternate estimate alorithm to the datasets (A) and (B) shown in Fi. 3, the same n ~ 1:T was obtained as before the autoreressive model was introduced. Table 3 shows transitions in estimates for the hyperparameters (µ ~ 1, µ ~ 2, a ~ 1:k, σ ~2 ). Althouh the lo likelihood was reduced to 3247 from the 3294 of Table 2, a hiher level is maintained than the 608 of Table 1. Additionally, the second raph from the bottom of Fi. 3 shows the obtained residual sequence w ~ def t = X n ~ t a ~ 1X n ~ t 1. The autocorrelation coefficients of this residual sequence are nearly the same as those in Fi. 8. The marinal smoothin distribution of this shown in Fi. 9 also converes better than the distribution in Fi. 6. Fi. 9. Fixed-interval smoothin distribution p(n t y 1:T ) (closeup view of t 520, modelin the differences e ~ 1:T as a first-order autoreressive process and applyin the alternate estimate alorithm). Now, we discuss the required time for the rid search. If the alternate estimate alorithm is repeated M times, the sum total of the required times is 100[G 1 G 2 G α + (M 1)G 1 G 2 ]. Since M = 2 for the calculations in Table 3, when all rid point counts were set to 5, the calculation time was reduced to approximately 4.2 hours from 17.4 hours. Note that M chaned from 2 to approximately 5. Therefore, the alternate estimate alorithm seems practical for obtainin n ~ 1:T from the taret data. 4. Discussion To verify the alternate estimate alorithm, we performed the followin simulation. (1) Establish fictional aue eometry and slippin or slidin to create data sequences for the supervised and trainin datasets (before the interpolation). Let n 1:T denote true values of the data points. (2) Independently create two white noise sequences as realization from the same normal distribution to simulate measurement noise and add them respectively to the two data sequences that were created in step (1). (3) Create Z 1:T, the sum of the supervised aue eometry and the noise. Then sequentially calculate Y 1 = Z 1, Y t = Z t + a * Y t 1 (2 t T) to create the supervised dataset Y 1:T. Also create the trainin dataset in a similar manner. (These sequences are shown in the top two raphs in Fi. 10.) Note that in this paper, a * = 0.8. (4) Estimate parameters (µ ~ 1, µ ~ 2, a ~ 1, σ ~2 ) and use those parameters to estimate n ~ 1:T. (End) The two results, obtained by the maximum likelihood estimation and by alternate estimate alorithm (see Table 4), are equal, and, therefore, the n ~ 1:T were also both the same. In addition, a ~ 1 and σ ~ 2 are both close to the true values. 68
9 We believe that this is because e 1:T has been modeled correctly. However, as Fi. 10 shows, n ~ 1:T n 1:T. The results of the parameter estimations usin actual measurement data differed for maximum likelihood estimation and the alternate estimate alorithm as shown in Tables 2 and 3. We believe that this is because e 1:T cannot be described accurately by a first-order autoreressive model for the followin reasons. That is, the main causes are that the measurement data contains irreular noise (for example, in the neihborhood of t = 2300) as is apparent from Fi. 3 and the filter that is actually applied is an analo low-pass filter that varies with time. 5. Conclusions Fi. 10. The result of adjustin samplin locations of the simulated datasets. The samplin locations of two measurement datasets of a track eometry obtained with different runs of a track inspection car can be approximately alined by dynamic prorammin if we model the wheel rotation and location detection pulse. Also, usin the fact that the optimal solution accordin to dynamic prorammin is a MAP estimate in the Bayesian framework enables unknown parameters included in the model to be estimated by usin the maximum likelihood method. When a low-pass filter is employed to represent the measurement datasets, its effect can be reduced if an autoreressive (AR) model is applied to the difference sequence. REFERENCES Table 4. Parameter estimation results for the simulation in Fi. 10 Maximum likelihood estimation Alternate estimate alorithm True value µ ~ µ ~ a ~ σ ~ LL Ibaraki T. Dou-teki Keikaku-hou, RisanSaitekikahou to Aruorizumu. Iwanami Shoten; p (in Japanese) 2. Akaike H. Likelihood and the Bayes procedure. In Bernardo JM, DeGroot MH, Lindley DV, Smith AFM (editors). Bayesian Statistics. University Press; p Godsill S, Doucet A, West M. Maximum a posteriori sequence estimation usin Monte Carlo particle filters. Ann Inst Statist Math 2001;53: Kitaawa G. Jikeiretsu-moderu no Yuudo-keisan to Parameta-suitei, FORTRAN77 Jikeiretsu-kaiseki Purouraminu. Iwanami Shoten; p (in Japanese) 5. Ishiuro M. Jiko-kaiki moderu no Atehame, Jikeiretsu-kaiseki no Houhou. Asakura Shoten; p (in Japanese) 69
10 AUTHORS (from left to riht) Masako Kamiyama completed her master s course in 1992 with a specialty in Resource Enineerin at Waseda University. Currently, she is employed in the Track Technoloy Division of the Railway Technical Research Institute, and she completed her doctoral course in 2004 in the Department of Statistical Science at the Graduate University for Advanced Studies (Sokendai). She is a member of the Japan Society for Industrial and Applied Mathematics and Japan Statistical Society. Tomoyuki Hiuchi (member) completed his doctoral course in the Department of Geophysics, in 1989 at the University of Tokyo. He holds a Ph.D. deree in Science. Currently, he is a professor in the Department of Statistical Modelin and Vice Director-General at the Institute of Statistical Mathematics, Research Oranization of Information and Systems. His specialty is statistical analysis, in particular, Bayesian modelin for knowlede discovery. He is a member of the Japan Statistical Society, the Society of Geomanetism and Earth, Planetary and Space Sciences, the American Geophysical Union, and the American Statistical Association. 70
PHY 133 Lab 1 - The Pendulum
3/20/2017 PHY 133 Lab 1 The Pendulum [Stony Brook Physics Laboratory Manuals] Stony Brook Physics Laboratory Manuals PHY 133 Lab 1 - The Pendulum The purpose of this lab is to measure the period of a simple
More informationNormed Distances and Their Applications in Optimal Circuit Design
Optimization and Enineerin, 4, 197 213, 2003 c 2003 Kluwer Academic Publishers. Manufactured in he Netherlands Normed Distances and heir Applications in Optimal Circuit Desin ABDEL-KARIM S.O. HASSAN Department
More informationEnergizing Math with Engineering Applications
Enerizin Math with Enineerin Applications Understandin the Math behind Launchin a Straw-Rocket throuh the use of Simulations. Activity created by Ira Rosenthal (rosenthi@palmbeachstate.edu) as part of
More informationIEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1
IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1 Intervention in Gene Reulatory Networks via a Stationary Mean-First-Passae-Time Control Policy Golnaz Vahedi, Student Member, IEEE, Babak Faryabi, Student
More informationEfficient method for obtaining parameters of stable pulse in grating compensated dispersion-managed communication systems
3 Conference on Information Sciences and Systems, The Johns Hopkins University, March 12 14, 3 Efficient method for obtainin parameters of stable pulse in ratin compensated dispersion-manaed communication
More informationStat260: Bayesian Modeling and Inference Lecture Date: March 10, 2010
Stat60: Bayesian Modelin and Inference Lecture Date: March 10, 010 Bayes Factors, -priors, and Model Selection for Reression Lecturer: Michael I. Jordan Scribe: Tamara Broderick The readin for this lecture
More informationMODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION
= = = MODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION Mauro Cettolo and Marcello Federico ITC-irst - Centro per la Ricerca Scientifica e Tecnoloica I-385 Povo, Trento, Italy ABSTRACT Robust acoustic
More informationA Multigrid-like Technique for Power Grid Analysis
A Multirid-like Technique for Power Grid Analysis Joseph N. Kozhaya, Sani R. Nassif, and Farid N. Najm 1 Abstract Modern sub-micron VLSI desins include hue power rids that are required to distribute lare
More informationTHE SIGNAL ESTIMATOR LIMIT SETTING METHOD
' THE SIGNAL ESTIMATOR LIMIT SETTING METHOD Shan Jin, Peter McNamara Department of Physics, University of Wisconsin Madison, Madison, WI 53706 Abstract A new method of backround subtraction is presented
More informationWire antenna model of the vertical grounding electrode
Boundary Elements and Other Mesh Reduction Methods XXXV 13 Wire antenna model of the vertical roundin electrode D. Poljak & S. Sesnic University of Split, FESB, Split, Croatia Abstract A straiht wire antenna
More informationEstimation of a Long-Term Variation of a Magnetic-Storm Index Using the Merging Particle Filter
1382 IEICE TRANS. INF. & SYST., VOL.E92 D, NO.7 JULY 2009 PAPER Special Section on Large Scale Algorithms for Learning and Optimization Estimation of a Long-Term Variation of a Magnetic-Storm Index Using
More informationLecture 5 Processing microarray data
Lecture 5 Processin microarray data (1)Transform the data into a scale suitable for analysis ()Remove the effects of systematic and obfuscatin sources of variation (3)Identify discrepant observations Preprocessin
More informationRESISTANCE STRAIN GAGES FILLAMENTS EFFECT
RESISTANCE STRAIN GAGES FILLAMENTS EFFECT Nashwan T. Younis, Younis@enr.ipfw.edu Department of Mechanical Enineerin, Indiana University-Purdue University Fort Wayne, USA Bonsu Kan, kan@enr.ipfw.edu Department
More informationV DD. M 1 M 2 V i2. V o2 R 1 R 2 C C
UNVERSTY OF CALFORNA Collee of Enineerin Department of Electrical Enineerin and Computer Sciences E. Alon Homework #3 Solutions EECS 40 P. Nuzzo Use the EECS40 90nm CMOS process in all home works and projects
More informationIT6502 DIGITAL SIGNAL PROCESSING Unit I - SIGNALS AND SYSTEMS Basic elements of DSP concepts of frequency in Analo and Diital Sinals samplin theorem Discrete time sinals, systems Analysis of discrete time
More informationLinearized optimal power flow
Linearized optimal power flow. Some introductory comments The advantae of the economic dispatch formulation to obtain minimum cost allocation of demand to the eneration units is that it is computationally
More informationRobust Semiparametric Optimal Testing Procedure for Multiple Normal Means
Veterinary Dianostic and Production Animal Medicine Publications Veterinary Dianostic and Production Animal Medicine 01 Robust Semiparametric Optimal Testin Procedure for Multiple Normal Means Pen Liu
More informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Imaes from Few Fourier Samples Gre Onie Department of Mathematics University of Iowa Iowa City, Iowa 54 Email: reory-onie@uiowa.edu Mathews Jacob Department of Electric and
More informationCausal Bayesian Networks
Causal Bayesian Networks () Ste7 (2) (3) Kss Fus3 Ste2 () Fiure : Simple Example While Bayesian networks should typically be viewed as acausal, it is possible to impose a causal interpretation on these
More informationSolution sheet 9. λ < g u = 0, λ = g α 0 : u = αλ. C(u, λ) = gλ gλ = 0
Advanced Finite Elements MA5337 - WS17/18 Solution sheet 9 On this final exercise sheet about variational inequalities we deal with nonlinear complementarity functions as a way to rewrite a variational
More informationTranslations on graphs with neighborhood preservation
1 Translations on raphs with neihborhood preservation Bastien Pasdeloup, Vincent Gripon, Nicolas Grelier, Jean-harles Vialatte, Dominique Pastor arxiv:1793859v1 [csdm] 12 Sep 217 Abstract In the field
More informationStochastic learning feedback hybrid automata for dynamic power management in embedded systems
Electrical and Computer Enineerin Faculty Publications Electrical & Computer Enineerin 2005 Stochastic learnin feedback hybrid automata for dynamic power manaement in embedded systems T. Erbes Eurecom
More informationStrong Interference and Spectrum Warfare
Stron Interference and Spectrum Warfare Otilia opescu and Christopher Rose WILAB Ruters University 73 Brett Rd., iscataway, J 8854-86 Email: {otilia,crose}@winlab.ruters.edu Dimitrie C. opescu Department
More information1 Introduction. Korea Atomic Energy Research Institute, South Korea.
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona, Spain, July 9-13, 2018 ICCFD10-0206 Multi-physics approach for nuclear reactor analysis usin thermal-hydraulics and
More informationDesign of Chevron Gusset Plates
017 SEAOC CONENTION PROCEEDINGS Desin of Chevron Gusset Plates Rafael Sali, Director of Seismic Desin Walter P Moore San Francisco, California Leih Arber, Senior Enineer American Institute of Steel Construction
More informationWhat Causes Image Intensity Changes?
Ede Detection l Why Detect Edes? u Information reduction l Replace imae by a cartoon in which objects and surface markins are outlined create line drawin description l These are the most informative parts
More informationLevenberg-Marquardt-based OBS Algorithm using Adaptive Pruning Interval for System Identification with Dynamic Neural Networks
Proceedins of the 29 IEEE International Conference on Systems, Man, and Cybernetics San Antonio, X, USA - October 29 evenber-marquardt-based OBS Alorithm usin Adaptive Prunin Interval for System Identification
More informationFLUID flow in a slightly inclined rectangular open
Numerical Solution of de St. Venant Equations with Controlled Global Boundaries Between Unsteady Subcritical States Aldrin P. Mendoza, Adrian Roy L. Valdez, and Carlene P. Arceo Abstract This paper aims
More informationTermination criteria in the Moore-Skelboe Algorithm for Global Optimization by Interval Arithmetic
To appear in series Nonconvex Global Optimization and Its Applications Kluwer Academic Publishers Termination criteria in the Moore-Skelboe Alorithm for Global Optimization by Interval Arithmetic M.H.
More informationRound-off Error Free Fixed-Point Design of Polynomial FIR Predictors
Round-off Error Free Fixed-Point Desin of Polynomial FIR Predictors Jarno. A. Tanskanen and Vassil S. Dimitrov Institute of Intellient Power Electronics Department of Electrical and Communications Enineerin
More informationModeling for control of a three degrees-of-freedom Magnetic. Levitation System
Modelin for control of a three derees-of-freedom Manetic evitation System Rafael Becerril-Arreola Dept. of Electrical and Computer En. University of Toronto Manfredi Maiore Dept. of Electrical and Computer
More informationA Mathematical Model for the Fire-extinguishing Rocket Flight in a Turbulent Atmosphere
A Mathematical Model for the Fire-extinuishin Rocket Fliht in a Turbulent Atmosphere CRISTINA MIHAILESCU Electromecanica Ploiesti SA Soseaua Ploiesti-Tiroviste, Km 8 ROMANIA crismihailescu@yahoo.com http://www.elmec.ro
More informationGraph Entropy, Network Coding and Guessing games
Graph Entropy, Network Codin and Guessin ames Søren Riis November 25, 2007 Abstract We introduce the (private) entropy of a directed raph (in a new network codin sense) as well as a number of related concepts.
More information6 Mole Concept. g mol. g mol. g mol ) + 1( g : mol ratios are the units of molar mass. It does not matter which unit is on the
What is a e? 6 Mole Concept The nature of chemistry is to chane one ecule into one or more new ecules in order to create new substances such as paints, fertilizers, food additives, medicines, etc. When
More informationANALYSIS OF POWER EFFICIENCY FOR FOUR-PHASE POSITIVE CHARGE PUMPS
ANALYSS OF POWER EFFCENCY FOR FOUR-PHASE POSTVE CHARGE PUMPS Chien-pin Hsu and Honchin Lin Department of Electrical Enineerin National Chun-Hsin University, Taichun, Taiwan e-mail:hclin@draon.nchu.edu.tw
More informationDesign and implementation of model predictive control for a three-tank hybrid system
International Journal of Computer cience and Electronics Enineerin (IJCEE) Volume 1, Issue 2 (2013) IN 2320 4028 (Online) Desin and implementation of model predictive control for a three-tank hybrid system
More informationBidirectional Clustering of Weights for Finding Succinct Multivariate Polynomials
IJCSNS International Journal of Computer Science and Network Security, VOL.8 No.5, May 28 85 Bidirectional Clusterin of Weihts for Findin Succinct Multivariate Polynomials Yusuke Tanahashi and Ryohei Nakano
More informationAutomatic Calibration Procedure for a Robotic Manipulator Force Observer
Automatic Calibration Procedure for a Robotic Manipulator Force Observer Gámez García, Javier; Robertsson, Anders; Gómez Ortea, Juan; Johansson, Rolf Published in: Robotics and Automation, 25. ICRA 25.
More informationHEAT TRANSFER IN EXHAUST SYSTEM OF A COLD START ENGINE AT LOW ENVIRONMENTAL TEMPERATURE
THERMAL SCIENCE: Vol. 14 (2010) Suppl. pp. S219 S232 219 HEAT TRANSFER IN EXHAUST SYSTEM OF A COLD START ENGINE AT LOW ENVIRONMENTAL TEMPERATURE by Snežana D. PETKOVIĆ a Radivoje B. PEŠIĆ b b and Jovanka
More informationTransactions of the VŠB Technical University of Ostrava, Mechanical Series No. 2, 2008, vol. LIV, article No. 1618
Transactions of the VŠB Technical University of Ostrava, Mechanical Series No., 008, vol. LIV, article No. 68 Dariusz JANECKI *, Lesze CEDRO ** SAMPLED SIGNALS DIFFERENTIATION USING A FREQUENCY CORRECTED
More information10log(1/MSE) log(1/MSE)
IROVED MATRI PENCIL METHODS Biao Lu, Don Wei 2, Brian L Evans, and Alan C Bovik Dept of Electrical and Computer Enineerin The University of Texas at Austin, Austin, T 7872-84 fblu,bevans,bovik@eceutexasedu
More informationTime Series Prediction by Kalman Smoother with Cross-Validated Noise Density
Time Series Prediction by Kalman Smoother with Cross-Validated Noise Density Simo Särkkä E-mail: simo.sarkka@hut.fi Aki Vehtari E-mail: aki.vehtari@hut.fi Jouko Lampinen E-mail: jouko.lampinen@hut.fi Abstract
More informationA Performance Comparison Study with Information Criteria for MaxEnt Distributions
A Performance Comparison Study with nformation Criteria for MaxEnt Distributions Ozer OZDEMR and Aslı KAYA Abstract n statistical modelin, the beinnin problem that has to be solved is the parameter estimation
More informationEquivalent rocking systems: Fundamental rocking parameters
Equivalent rockin systems: Fundamental rockin parameters M.J. DeJon University of Cambride, United Kindom E.G. Dimitrakopoulos The Hon Kon University of Science and Technoloy SUMMARY Early analytical investiations
More informationOPTIMAL BEAMFORMING AS A TIME DOMAIN EQUALIZATION PROBLEM WITH APPLICATION TO ROOM ACOUSTICS
OPTIMAL BEAMFORMING AS A TIME DOMAIN EQUALIZATION PROBLEM WITH APPLICATION TO ROOM ACOUSTICS Mark R P Thomas, Ivan J Tashev Microsoft Research Redmond, WA 985, USA {markth, ivantash}@microsoftcom Felicia
More informationGeneralized Distance Metric as a Robust Similarity Measure for Mobile Object Trajectories
Generalized Distance Metric as a Robust Similarity Measure for Mobile Object rajectories Garima Pathak, Sanjay Madria Department of Computer Science University Of Missouri-Rolla Missouri-6541, USA {madrias}@umr.edu
More informationContinuous Digital Simulation of the Second- Order Slowly Varying Wave Drift Force
Diital Commons@ Loyola Marymount University and Loyola Law School Systems Enineerin Faculty Works Systems Enineerin 9-1-1980 Continuous Diital Simulation of the Second- Order Slowly Varyin Wave Drift Force
More informationSimple Optimization (SOPT) for Nonlinear Constrained Optimization Problem
(ISSN 4-6) Journal of Science & Enineerin Education (ISSN 4-6) Vol.,, Pae-3-39, Year-7 Simple Optimization (SOPT) for Nonlinear Constrained Optimization Vivek Kumar Chouhan *, Joji Thomas **, S. S. Mahapatra
More informationActive filter synthesis based on nodal admittance matrix expansion
Tan et al. EURASIP Journal on Wireless Communications and Networkin (1) 1:9 DOI 1.118/s18-1-881-8 RESEARCH Active filter synthesis based on nodal admittance matrix expansion Linlin Tan, Yunpen Wan and
More informationThe Compact Muon Solenoid Experiment. Conference Report. Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland
Available on CMS information server CMS CR-28/5 he Compact Muon Solenoid Experiment Conference Report Mailin address: CMS CERN, CH-1211 GENEVA 23, Switzerland 14 February 28 Observin Heavy op Quarks with
More informationA Social Interactions Model with Endogenous Friendship Formation and Selectivity*
A Social Interactions Model with Endoenous Friendship Formation and Selectivity* Chih-Shen Hsieh and Lun-Fei Lee Department of Economics Ohio State University Columbus, Ohio, USA Date: Auust, 2011 Abstract
More informationEvaluation of the SONAR Meter in Wet Gas Flow for an Offshore Field Development
Evaluation of the SONAR Meter in Wet Gas Flow for an Offshore Field Development Anela Floyd, BP Siddesh Sridhar and Gabriel Dranea, Expro Meters 1 INTRODUCTION The ABC project is a hih pressure as condensate
More informationChapter 8 Applications of Newton s Second Law
81 Force Laws 2 Chapter 8 Applications of Newton s Second Law 811 Hooke s Law 2 822 Principle of Equivalence: 6 823 Gravitational Force near the Surface of the Earth 7 824 Electric Chare and Coulomb s
More information(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2
11 th Physics - Unit 2 Kinematics Solutions for the Textbook Problems One Marks 1. Which one of the followin Cartesian coordinate system is not followed in physics? 5. If a particle has neative velocity
More informationCode_Aster. Element CABLE_GAINE
Titre : Élément CABLE_GAINE Date : 28/07/2015 Pae : 1/12 Element CABLE_GAINE Summary: The element CABLE_GAINE presented in this document is to model cables of prestressin bein able to rub or slip into
More informationONLINE: MATHEMATICS EXTENSION 2 Topic 6 MECHANICS 6.3 HARMONIC MOTION
ONINE: MATHEMATICS EXTENSION Topic 6 MECHANICS 6.3 HARMONIC MOTION Vibrations or oscillations are motions that repeated more or less reularly in time. The topic is very broad and diverse and covers phenomena
More informationJoint multi-target detection and localization with a noncoherent statistical MIMO radar
Joint multi-taret detection and localization with a noncoherent statistical MIMO radar Yue Ai*, Wei Yi*, Mark R. Morelande and Linjian Kon* *School of Electronic Enineerin, University of Electronic Science
More informationAn Optimization Problem With a Surprisingly Simple Solution
An Optimization Problem With a Surprisinly Simple Solution D. Drinen, K. G. Kennedy, and W. M. Priestley 1 INTRODUCTION. Suppose you and n of your friends play the followin ame. A random number from the
More informationDavid R. Lovell, Tom Downs & Ah Chung Tsoi. given in [9].) Hildebrandt's method for adjusting selectivities is briey
An Evaluation of The Neoconitron David R Lovell, Tom Downs & Ah Chun Tsoi Abstract We describe a sequence of experiments investiatin the strenths and limitations of Fukushima's neoconitron as a handwritten
More informationGenetic Algorithm Approach to Nonlinear Blind Source Separation
Genetic Alorithm Approach to Nonlinear Blind Source Separation F. Roas, C.G. Puntonet, I.Roas, J. Ortea, A. Prieto Dept. of Computer Architecture and Technoloy, University of Granada. fernando@atc.ur.es,
More informationAcoustic Noise Reduction of Switched Reluctance Motor Drives
Middle-East Journal of Scientific Research 8 (6): 118-16, 11 ISSN 199-933 IDOSI Publications, 11 Acoustic Noise Reduction of Switched Reluctance Motor Drives M. Divandari and M.M. Kabir Department of Enineerin,
More informationOptimization of Mechanical Design Problems Using Improved Differential Evolution Algorithm
International Journal of Recent Trends in Enineerin Vol. No. 5 May 009 Optimization of Mechanical Desin Problems Usin Improved Differential Evolution Alorithm Millie Pant Radha Thanaraj and V. P. Sinh
More informationSynthesis-Gas Production from CH 4 -CO 2 -Ar via. Microwave Plasma Torch
Adv. Studies Theor. Phys., Vol. 6, 2012, no. 26, 1273-1287 Synthesis-Gas Production from CH 4 -CO 2 -Ar via Microwave Plasma Torch Neda H. Khiabani Laser and Plasma Research Institute, Shahid Beheshti
More informationAMERICAN INSTITUTES FOR RESEARCH
AMERICAN INSTITUTES FOR RESEARCH LINKING RASCH SCALES ACROSS GRADES IN CLUSTERED SAMPLES Jon Cohen, Mary Seburn, Tamas Antal, and Matthew Gushta American Institutes for Research May 23, 2005 000 THOMAS
More informationThe Bionic Lightweight Design of the Mid-rail Box Girder Based on the Bamboo Structure
Wenmin Chen 1, Weian Fu 1, Zeqian Zhan 1, Min Zhan 2 Research Institute of Mechanical Enineerin, Southwest Jiaoton University, Chendu, Sichuan, China (1), Department of Industrial and System Enineerin,
More informationNonlinear Model Reduction of Differential Algebraic Equation (DAE) Systems
Nonlinear Model Reduction of Differential Alebraic Equation DAE Systems Chuili Sun and Jueren Hahn Department of Chemical Enineerin eas A&M University Collee Station X 77843-3 hahn@tamu.edu repared for
More information2.2 Differentiation and Integration of Vector-Valued Functions
.. DIFFERENTIATION AND INTEGRATION OF VECTOR-VALUED FUNCTIONS133. Differentiation and Interation of Vector-Valued Functions Simply put, we differentiate and interate vector functions by differentiatin
More informationNing Wu Institute for Traffic Engineering Ruhr University Bochum, Germany Tel: ; Fax: ;
MODELLING THE IMPACT OF SIDE-STREET TRAFFIC VOLUME ON MAJOR- STREET GREEN TIME AT ISOLATED SEMI-ACTUATED INTERSECTIONS FOR SIGNAL COORDINATION DECISIONS Donmei Lin, Correspondin Author Center for Advanced
More informationSUSPENSION ANALYSIS EQUIPMENT BASED ON A HYDRAULIC SERVO-CYLINDER. P. M. Talaia, A.M. Rocha, and J. A. Ferreira
SUSPENSION ANALYSIS EQUIPMENT BASED ON A HYDRAULIC SERVO-CYLINDER P. M. Talaia, A.M. Rocha, and J. A. Ferreira Department of Mechanical Enineerin, P-38-93 Aveiro, Portual Abstract: This work describes
More informationSlip-Flow and Heat Transfer in Isoflux Rectangular Microchannels with Thermal Creep Effects
Journal of Applied Fluid Mechanics, Vol. 3, No. 2, pp. 33-4, 200. Available online at www.jafmonline.net, ISSN 735-3645. Slip-Flow and Heat Transfer in Isoflux Rectanular Microchannels with Thermal Creep
More informationLP Rounding and Combinatorial Algorithms for Minimizing Active and Busy Time
LP Roundin and Combinatorial Alorithms for Minimizin Active and Busy Time Jessica Chan, Samir Khuller, and Koyel Mukherjee University of Maryland, Collee Park {jschan,samir,koyelm}@cs.umd.edu Abstract.
More informationConical Pendulum: Part 2 A Detailed Theoretical and Computational Analysis of the Period, Tension and Centripetal Forces
European J of Physics Education Volume 8 Issue 1 1309-70 Dean onical Pendulum: Part A Detailed heoretical and omputational Analysis of the Period, ension and entripetal orces Kevin Dean Physics Department,
More informationPath Loss Prediction in Urban Environment Using Learning Machines and Dimensionality Reduction Techniques
Path Loss Prediction in Urban Environment Usin Learnin Machines and Dimensionality Reduction Techniques Mauro Piacentini Francesco Rinaldi Technical Report n. 11, 2009 Path Loss Prediction in Urban Environment
More informationMatrix multiplication: a group-theoretic approach
CSG399: Gems of Theoretical Computer Science. Lec. 21-23. Mar. 27-Apr. 3, 2009. Instructor: Emanuele Viola Scribe: Ravi Sundaram Matrix multiplication: a roup-theoretic approach Given two n n matrices
More informationThis is a repository copy of The effects of gas-phase and in-depth radiation absorption on ignition and steady burning rate of PMMA.
This is a repository copy of The effects of as-phase and in-depth radiation absorption on inition and steady burnin rate of PMMA. White ose esearch Online UL for this paper: http://eprints.whiterose.ac.uk/9515/
More informationInvestigation of ternary systems
Investiation of ternary systems Introduction The three component or ternary systems raise not only interestin theoretical issues, but also have reat practical sinificance, such as metallury, plastic industry
More informationAnalysis of visibility level in road lighting using image processing techniques
Scientific Research and Essays Vol. 5 (18), pp. 2779-2785, 18 September, 2010 Available online at http://www.academicjournals.or/sre ISSN 1992-2248 2010 Academic Journals Full enth Research Paper Analysis
More informationarxiv: v2 [math.oc] 11 Jun 2018
: Tiht Automated Converence Guarantees Adrien Taylor * Bryan Van Scoy * 2 Laurent Lessard * 2 3 arxiv:83.673v2 [math.oc] Jun 28 Abstract We present a novel way of eneratin Lyapunov functions for provin
More informationarxiv: v1 [cs.ai] 15 Nov 2013
Inferrin Multilateral Relations from Dynamic Pairwise Interactions arxiv:1311.3982v1 [cs.ai] 15 Nov 2013 Aaron Schein, Juston Moore, Hanna Wallach School of Computer Science University of Massachusetts
More informationGeneration of random waves in time-dependent extended mild-slope. equations using a source function method
Generation of random waves in time-dependent etended mild-slope equations usin a source function method Gunwoo Kim a, Chanhoon Lee b*, Kyun-Duck Suh c a School of Civil, Urban, and Geosystem Enineerin,
More informationStandard Test Method for Sulfur in Automotive Fuels by Polarization X-ray Fluorescence Spectrometry 1
Desination: D 7220 06 An American National Standard Standard Test Method for Sulfur in Automotive Fuels by Polarization X-ray Fluorescence Spectrometry 1 This standard is issued under the fixed desination
More informationImpact of Time Delays on Stability of Inertia-less Power Systems
P L Power Systems Laboratory Jialun Zhan Impact of Time Delays on Stability of Inertia-less Power Systems Semester Thesis PSL 186 EEH Power Systems Laboratory Swiss Federal Institute of Technoloy (ETH)
More informationARTICLE IN PRESS. Nuclear Instruments and Methods in Physics Research A
Nuclear Instruments and Methods in Physics Research A 66 (29) 517 522 Contents lists available at ScienceDirect Nuclear Instruments and Methods in Physics Research A journal homepae: www.elsevier.com/locate/nima
More informationExtended Target Poisson Multi-Bernoulli Filter
1 Extended Taret Poisson Multi-Bernoulli Filter Yuxuan Xia, Karl Granström, Member, IEEE, Lennart Svensson, Senior Member, IEEE, and Maryam Fatemi arxiv:1801.01353v1 [eess.sp] 4 Jan 2018 Abstract In this
More informationRESIDUAL 11-YEAR SIGNAL IN COEFFICIENTS OF MAIN FIELD MODELS. Several recently developed main geomagnetic field models, based on both observatory and
RESIDUAL -YEAR SIGNAL IN COEFFICIENTS OF MAIN FIELD MODELS CRISTIANA ŞTEFAN, VENERA DOBRICĂ, CRIŞAN DEMETRESCU Sabba S. Ştefănescu Institute of Geodynamics of te Romanian Academy, 9-, Jean-Louis Calderon
More informationTHE USE OF POWER FOR THE CHARACTERIZATION OF EARTHQUAKE GROUND MOTIONS
13 th World Conference on Earthquake Enineerin Vancouver, B.C., Canada Auust 1-6, 004 Paper No. 574 THE UE OF POWER FOR THE CHARACTERIZATION OF EARTHQUAKE GROUND MOTION Kenneth ELWOOD 1 and Erik NIIT UMMARY
More informationDesign of Impeller Blade by Varying Blades and Type of Blades Using Analytical
Desin of Impeller Blade by Varyin Blades and Type of Blades Usin Analytical Mr. C.V.S.Rajesh M.Tech Student, Department of Mechanical Enineerin, Avanthi Institute of Enineerin & Technoloy, Thaarapalasa,
More informationOpen boundary conditions for barotropic waves
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. C5, 3168, doi:10.109/00jc00159, 003 Open boundary conditions for barotropic waves J. Nycander and K. Döös Department of Meteoroloy, Stockholm University,
More informationMODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM. Apurba Sakti EGEE 520, Mathematical Modeling of EGEE systems Spring 2007
MODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM Apurba Sakti EGEE 520, Mathematical Modelin of EGEE systems Sprin 2007 Table of Contents Abstract Introduction Governin Equations Enery balance Momentum
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationDielectric characteristics of glass fibre reinforced plastics and their components
Plasticheskie Massy, No. 11, 2004, pp. 20 22 Dielectric characteristics of lass fibre reinforced plastics and their components V. I. Sokolov, S. I. Shalunov, I. G. Gurtovnik, L. G. Mikheeva, and I. D.
More informationComputer Vision Lecture 3
Demo Haribo Classification Computer Vision Lecture 3 Linear Filters 23..24 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Code available on the class website...
More informationSOFTWARE FOR THE GALE TRANSFORM OF FEWNOMIAL SYSTEMS AND A DESCARTES RULE FOR FEWNOMIALS
SOFTWARE FOR THE GALE TRANSFORM OF FEWNOMIAL SYSTEMS AND A DESCARTES RULE FOR FEWNOMIALS DANIEL J. BATES, JONATHAN D. HAUENSTEIN, MATTHEW E. NIEMERG, AND FRANK SOTTILE Abstract. We ive a Descartes -like
More informationEnergy-Efficient Resource Allocation for Ultra-reliable and Low-latency Communications
Enery-Efficient Resource Allocation for Ultra-reliable and Low-latency Communications Chenjian Sun, Chanyan She and Chenyan Yan School of Electronics and Information Enineerin, Beihan University, Beijin,
More informationThermal Conductivity of VIPs as a Function of Internal Pressure
Thermal Conductivity of VIPs as a Function of Internal Pressure The insulation properties of a VIP panel are determined by its effective thermal conductivity, as described in Eq.1. The lower the thermal
More informationAdaptive and predictive controllers applied to onshore wind energy conversion system
http://www.aimspress.com/journal/enery AIMS Enery, 6(4): 65 63. DOI:.3934/enery.8.4.65 Received: 3 May 8 Accepted: 3 Auust 8 Published: 8 Auust 8 Research article Adaptive and predictive controllers applied
More informationREPORT DOCUMENTATION PAGE
REPORT DOCMENTATION PAGE Form Approved OMB No. 74-88 Public reportin burden for this collection of information is estimated to averae hour per response, includin the time for reviewin instructions, searchin
More informationDo We Really Have to Consider Covariance Matrices for Image Feature Points?
Electronics and Communications in Japan, Part 3, Vol. 86, No. 1, 2003 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J85-A, No. 2, February 2002, pp. 231 239 Do We Really Have to Consider Covariance
More informationDecomposing compositional data: minimum chi-squared reduced-rank approximations on the simplex
Decomposin compositional data: minimum chi-squared reduced-rank approximations on the simplex Gert Jan Welte Department of Applied Earth Sciences Delft University of Technoloy PO Box 508 NL-600 GA Delft
More informationViscoData & ViscoShift
Introductory Theory Manual ViscoData & ViscoShift (Preliminary Version) Contents ) Preface 2) Introduction 3) Dynamic mechanical properties 4) sin ViscoShift and ViscoData 5) Concludin remarks 6) Refereces
More information