NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS
|
|
- Sheila Douglas
- 5 years ago
- Views:
Transcription
1 NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs when an unstabe deterministic system of ordinary differentia equations is stabiized by the addition of romness into the system. Noise-induced stabiization is quite an intriguing surprising phenomenon as one s first intuition is often that noise wi ony serve to further destabiize the system. In this paper, we investigate under what conditions one-dimensiona, autonomous stochastic differentia equations are stabe, where we take the notion of stabiity to be that of goba stochastic boundedness. Specificay, we find the minimum amount of noise necessary for noiseinduced stabiization to occur when the drift noise coefficients are power, onentia, or ogarithmic functions.. Introduction This work is motivated by the intriguing phenomenon of noise-induced stabiization. Noise-induced stabiization occurs when the addition of romness to an unstabe deterministic system of ordinary differentia equations (ODEs) resuts in a stabe system of stochastic differentia equations (SDEs). We are particuary motivated to find the minimum amount of noise required for the phenomenon of noise-induced stabiization to occur for one-dimensiona, autonomous SDEs of the form dx(t) = b(x(t))dt + σ(x(t))db(t). Here b(x) is the drift coefficient, which pushes the soution deterministicay in some direction, σ(x) is the noise coefficient, which contros the strength of the noise, B(t) is Brownian motion, a stochastic process where for each fixed t, B(t) is a normay distributed rom variabe with mean 0 variance t. We assume that b(x) σ(x) are West Virginia University, Morgantown, WV 6506 Mercyhurst University, Erie, PA Bethany Lutheran Coege, Mankato, MN Department of Mathematics Statistics, Vaparaiso University, Vaparaiso, IN 46383
2 ALLEN, GEBHARDT, KLUBALL, KOLBA continuous functions that there exists 0 0 such that σ(x) 0 for a x 0. Our sense of stabiity comes from that of goba stochastic boundedness, which is defined more formay beow 3. Definition. X(t) is stabe if for a initia conditions a ɛ > 0, there exists some bound R such that for a t 0. P ( X(t) R) > ɛ In the deterministic setting where X(t) is a soution to an ODE, the definition of stabiity reduces to that of boundedness. Hence, the unstabe deterministic systems that we consider either bow up in finite time or wer off to infinity. We say that noise-induced stabiization occurs when the addition of noise to an unstabe ODE resuts in a stabe SDE. Work by Scheutzow 5 has shown sufficient conditions for the occurrence of noise-induced stabiization. In this paper, we find necessary sufficient conditions for noise-induced stabiization to occur when the drift noise coefficients are restricted to certain forms. Section discusses usefu background information on the techniques that we use to prove noise-induced stabiization. Sections 3, 4, 5, 6 present prove our resuts concerning noise-induced stabiization for when the drift noise coefficients are genera power functions, poynomias, onentia functions, ogarithmic functions, respectivey. Section 7 discusses the probem of the minimum amount of noise required to stabiize an unstabe system when the noise coefficient is not restricted to a particuar form.. Background In this section, we discuss preiminary information on the methods used to find our resuts. In particuar, our work uses a we-known resut from to determine the stabiity of SDEs, specificay whether noise-induced stabiization occurs. For ease, we wi refer to this resut as the Stochastic Stabiity Theorem. Stochastic Stabiity Theorem. Consider the genera form for autonomous first-order SDEs: dx(t) = b(x(t))dt + σ(x(t))db(t) where b(x) σ(x) are continuous functions there exists 0 0 such that σ(x) 0 for a x 0.
3 NOISE-INDUCED STABILIZATION 3 Define the foowing quantities: x b(z) s(x) = σ (z) dz when x > x b(z) σ (z) dz when x < x s(y)dy for x > S(x) = x s(y)dy for x < m(x) = s(x)σ (x) x m(y)dy M(x) = x m(y)dy for x > for x < The SDE is stabe if any ony if there exists 0 such that S( ) =, S( ) =, M( ) <, M( ) <. Note that since the ower imit of integration does not affect the convergence or divergence of these integras, can be any rea number greater than or equa to 0. Due to the terms defined in the Stochastic Stabiity Theorem, there are severa reoccurring, non-trivia integras in our proofs. Hence, we prove beow the convergence or divergence of these types of integras to use as a reference. Lemma. Consider the integra where n is any rea number. e cxn dx () If c 0, then the integra is infinite. () If c < 0 n 0, then the integra is infinite. (3) If c < 0 n > 0, then the integra is finite Proof. () If c 0, then cx n 0 for a x, which impies e cxn for a x. Thus, e cxn dx dx =.
4 4 ALLEN, GEBHARDT, KLUBALL, KOLBA () If c < 0 n 0, then x n for a x, which impies cx n c for a x, which in turn impies that e cxn e c for a x. Thus e cxn dx e c dx =. (3) By the Tayor series ansion for e y, e y yk for any y > 0 k! any nonnegative integer k. Hence, k! for any y > 0 e y y k any nonnegative integer k. If c < 0 n > 0, then e cxn = e c xn where c x n > 0 for a x > 0. Hence, e cxn k! = k! for ( c x n ) k c k x nk any nonnegative integer k for a x > 0. Since it hods for any nonnegative integer k n > 0, we can choose k 0 such that nk 0 >. Thus, k e cxn 0! dx c k 0 x nk 0 dx k 0! = c k 0 (nk0 )x nk 0 Consider the integra = k 0! c k 0 (nk0 ) <. e cxn dx. This integra is we-defined ony when n is an integer, so we wi now restrict to this case. Coroary. Consider the integra e cxn dx, where n is an even integer. () If c 0, then the integra is infinite for any n. () If c < 0 n 0, then the integra is infinite. (3) If c < 0 n > 0, then the integra is finite Proof. If n is an even integer, then e cxn = e c( x)n e cxn dx = e cxn dx.
5 NOISE-INDUCED STABILIZATION 5 Thus, we get the same cassification of the integra being finite versus infinite as in the previous emma. Coroary. Consider the integra e cxn dx, where n is an even integer. () If c 0, then the integra is infinite for any n. () If c > 0 n 0, then the integra is infinite. (3) If c > 0 n > 0, then the integra is finite Proof. If n is an odd integer, then e cxn = e c( x)n e cxn dx = e cxn dx. Thus, we get amost the same cassification of the integra being finite versus infinite as in the previous emma, with just the sign of c reversed. 3. Power Function Stabiization In this section, we consider ODE s where the drift coefficient is a genera power function i.e. dx(t) = r X(t) q dt. These ODE s are unstabe for any r 0. Theorem. Consider the SDE dx(t) = b(x(t))dt + σ(x(t))b(t) where b(x(t)) = a X(t) q if q 0 σ(x(t)) = r X(t) p if p 0. If q < 0 then the drift coefficient takes the foowing form: { r X(t) q for X(t) b(x(t)) = r for X(t) <. Likewise, if p < 0, then the noise coefficient takes the foowing form: { a X(t) p for X(t) σ(x(t)) = a for X(t) <. Aso, r, a, q, p are any rea numbers. Noise-induced stabiization occurs if ony if r 0 one of the foowing sets of conditions is met: { p > max, q + } for any a 0 or
6 6 ALLEN, GEBHARDT, KLUBALL, KOLBA Figure. This image depicts an unstabe power ODE that is stabiized with the sufficient amount of noise. p = q + { a q > r for 0 < q a r for q >. Figure shows three separate graphs depicting the phenomenon of noise-induced stabiization. The graph on the far eft shows an ODE that diverges off to infinity for a positive initia vaues; therefore, the ODE is unstabe. The midde graph depicts the noise vaue of X(t)dB(t) added to the previousy unstabe ODE. This specific amount of noise added does not stabiize the system as a < r. The fina image depicts 3X(t)dB(t) added to the unstabe ODE. It can be seen that the SDE converges. This occurs since a > r. Proof. The Stochastic Stabiity Theorem wi be used to show under what conditions noise-induced stabiization occurs. With = evauation of the s(x) term gives x r z q dz when x > a s(x) = z p x r z q dz when x < a z p { r x a = zq p dz when x > x ( z)q p dz when x < r a where the integration depends on the vaue of q p. The non-critica q+ case where p wi be considered first.
7 NOISE-INDUCED STABILIZATION 7 Simpification of Proof: Y (t) = X(t) must have the same exact stabiity as X(t) since they have the same magnitude. dy (t) = dx(t) = r X(t) q dt a X(t) p db(t) = r Y (t) q dt a Y (t) p db(t) = r Y (t) q dt a Y (t) p db(t) Hence, the stabiity of X(t) with r, q, a, p must be equivaent to the stabiity with r, q, p, a. Thus, when proving Theorem 3, it suffices to just prove the case with r > 0. Non-Critica Case: Assume q p et c = n = q p +. Then { c(x n ) when x > s(x) = c(( x) n ) when x < { c cx n when x > = c c( x) n when x < r a (q p+) therefore S(x) = { c x cyn dy when x > c x c( y)n dy when x <. Evauation of the m(x) term gives then m(x) = M(x) = { { Letting x = gives c x p cx n a when x > c ( x) p c( x) n a when x < c x a y p cy n dy when x > c x a ( y) p c( y) n dy when x < S( ) = c cy n dy () M( ) = a c y p cy n dy ()
8 8 ALLEN, GEBHARDT, KLUBALL, KOLBA etting x = gives S( ) = c = c = c c( y) n dy c( y) n dy cy n dy (3) M( ) = a c ( y) p c( y) n dy = a c ( y) p c( y) n dy = a c y p cy n dy (4) Non-Critica Case Assume r > 0, p > q+, p > which means r n = q p + < 0 c = > 0. It foows that c a (q p+) cyn > 0 therefore c cy n > on the interva from to. The ower bound of cy n aows for the foowing comparison of (): S( ) > c dy =. Therefore, S( ) =. The upper bound of cy n aows for the foowing comparison of (4): M( ) a c y p cdy = y p dy a > After canceation of the onentia terms, the integra takes the form y p dy. This integra converges if p < which is the case since it was assumed that p >. Therefore, M( ) <. The constants n c are both negative so 0 > cy n c on the interva. Then > cy n c. The ower bound of
9 cy n aows for the comparison S( ) c NOISE-INDUCED STABILIZATION 9 cdy = dy = of () the concusion that S( ) =. Therefore, S( ) =. The upper bound of cy n enabes the comparison M( ) < a c dy < yp of (3). Once again, the integra y p dy converges since p >. Therefore, M( ) < 0. Through appication of the Stochastic Stabiity Theorem, it has been shown that noise-induced stabiization occurs when r > 0, p > q+, p >. Recaing the proof simpification, it can aso be concuded that noise-induced stabiization occurs when r > 0, p > q+, p >. Non-Critica Case : Assume r > 0 p < q+. Then n = q p+ > r 0, c = < 0, c > 0. Therefore, the integra contained a (q p+) in () is of the form cy n dy where c is negative n is positive. Lemma.3.4 indicates that integras of this form con- verge. Therefore, S( ) <. Appication of the Stochastic Stabiity Theorem shows that noise-induced stabiization does not occur when r > 0 p < q+. The proof simpification indicates that noiseinduced stabiization aso does not occur when r > 0 p < q+ Non-Critica Case 3: Assume r > 0 > p > q+. Then n is a negative constant c is a positive constant. Therefore, cy n > on the interva from to. This aows for the foowing comparison of (4): M( ) < a c y p dy. Since p <, the integra y p dy must diverge. Therefore, M( ) = so noise-induced stabiization does not occur when r > 0 > p > q+. The simpification of the proof indicates that stabiization aso does not occur when r > 0 > p > q+ Critica Case: Assume q p =. Then the s(x) term takes the..
10 0 ALLEN, GEBHARDT, KLUBALL, KOLBA form { r x a s(x) = z dz when x > r x a ( z) dz when x < { r n(x) when x > = a r n( x) when x < a {x r a for x > = ( x) r a for x < { x s(y)dy for x > S(x) = x s(y)dy for x < { x = y r a dy for x > r ( y) a dy for x <. x Evauation of the m(x) term gives { x r m(x) = a a p for x > ( x) r a a p for x < M(x) = { a x y r a x a p dy for x > r ( y) a p dy for x <. Evauation of the terms at x = x = gives S( ) = S( ) = = = M( ) = a y r a dy, (5) ( y) r a dy ( y) r a dy y r a dy, (6) y r a p dy, (7)
11 NOISE-INDUCED STABILIZATION M( ) = ( y) r a = ( y) r a = a a p dy a p dy y r a p dy. (8) Critica Case : Assume r > 0, p = q+ >, a r a q > r. It foows from the assumptions that 0 < r for q >. Therefore, a r < r which means the integras in both (5) a a (6) S( ) = S( ) = y r a dy y r a dy must diverge. Likewise, if q, it must be the case that q < r < r < q so that the integras must sti diverge. Therefore, a a S( ) = S( ) =. Note, that in the case where r is negative, r > r. This is why the magnitude of r must be considered in a a the stabiization conditions. The assumption that r < a q p = q+ is equivaent to the assumption that r p < for q. Since r is positive it can further a be impied that r p <. The fact that both of these terms are a stricty ess than negative means the integras in both (7) (8) M( ) = a M( ) = a y r a p dy y r a p dy must converge. For q >, it must be the case that q > r which a impies r p < r p < so the integras once again converge. a a Therefore, M( ) M( ) are both finite. Through appication of the Stochastic Stabiity Theorem, it has been
12 ALLEN, GEBHARDT, KLUBALL, KOLBA shown that noise-induced stabiization occurs when r > 0, p = q+ >, a r a q > r. Recoection of the proof simpification motivates the concusion that stabiization aso occurs when r > 0, p = q+ >, a r a q > r. Critica Case : Assume r > 0, p = q+, a < r. Then it must be the case that r < which impies the integra contained a in (5) must converge. Therefore, S( ) < so noise-induced stabiization does not occur. Likewise, if r < 0, p = q+ a < r, then r > r < so the integra contained in (6) a a must converge. Therefore, S( ) > noise-induced stabiization does not occur when r > 0, p = q+, a < r. The proof simpification aows for the additiona concusion that stabiization does not occur when r > 0, p = q+, a < r. Critica Case 3: Assume r > 0, p = q+ a q r. Then it must be the case that r p which impies the integra contained in (7) a diverges. Therefore, M( ) = noise-induced stabiization does not occur. Now assume r < 0, p = q+ a q r. Then it foows that r p which impies the integra contained in (8) diverges. a Therefore, M( ) = noise-induced stabiization does not occur when r > 0, p = q+ a q r. The proof simpification indicates that stabiization aso does not occur when r > 0, p = q+ a q r. Critica Case 4: Assume p = q+ assume a q > r. These assumptions impy r < q 0. However, this is not possibe. Therefore, a there are no cases where both p = q+ a q > r. 4. Poynomia Stabiization In this section, we investigate stabiizing dx(t) = b(x(t))dt, where b(x) is any poynomia of degree q. If r is the coefficient of the highest degree term, this ODE is unstabe when q is even for any r 0 when q is odd for any r > 0. The emma beow, as we as the foowing coroaries, describe genera resuts about poynomias which wi be used in the proof of noiseinduced stabiization for when the drift noise coefficients are poynomias.
13 NOISE-INDUCED STABILIZATION 3 Lemma. Let f(x) be any poynomia of degree q, where the coefficient of the highest degree term is r > 0. Then, there exists k 0 such that for a x k. r(x k) q f(x) r(x + k) q Proof. If q = 0, then f(x) = r the caim is true for k = 0. Hence, assume q. Since f(x) is a poynomia of degree q, it has the foowing form: f(x) = rx q + r x q + r x q r q x + r q ( = r x q + r r xq + r r xq r q r x + r ) q r where r, r,..., r q, r q are any rea constants. Let r max = max{ r r, r r,..., r q r, r q r } k = max{q r max, }. Then for x k, factoring gives Hence for x k, rx q kx q f(x) rx q + kx q rx q (x k) f(x) rx q (x + k). r(x k) q f(x) r(x + k) q. Coroary 3. Let f(x) be any poynomia of degree q, where the coefficient of the highest degree term is r < 0. Then, there exists k 0 such that r(x + k) q f(x) r(x k) q for a x k. Coroary 4. Let f(x) be any poynomia of degree q, where the coefficient of the highest degree term is r > 0. Then, there exists k 0 such that if q is even, if q is odd, for a x k. r(x + k) q f(x) r(x k) q r(x k) q f(x) r(x + k) q
14 4 ALLEN, GEBHARDT, KLUBALL, KOLBA Figure. This image depicts an unstabe poynomia ODE that is stabiized with the sufficient amount of noise. Coroary 5. Let f(x) be any poynomia of degree q, where the coefficient of the highest degree term is r < 0. Then, there exists k 0 such that if q is even, if q is odd, for a x k. r(x k) q f(x) r(x + k) q r(x + k) q f(x) r(x k) q Theorem. Let b(x) be any poynomia of degree q, where the coefficient of the highest degree term is r, et σ(x) be any poynomia of degree p, where the coefficient of the highest degree term is a 0. Then noise-induced stabiization of the SDE dx(t) = b(x(t))dt + σ(x(t))db(t) occurs if ony if the corresponding ODE is unstabe one of the foowing conditions is met: p > q+ or p = q+ { a > r for q = a r for q 3. Figure shows three separate graphs depicting the phenomenon of noise-induced stabiization where the drift noise coefficients are poynomias. The graph on the far eft shows an ODE that diverges off to infinity for a positive initia vaues; therefore, the ODE is unstabe. The midde graph depicts the noise vaue of (X(t) + )db(t) added to the previousy unstabe ODE. This specific amount of noise added does not stabiize the system as p < q+. The fina image depicts
15 NOISE-INDUCED STABILIZATION 5 (X(t) + ) db(t) added to the unstabe ODE. It can be seen that the SDE converges. This occurs since p > q+. Proof. Suppose r > 0. By Lemma Coroary 4, there exist some k b k σ such that x s(x) x x s(x) x r(z+k b ) q a (z k σ) p dz r(z k b ) q a (z+k σ) p dz for even q, x s(x) x x s(x) x r(z k b ) q a (z+k σ) dz p r(z+k b ) q a (z k σ) p dz r(z+k b ) q a (z k σ) dz p r(z+k b ) q a (z k σ) p dz r(x k b ) q a (x+k σ) dz p r(z k b ) q a (z+k σ) p dz when x > max (k b, k σ ) when x < max (k b, k σ ) when x > max (k b, k σ ) when x < max (k b, k σ ) when x > max (k b, k σ ) when x < max (k b, k σ ) when x > max (k b, k σ ) when x < max (k b, k σ ) (9) (0) () () for odd q. These bounds aow for simper anaysis of conditions set by the Stochastic Stabiity Theorem. First, consider the case where q is even. Suppose p > q+. By (9), x S( ) r(z + k b ) q dz dx. a (z k σ ) p Let y = z k σ. Then, by Lemma, there exists a k > 0 such that x r(y q + S( ) k y q ) dy dx. a y p By dividing through integrating, it can be found that r x q p+ r k x q p S( ) c c dx a q p + a q p r where c = q p+ r k c a q p+ = q p. Since q p <, a q p r k therefore > 0, the smaest r k x q p can be is. a (q p) a q p
16 6 ALLEN, GEBHARDT, KLUBALL, KOLBA Therefore, S( ) c c r x q p+ dx. a q p + Since q p+ < 0, this integra diverges by Lemma. By comparison, S( ) =. With these bounds on s(x), S( ) x r(z + k b ) q dz dx. a (z k σ ) p Let y = z k σ. Then, by Coroary 4, there exists a k > 0 such that x r(y q + S( ) k y q ) dy dx. a y p By methods simiar to those used previousy, it can be found that r x q p+ S( ) c dx a q p + where c = c = r ( ) q p+ a q p+ r k a ( ) q p q p, this integra diverges to. Thus, S( ) =. Consider M( ). By Lemma, x M( ) a (x k σ ) p. By Coroary r(z q + k z q ) dz dx. a z p There argest the onentia term can be is, therefore where c = M( ) c c a (x k σ ) r q p+ c a q p+ = r k a p dx q p q p. This integra converges for p > q+, which is aways true since p > q 0. Thus, M( ) <. Proof that M( ) < foows simiary. By Coroary 4, x r(z q + M( ) k z q ) dz dx a (x + k σ ) p a z p. The argest the onentia term can be is c where c = Thus, M( ) c dx a (x + k σ ) p, r a ( ) q p+ q p+.
17 where c = NOISE-INDUCED STABILIZATION 7 r k ( ) q p. This integra converges for p >, a q p which again is aways true. Therefore M( ) <, noiseinduced stabiization does occur when p > q+, r > 0, q is even. Now, assume q is odd. The proofs that S( ) = M( ) < are identica to the case where q is even. Now consider S( ). By (), S( ) x r(z k b ) q dz dx. a (z + k σ ) p Let y = z + k σ. Then, by Coroary 4, there exists a k > 0 such that x r(y q S( ) k y q ) dy dx. a y p By dividing through integrating, it can be found that r x q p+ r k x q p S( ) c c dx a q p + a q p r ( ) where c = q p+ r k c a q p+ = ( ) q p a q p is odd,. Therefore, r k a x q p q p S( ) c c r x q p+ dx. a q p +. Since q p By Coroary, this integra diverges to negative infinity, therefore S( ) =. Now consider M( ). By Coroary 4, x r(z q M( ) k z q ) dz dx a (x + k σ ) p a z p. The argest the onentia term can be is. Therefore, M( ) c c dx a (x + k σ ) p, r ( ) where c = q p+ r k c a q p+ = ( ) q p. Again, this a q p integra converges therefore M( ) <. Now, it is sufficient to say that noise-induced stabiization occurs when p > q+, r > 0, q is odd.
18 8 ALLEN, GEBHARDT, KLUBALL, KOLBA Suppose p = q. By (0), S( ) can be bounded such that x r(z k b ) q S( ) dz dx. a (z + k σ ) p Let y = z + k σ,, by Lemma, there exists a k > 0 such that x r(y q S( ) k y q ) dy dx. a y p By dividing through integrating, it can be found that x r x S( ) a dz r k a z dz dx rx = c c x r k a dx, a where c = r a c = r k a. This integra converges, therefore S( ) <. Suppose p < q. By the same bounds on s(x), S( ) x r(z k b ) q dz a (z + k σ ) p dx. Let y = z + k σ. By Lemma, there exists a k > 0 such that x r(y q S( ) k y q ) dy dx a y p r x q p+ r k x q p = c c dx, a q p + a q p r where c = q p+ r k c a q p+ = q p. For sufficienty a q p arge x, r x q p+ r x q p+ S( ) c c dx a q p + a q p + r x q p+ = c c dx. a q p + By Lemma this integra converges, thus S( ) <. Therefore noise-induced stabiization does not occur when p < q+ r > 0.
19 NOISE-INDUCED STABILIZATION 9 Suppose p = q+, suppose a > r when q = a r when q 3. By (9), S( ) x r(z + k b ) q dz dx. a (z k σ ) p Let y = z k σ. By Lemma, there exists a k > 0 such that x r(y q + S( ) k z q ) dy dx a y p = c c x r r k a dx, a x where c = r a c = r k. This integra diverges for a a r, therefore S( ) =. Now, consider S( ). Note that since p = q+, since p is an integer, q must be odd. By (), S( ) x r(z k b ) q dz dx. a (z + k σ ) p Let y = z + k σ. By Coroary 4, there exists a k > 0 such that x r(z q S( ) k z q ) dz dx a z p c c x r r k a dx, a x where c = r a c = r k. This integra diverges to a negative infinity for a r, therefore S( ) =. Consider M( ). By Lemma, M( ) x a (x k σ ) p r(z q + k z q ) dz dx. a z p The onentia term can be partiay bounded such that, M( ) c r x a a (x k σ ) p dx,
20 0 ALLEN, GEBHARDT, KLUBALL, KOLBA where c = r a. For sufficienty arge x, a (x k σ) p p a x p. Therefore, M( ) c = c p r x a a xp dx p a x r a p dx. This integra converges for a q > r, therefore M( ) <. Consider M( ). By Coroary 4, x r(z q M( ) k z q ) dz dx a (x + k σ ) p a z p = c r r k c x a a (x + k σ ) p dx a x, where c = r a r k c = r k. The argest can a a x be is. Aso, for sufficienty arge x, M( ) c p r c x a a xp dx = c p r c x a a p dx This integra converges, therefore M( ) <. Now, suppose a r when q = a < r when q 3. Consider S( ). By the bounds on s(x), x r(z k b ) q S( ) dz dx. a (z + k σ ) p Let y = z + k σ. By Lemma, there exists a k > 0 such that x r(y q S( ) k z q ) dy dx a y p = c c x r r k a dx, a x where c = r a c = r k. Since r k, a a x S( ) c c x r a dx. This integra converges when a r, therefore S( ) <. Thus, when p = q+, r > 0, noise-induced stabiization occurs if ony
21 NOISE-INDUCED STABILIZATION if a > r when q = a r when q 3. Let Y (t) = X(t), suppose q is even. Since Y (t) is bounded if ony if X(t) is bounded, Y (t) must have the exact same stabiity as X(t). It foows that dy (t) = dx(t) = b(x(t)dt + σ(x(t))db(t) = b( Y (t)dt + σ( Y (t))db(t). Since the highest degree term of the drift coefficient determines the stabiity, since q is even, the stabiity of X(t) with eading coefficient r is equivaent to its stabiity with eading coefficient r. Since the noise coefficient, σ(x), is aways squared in the terms defined by the Stochastic Stabiity Theorem, the sign of the coefficient pays no roe in the stabiity of the SDE. Note that this ony hods for even q since the origina ODE is stabe for r < 0 q odd. 5. Exponentia Function Stabiization This section considers ODE s where the drift coefficient is an onentia function, i.e. dx(t) = r(x(t)) q dt. where r q are any rea number. The ODE is unstabe for any r Genera Exponentia Stabiization. Theorem 3. Consider the SDE dx(t) = r(x(t)) q dt + a(x(t)) p db(t) where r, a, q, p are rea numbers. Noise-induced stabiization occurs if ony if r 0, a 0 one of the foowing sets of conditions is met: r > 0 p > max{0, q } or r < 0 p > max{0, q }. Figure 3 shows three separate graphs depicting the phenomenon of noise-induced stabiization. The graph on the far eft shows an ODE that diverges off to infinity for a positive initia vaues; therefore, the ODE is unstabe. The midde graph depicts the noise vaue of 0 X(t) db(t) added to the previousy unstabe ODE. This specific amount of noise added does not stabiize the system as p <. The 4
22 ALLEN, GEBHARDT, KLUBALL, KOLBA Figure 3. This image depicts an unstabe onentia ODE that is stabiized with the sufficient amount of noise. fina image depicts X(t) db(t) added to the unstabe ODE. It can be seen that the SDE converges. This occurs since p >. Proof. The Stochastic Stabiity Theorem wi be used to show under what conditions noise-induced stabiization occurs. Simpification of Proof: Y (t) = X(t) must have the same exact stabiity as X(t) since they have the same magnitude. dy (t) = dx(t) = r(x(t)) q dt a(x(t)) p db(t) = r( Y (t)) q dt a( Y (t)) p db(t) = r(y (t)) q dt a(y (t)) p db(t) Hence, the stabiity of X(t) with r must be equivaent to the stabiity with r, but with q, p, a substituted for q, p, a. Thus, when proving Theorem 3, it suffices to just prove the case with r > 0.
23 NOISE-INDUCED STABILIZATION 3 With =, the s(x) term takes the form x r qz a s(x) = dz for x > pz x r qz dz for x < a pz { r x a = (q p)zdz for x > r x (q p)zdz for x < a { r x a = kzdz for x > r x kzdz for x < a where k = q p. The integration depends on the vaue of k. The non-critica case where k 0 wi first be considered. Non-Critica Case: Assume k 0. { r s(x) = (kx k) for x > a k r(kx k) for x < a k r k r a = kx for x > k a k r k r kx for x < a k a k S(x) = = where c = m(x) term gives m(x) = M(x) = = r k a k r k a k { x c r x c r r k a k { x r ky dy for x > a k x r ky dy for x < a k ky dy for x > a k ky dy for x < a k c = r k a k. Moving on to the a c r kx px for x > a k a c r kx px for x < a k { x a c { x a c r ky pydy for x > a k x r ky a c pydy for x < a k r ky py dy for x > a k x r ky py a c dy for x <. a k
24 4 ALLEN, GEBHARDT, KLUBALL, KOLBA Setting x equa to positive infinity gives r S( ) = c a k ky dy (3) M( ) = a c r ky py dy. (4) a k Setting x equa to negative infinity gives r S( ) = c a k ky dy r = c a k ky dy r = c a k ky dy (5) M( ) = a c = a c = a c = a c r a k ky py dy r ky py dy a k r ky + py dy a k ( ) r ky + a kp a k r y dy. (6) Non-Critica Case : Assume r > 0 p > max{0, q }. Since p > max{0, q }, it must be the case that k = q p < 0 p > 0. Since k is negative r is positive, it foows that the constant r a k is positive therefore the onentia terms r ky a k r ky must be greater than on the interva from to. a k This aows for the comparison S( ) c dy = S( ) c dy = the concusion that (3) (5) diverge. Therefore S( ) = S( ) =. The onentia decay term ky bounds
25 NOISE-INDUCED STABILIZATION 5 the onentia term r ky between r a k a k (ower bound) (upper bound). This aows for the foowing comparison of (4) to a smaer integra: M( ) a c < py dy The smaer integra has the form pydy with where the constant p is negative. Then Lemma indicates that the integra converges M( ) <. Since k is positive, the Tayor series ansion indicates ky > ky on the interva from to. Now, assuming k + a pk > 0, the r foowing comparison of (6) to a smaer magnitude integra can be made: M( ) a c = a c >. r a k r a k ( ky + a pk r ( k + a pk r ) y ) y Since it was assumed that k + a pk > 0, the smaer magnitude integra takes the form r cy n dy where c is negative n =. Then Lemma indicates that integras of this form converge so M( ) < if k + a pk > 0. r Now, assume k + a pk 0 which is equivaent to the assumption that r k a pk. This inequaity aows for the comparison ky r a pk the foowing comparison of (6) to a smaer magnitude r integra: M( ) a c ( r a pk a k r y ) + a pk r y. Recaing that a pk > 0 considering the Tayor Series Expansion, the comparison a pk r ( ) ( ) y > a pk y + a pk r r r y can
26 6 ALLEN, GEBHARDT, KLUBALL, KOLBA be made. This aows for the foowing comparison of the smaer magnitude integra: M( ) a c ( r a pk a k r y + = ( r a pk a c a k r >. ( a pk r ) y ) ) y + a pk r y The fina integra takes the form cy n dy where c = r a k ( a pk r is a negative constant n = is a positive constant. Then Lemma indicates that integras of this form converge so then M( ) <. Appication of the Stochastic Stabiity Theorem indicates that noiseinduced stabiization occurs when r > 0 p > max{0, q }. The simpification of the proof indicates that noise-induced stabiization aso occurs when r > 0 p > max{0, q }. Non-Critica Case : Assume r > 0 q > p. Since q > p, it must be the case that k = q p > 0 r < 0. The Tayor series a k ansion indicates that ky > ky on the interva which aows for the foowing comparison of (3) to a smaer magnitude integra: r S( ) c a k ky dy < The smaer magnitude integra takes the form cy n dy where c = r is a negative constant n =. Then Lemma indicates a that integras of this form converge so S( ) <. Appication of the Stochastic Stabiity Theorem indicates that noiseinduced stabiization does not occur when r > 0 p < q. The simpification of proofs indicates that noise-induced stabiization aso does not occur when r > 0 p < q. Non-Critica Case 3: Assume r > 0 q < p 0. It foows from the assumptions that k = q p < 0 which impies r < 0. The onentia decay term ky bounds r ky between a k r a k a k )
27 NOISE-INDUCED STABILIZATION 7 which aows for the foowing comparison of (4) to a smaer magnitude integra: M( ) a c r pydy. a k Since p is a non-negative constant, it must be the case that py on the interva which aows for the second comparison M( ) r dy a c a k where (4) is being compared to an even smaer magnitude integra which diverges. Therefore, M( ) =. Appication of the Stochastic Stabiity Theorem indicates that noiseinduced stabiization does not occur when r > 0 q < p 0. The simpification of proof indicates that noise-induced stabiization does not occur when r > 0 q < p 0 which can aso be stated as when r < 0 q > p 0. Critica Case: Assume q = p so k = 0 consider when x >. Then the s(x) term takes the form r s(x) = a r = x 0dz (x ) a r rx = a so then the S(x) term takes the form r x ry S(x) = dy a a ( ) r a rx r = a r a a = a r rx r + a a a r a
28 8 ALLEN, GEBHARDT, KLUBALL, KOLBA Considering when r < 0 setting x = gives S( ) = a r = a r <. r a r a + a r The Stochastic Stabiity Theorem indicates that noise-induced stabiization does not occur when r < 0 p = q. The simpification of proof indicates that noise-induced stabiization aso does not occur when r < 0 p = q. The genera onentia stabiization theorem shows that sufficienty arge onentia growth can stabiize any onentia drift coefficient. Recognizing the difference between the drift coefficient eriencing onentia growth as opposed to onentia decay, the foowing sections consider them separatey. 5.. Exponentia Growth Stabiization. In this section, stabiization of onentia growth in the drift coefficient is considered which means r q have the same sign. Theorem 4. Consider the SDE dx(t) = r qx(t) dt+σ(x(t))db(t) where a X(t) m px(t) for sgn(r)x σ(x(t)) = a px(t) for 0 sgn(r)x < a nx(t) for sgn(r)x < 0 r, a, q, p, m, n R. More specificay, consider when r, q n have the same sign p = q. Noise-induced stabiization occurs if ony if r 0, a 0, one of the foowing sets of conditions is met: m >, n > q, m =, n > q, a r, m >, n = q a n < r or m =, n = q <, a r.
29 NOISE-INDUCED STABILIZATION 9 Simpification of Proof: Y (t) = X(t) must have the same exact stabiity as X(t) since they have the same magnitude. dy (t) = dx(t) r(x(t)) q dt a X(t) m (X(t)) p db(t) for sgn(r)x(t) > r(x(t)) q dt a(x(t)) p db(t) = for 0 sgn(r)x(t) < r(x(t)) q dt a(x(t)) n db(t) for sgn(r)x(t) < r( Y (t)) q dt a Y (t) m ( Y (t)) p db(t) for sgn(r) Y (t) > r( Y (t)) q dt a( Y (t)) p db(t) = for 0 sgn(r) Y (t) < r( Y (t)) q dt a( Y (t)) n db(t) for sgn(r) Y (t) < r(y (t)) q dt a Y (t) m (Y (t)) p db(t) for sgn( r)y (t) > r(y (t)) q dt a(y (t)) p db(t) = for 0 sgn( r)y (t) < r(y (t)) q dt a(y (t)) n db(t) for sgn( r)y (t) < Hence, the stabiity of X(t) with r q, m, a, p, n must be equivaent to the stabiity with r, q, a, m, p, n. Thus, when proving Theorem 3, it suffices to just prove the case with r > 0. Proof. The Stochastic Stabiity Theorem wi be used to show the conditions for noise-induced stabiization. With = r > 0, the s(x) term takes the form { r x a s(x) = z m dz for x > r x (q n)z dz for x < a where the integration depends on the vaue of m.
30 30 ALLEN, GEBHARDT, KLUBALL, KOLBA Non-Critica Case: Assuming m n q gives r a s(x) = ( m) (x m ) for x > r ((q n)x n q for x < a (q n) r r a = ( m) a ( m) x k for x > r n q r (q n)x for x < a (q n) a (q n) r c a = ( m) x m for x > r c (q n)x for x < a (q n) x c S(x) = r y m dy for x > a( m) x c r (q n)y dy for x < a (q n) where c = the m(x) term gives r a ( m) c = c m(x) = x m r a a ( m) x m px r c (q n)x nx a a (q n) r n q. Evauation of a (q n for x > for x < x c M(x) = a y m r a ( m) y m py dy for x > x r (q n)y ny c dy for x <. a a (n q) Letting x = gives r S( ) = c a ( m) y m dy (7) M( ) = c a y m r a ( m) y m py dy. (8)
31 NOISE-INDUCED STABILIZATION 3 Likewise, etting x = gives r S( ) = c a (q n) r = c a (n q) = c M( ) = c a = c a = c a (q n)y dy (q n)y r (n q)y a (n q) r a (n q) r a (n q) dy dy (9) (q n)y ny dy (q n)y ny r (n q)y + ny a (n q) dy dy. (0) Non-Critica Case : Assume r > 0, q > 0, n > 0, m > n > q. The integra in (7) takes the form cy n dy where c = r is a positive constant so the integra must diverge by Lemma a ( m). Therefore, S( ) =. Since m >, then m < 0 so the term y m r bounds a ( m) y m r between a ( m) y m. This aows for the foowing comparison of (8) to the smaer magnitude integra: M( ) c a r y m py. a ( m) The onentia decay term py aows for the further comparison of (8) to an even smaer magnitude integra: M( ) c a r y m dy. a ( m) Once again, since m >, it must be the case that this fina integra converges. Therefore, M( ) <. The conditions for x < are the same as for Case of the genera onentia theorem. Therefore, S( ) = M( ) <. Appication of the Stochastic Stabiity Theorem indicates that noise-induced stabiization occurs when r > 0, q > 0, m > n > q. The proof simpification indicates that stabiization aso occurs when r > 0, q > 0, m >, n > q.
32 3 ALLEN, GEBHARDT, KLUBALL, KOLBA Non-Critica Case : Assume r > 0, q > 0, m <. Then the integra in (7) takes the form cy n r dy where c = is a ( m) a negative constant n = m is a positive constant. Then Lemma indicates that integras of this form converge so S( ) < noise-induced stabiization does not occur when r > 0, q > 0, m <. Then the proof simpification aows for the concusion that noise-induced stabiization aso does not occur when r > 0, q > 0, m <. Non-Critica Case 3: Assume r > 0 0 < n < q. It foows from the assumption that n q is a negative constant. Therefore, the onentia decay term (n q)y bounds the onentia term r (n q)y a (n q) between (upper bound) ower bound on the interva r a (n q) from to. This fact aows for the foowing comparison of (0) M( ) c a ny dy The integra (0) is compared to has the form cy n dy where c = n n =. Then Lemma indicates that integras of this form diverge so M( ) = noise-induced stabiization does not occur when r > 0 0 < n < q. Additionay, the proof simpification indicates that noise-induced stabiization aso does not occur when r > 0 0 < n < q. Critica Case: Now, considering when m = x, the s(x) term takes the form r x s(x) = z dz a r = n(x) then = x r a S(x) = a x y r a dy.
33 Evauation of m(x) gives Letting x = gives NOISE-INDUCED STABILIZATION 33 m(x) = a x r a px M(x) = a x S( ) = M( ) = a y r a pydy. y r a dy () y r a pydy. () Critica Case : Assume r > 0, m =, n > q > 0, a r. Once again, the conditions for x < are the same as for Case of the previous theorem. Therefore, S( ) = M( ) <. The assumption that a r is equivaent to the assumption that r. Then the integra in () diverges S( ) =. The a assumptions aso impy that r 0. Then y r a a on the interva which aows for the foowing comparison of () to a arger magnitude integra. M( ) pydy a The smaer magnitude integra has the form cy n dy where c = p is a negative constant n =. Then Lemma indicates that integras of this form converge so M( ) <. Appication of the Stochastic Stabiity Theorem indicates that noise-induced stabiization occurs when r > 0, m =, n > q > 0, a r. The proof simpification further concudes that stabiization aso occurs r > 0, m = q, n > > 0, a r. Critica Case : Assume r > 0, m =, a < r. It foows from the assumptions that r <. Therefore, the integra contained in a () converges. So S( ) = noise-induced stabiization does not occur when r > 0, m =, a < r. The proof simpification indicates that noise-induced stabiization aso does not occur when r > 0, m =, a < r.
34 34 ALLEN, GEBHARDT, KLUBALL, KOLBA Semi-Critica Case: Consider when r > 0 n = q. Then for x <, the s(x) term takes the form r x s(x) = dz a a r = r r x r S(x) = a a Evauation of the m(x) term gives m(x) = a r a a x, y dy. ( ) r a n x M(x) = r x ( ) r a a a n y dy. Setting x = gives r S( ) = a r a y dy. (3) M( ) = r ( a n r ) y dy. (4) a a Semi-Critica Case : Assume r > 0, m >, n = q > 0 a n < r. It was proven in Non-Critica Case that S( ) = M( ) < under these conditions. The integra in (3) has the form cy n dy where c = r is positive n =. Then Lemma indicates that integras of this form diverge. Therefore, S( ) =. The assumption a that a n < r is equivaent to the assumption that n r < 0. Therefore, the integra contained in (4) takes the form a cy n dy where c = n r is a negative constant n =. Then Lemma indicates a that integras of this form converge so M( ) <. Appication of the Stochastic Stabiity Theorem indicates that noise-induced stabiization occurs when r > 0, m >, n = q > 0 a n < r. The proof simpification aows for the concusion that stabiization aso occurs when r > 0, m > q, n > > 0, a n < r. Semi-Critica Case : Assume r > 0, m =, 0 < n = q <, a r. It was proven in Critica Case that S( ) =
35 NOISE-INDUCED STABILIZATION 35 M( ) < under these conditions. It was proven in Semi-Critica Case that S( ) = under these conditions. Since a r, it must be the case that r. Then the assumption that n < impies a n r < 0. The integra contained in (4) has the form cy n dy a where c = n r is negative n. Then Lemma indicates that a integras of this form converge. Therefore M( ) < the Stochastic Stabiity Theorem indicates that noise-induced stabiization occurs when r > 0, m =, 0 < n = q <, a r. The proof simpification aows for the concusion that stabiization aso occurs when r > 0, m = q, 0 < n = <, a r Semi-Critica Case 3: Assume r > 0, n = q > 0 a n r. Then n r 0. The integra contained in (4) takes the form a cy n dy where c = n r 0. Then Lemma indicates that a integras of this form diverge so M( ) = noise-induced stabiization does not occur when r > 0, n = q > 0 a n r. Then the proof simpification indicates that stabiization aso does not occur when r > 0, n = q > 0, a n r. Semi-Critica Case 4: Assume r > 0, m =, a r, n = q r. It foows that since n, it must be the case a that n r 0. Since (4) takes the form cy n dy where a c = n r 0, the integra must diverge by Lemma. Therefore, a noise-induced stabiization must not occur when r > 0, m =, a r, n = q. Then stabiization aso does not occur when r > 0, m =, a r, q n =. The onentia growth theorem s on the genera onentia theorem by defining the noise coefficient piece-wise therefore aowing for the critica case where p = q Exponentia Decay Stabiization. In this section, stabiization of onentia decay in the drift coefficient is considered which means r q have the opposite sign. Theorem 5. Consider the SDE dx(t) = r qx(t) dt+a X(t) p db(t) where r q have opposite signs. Noise-induced stabiization occurs if ony if r 0, a 0, p >. Proof. The Stochastic Stabiity Theorem wi be used to show under what conditions noise-induced stabiization occurs.
36 36 ALLEN, GEBHARDT, KLUBALL, KOLBA Simpification of Proof: Y (t) = X(t) must have the same exact stabiity as X(t) since they have the same magnitude. dy (t) = dx(t) = r(x(t)) q dt a X(t) p db(t) = r( Y (t)) q dt a Y (t) p db(t) = r(y (t)) q dt a Y (t) p db(t) Hence, the stabiity of X(t) with r, q, a, p must be equivaent to the stabiity with r, q, p, a. Thus, when proving Theorem 3, it suffices to just prove the case with r > 0. The s(x) term takes the form { r x a s(x) = z p qzdz for x > x ( z) p qzdz for x < r a { x S(x) = r a x r a y z p qzdz dy for x > y ( z) p qzdz dy for x <. The m(x) term then takes the form { x p r x a m(x) = a z p qzdz for x > ( x) p r x a a ( z) p qzdz for x < { x a M(x) = y p r y a z p qzdz dy for x > x a ( y) p r y a ( z) p qzdz dy for x <. If x =, then S( ) = M( ) = a r y z p qzdz dy (5) a r y y p z p qzdz dy. (6) a
37 If x =, then S( ) = = NOISE-INDUCED STABILIZATION 37 r a r a y y M( ) = r ( y) p a a = r y p a a ( z) p qzdz dy z p qzdz y y dy (7) ( z) p qzdz dy z p qzdz dy (8) Case : Assume r > 0, q < 0, p > r. Then is a negative constant a since qz is an onentia decay term, the foowing comparison of (5) can be made: r y S( ) z p dz dy a With the assumption that p >, the comparison takes the form r r S( ) a ( p) a ( p) y p dy. Now the integra takes the form cy n dy where c = r is a a ( p) positive constant. Then Lemma indicates that integras of this form diverge so S( ) =. Since r > 0 qy is an onentia a growth term, the foowing comparison hods for (7): r y S( ) z p dz dy a r r = a ( p) a ( p) y p dy The fina comparison integra takes the form cy n dy where c = r is a negative constant n = p is a negative constant. a ( p) Then Lemma indicates that integras of this form diverge so S( ) =. Using the fact that r > 0 qz is an onentia decay a term, the foowing comparison of (6) can be made: M( ) r y y p z p dz dy a a = a r y p r a ( p) a ( p) y p dy.
38 38 ALLEN, GEBHARDT, KLUBALL, KOLBA The onentia term r a ( p) y p must be ess than or equa to on the interva which aows for the foowing comparison of (6): M( ) a r y p dy. a ( p) This fina integra must converge since p must be ess than. Therefore, M( ) <. Once again using the fact that r is a a negative constant qz is an onentia growth term, the foowing comparison of (8) can be made: M( ) r y y p z p dz dy a a = a r y p r a ( p) a ( p) y p dy. Now, the onentia term r a ( p) y p r is bounded between (upper bound) (ower bound) a ( p) due to p being a negative power. This aows for the foowing comparison: M( ) a r y p r dy a ( p) a ( p) = y p dy a >. Once again, this fina integra must converge since p is stricty ess than. Therefore, M( ) <. Appication of the Stochastic Stabiity Theorem indicates that noise-induced stabiization occurs when r > 0, q < 0, p >. Then stabiization aso occurs when r > 0, q < 0, p >. Case : Assume r > 0, q < 0, p r. Since is a positive constant, a qz is an onentia decay term, the foowing comparison of
39 (6) hods: M( ) r y p a a = a r a ( p) NOISE-INDUCED STABILIZATION 39 y z p dz dy y p r a ( p) y p dy. Since p r, the constant must be positive so that the onentia term a ( p) r a ( p) y p is greater than or equa to on the interva. This fact aows for the additiona comparison of (6) to an even smaer magnitude integra: M( ) a r y p dy. a ( p) Since p, this fina integra must diverge. Therefore, M( ) = noise-induced stabiization does not occur when r > 0, q < 0, p. Then it aso does not occur when r > 0, q < 0, p. The onentia decay theorem shows that an onentia noisecoefficient is not required for stabiization of a onentia drift coefficients. In this theorem, a power noise coefficient was used to stabiize the onentia decay term. 6. Logarithmic Function Stabiization In this section we wi investigate the stabiization of the ODE dx(t) = rn( X(t) ) q dt, where r q are rea numbers with q 0. This ODE is unstabe if ony if r 0. Theorem 6. Consider the SDE { dx(t) = rn( X(t) ) q a X(t) m n( X(t) ) p db(t) for x > dt + a( m )n() p db(t) for x where r, a, q, p, m are rea numbers q, p, m 0. Noise-induced stabiization occurs if ony if the ODE is unstabe, one of the foowing conditions is met: m >, m = q+ p >, or m = q+ p = a q > r.
40 40 ALLEN, GEBHARDT, KLUBALL, KOLBA Figure 4. This image depicts an unstabe ogarithmic ODE that is stabiized with the sufficient amount of noise. Figure 4 shows three separate graphs depicting the phenomenon of noise-induced stabiization. The graph on the far eft shows an ODE that diverges off to infinity for a positive initia vaues; therefore, the ODE is unstabe. The midde graph depicts the noise vaue of X(t) 4 n( X(t) +)db(t) added to the previousy unstabe ODE. This specific amount of noise added does not stabiize the system as m <. The fina image depicts X(t) n( X(t) +) added to the unstabe ODE. It can be seen that the SDE converges. This occurs since m >. Proof. For this SDE we are ooking at b(x) = r n( x ) q σ(x) = a x m n( x ) p where r > 0. Assume there exists some bound such that σ(x) is not undefined, then by the Stochastic Stabiity Theorem x s(x) = x After further simpification, x s(x) = x r n( z ) q a z m n( z ) p dz r n( z ) q a z m n( z ) p dz r n( z ) q p dz a z m r n( z ) q p dz a z m for x > for x <. for x > for x <. (9) In order to find where noise-induced stabiization wi occur, we must separate the vaue of m into three different possibiities: m >, m <, m =. Case : m >
41 NOISE-INDUCED STABILIZATION 4 To show that when m > a vaues of r, a, q, p, m there is noiseinduced stabiization we want s(x) to be greater than or equa to some integra that diverges. Since n( z ) n() for z. Thus, x r n() (q p) dz for x > a s(x) z m x for x <. After integrating, s(x) where c = c = c m(x) r n() (q p) dz a z m r n() (q p) z ( m+) a ( m+) x r n() (q p) z ( m+) a ( m+) r n() (q p) x ( m+) a ( m+) r n() (q p) x ( m+) a ( m+) x r n() (q p) x ( m+) a ( m+) c a x m n( x ) p r n() (q p) x ( m+) a ( m+) c a x m n( x ) p r n() (q p) ( m+) a ( m+) c = for x > for x < for x > for x < for x > for x < First assume q p < 0. Using (30) (3) we know r n() (q p) y ( m+) S( ) c dy, a ( m + ) M( ) S( ) M( ) r n() (q p) y ( m+) a ( m+) c a y m n( y ) p dy, (30) (3) r n() (q p) ( m+) a ( m+) r n() (q p) y ( m+) c dy, a ( m + ) r n() (q p) y ( m+) a ( m+) c a y m n( y ) p Looking at S( ) using Lemma aows us to concude that S( ). For M( ) the onentia wi converge since the coefficient is negative; therefore, as a whoe M( ) <. For S( ) the argest vaue dy..
42 4 ALLEN, GEBHARDT, KLUBALL, KOLBA the onentia function can be is ; therefore, S( ) c dy, which equas. From this it can be concuded that S( ) is infinite. Then M( ) < since the onentia function is converging from its coefficient being negative. Now assume q p > 0. To show s(x) is greater than something that diverges we use the comparison n( z ) q p z m. Then, x s(x) x r z (m ) dz a z m r z (m ) dz a z m for x > for x <, which simpifies to x s(x) x r dz a z (m+ ) r dz a z (m+ ) for x > for x <. Then, s(x) r z ( m+ ) a ( m+ ) r z ( m+ ) a ( m+ ) x x for x > for x < c = c r x ( m+ ) a ( m+ ) r x ( m+ ) a ( m+ ) for x > for x < (3) m(x) r x ( m+ ) a ( m+ ) c a x m n( x ) p r x ( m+ ) a ( m+ ) c a x m n( x ) p for x > for x < (33)
43 NOISE-INDUCED STABILIZATION 43 r where c = ( m+ ) r c a ( m+ ) = ( m+ ). Using (3) a ( m+ ) (33), the Stochastic Stabiity Theorem it can be proven that r y ( m+ ) S( ) c a ( m + ) dy, M( ) S( ) M( ) r y ( m+ ) a ( m+ ) dy, c a y m n( y ) p c r y ( m+ ) a ( m + ) dy, r y ( m+ ) a ( m+ ) dy. c a y m n( y ) p When evauating S( ) the coefficient in the onentia is positive; therefore, by the integra comparison test it can be concuded that S( ) diverges. For M( ) the coefficient of the onentia function is negative, which causes it to converge. Since the onentia is converging, m p 0, then the integra as a whoe is aso converging. In concusion, M( ) <. As for S( ), the smaest vaue the onentia function can be is ; therefore S( ) c dy. This integra diverges to, thus S( ). For simiar reasons as M( ), M( ) converges as we. Finay, assume q p = 0 when m >. Then, x s(x) = x r dz a z m r dz a z m After integrating, it can be said that s(x) = r z ( m+) a ( m+) r z ( m+) a ( m+) x x for x > for x < for x > for x <
44 44 ALLEN, GEBHARDT, KLUBALL, KOLBA r x ( m+) c = c a ( m+) r x ( m+) a ( m+) for x > for x < (34) r x ( m+) a ( m+) c m(x) = for x > a x m n( x ) p (35) r x ( m+) a ( m+) c for x <, a x m n( x ) p where c = r ( m+) c a ( m+) = r ( m+). Again, a ( m+) using (34) (35), the Stochastic Stabiity Theorem it can be determined that r y ( m+) S( ) = c dy, a ( m + ) M( ) = S( ) = M( ) = r y ( m+) a ( m+) dy, c a y m n( y ) p r y ( m+) c dy, a ( m + ) r y ( m+) a ( m+) dy. c a y m n( y ) p When evauating the integra of S( ), Lemma says since the coefficient of the onentia function is positive, S( ) is infinite. The coefficient of the onentia function for M( ) is negative; therefore, the onentia function wi converge. Since m p are positive, the onentia function is converging, the integra of those functions wi converge as we. Thus, M( ) <. As for S( ), the ry c ( m+) dy = c ry a ( m+) ( m+) dy, which diverges to negative infinity. Therefore, S( ) =. From M( ) a ( m+) it is aready known that the function in M( ) is converging. Thus, M( ) <. In summary, when m > aways occurs. we found that noise-induced stabiization
45 NOISE-INDUCED STABILIZATION 45 Case : m < Looking back at (9), in order to show that noise-induced stabiization does not occur for m <, s(x) must be ess than or equa to some function that converges. In this case, we use the comparison n() q p n( z ) q p for z to show x r n() (q p) dz for x > a s(x) z m x for x <. After integration, where c = s(x) c c m(x) r n() (q p) dz a z m r n() (q p) z ( m+) a ( m+) r n() (q p) z ( m+) a ( m+) r n() (q p) x ( m+) a ( m+) r n() (q p) x ( m+) a ( m+) r n() (q p) x ( m+) a ( m+) c a x m n( x ) p r n() (q p) x ( m+) a ( m+) c a x m n( x ) p r n() (q p) ( m+) a ( m+) x x for x > for x < for x > for x < for x > for x < c = (36) r n() (q p) ( m+) a ( m+) First suppose q p > 0, then referring back to (36) can be said that r n() (q p) y ( m+) S( ) c dy. a ( m + ) Using Lemma, S( ) is ess than or equa to a converging integra; therefore, noise-induced stabiization does not occur for m < when q p > 0. Now assume q p < 0. Using the comparison n( z ) q p z m, x s(x) x r z (m ) dz a z m r z (m ) dz a z m for x > for x <..
Wave Equation Dirichlet Boundary Conditions
Wave Equation Dirichet Boundary Conditions u tt x, t = c u xx x, t, < x 1 u, t =, u, t = ux, = fx u t x, = gx Look for simpe soutions in the form ux, t = XxT t Substituting into 13 and dividing
More informationWeek 6 Lectures, Math 6451, Tanveer
Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationMat 1501 lecture notes, penultimate installment
Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationSome Measures for Asymmetry of Distributions
Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationCourse 2BA1, Section 11: Periodic Functions and Fourier Series
Course BA, 8 9 Section : Periodic Functions and Fourier Series David R. Wikins Copyright c David R. Wikins 9 Contents Periodic Functions and Fourier Series 74. Fourier Series of Even and Odd Functions...........
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationTHE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES
THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia
More informationB. Brown, M. Griebel, F.Y. Kuo and I.H. Sloan
Wegeerstraße 6 53115 Bonn Germany phone +49 8 73-347 fax +49 8 73-757 www.ins.uni-bonn.de B. Brown, M. Griebe, F.Y. Kuo and I.H. Soan On the expected uniform error of geometric Brownian motion approximated
More informationu(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0
Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,
More informationarxiv: v1 [math.co] 17 Dec 2018
On the Extrema Maximum Agreement Subtree Probem arxiv:1812.06951v1 [math.o] 17 Dec 2018 Aexey Markin Department of omputer Science, Iowa State University, USA amarkin@iastate.edu Abstract Given two phyogenetic
More informationResearch Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation
Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationA NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC
(January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationTHINKING IN PYRAMIDS
ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think
More informationMARKOV CHAINS AND MARKOV DECISION THEORY. Contents
MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After
More informationBASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set
BASIC NOTIONS AND RESULTS IN TOPOLOGY 1. Metric spaces A metric on a set X is a map d : X X R + with the properties: d(x, y) 0 and d(x, y) = 0 x = y, d(x, y) = d(y, x), d(x, y) d(x, z) + d(z, y), for a
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationMath 220B - Summer 2003 Homework 1 Solutions
Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary
More informationGeneralized Bell polynomials and the combinatorics of Poisson central moments
Generaized Be poynomias and the combinatorics of Poisson centra moments Nicoas Privaut Division of Mathematica Sciences Schoo of Physica and Mathematica Sciences Nanyang Technoogica University SPMS-MAS-05-43,
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationFourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form
Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationarxiv: v1 [math.fa] 23 Aug 2018
An Exact Upper Bound on the L p Lebesgue Constant and The -Rényi Entropy Power Inequaity for Integer Vaued Random Variabes arxiv:808.0773v [math.fa] 3 Aug 08 Peng Xu, Mokshay Madiman, James Mebourne Abstract
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationHomework 5 Solutions
Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity
More informationDo Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix
VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix
More informationOn the Goal Value of a Boolean Function
On the Goa Vaue of a Booean Function Eric Bach Dept. of CS University of Wisconsin 1210 W. Dayton St. Madison, WI 53706 Lisa Heerstein Dept of CSE NYU Schoo of Engineering 2 Metrotech Center, 10th Foor
More informationAnalysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case
Technica Report PC-04-00 Anaysis of Emerson s Mutipe Mode Interpoation Estimation Agorithms: The MIMO Case João P. Hespanha Dae E. Seborg University of Caifornia, Santa Barbara February 0, 004 Anaysis
More informationGeneral Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping
Journa of Mathematica Research with Appications Jan.,, Vo. 3, No., pp. 53 6 DOI:.377/j.issn:95-65...7 Http://jmre.dut.edu.cn Genera Decay of Soutions in a Viscoeastic Equation with Noninear Locaized Damping
More informationPattern Frequency Sequences and Internal Zeros
Advances in Appied Mathematics 28, 395 420 (2002 doi:10.1006/aama.2001.0789, avaiabe onine at http://www.ideaibrary.com on Pattern Frequency Sequences and Interna Zeros Mikós Bóna Department of Mathematics,
More information1 Heat Equation Dirichlet Boundary Conditions
Chapter 3 Heat Exampes in Rectanges Heat Equation Dirichet Boundary Conditions u t (x, t) = ku xx (x, t), < x (.) u(, t) =, u(, t) = u(x, ) = f(x). Separate Variabes Look for simpe soutions in the
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationThe distribution of the number of nodes in the relative interior of the typical I-segment in homogeneous planar anisotropic STIT Tessellations
Comment.Math.Univ.Caroin. 51,3(21) 53 512 53 The distribution of the number of nodes in the reative interior of the typica I-segment in homogeneous panar anisotropic STIT Tesseations Christoph Thäe Abstract.
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationNearly Optimal Constructions of PIR and Batch Codes
arxiv:700706v [csit] 5 Jun 07 Neary Optima Constructions of PIR and Batch Codes Hia Asi Technion - Israe Institute of Technoogy Haifa 3000, Israe shea@cstechnionaci Abstract In this work we study two famiies
More informationTheory of Generalized k-difference Operator and Its Application in Number Theory
Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
SIAM J. NUMER. ANAL. Vo. 0, No. 0, pp. 000 000 c 200X Society for Industria and Appied Mathematics VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW
More informationInput-to-state stability for a class of Lurie systems
Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationStat 155 Game theory, Yuval Peres Fall Lectures 4,5,6
Stat 155 Game theory, Yuva Peres Fa 2004 Lectures 4,5,6 In the ast ecture, we defined N and P positions for a combinatoria game. We wi now show more formay that each starting position in a combinatoria
More informationPSEUDO-SPLINES, WAVELETS AND FRAMELETS
PSEUDO-SPLINES, WAVELETS AND FRAMELETS BIN DONG AND ZUOWEI SHEN Abstract The first type of pseudo-spines were introduced in [1, ] to construct tight frameets with desired approximation orders via the unitary
More informationOn Some Basic Properties of Geometric Real Sequences
On Some Basic Properties of eometric Rea Sequences Khirod Boruah Research Schoar, Department of Mathematics, Rajiv andhi University Rono His, Doimukh-791112, Arunacha Pradesh, India Abstract Objective
More informationAsymptotic Properties of a Generalized Cross Entropy Optimization Algorithm
1 Asymptotic Properties of a Generaized Cross Entropy Optimization Agorithm Zijun Wu, Michae Koonko, Institute for Appied Stochastics and Operations Research, Caustha Technica University Abstract The discrete
More informationarxiv: v1 [math.ap] 8 Nov 2014
arxiv:1411.116v1 [math.ap] 8 Nov 014 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS FOR m-component REACTION-DIFFUSION SYSTEMS WITH A TRIDIAGONAL SYMMETRIC TOEPLITZ MATRIX OF DIFFUSION COEFFICIENTS SALEM ABDELMALEK
More informationAnother Class of Admissible Perturbations of Special Expressions
Int. Journa of Math. Anaysis, Vo. 8, 014, no. 1, 1-8 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.014.31187 Another Cass of Admissibe Perturbations of Specia Expressions Jerico B. Bacani
More informationFourier series. Part - A
Fourier series Part - A 1.Define Dirichet s conditions Ans: A function defined in c x c + can be expanded as an infinite trigonometric series of the form a + a n cos nx n 1 + b n sin nx, provided i) f
More informationConvergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems
Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,
More informationJENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF SEVERAL VARIABLES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Voume 128, Number 7, Pages 2075 2084 S 0002-99390005371-5 Artice eectronicay pubished on February 16, 2000 JENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF
More informationApproximated MLC shape matrix decomposition with interleaf collision constraint
Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationEquilibrium of Heterogeneous Congestion Control Protocols
Equiibrium of Heterogeneous Congestion Contro Protocos Ao Tang Jiantao Wang Steven H. Low EAS Division, Caifornia Institute of Technoogy Mung Chiang EE Department, Princeton University Abstract When heterogeneous
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationAkaike Information Criterion for ANOVA Model with a Simple Order Restriction
Akaike Information Criterion for ANOVA Mode with a Simpe Order Restriction Yu Inatsu * Department of Mathematics, Graduate Schoo of Science, Hiroshima University ABSTRACT In this paper, we consider Akaike
More information2M2. Fourier Series Prof Bill Lionheart
M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier
More informationControl with Minimum Communication Cost per Symbol
Contro with Minimum Communication Cost per ymbo Justin Pearson, João P. Hespanha, Danie Liberzon Abstract We address the probem of stabiizing a continuous-time inear time-invariant process under communication
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized
More informationLemma 1. Suppose K S is a compact subset and I α is a covering of K. There is a finite subcollection {I j } such that
2 Singuar Integras We start with a very usefu covering emma. Lemma. Suppose K S is a compact subset and I α is a covering of K. There is a finite subcoection {I j } such that. {I j } are disjoint. 2. The
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationA. Distribution of the test statistic
A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch
More informationNEW FIBONACCI-LIKE WILD ATTRACTORS FOR UNIMODAL INTERVAL MAPS ZHANG RONG. (B.Sc., Nanjing University, China)
NEW FIBONACCI-LIKE WILD ATTRACTORS FOR UNIMODAL INTERVAL MAPS ZHANG RONG B.Sc., Nanjing University, China) A THESIS SUBMITTED FOR THE DEGREE OF DOCTOR OF PHILOSOPHY DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY
More informationCoded Caching for Files with Distinct File Sizes
Coded Caching for Fies with Distinct Fie Sizes Jinbei Zhang iaojun Lin Chih-Chun Wang inbing Wang Department of Eectronic Engineering Shanghai Jiao ong University China Schoo of Eectrica and Computer Engineering
More informationDavid Eigen. MA112 Final Paper. May 10, 2002
David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationT.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA
ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network
More informationValidated Numerical Bounds on the Global Error for Initial Value Problems for Stiff Ordinary Differential Equations. Chao Yu
Vaidated Numerica Bounds on the Goba Error for Initia Vaue Probems for Stiff Ordinary Differentia Equations by Chao Yu A thesis submitted in conformity with the requirements for the degree of Master of
More information(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].
PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform
More informationAbsolute Value Preconditioning for Symmetric Indefinite Linear Systems
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce
More informationInvestigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l
Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA
More informationImproving the Reliability of a Series-Parallel System Using Modified Weibull Distribution
Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution
More informationHAMILTON DECOMPOSITIONS OF ONE-ENDED CAYLEY GRAPHS
HAMILTON DECOMPOSITIONS OF ONE-ENDED CAYLEY GRAPHS JOSHUA ERDE, FLORIAN LEHNER, AND MAX PITZ Abstract. We prove that any one-ended, ocay finite Cayey graph with non-torsion generators admits a decomposition
More informationRestricted weak type on maximal linear and multilinear integral maps.
Restricted weak type on maxima inear and mutiinear integra maps. Oscar Basco Abstract It is shown that mutiinear operators of the form T (f 1,..., f k )(x) = R K(x, y n 1,..., y k )f 1 (y 1 )...f k (y
More informationKing Fahd University of Petroleum & Minerals
King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161
More informationConvergence of the spectral measure of non normal matrices
Convergence of the spectra measure of non norma matrices Aice Guionnet Phiip Wood Ofer Zeitouni October, 2 Abstract We discuss reguarization by noise of the spectrum of arge random non- Norma matrices.
More informationRate-Distortion Theory of Finite Point Processes
Rate-Distortion Theory of Finite Point Processes Günther Koiander, Dominic Schuhmacher, and Franz Hawatsch, Feow, IEEE Abstract We study the compression of data in the case where the usefu information
More informationUniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete
Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity
More informationMIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI
MIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI KLAUS SCHMIDT AND TOM WARD Abstract. We prove that every mixing Z d -action by automorphisms of a compact, connected, abeian group is
More informationPricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications
Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA
More informationLECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL
More informationSymbolic models for nonlinear control systems using approximate bisimulation
Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationThroughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay
Throughput Optima Scheduing for Wireess Downinks with Reconfiguration Deay Vineeth Baa Sukumaran vineethbs@gmai.com Department of Avionics Indian Institute of Space Science and Technoogy. Abstract We consider
More information