Stochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences
|
|
- Jeffery Carter
- 5 years ago
- Views:
Transcription
1 Applied athematics and Computation 6 (27) Stochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences. Rahman a, *,. Jackiewicz b, B.D. Welfert b a Department of athematics and Statistics, University of North Florida, 4567 St. Johns Bluff Rd. South, Jacksonville, FL 32224, USA b Department of athematics and Statistics, P.O. Box 527, Tempe, A 527-4, USA Abstract This paper extends the results of synaptically generated wave propagation through a network of connected excitatory neurons to a continuous model, defined by a Fredholm Volterra integro-differential equation (FVIDE), which includes memory effects of the past in the propagation. Stochastic approximation and numerical simulations are discussed. Ó 26 Elsevier Inc. All rights reserved. Keywords: Fredholm Volterra integro-differential equation; Euler Hermite method; Stochastic approximation; Error; Neural network. Introduction Noisy systems can be modelled in several different ways: For example, Langevin s equation describes a linear physical system to which white noise is added, and the linear theory for it has been extended to nonlinear stochastic differential equations with additive white noise, see []. This approach is based on additive white noise. Another approach was derived from work of Bogoliubov, see [9,] on averaging nonlinear oscillatory systems. In this system parameters are allowed to be random processes, and methods based on averaging and ergodic theory provide useful predictions from the model. This approach is referred to as being based on parametric noise. A third approach is to derive models, such as arkov Chains, for the systems state variable as being random variables, and then use the method of probability theory; for analysis, see []. Nonlinear Fredholm integro differential equation arises whenever the connections of excitable neurons are distributed according to the normal probability distribution and it is appropriate to take kernel as a Gaussian, or normal distribution of influence from any site. Gaussian functions, or normal distributions, arise naturally in neuroscience model, see [5] and [6]. It is known from the central limit theorem of probability theory that most random processes, when correctly scaled have (approximately) normal distributions. If we include some * Corresponding author. addresses: mrahman@unf.edu (. Rahman), jackiewi@math.asu.edu (. Jackiewicz), bdw@asu.edu (B.D. Welfert) /$ - see front matter Ó 26 Elsevier Inc. All rights reserved. doi:.6/j.amc
2 74. Rahman et al. / Applied athematics and Computation 6 (27) 73 2 memory effects of the past and then Fredholm integro-differential equation becomes Fredholm Volterra integro-differential equation. The purpose of this paper is to apply a random perturbation method to Fredholm Volterra integro differential equation that has a parametric noise, and demonstrate that for relatively small variance of noise the propagation sustains. We have shown this result using perturbation method in Fredholm Volterra integro differential equation described in [2]. 2. odel Consider the initial-value problem for the Fredholm Volterra integro differential equation of convolution type _h a ðs; tþ ¼ þ cos hðs; tþþfðtþ þl R t Kðt t ; yðt =ÞÞ R kðs s Þ h _ a ðs ; t Þds dt ðþ f a ðs; t; yðt=þ; h a ; aþ h a ðs; Þ ¼gðsÞ; 6 s 6,t P with a given smooth initial function g and _ h a ðs ; t Þ d dt h a ðs ; t Þ. n is vector of random processes that satisfies certain natural conditions to be listed below. The kernel K(t, n) is defined by Kðt; nþ ¼AðnÞ a e a t ; ð2þ where A(n) = n with a >. In applications n = y(t/) and is a small parameter. It measures the ratio of time scales between noise (fast) and the system response. The strength of the connection is defined by kðsþ ¼pffiffiffiffiffiffiffiffiffiffi 2pr 2 e s2 r 2 with r >. F(t) is the forcing function. < l is small, represents dimensionless parameter, reflecting the strength of connection between neurons. n takes values in a measurable space Y R N. We suppose that n satisfies the following conditions: It is a stationary arkov processes. It is ergodic process (in the sense of Birkhoff) having ergodic measure, say q, iny. It satisfies strong mixing conditions. k is a measurable function of Y with respect to q. When these conditions are satisfied we can obtain a useful approximations to the solution h a ðs; tþ ¼ p h a ðs; tþþ ffiffi ~ha ðs; tþþerror and d ~ h a ðs; tþ ¼f ha ðt; s; h a ; aþ ~ h a ðs; tþdt þ Bðs; tþdw; ð5þ p where derivation B(s,t) is shown in Section 4 and dw ¼ ffiffiffiffi dt Nð; Þ and ha ðs; tþ solves corresponding averaged system is _h a ðs; tþ ¼ þ cos hðs; tþþfðtþ þl R t Kðt t Þ R kðs s Þ _ ha ðs ; t Þds dt ð6þ f a ðs; t; h a ; aþ h a ðs; Þ ¼gðsÞ; 6 s 6,t P with initial function g and _ ha ðs ; t Þ d dt h a ðs ; t Þ. ð3þ ð4þ
3 . Rahman et al. / Applied athematics and Computation 6 (27) KðtÞ ¼ Kðt; yþqðdyþ q Kðt; y Þþq 2 Kðt; y 2 Þ: Y ð7þ The problem () is a continuous analog of a discrete Voltage Controlled Oscillator Neuron (VCON) model with perturbed kernel of transmission line in neural networks discussed in [5,6]. The angle h a (s,t) represents the phase at time t associated with a neuron located at the point s, and the integral appearing in () corresponds to the influence of the neighboring neuron along with some some memory effects of the past. It is shown by numerical experimentation that solution can be expanded using Central Limit Theorem for arkov processes (see [2,3]): In particular h a ¼ p h a þ ffiffi ~ha þ error; ðþ where d ~ h a ¼ f ha ðt; s; h a ; aþ ~ h a dt þ BðtÞdw; ð9þ p where, B(t) shown in Section 4 and dw ¼ ffiffiffiffi dt Nð; Þ. A solution strategy for () is based on the use of Gaussian rules and interpolation to a regular grid on bounded interval for some final time as we used in our earlier paper [4] and we also used random path described by 2-state transition probability matrix. In Section 3, we gave numerical method to solve the perturbed model. In Section 4, we gave experimental results. Finally, we investigate the asymptotic behavior of the model as!. 3. Simulating stochastic process Noise can arise in the network in a variety of ways: The cells have irregularities in their metabolism and in their membrane properties. The chemical bath in which they reside fluctuates. Input signal arriving from other cells or from outside the body are irregular, etc. We will use two state arkov process to generate this noise. Thus consider the case where y(t) = y(t/) takes values in a countable set Y ¼fy ; y 2 g; ðþ where we choose y =, y 2 =. The transition probability matrix P(,y i,y j ), where i,j =,2 can be represented by the 2 2 matrix P such that P i,j = P(,y i,y j ) denotes the probability y() =y j when starting with y() = y i. Let P ¼ P! P! ðþ P! P! ¼ a a b b ; ð2þ where < a,b < and << 2. The matrix P satisfies the following properties (see [5]). P is nonnegative (its entries are probabilities). P is stochastic, i.e., Pe = e (for sure one of the y j is the outcome of the transition from y i. P is irreducible, i.e., any state y j can be eventually be reached in a finite number of steps with nonzero probabilities. P is aperiodic (or acyclic, e.g., not a permutation). This implies that jk 2 j <.
4 76. Rahman et al. / Applied athematics and Computation 6 (27) 73 2 The above conditions guarantee the existence of a unit vector q > such that q T e ¼ ð3þ and lim k! P k = eq T. Thus q is is the unique positive left eigen vector associated to k = satisfying (3). Itis natural to consider the spectral decomposition (see [3]) and thus!! b a a a P ¼ þð a bþ ¼ P þð a bþp 2 b a b with P 2 i ¼ P i, P i P j =. Then P k ¼ðk P þ k 2 P 2 Þ k ¼ k k P k þ kk 2 P k 2 ¼ P þð a bþ k P 2! b a ¼ eq T : Thus q ¼ b and q 2 ¼ a and in state with probability 4. An Euler Hermite approach a b. This tells us such a system, in the long run, will be in state with probability, irrespective of the initial state in which the system started. a In this section we will describe numerical method for where integral term R is approximated by the Gauss Hermite quadrature rule for the kernel k defined by (). Then using () in () we have _h a ðs; tþ ¼ þ cos h a ðs; tþþfðtþ þl Kðt t ; yðt l =ÞÞ pffiffiffiffiffi 2p r h a ðs; Þ ¼gðsÞ; e ðs s Þ2 r 2 _ h a ðs ; t Þds dt ð4þ 6 s 6,t P. Using change of variables s s ¼ g, r s ¼ s rg, ds ¼ rdg, (4) takes the form, _h a ðs; tþ ¼ þ cos h a ðs; tþþfðtþ þ p l ffiffiffiffiffi Kðt t ; yðt =ÞÞ e g2 ha _ ðs rg; t Þdg dt 2p h a ðs; Þ ¼gðsÞ; ð5þ 6 s 6,t P. To approximate (5) we restrict the range of the variable s to the symmetric interval [ A,A], where A is large enough, and and introduce the grid s j = A +(j )Ds,j =,2,3,...,2N +. Then s = A,s N+ =,s 2N+ = A. Then (5) for s = s j, j =,...,2N +, and l ¼ p l ffiffiffi takes the form, 2p _h a ðs j ; tþ ¼ þ cos h a ðs j ; tþþfðtþ þl Kðt t ; yðt =ÞÞ e g2 ha _ ðs j rg; t Þdg dt ð6þ h a ðs j ; Þ ¼gðs j Þ; j =,...,2N +. We discretize by the Gauss Hermite quadrature rule e x2 dx ¼ Xn w k f ðx k Þþ n! pffiffiffi p 2 n ð2nþ! f ð2nþ ðnþ ð7þ k¼ 6 n 6. The rule (7) is convergent if the function f satisfies the inequality jf ðxþj 6 ex2 jxj þq
5 for some q >, compare [7]. oreover, if f has bounded derivatives of any order then the errors of (7) tend to zero as n!with spectral rate. Then,! t _h a ðs j ; tþ ¼ þ cos h a ðs j ; tþþfðtþþl Kðt t ; yðt =ÞÞ XNH w k ha _ ðs j rg k ; t Þ dt k¼ ðþ h a ðs j ; Þ ¼gðs j Þ; j =,2,...,2N +. Here w,...,w NH and g,g 2,...,g NH are the weights and abscissas of the Gauss Hermite quadrature rules, which can be computed using eigen values of the tridiagonal matrices (see [7]). Now introducing the grid t i ¼ði ÞDt; Dt ¼ T ; i ¼ ; 2;...; þ on the time interval [,T]. Then for t = t i, () takes the form! ti _h a ðs j ; t i Þ¼þcos h a ðs j ; t i ÞþFðtðiÞÞ þ l Kðt i t ; yðt =ÞÞ XNH w k ha _ ðs j rg k ; t Þ dt k¼ ð9þ h a ðs j ; t Þ¼gðs j Þ; j =,2,...,2N +. This implies _h a ðs j ; t i Þ¼þcos h a ðs j ; t i ÞþFðtðiÞÞ þ l Pi h a ðs j ; t Þ¼gðs j Þ; l¼ lþ t l Kðt i t ; yðt =ÞÞ XNH k¼! w k ha _ ðs j rg k ; t Þ dt where j =,2,...,2N +,i =,2,..., +. We will approximate h _ a ðs j rg k ; t Þ by h _ a ðs j rg k ; t l Þ, if t 2 [t l,t l+ ]! _h a ðs j ; t i Þ¼þcos h a ðs j ; t i ÞþFðtðiÞÞ þ l Pi tlþ Kðt i t ; yðt =ÞÞ XNH w k ha _ ðs j rg k ; t l Þ dt l¼ t l k¼ ð2þ h a ðs j ; t Þ¼gðs j Þ; j =,2,...,2N +,i =,2,..., +. Now introducing mini-grid t l n ¼ t l þðn Þ Dt ; n ¼ ; 2;...; þ ; ¼ Then tl ¼ t l; t l ¼ t þ lþ. Then (2) takes the form _h a ðs j ; t i Þ¼þcos h a ðs j ; t i ÞþFðtðiÞÞ þ l Pi P Dt Kðt i t l n ; PNH yðtl n =ÞÞ w k ha _ ðs j rg k ; t l Þ dt l¼ n¼ k¼ h a ðs j ; t Þ¼gðs j Þ; ð22þ j =,2,...,2N +,i =,2,..., +. Discretizing h _ a ðs j ; t i Þ and h _ a ðs j rg k ; t l Þ appearing in (22) takes the form h a ðs j ;t iþ Þ h a ðs j ;t i Þ P ¼ þ cos h Dt a ðs j ; t i ÞþFðtðiÞÞ þ l i P ðkðt i t l n ; yðtl n =ÞÞ h a ðs j ; t Þ¼gðs j Þ;. Rahman et al. / Applied athematics and Computation 6 (27) PNH k¼ l¼ n¼ w k ðh a ðs j rg k ; t lþ Þ h a ðs j rg k ; t l ÞÞÞdt j =,2,...,2N +,i =,2,..., +. Then for n =,2,..., +, t l n ¼ t l n ¼ t l þðn Þ Dt ¼ ðl ÞDt þðn Þ Dt ðl Þ þ n ADt fflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflffl} index ð2þ ð23þ ð24þ
6 7. Rahman et al. / Applied athematics and Computation 6 (27) 73 2 so that y tl n ¼ yððl Þ þ nþþ and thus (24) takes the form h a ðs j ; t iþ Þ¼h a ðs j ; t i ÞþDtðþcos h a ðs j ; t i ÞþFðtðiÞÞÞ P þ ldt i P Kðt i t l n ; yððl Þ þ nþþ PNH w k ðh a ðs j rg k ; t lþ Þ h a ðs j rg k ; t l Þ dt l¼ n¼ k¼ h a ðs j ; t Þ¼gðs j Þ; ð25þ j =,2,...,2N +,i =,2,..., +. Thus using (25) Euler Hermite method takes the form P h a ðs j ; t iþ Þ¼h a ðs j ; t i ÞþDt þ cos h a ðs j ; t i ÞþFðtðiÞ Þþl Dt i P Kðt i t l n ; yðn þ ðl ÞÞÞ l¼ n¼l PNH ð26þ w k ðh a ðs j rg k ;;t lþ Þ h a ðs j rg k ; t l ÞÞ k¼ h a ðs j ; t Þ¼gðsðjÞÞ; j ¼ ; 2;...; 2N þ ; i ¼ ; 2;...; : We have to approximate h a (s j rg k,t l+ ) and h a (s j rg k,t l ) using the values h a (s j,t l ) and h a (s j,t l+ ) defined at the grid points s j.if s j rg k 6 s then for any t we use h a ðs j rg k ; tþ ¼h a ðs j ; tþ: If sðþ < s j g k <6 s 2Nþ then we use piecewise linear interpolation h a ðs j rg k ; tþ ¼hðs j ; l þ Þþðh a ðs j ; l þ Þ h a ðs j ; lþþ s j rg k s j ; Ds where s j rg k 2 (x q,s q ]. This corresponds to the index q defined by s j rg q ¼ ceil k s þ : Dt If s j rg k P s 2N+ = A, we choose h a ðs j rg k ; tþ ¼h a ðs 2Nþ ; tþ: We also generate random sequence y using random number in ATLAB generated function rand for two state transition matrix (see Section 2) given by () with initial state y() =. Since t i = (i )Dt,i =,2,..., +, +2,..., +. Then t þ ¼ Dt ¼ T ¼ T. Thus for i =,2,..., *, if rand < a and y(i) = then y(i + ) = ; if rand 6 a and y(i) = then y(i + ) = ; if rand < b and y(i) = then y(i + ) = ; if rand 6 b and y(i) = then y(i +)=. 5. Stochastic approximation In this section we derived explicit expression for Stochastic Correction Equation defined by (9) to estimate the asymptotic behavior of h a ð p h a þ ffiffi ~ha Þ ¼ oðþ ð27þ as!, where h a and h a are the solution of () and (6) respectively and ~ h a is the solution of the Stochastic Correction Eq. (9), for more detail see [3], is given by B 2 ðtþ ¼ b aða þ bþ gðt; s; h a; y ; aþ 2 ½2 a bš; ð2þ
7 . Rahman et al. / Applied athematics and Computation 6 (27) where derivation of g(t,s,h a,y,a) is obtained using f ðt; s; h a ; y ; aþ, f ðt; s; h a ; y 2 ; aþ and f ðt; s; h a ; aþ, such that q g(t,s,h a,y,a)+q 2 g(t,s,h a,y 2,a) =, with q + q 2 =, where and gðt; s; h a ; y ; aþ f ðt; s; h a ; aþ f ðt; s; h a ; y ; aþ ð29þ gðt; s; h a ; y 2 ; aþ f ðt; s; h a ; aþ fðt; s; h a ; y 2 ; aþ: Thus using (7), (6) takes the form f ðt; s; h a ; aþ ¼ þ cos h a ðs; tþþfðtþ þ l ðq Kðt t ; y Þþq 2 Kðt t ; y 2 ÞÞ kðs s Þ hðs _ ; t Þds dt : ð3þ ð3þ Thus again using (2), we obtain f ðt; s; h a ; aþ ¼ þ cos h a ðs; tþþfðtþ þ l q y a e a ðt t Þ þ q 2 y 2 a e a ðt t Þ kðs s Þ _ ha ðs ; t Þds dt : Then using y = y and y = y 2,in() we have f ðt; s; h a ; y ; aþ ¼ þ cos t h a ðs; tþþfðtþþl Kðt t ; y Þ kðs s Þ _ ha ðs ; t Þds dt and ¼ þ cos h a ðs; tþþf ðtþ f ðt; s; h a ; y 2 ; aþ ¼ þ cos h a ðs; tþþf ðtþþl ¼ þ cos h a ðs; tþþf ðtþþl Thus using (33) and (3), (29) becomes, Kðt t ; y 2 Þ q y 2 a e a ðt t Þ gðt; s; h a ; y ; aþ ¼ þ cos h a ðs; tþþfðtþ þ l q y a e a ðt t Þ þ q 2 y 2 a e a ðt t Þ q 2 y 2 a e a ðt t Þ kðs s Þ _ ha ðs ; t Þds kðs s Þ _ ha ðs ; t Þds kðs s Þ _ ha ðs ; t Þds cos h a ðs; tþ FðtÞ l y a e a ðt t Þ kðs s Þ _ ha ðs ; t Þds ¼ l kðs s Þ _ ha ðs ; t Þds dt dt dt : dt dt ð32þ ð33þ ð34þ ð35þ and thus using (34), (3), and (3) becomes, gðt; s; h a ; y 2 ; aþ ¼ þ cos h a ðs; tþþfðtþ þ l q y a e a ðt t Þ þ q 2 y 2 a e a ðt t Þ ¼ l cos h a ðs; tþ F ðtþ l ðq 2 Þy 2 a e a ðt t Þ y 2 a e a ðt t Þ kðs s Þ _ ha ðs ; t Þds kðs s Þ _ ha ðs ; t Þds kðs s Þ _ ha ðs ; t Þds dt : dt dt ð36þ
8 . Rahman et al. / Applied athematics and Computation 6 (27) 73 2 Using (29), (2), takes the form B 2 bð2 a bþ t ðs; tþ ¼ l q aða þ bþ 2 y 2 a e a ðt t Þ kðs s Þ _ 2 ha ðs ; t Þds dt : ð37þ Thus sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi bð2 a bþ t Bðs; tþ ¼ l q aða þ bþ 2 y 2 a e a ðt t Þ kðs s Þ _ ha ðs ; t Þds and thus the Stochastic correction equation becomes d h a ¼½ sin h a ðs; tþš ~ h a dt þ Bðs; tþdw: dt ð3þ ð39þ 6. Numerical results We solved Eq. () numerically over time interval [, ] using. Euler Hermite with smaller time scale. 2. Evaluate p h a ðs; tþþ ffiffi ~ha ðs; tþ, where h a ðs; tþ and ~ hðs; tþ are computed using Eqs. (6) and (39) respectively. 3. Compare 2 with where h a ðs; tþ is the solution of (6) using large time scale, H and B 2 (s,t) is also computed using the solution of (3). Fig.. Contour plot of cosh a (s,t) for the Euler Hermite quadrature method with piecewise linear interpolation of FVIDE. Fig. 2. Contour plot of cosh a (s,t) for Euler Hermite quadrature method with piecewise linear interpolation of average FVIDE.
9 . Rahman et al. / Applied athematics and Computation 6 (27) 73 2 we using Bt w/ euler 2 Using E h a ðs; tþ ð p h a ðs; tþþ ffiffi ~ha ðs; tþþ, where E denotes the expected value, leads the concept of strong convergence. Thus numerical experimentations suggests that h a ðs; tþ ¼ p h a ðs; tþþ ffiffi ~ha ðs; tþ for! for perturbed Fredholm Volterra integro-differential equation. Thus the asymptotic result for Volterra integral equations given in, [2] can extended to more general FVIDE (see Figs. 3). 7. Concluding remarks In simulation we observed, the elapsed time for solving the averaged system approximately. times that for solving the perturbed system. While experimenting the error, we implicitly assumed that number of other sources of error are negligible, including the following, for more detail, (see []): Sampling error: error arising from approximating an expected value by sample mean. Random number bias: inherent error in the random generator. Rounding error: floating point roundoff errors. For a typical computation the sampling error is likely to be the most significant of these three. In preparing the programs in these simulations we found that some experimentation is required to make the number of samples sufficiently large and the time step is sufficiently small for the predicted order of convergence to be observable. Experimentation indicates that as step size decreases, the lack of independence in the samples from a random generator typically degrades the computation before rounding errors becomes significant. Future generalization of this work can be extended in the networks of excitatory and inhibitory neurons with sparse and random connectivity. References 3 2 Fig. 3. Strong error plots using Euler Hermite and dashed line is the line with slope. [] D.J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, Siam Rev. 43 (3) (22) [2] A.V. Skorokhod, F.C. Hoppensteadt, Random perturbations methods with applications in science and engineering, Appl. ath. Sci. 5 (997). [3] F.C. Hoppensteadt,. Rahman, B.D. Welfert, Central limit theorems for arkov processes with applications to circular processes, Adv. Appl. Prob., submitted for publication. [4]. Jackiewicz,. Rahman, B.D. Welfert, Numerical solution of Fredholm integro-differential equation modelling neural networks, Appl. Numer. ath. 56 (26) [5] F.C. Hoppensteadt, An Introduction to the athematical Neurons, second ed., Cambridge University Press, 997.
10 2. Rahman et al. / Applied athematics and Computation 6 (27) 73 2 [6] F.C. Hoppensteadt, _ h-networks, unpublished manuscript, Arizona State University, 23. [7] P.J. Davis, P. Rabinowitz, ethods of Numerical Integration, second ed., Academic Press, New York, 94. [] B. Oksendal, Stochastic differential equations, Springer-Verlag, Berlin, 2. [9] I.I. Gikhmann, Differential equations with random functions, AS transl. 2 (973). [] R.. Kahsminski, On random perturbations determined by differential equations with a small parameter, Theo. Prob. Appl. (966) [] S. Karlin, A First Course in Stochastic Processes, Academic press, 966.
Numerical solution of a Fredholm integro-differential equation modelling _ h-neural networks
Available online at www.sciencedirect.com Applied Mathematics and Computation 195 (2008) 523 536 www.elsevier.com/locate/amc Numerical solution of a Fredholm integro-differential equation modelling _ h-neural
More informationHaar wavelet method for nonlinear integro-differential equations
Applied Mathematics and Computation 176 (6) 34 333 www.elsevier.com/locate/amc Haar wavelet method for nonlinear integro-differential equations Ülo Lepik Institute of Applied Mathematics, University of
More informationRounding Error in Numerical Solution of Stochastic Differential Equations
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 21, No. 2, pp. 281 300, 2003 Rounding Error in Numerical Solution of Stochastic Differential Equations Armando Arciniega* and Edward Allen Department of Mathematics
More informationSynchronization of self-sustained oscillators by common white noise
Physica A 351 (5) 16 13 www.elsevier.com/locate/physa Synchronization of self-sustained oscillators by common white noise D.S. Goldobin a,b, A.S. Pikovsky a, a Department of Physics, University of Potsdam,
More informationDifferentiation matrices in polynomial bases
Math Sci () 5 DOI /s9---x ORIGINAL RESEARCH Differentiation matrices in polynomial bases A Amiraslani Received January 5 / Accepted April / Published online April The Author(s) This article is published
More informationExistence of positive periodic solutions for a periodic logistic equation
Applied Mathematics and Computation 139 (23) 311 321 www.elsevier.com/locate/amc Existence of positive periodic solutions for a periodic logistic equation Guihong Fan, Yongkun Li * Department of Mathematics,
More informationA numerical analysis of a model of growth tumor
Applied Mathematics and Computation 67 (2005) 345 354 www.elsevier.com/locate/amc A numerical analysis of a model of growth tumor Andrés Barrea a, *, Cristina Turner b a CIEM, Universidad Nacional de Córdoba,
More informationTwo numerical methods for solving a backward heat conduction problem
pplied Mathematics and Computation 79 (6) 37 377 wwwelseviercom/locate/amc Two numerical methods for solving a backward heat conduction problem Xiang-Tuan Xiong, Chu-Li Fu *, hi Qian epartment of Mathematics,
More informationTHE ALGORITHM TO CALCULATE THE PERIOD MATRIX OF THE CURVE x m þ y n ¼ 1
TSUKUA J MATH Vol 26 No (22), 5 37 THE ALGORITHM TO ALULATE THE PERIOD MATRIX OF THE URVE x m þ y n ¼ y Abstract We show how to take a canonical homology basis and a basis of the space of holomorphic -forms
More informationChapter 2 Linear Systems
Chapter 2 Linear Systems This chapter deals with linear systems of ordinary differential equations (ODEs, both homogeneous and nonhomogeneous equations Linear systems are etremely useful for analyzing
More informationInitial inverse problem in heat equation with Bessel operator
International Journal of Heat and Mass Transfer 45 (22) 2959 2965 www.elsevier.com/locate/ijhmt Initial inverse problem in heat equation with Bessel operator Khalid Masood *, Salim Messaoudi, F.D. Zaman
More informationJournal of Computational Physics
Journal of Computational Physics 8 (009) 8 93 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A mixed implicit explicit finite difference
More informationA New Elliptic Mean. Peter Kahlig. (Vorgelegt in der Sitzung der math.-nat. Klasse am 21. März 2002 durch das w. M. Ludwig Reich)
Sitzungsber. Abt. II (00) 11: 137 14 A New Elliptic Mean By Peter Kahlig (Vorgelegt in der Sitzung der math.-nat. Klasse am 1. März 00 durch das w. M. Ludwig Reich) Dedicated to Professor Janusz Matkowski
More informationNONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM. PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO.
NONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM BY PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO. 24 COWLES FOUNDATION FOR RESEARCH IN ECONOMICS YALE UNIVERSITY Box 2828
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More informationCOMPEL 30,2. Ramazan Rezaeyan Department of Mathematics, Sciences and Research Branch, Islamic Azad University, Tehran, Islamic Republic of Iran, and
The current issue and full text archive of this journal is available at www.emeraldinsight.com/332-1649.htm COMPEL 3,2 812 A stochastic perspective of RL electrical circuit using different noise terms
More informationLong-term behavior of polynomial chaos in stochastic flow simulations
Comput. Methods Appl. Mech. Engrg. 195 (26) 5582 5596 www.elsevier.com/locate/cma Long-term behavior of polynomial chaos in stochastic flow simulations Xiaoliang Wan, George Em Karniadakis * Division of
More informationShooting Methods for Numerical Solution of Stochastic Boundary-Value Problems
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 22, No. 5, pp. 1295 1314, 24 Shooting Methods for Numerical Solution of Stochastic Boundary-Value Problems Armando Arciniega and Edward Allen* Department of Mathematics
More informationA two level finite difference scheme for one dimensional PennesÕ bioheat equation
Applied Mathematics and Computation 171 (5) 3 331 wwwelseviercom/locate/amc A two level finite difference scheme for one dimensional PennesÕ bioheat equation Jennifer J Zhao a,1, Jun Zhang b, *,, Ning
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationChapter 2 Torsion Stresses in Thin-Walled Multi-Cell Box-Girders
Chapter Torsion Stresses in Thin-Walled Multi-Cell Box-Girders. Torsion of Uniform Thin-Walled Two-Cell Box-Girders The thin-walled box section with uniform thickness t as shown in Fig.., is subjected
More informationApplied Mathematics and Computation
Applied Mathematics and Computation 245 (2014) 86 107 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc An analysis of a new family
More informationLinearization techniques for singular initial-value problems of ordinary differential equations
Applied Mathematics and Computation 161 (25) 52542 www.elsevier.com/locate/amc Linearization techniques for singular initial-value problems of ordinary differential equations J.I. Ramos Room I-32-D, E.T.S.
More informationApproximate step response of a nonlinear hydraulic mount using a simplified linear model
Journal of Sound and Vibration 99 (007) 656 663 Short Communication JOURNAL OF SOUND AND VIBRATION Approximate step response of a nonlinear hydraulic mount using a simplified linear model Song He, Rajendra
More informationMathematics MPC1 (JUN15MPC101) General Certificate of Education Advanced Subsidiary Examination June Unit Pure Core TOTAL
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 1 Wednesday 13 May 2015 General Certificate of Education Advanced
More informationExperimental designs for precise parameter estimation for non-linear models
Minerals Engineering 17 (2004) 431 436 This article is also available online at: www.elsevier.com/locate/mineng Experimental designs for precise parameter estimation for non-linear models Z. Xiao a, *,
More informationFractional Ornstein Uhlenbeck processes. Joseph effect in models with infinite variance
Physica A 387 (28) 123 133 www.elsevier.com/locate/physa Fractional Ornstein Uhlenbeck processes. Joseph effect in models with infinite variance Marcin Magdziarz Institute of Mathematics and Computer Science,
More informationJournal of Computational Physics
Journal of Computational Physics 242 (213) 561 58 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Higher-order time integration
More informationPreservation of local dynamics when applying central difference methods: application to SIR model
Journal of Difference Equations and Applications, Vol., No. 4, April 2007, 40 Preservation of local dynamics when applying central difference methods application to SIR model LIH-ING W. ROEGER* and ROGER
More informationDynamics of four coupled phase-only oscillators
Communications in Nonlinear Science and Numerical Simulation 13 (2008) 501 507 www.elsevier.com/locate/cnsns Short communication Dynamics of four coupled phase-only oscillators Richard Rand *, Jeffrey
More informationQuantum damped harmonic oscillator on non-commuting plane
Physica A 385 (7) 4 36 www.elsevier.com/locate/physa Quantum damped harmonic oscillator on non-commuting plane Antony Streklas Department of Mathematics, University of Patras, 65 Patras, Greece Received
More informationSequences, Series, and the Binomial Formula
CHAPTER Sequences, Series, nd the Binomil Formul. SEQUENCES. ; ; ; ; 6 ; 6 6. ðþ ; ðþ ; ð Þ 6; ðþ ; ðþ 6; 6 ð6þ. ðþ ; ðþ : ðþ ; ðþ ; ðþ ; 6 ðþ 6 6 6. ; ; ; ; ; 6 6 6. 0 ; ; ; 8 ; 6 8 ; 6. 0; ; 6 ; ; 6
More informationBootstrap tests of multiple inequality restrictions on variance ratios
Economics Letters 91 (2006) 343 348 www.elsevier.com/locate/econbase Bootstrap tests of multiple inequality restrictions on variance ratios Jeff Fleming a, Chris Kirby b, *, Barbara Ostdiek a a Jones Graduate
More informationChebyshev pseudospectral method for computing numerical solution of convection diffusion equation
Available online at wwwsciencedirectcom Applied Mathematics and Computation (8) 537 546 wwwelseviercom/locate/amc Chebyshev pseudospectral method for computing numerical solution of convection diffusion
More informationArbitrary order Krylov deferred correction methods for differential algebraic equations
Journal of Computational Physics 221 (2007) 739 760 www.elsevier.com/locate/jcp Arbitrary order Krylov deferred correction methods for differential algebraic equations Jingfang Huang *,1, Jun Jia 1, Michael
More informationNewton-homotopy analysis method for nonlinear equations
Applied Mathematics and Computation 188 (2007) 1794 1800 www.elsevier.com/locate/amc Newton-homotopy analysis method for nonlinear equations S. Abbasbandy a, *, Y. Tan b, S.J. Liao b a Department of Mathematics,
More informationFrom graph to manifold Laplacian: The convergence rate
Appl. Comput. Harmon. Anal. 2 (2006) 28 34 www.elsevier.com/locate/acha Letter to the Editor From graph to manifold Laplacian: The convergence rate A. Singer Department of athematics, Yale University,
More informationApplied Mathematical Modelling
Applied Mathematical Modelling 33 (29) 2293 2297 Contents lists available at ScienceDirect Applied Mathematical Modelling journal homepage: www.elsevier.com/locate/apm Stability criteria for a nonlinear
More informationARTICLE IN PRESS. Journal of Sound and Vibration
Journal of Sound and Vibration 9 (1) 18 Contents lists available at ScienceDirect Journal of Sound and Vibration journal homepage: www.elsevier.com/locate/jsvi The decoupling of damped linear systems in
More informationAN INVESTIGATION OF STOKES SECOND PROBLEM FOR NON-NEWTONIAN FLUIDS
Numerical Heat Transfer, Part A, 47: 955 980, 2005 Copyright # Taylor & Francis Inc. ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407780590926390 AN INVESTIGATION OF STOKES SECOND PROBLEM FOR
More informationImage thresholding using Tsallis entropy
Pattern Recognition Letters 25 (2004) 1059 1065 www.elsevier.com/locate/patrec Image thresholding using Tsallis entropy M. Portes de Albuquerque a, *, I.A. Esquef b, A.R. Gesualdi Mello a, M. Portes de
More informationExponential and Logarithmic Functions
CHAPTER Eponential and Logarithmic Functions. EXPONENTIAL FUNCTIONS p. ffiffi ¼ :. ¼ :00. p pffiffi ffiffi ¼ :.. 0 0 f ðþ 0 0. 0 f ðþ 0 0 hðþ : : : :0 : 0 0 CHAPTER EXPONENTIAL AND LOGARITHMIC FUNCTIONS..
More informationStabilization and destabilization effects of the electric field on stochastic precipitate pattern formation
Chemical Physics 303 (2004) 151 155 www.elsevier.com/locate/chemphys Stabilization and destabilization effects of the electric field on stochastic precipitate pattern formation Istvan Lagzi a, Ferenc Izsak
More informationNew design of estimators using covariance information with uncertain observations in linear discrete-time systems
Applied Mathematics and Computation 135 (2003) 429 441 www.elsevier.com/locate/amc New design of estimators using covariance information with uncertain observations in linear discrete-time systems Seiichi
More informationA connection between time domain model order reduction and moment matching for LTI systems
Mathematical and Computer Modelling of Dynamical Systems Methods, Tools and Applications in Engineering and Related Sciences ISSN: 1387-3954 (Print) 1744-551 (Online) Journal homepage: http://www.tandfonline.com/loi/nmcm
More informationLocal spectral time splitting method for first- and second-order partial differential equations
Journal of Computational Physics 26 (25) 727 78 www.elsevier.com/locate/jcp Local spectral time splitting method for first- and second-order partial differential equations Sining Yu a, Shan Zhao a, G.W.
More informationA multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system
Chaos, Solitons and Fractals 30 (006) 197 03 www.elsevier.com/locate/chaos A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Qi Wang a,c, *,
More informationGeneralized diffusion equation
Physica A 368 (006) 55 6 www.elsevier.com/locate/hysa Generalized diffusion equation Jean Pierre Boon, James F. Lutsko Center for Nonlinear Phenomena and Comlex Systems, Université Libre de Bruxelles,
More informationSynchronization of an uncertain unified chaotic system via adaptive control
Chaos, Solitons and Fractals 14 (22) 643 647 www.elsevier.com/locate/chaos Synchronization of an uncertain unified chaotic system via adaptive control Shihua Chen a, Jinhu L u b, * a School of Mathematical
More informationMath 132, Spring 2003 Exam 2: Solutions
Math 132, Spring 2003 Exam 2: Solutions No calculators with a CAS are allowed. e sure your calculator is set for radians, not degrees if you do any calculus computations with trig functions. Part I, Multiple
More informationA new modified Halley method without second derivatives for nonlinear equation
Applied Mathematics and Computation 189 (2007) 1268 1273 www.elsevier.com/locate/amc A new modified Halley method without second derivatives for nonlinear equation Muhammad Aslam Noor *, Waseem Asghar
More informationSynchronization of two chaotic oscillators via a negative feedback mechanism
International Journal of Solids and Structures 40 (2003) 5175 5185 www.elsevier.com/locate/ijsolstr Synchronization of two chaotic oscillators via a negative feedback mechanism Andrzej Stefanski *, Tomasz
More informationParameter Estimation for Partially Complete Time and Type of Failure Data
Biometrical Journal 46 (004), 65 79 DOI 0.00/bimj.00004 arameter Estimation for artially Complete Time and Type of Failure Data Debasis Kundu Department of Mathematics, Indian Institute of Technology Kanpur,
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationASYMPTOTICS FOR STATIONARY VERY NEARLY UNIT ROOT PROCESSES
doi:1.1111/j.1467-9892.27.552.x ASYMPTOTICS FOR STATIONARY VERY NEARLY UNIT ROOT PROCESSES By Donald W. K. Andrews and Patrik Guggenberger Cowles Foundation, Yale University and Department of Economics,
More informationJournal of Computational Physics
Journal of Computational Physics 35 (3) 43 445 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics ournal homepage: www.elsevier.com/locate/cp Numerical methods and comparison
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More informationDeterminants and polynomial root structure
International Journal of Mathematical Education in Science and Technology, Vol 36, No 5, 2005, 469 481 Determinants and polynomial root structure L G DE PILLIS Department of Mathematics, Harvey Mudd College,
More informationMapping from a fragile glass-forming system to a simpler one near their glass transitions
Physica A 385 (7) 439 455 www.elsevier.com/locate/physa Mapping from a fragile glass-forming system to a simpler one near their glass transitions Michio Tokuyama a,, Takayuki Narumi b, Eri Kohira b a Institute
More information1 Basic tools. 1.1 The Fourier transform. Synopsis. Evenness and oddness. The Fourier transform of f(x) is defined as Z 1
1 Basic tools 1.1 The Fourier transform 1 Synopsis The Fourier transform of f(x) is defined as FðsÞ ¼ f ðxþe i2pxs dx The inverse Fourier transform is given by f ðxþ ¼ FðsÞe þi2pxs ds Evenness and oddness
More informationPossible numbers of ones in 0 1 matrices with a given rank
Linear and Multilinear Algebra, Vol, No, 00, Possible numbers of ones in 0 1 matrices with a given rank QI HU, YAQIN LI and XINGZHI ZHAN* Department of Mathematics, East China Normal University, Shanghai
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationGeneralised projections in finite state automata and decidability of state determinacy
International Journal of Control Vol. 81, No. 10, October 2008, 1626 1644 Generalised projections in finite state automata and decidability of state determinacy Ishanu Chattopadhyay and Asok Ray* Pennsylvania
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationYimin Wei a,b,,1, Xiezhang Li c,2, Fanbin Bu d, Fuzhen Zhang e. Abstract
Linear Algebra and its Applications 49 (006) 765 77 wwwelseviercom/locate/laa Relative perturbation bounds for the eigenvalues of diagonalizable and singular matrices Application of perturbation theory
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
More information1. INTRODUCTION 2. THE DATA
The Astrophysical Journal, 591:695 713, 003 July 10 # 003. The American Astronomical Society. All rights reserved. Printed in U.S.A. E THE CROSS-WAVELET TRANSFORM AND ANALYSIS OF QUASI-PERIODIC BEHAVIOR
More informationApproximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach
Approximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach Henrik T. Sykora, Walter V. Wedig, Daniel Bachrathy and Gabor Stepan Department of Applied Mechanics,
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationAnalysis of the Global Relation for the Nonlinear Schro dinger Equation on the Half-line
Letters in Mathematical Physics 65: 199 212, 23. 199 # 23 Kluwer Academic Publishers. Printed in the Netherlands. Analysis of the Global Relation for the Nonlinear Schro dinger Equation on the Half-line
More informationSubthreshold cross-correlations between cortical neurons: Areference model with static synapses
Neurocomputing 65 66 (25) 685 69 www.elsevier.com/locate/neucom Subthreshold cross-correlations between cortical neurons: Areference model with static synapses Ofer Melamed a,b, Gilad Silberberg b, Henry
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationHow do our representations change if we select another basis?
CHAPTER 6 Linear Mappings and Matrices 99 THEOREM 6.: For any linear operators F; G AðV Þ, 6. Change of Basis mðg FÞ ¼ mðgþmðfþ or ½G FŠ ¼ ½GŠ½FŠ (Here G F denotes the composition of the maps G and F.)
More informationLocalized Excitations in Networks of Spiking Neurons
Localized Excitations in Networks of Spiking Neurons Hecke Schrobsdorff Bernstein Center for Computational Neuroscience Göttingen Max Planck Institute for Dynamics and Self-Organization Seminar: Irreversible
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationGlobal attractivity in a rational delay difference equation with quadratic terms
Journal of Difference Equations and Applications ISSN: 1023-6198 (Print) 1563-5120 (Online) Journal homepage: http://www.tandfonline.com/loi/gdea20 Global attractivity in a rational delay difference equation
More informationOPTIMAL SIGNAL PROCESSING IN CAVITY RING-DOWN SPECTROSCOPY
CHAPTER 18 OPTIMAL SIGNAL PROCESSING IN CAVITY RING-DOWN SPECTROSCOPY Kevin K. Lehmann and Haifeng Huang Contents 1. Introduction 64. The Model 65 3. Generalized Least Squares Fit with Correlated Data
More informationRobust synchronization of Sprott circuits using sliding mode control
Chaos, Solitons and Fractals 3 (6) 11 18 www.elsevier.com/locate/chaos Robust synchronization of Sprott circuits using sliding mode control David I. Rosas Almeida *, Joaquín Alvarez, Juan Gonzalo Barajas
More informationChapter 2 Reliability and Confidence Levels of Fatigue Life
Chapter Reliability and Confidence Levels of Fatigue Life.1 Introduction As is well known, fatigue lives of nominally identical specimens subjected to the same nominal cyclic stress display scatter as
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationFuzzy age-dependent replacement policy and SPSA algorithm based-on fuzzy simulation
Available online at wwwsciencedirectcom Information Sciences 178 (2008) 573 583 wwwelseviercom/locate/ins Fuzzy age-dependent replacement policy and SPSA algorithm based-on fuzzy simulation Jiashun Zhang,
More informationSeries solutions of non-linear Riccati differential equations with fractional order
Available online at www.sciencedirect.com Chaos, Solitons and Fractals 40 (2009) 1 9 www.elsevier.com/locate/chaos Series solutions of non-linear Riccati differential equations with fractional order Jie
More informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
More informationNumber of Complete N-ary Subtrees on Galton-Watson Family Trees
Methodol Comput Appl Probab (2006) 8: 223 233 DOI: 10.1007/s11009-006-8549-6 Number of Complete N-ary Subtrees on Galton-Watson Family Trees George P. Yanev & Ljuben Mutafchiev Received: 5 May 2005 / Revised:
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx ) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: wwwelseviercom/locate/cnsns A note on the use of Adomian
More informationApplied Mathematics and Computation
Alied Mathematics and Comutation 18 (011) 944 948 Contents lists available at ScienceDirect Alied Mathematics and Comutation journal homeage: www.elsevier.com/locate/amc A generalized Birkhoff Young Chebyshev
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationEstimation of cusp in nonregular nonlinear regression models
Journal of Multivariate Analysis 88 (2004) 243 251 Estimation of cusp in nonregular nonlinear regression models B.L.S. Prakasa Rao Indian Statistical Institute, 7, SJS Sansanwal Marg, New Delhi, 110016,
More informationCharacteristics of forced convection heat transfer in vertical internally finned tube B
International Communications in Heat and Mass Transfer 32 (2005) 557 564 www.elsevier.com/locate/ichmt Characteristics of forced convection heat transfer in vertical internally finned tube B A. Al-Sarkhi*,
More informationTail probability of random variable and Laplace transform
Applicable Analysis Vol. 84, No. 5, May 5, 499 5 Tail probability of random variable and Laplace transform KENJI NAKAGAWA* Department of Electrical Engineering, Nagaoka University of Technology, Nagaoka,
More informationIrreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1
Irreducibility Irreducible every state can be reached from every other state For any i,j, exist an m 0, such that i,j are communicate, if the above condition is valid Irreducible: all states are communicate
More informationParameter Selection Techniques and Surrogate Models
Parameter Selection Techniques and Surrogate Models Model Reduction: Will discuss two forms Parameter space reduction Surrogate models to reduce model complexity Input Representation Local Sensitivity
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 4 (9) 6 96 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A constructive proof on the existence
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 7 (0) 3499 3507 Contents lists available at SciVerse ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Application of the
More informationComputing with inter-spike interval codes in networks ofintegrate and fire neurons
Neurocomputing 65 66 (2005) 415 420 www.elsevier.com/locate/neucom Computing with inter-spike interval codes in networks ofintegrate and fire neurons Dileep George a,b,, Friedrich T. Sommer b a Department
More informationAcommon practice for computing elementary transcendental
678 IEEE TRANSACTIONS ON COMPUTERS, VOL. 53, NO. 6, JUNE 004 Near Optimality of Chebyshev Interpolation for Elementary Function Computations Ren-Cang Li Abstract A common practice for computing an elementary
More informationRecurrent neural networks with trainable amplitude of activation functions
Neural Networks 16 (2003) 1095 1100 www.elsevier.com/locate/neunet Neural Networks letter Recurrent neural networks with trainable amplitude of activation functions Su Lee Goh*, Danilo P. Mandic Imperial
More informationA peak factor for non-gaussian response analysis of wind turbine tower
Journal of Wind Engineering and Industrial Aerodynamics 96 (28) 227 2227 www.elsevier.com/locate/jweia A peak factor for non-gaussian response analysis of wind turbine tower Luong Van Binh a,, Takeshi
More informationOn the representation of hysteresis operators of Preisach type
Physica B 372 (2006) 40 44 www.elsevier.com/locate/physb On the representation of hysteresis operators of Preisach type R.V. Iyer, R. Paige Department of Mathematics and Statistics, Room 201, Broadway
More information