RISK-SENSITIVE ASSET MANAGEMENT WITH GENERAL FACTOR MODELS. ( ) : dst 0 = r(y t )St 0 dt, S0 0 = s 0 0 i(i = 1,, m) : ds i t = S i t : { µ i (Y t )dt
|
|
- なぎさ ののした
- 4 years ago
- Views:
Transcription
1 RISK-SENSITIVE ASSET MANAGEMENT WITH GENERAL FACTOR MODELS. ( ) : dst = r(y t )St dt, S = s i(i =,, m) : ds i t = S i t : { µ i (Y t )dt + n+m k= dy t = g(y t )dt + λ(y t )dw t, σ i k(y t )dw k t }, S i = s i Y () = y R n W t = (W k t ) k=,,(n+m) n + m σ, λ, µ, g, r R m (n+m), R n (n+m), R m, R n, R- (A) λ, g, σ, µ, r, (A) x, η R n, ξ R m, µ ξ ξ σσ (x)ξ µ ξ, µ η η λλ (x)η µ η, µ, µ > ( ) (A3) r πt i i ( ), = (,, ) Xt π dx π t X π t = ( π t ) ds t S t + m i= π i t dst i, X π St i = (FRS) Γ T (γ) := sup π A T γ log E [(Xπ T ) γ ] γ (, ) (, ) A T { γ (, ) ( ) Nagai(3) γ (, ) ( ) () ( )HJB (.) ( HJB (.) sup ] sup ˇπ ) π R m [
2 () HJB (.) ˇπ (.)( ) (3) (.) Verification Theorem ( ˇπ ) (FRS) HJB (.) v t + tr(λ(y)λ(y) D v) + γ (Dv) λ(y)λ(y) Dv + g(y) Dv + r(y) [ + sup γ ] π σ(y)σ(y) π + π {µ(y) r(y) + γσ(y)λ(y) Dv} π R m v(t, y) = v t + tr(λ(y)λ(y) D v) + B(y) Dv + (Dv) Q(y)Dv + U(y) =, (.) v(t, y) = Q, B, U { Q(y) := γλ(y) I + γ } γ σ(y) (σ(y)σ(y) ) σ(y) λ(y), B(y) := g(y) + γ γ λ(y)σ(y) (σ(y)σ(y) ) (µ(y) r(y)), U(y) := ( γ) (µ(y) r(y)) (σ(y)σ(y) ) (µ(y) r(y)) + r(y) Theorem.. (A) (A3) (A4) Dξ (y) C ( y + ), C > tr(λ(y)λ(y) D ξ ) + B(y) Dξ + (Dξ ) Q(y)Dξ + U(y), y ξ (A5) (µ(y) r(y)) (σ(y)σ(y) ) (µ(y) r(y)) + r(y), y. γ D := {γ (, ) : (A4) is satisfied} (.) v v v, t, D i v, D ij v L p (, T ; L p loc (Rn )), < p <.. γ D π(t, Y t ) := γ (σ(y t)σ(y t ) ) {µ(y t ) r(y t ) + γσ(y t )λ(y t ) D v(t, Y t )} (FRS) Γ T (γ) = v(, y) [] H. Hata (5) Risk-sensitive asset management with general factor models, preprint. [] H. Nagai (3) Optimal strategies for risk-sensitive portfolio optimization problems for general factor models, SIAM J. Cont. Optim. 4, 779 8, MR =,
3 Local risk-minimization for Barndorff-Nielsen and Shephard models ( ) stochastic volatility Barndorff-Nielsen and Shephard (BNS) locally risk-minimizing (LRM) BNS T > S t : { t ( ) t } S t = S exp µ + βσs ds + σ s dw s. S >, µ, β R, ρ W dσ t = λσ t dt + dh λt, σ >, volatility σ λ >, H subordinator( Lévy ) LRM X c R ξ T X = c + ξ t ds t c F ξ F selffinancing( ) self-financing L - t η t, ξ t V t = η t + ξ t S t t C t C t = V t t ξ sds s C = V C t F ( ) ( )
4 F = V T = C T + T ξ sds s C t C t η ξ φ = (η, ξ) C t (φ) R t (φ) R t (φ) := E [(C T (φ) C t (φ)) ] Ft φ R t (φ) R t ( φ) a.s. for every t [, T], risk-minimization LRM K > X = (S T K) + BNS LRM. β = ([3]). ρ = ([]) Föllmer-Schweizer Lévy Malliavin [4] Clark-Ocone ([5]) Malliavin LRM minimal martingale measure(mmm) equivalent martingale measure density Malliavin S t MMM Malliavin [] References [] Arai, T. (5) Local risk-minimization for Barndorff-Nielsen and Shephard models with volatility risk premium. Available at [] Arai, T., Imai, Y. and Suzuki, R. (5) Numerical analysis on local riskminimization for exponential Lévy models, to appear in International Journal of Theoretical and Applied Finance. [3] Arai, T., Imai, Y. and Suzuki, R. (5) Local risk minimization for Barndorff-Nielsen and Shephard models. Available at [4] Arai, T and Suzuki, R. (5) Local risk-minimization for Lévy markets, International Journal of Financial Engineering, Vol., 555. [5] Suzuki, R. (3) A Clark-Ocone type formula under change of measure for Lévy processes with L -Lévy measure, Commun. Stoch. Anal., Vol.7,
5 On the Euler-Maruyama approximation for one-dimensional SDEs with irregular coefficients Dai Taguchi (Ritsumeikan University) joint work with Hoang-Long Ngo (Hanoi National University of Education) In this talk, we consider the Euler-Maruyama approximation for one-dimensional stochastic differential equations. We provide the strong rate of convergence when the drift coefficient is the sum of a bounded variation function on compact sets and a Hölder continuous function, and the diffusion coefficient is a Hölder continuous function. Let X = (X t ) t T be a one-dimensional stochastic differential equation (SDE) X t = x + t b(x s )ds + t σ(x s )dw s, x R, t [, T ], () where W := (W t ) t T is a standard one-dimensional Brownian motion on a probability space (Ω, F, P). Veretennikov [6] show that if the drift coefficient b is bounded measurable and the diffusion coefficient σ is bounded, uniformly elliptic and / + α-hölder continuous for some α [, /], then the equation () has a unique strong solution. One often approximates X by using the Euler-Maruyama scheme which is defined by X (n) t = x + t ( ) b X (n) η n(s) ds + t ( ) σ X (n) η n(s) dw s, t [, T ], where η n (s) = kt/n if s [kt/n, (k + )T/n). It is well-known that if the coefficients b and σ are Lipschitz continuous functions, then the Euler-Maruyama scheme has strong rate of convergence /, that is for any p >, there exists C > such that E[ sup X t X (n) t p ] C t T n. p/ Recently, Yan [7] and Gyöngy and Rásonyi [] have shown that the strong rate in onedimensional setting with Hölder continuous diffusion coefficient. Ngo and Taguchi [4] extended the results in [, 7] for multi-dimensional SDEs with discontinuous drift coefficient (but it is onesided Lipschitz function). Halidias and Kloeden [] prove that if b is increasing, continuous from bellow and σ is a Lipschitz continuous, then X (n) converges to X in L -norm. Their proof is based on upper and lower solutions of the SDE and its Euler-Maruyama approximation, so it is difficult to get any
6 rate of convergence by using their approach. Leobacher and Szölgyenyi [3] introduce a clever way to transfer equation () with piecewise Lipschitz drift coefficient with finite number of discontinuous points, and Lipschitz continuous diffusion coefficient. Using their transformation technique, the equation () is equivalent to SDE with Lipschitz continuous coefficients. Therefore, the new equation can be approximated by its Euler-Maruyama scheme with the standard strong rate of convergence /. In this talk, we provide the rates of strong convergence of the Euler-Maruyama approximation for SDE () when the coefficients b and σ may have a very low regularity which includes every cases of [,, 3, 4, 7]. More preciously, we suppose that the drift coefficient b = b A + b H L (R) is bounded measurable, b A is a function of bounded variation on compact sets and b H is a Hölder continuous with β (, ], and the diffusion coefficient σ is bounded, uniformly elliptic and /+α- Hölder continuous with α [, /]. Under these assumptions for the coefficients b and σ, we obtain that the L -convergence rate for the Euler-Maruyama scheme is β α if α (, /] and log n if α = (Theorem.3 of [5]). For the case b / L (R), we also get the L -convergence rate, by using localization arguments, (Theorem.4 of [5]). Moerover, we obtain the L p -sup convergence rate for any p (Theorem.6,.7,.8 of [5]). The idea of proof is to use the Yamada-Watanabe approximation technique, the removal drift method and the Gaussian upper bounded for the density of the Euler-Maruyama approximation. References [] Gyöngy, I. and Rásonyi, M.: A note on Euler approximations for SDEs with Hölder continuous diffusion coefficients. Stochastic. Process. Appl., 89 (). [] Halidias, N. and Kloeden, P.E.: A note on the Euler-Maruyama scheme for stochastic differential equations with a discontinuous monotone drift coefficient. BIT 48() 5 59 (8). [3] Leobacher, G. and Szölgyenyi M.: A numerical method for SDEs with discontinuous drift. BIT Numer. Math. (5). [4] Ngo, H-L., and Taguchi, D.: Strong rate of convergence for the Euler-Maruyama approximation of stochastic differential equations with irregular coefficients. To appear in Mathematics of Computation. [5] Ngo, H-L., and Taguchi, D.: On the Euler-Maruyama approximation for one-dimensional stochastic differential equations with irregular coefficients. Preprint, arxiv: (5). [6] Veretennikov, A.Yu.: On strong solution and explicit formulas for solutions of stochastic integral equations. Math. USSR Sb. 39, (98). [7] Yan, B. L.: The Euler scheme with irregular coefficients. Ann. Probab. 3, no. 3, 7 94 ().
7 Martingale problems for diffusions on metric graphs TOMOYUKI ICHIBA We shall extend martingale problems for the Walsh Brownian motion ([], [3], [4]) on a star graph to those for diffusions on a metric graph (G, d) with a finite number of vertices (cf. []). Here the metric graph (G, d) is a collection of finite or semi-infinite edges {e E} with N vertices V := {v k, k =,..., N} for some N < + and metric d. We assume that the graph is imbedded in the Euclidian space R and that any two edges can only meet at a vertex. The metric d is defined in the canonical way as the length of a shortest path between two points in G along the edges and the length along each edge is measured with the usual Euclidian metric. The semi-infinite edges isomorphic to R + = [, ) are called external, while the finite edges isomorphic to a finite open interval are called internal. We assign the initial vertex and the terminal vertex for each edge e E via a map δ. The map δ associates each internal edge e to a pair (δ (e), δ (e)) V V with its initial vertex δ (e) and its terminal vertex δ (e). An external edge e has only the initial vertex δ (e) V and we define δ (e) = + formally. We call δ i (e), i =, the end points of e E. We write v e, if v V is one of the end points of e E, and define the set E(v) := {e E : v e} for each v V. We also assume that for each point x G \ V there exist unique e E, v V and r R +, such that x belongs to the edge e(x) with the initial vertex v(x), the terminal vertex δ (e(x)), and the length between x and the initial vertex is r(x) defined by v(x) := δ (e(x)), r(x) := d(v(x), x) ; x G. () Thus we identify the point x G in terms of such triplet (e, v, r) E V R +, and define each function g : G R with this coordinate, i.e., g(x) g(e, v, r), x G. We consider the class D of BOREL-measurable functions g : G R with the following properties: (i) they are continuous in the topology induced by the metric d ; (ii) for every e, v, the function r g e,v (r) := g(e, v, r) is twice continuously differentiable on (, ) and has finite first and second right-derivatives at the origin; (iii) the resulting functions (e, v, r) g e,v(r) and (e, v, r) g e,v(r) are BOREL measurable; and (iv) sup <r<k,e E(v),v V g e,v(r) + g e,v(r) < + holds for all finite K. For notational simplicity we write G (x) := g e,v(r), G (x) := g e,v(r) for x = (e, v, r) with r > () and for each vertex v V with e E(v) we define G (e, v, +) := lim x v,x E(v) ( )i+ G (x), if v = δ i (e), i =,. (3) Let us consider the canonical space Ω := C([, ); R ) of R -valued continuous functions on [, ) endowed with the σ -algebra F of its BOREL sets. We consider also its coordinate mapping ω( ) and the natural filtration F := {F(t), t < } with F(t) := σ(ω(s), s t), t <. For each vertex v V we shall consider a probability measure ν v (de), e E(v) on the set E(v) of edges. For example, if the total number of edges at v is finite, the measure ν v ( ) is a discrete probability measure (cf. []). In the following we shall construct the diffusion on G with drifts b := {b(x), x G}, dispersions σ := {σ(x), x G} and the vertex measures ν := {ν v ( ), v V} via the martingale problem associated with (b, σ, ν). Department of Statistics and Applied Probability, South Hall, University of California, Santa Barbara, CA 936
8 Given the vertex measures ν and BOREL measurable functions b : G R +, σ : G R \ {} we define for every g D the process M g ( ; ω) := g(ω( )) g(ω ()) Lg(ω(t)) {ω(t) E} dt, (4) where with a( ) := σ ( ) and the derivatives in () we define the infinitesimal generator Lg(x) := b(x)g (x) + a(x)g (x) ; x G. Local Martingale Problem: For every fixed x G to find a probability measure P on the canonical space (Ω, F), such that ω() = x holds P -a.e.; (ii) {ω(t) V} dt holds P - a.e.; (iii) for every function g D + (respectively, D ), the process M g ( ) in (4) is a continuous local submartingale (resp., martingale) with respect to the filtration F := {F(t+), t < }, where { } D + (resp., D ) := g D : G (e, v, +)ν v (de) (resp., ), v V. E(v) Proposition. Suppose that the drifts are identically zero and the reciprocal of the dispersion coefficient r σ e,v (r) := σ(e, v, r) = σ(x) is locally square integrable, i.e., K ( / σ e,v(r))dr < + for every compact subset K of e := e {δ (e)} {δ (e)}, e E(v), v V. Then the local martingale problem associated with the triplet (, σ, ν) is well-posed. More generally, under appropriate conditions, by the method of time-change and the scale functions we may construct the diffusion on the graph G associated with the triplet (b, σ, ν). For a fixed initial vertex v V and its neighborhood B (as a subset of E(v ) ), the local behavior of the resulting coordinate process X(, ω) := ω( ) in B is characterized by the distance R( ) := r(x( )) (defined as in ()) from the vertex v which satisfies dr(t) = b(x(t))dt + σ(x(t))dβ(t) + dl R (t) ; t (5) locally for some Brownian motion β( ) and the local time growths dl RA (t) = ν v (A) dl R (t) ; t, (6) where L ξ ( ) is the semimartingale local time at the origin for the generic semimartingale ξ( ) and R A ( ) := R( ) {e(x( )) A} for every BOREL subset A of E(v ). This is an extension of the results in [] from the graph with a finite number of edges to that with possibly uncountable number of edges but still with a finite number N of vertices. We shall also discuss some examples (including the WALSH diffusions with N = ) and further extensions. References [] BARLOW, M., PITMAN, J. & YOR, M. (989) On Walsh s Brownian motions. Séminaire de Probabilités XXIII. Lecture Notes in Mathematics 37, [] FREIDLIN, M. & SHEU, S. () Diffusion processes on graphs: stochastic differential equations, large deviation principle. Probab. Theory and Relat. Fields [3] ICHIBA, T., KARATZAS, I., PROKAJ, V. & YAN, M. (5) Stochastic integral equations for Walsh semimartingales. Preprint available at arxiv: [4] WALSH, J. (978) A diffusion with a discontinuous local time. Astérisque 5-53,
9 ( ), ( ), ( ), Poisson {X(t)} ( Surplus ) : (.) X(t) = x + κ t t dn(s) Z N(s), x., κ >, {N(t), t } γ Poisson, {Z k, k =,, } F i.i.d. r.v. s {N(t)}. {X(t)}, T. Lundberg model,., T X(T ) κ, x, γ. F, Laplace Fourier v(x, α, β) E x [ exp{ α T + i β X T } ], (α >, β R ) ( v ). δ (cf.,.3, D ),., D dense,, F. Surplus (.) Ito,.. ( ). α >, β R, v(x, α, β) = E x [exp{ α T + i β X(T )}] : F (λ) df (y) e λ y [ x α v(x, α, β) = κ x v(x, α, β) + γ df (y) { v(x y, α, β) e i β (x y)} + + e i β x F (i β) v(x, α, β) ], x, v(x, α, β) = e i β x, x <., x Laplace. ρ γ /κ, (.) v(λ, α, β) = v(, α, β) ρ { } F (i β) F (λ) /(λ i β) { }. λ ρ F (λ) α/κ, (.), v(, α, β).
10 {X(t)} random walk, Feller ([], Ch 8, Lemma ), v(, α, β) :.. F Fourier F (ξ) = E[e i ξ Z ], (.) log v(, α, β) = π h(ξ) iξ κ γ + γ F ( ξ). lim R α ds dx e iβx R R dξ e iξ x s h(ξ)..3. F δ, M >, < a < a <... < a M,, p k <, k =,,..., M (p + + p M = ), M (.3) df (x) = p k δ(x, a k ) dx.. δ(x, a) a δ. k= (.) (.),.,.4 ( ). F δ (.3), ρ γ κ, ζ ρ + α κ { v(, α, β) = exp k Z M + v(x, α, β) = E x [ exp{ α T + i β X(T ) } ] = v(, α, β) + M l= k Z M + k Z M +, v(x, α, β) p k k! ρ k e i β k,a k,a dy e (i β ζ) y y k }, p k k! ( ρ ) k ( x k, a ) k e ζ (x k,a ) I {x> k,a } p k k! ( ρ ) k + p l dy e ζ (y k,a )+i β (x y a l) (y k, a ) k I {y> k,a } I {<x y<al }. k Z M + ( ).. [] Cramér, H., Historical review of Fillip Lundberg s work on risk theory, Skand. Aktua. (Suppl.), 5 (969), 6-. [] Feller, W., An Introduction to Probability Theory and its Aplications, Vol. II, John Wiley & Sons, 97. [3] Lundberg, F., Über die Wahrscheinlichkeitsfunktion einer Risikenmasse. Skand. Aktua., 3 (93), -83. [4] Nishioka, K. and Igarashi, T., Non-ruin probability of an insurer under the Lundberg model,,, 93 (4), 4 47.
11 *,, 3 Ryo Oizumi*, Toshikazu Kuniya, Yoichi Enatsu,, 3. oizumi@ms.u-tokyo.ac.jp A R d a (, α) X a A X = x A Ito { dx j a (t) = g j (X a (t), v a, Γ t,x ) da + N k= σ jk (X a (t), v a, Γ t,x ) db k a, j d X j = xj () dx j a (t) := X j a+ε (t + ε) X j a (t) v = v a := (v (a), v (a),, v l (a))
12 [] t Γ t,x := ( ) Γ t,x, Γ t,x,, Γ m t,x,, Γ M t,x t a y A P t (a, x y) P t+ε (a + ε, x y) = A dξ K ε,t (ξ y) P t (a, x ξ) P t (, x y) = n t (x) δ d (x y) }{{} α. () n t (x) = dady F (y, Γ t,x ) P t (, x y) A }{{} t a ε Ito nonlocal Eq.(). Hamilton Jacobi Bellman a w {[ λ,a (x, Γ) inf Hv x (Γ) + λ ] w λ,a (x, Γ) } = v V w λ,α (x, Γ) = F (x, Γ) d H x v (Γ) := g j (x, v, Γ) x j d c jj (x, v, Γ) Γ R M +, j= ψ λ (x, Γ) = α j,j = x j x j + µ (x, v, Γ) (3) da w λ,a (x, Γ) =, (4) [] R. Oizumi, Unification theory of optimal life histories and linear demographic models in internal stochasticity, PLOS ONE 9 (6) (4) e98746.
13 (cf. [, ]) sharp Varadhan Z Z d Bernoulli Ginzburg-Landau [] C. Kipnis and C. Landim, Scaling Limits of Interacting Particle Systems, 999, Springer. [] S.R.S. Varadhan, Nonlinear diffusion limit for a system with nearest neighbor interactions II, in Asymptotic Problems in Probability Theory, Stochastic Models and Diffusions on Fractals, Pitman Res. Notes Math. Ser., 83 (994), 75-8., sasada@ms.u-tokyo.ac.jp
14 On fluctuations of eigenvalues of Gaussian beta ensembles at high temperature Trinh Khanh Duy Institute of Mathematics for Industry, Kyushu University Introduction Gaussian beta ensembles, as generalizations of Gaussian orthogonal ensembles, Gaussian unitary ensembles and Gaussian symplectic ensembles, were initially defined as eigenvalues ensembles with the joint density, (λ,..., λ n ) (λ) β e β ( n n )) i= V (λi) dλ = exp β ln λ j λ i V (λ i dλ, where V (λ) = λ /4 and (λ) = i<j (λ j λ i ) denotes the Vandermonde determinant. They can be also viewed as the equilibrium measure of a one dimensional Coulomb log-gas at the inverse temperature β. Let L n,β = n δ n λj / n j= be the empirical distributions of the scaled Gaussian beta ensembles. Then for fixed β, it is well known that the empirical measures L n,β converge weakly to the semicircle distribution, where the semicircle distribution, denoted by sc, is a probability measure supported on [, ] with the density 4 x /(π). This means that for any bounded continuous function f, L n,β, f = n f(λ j / n) sc, f = f(x) 4 x n π dx a.s. as n. j= The fluctuations around the semicircle distribution were also investigated. More precisely, it was shown in [7] that for a sufficiently nice function f with sc, f =, n j= i<j f(λ j / n) N d (, σf ) as n. d Here denotes the convergence in distribution. Note that beta ensembles with general potential V were considered in [7]. Matrix models for Gaussian beta ensembles were introduced by Dumitriu and Edelman [4]. They are symmetric tridiagonal matrices, also called Jacobi matrices, whose components are independent and are distributed as N (, ) χ (n )β χ (n )β N (, ) χ (n )β T n,β = β χ β N (, ) i=
15 Here χ k, for k >, denotes the (/ )-chi distribution with k degrees of freedom or the square root of the gamma distribution Γ(k, ). By those matrix models, Dumitriu and Edelman [5] established the above central limit theorem for polynomials f by different methods. The limiting behaviour of Gaussian beta ensembles in the regime where n with β = α/n (α being a positive constant) has been considered recently. It was shown first in [, 6] that the mean empirical distribution (also the mean spectral measure) converges weakly to a deterministic distribution µ α, and then in [3] that the empirical distributions L n,β themselves converge weakly to the same limit. Namely, for any polynomial p, as n with nβ = α, L n,β, p µ α, p in probability. The probability measure µ α here is the (scaled) measure of associated Hermite polynomial []. Note that the support µ α is the whole real line and the measure µ α is determined by its moments. It is the purpose of this talk to investigate the fluctuations of the empirical distributions around µ α. The main result is as follows. Theorem.. For a polynomial p of positive degree, as n with nβ = α, L n, p E[ L n, p ] n d N (, ˆσ p). The main tool used here is the martingale difference central limit theorem. It is expected that the central limit theorem should hold for a larger class of nice function f. References [] R. Allez, J-P. Bouchaud, and A. Guionnet: Invariant beta ensembles and the Gauss- Wigner crossover, Physical review letters 9 (), no. 9, 94. [] R. Askey and J. Wimp: Associated Laguerre and Hermite polynomials, Proc. Roy. Soc. Edinburgh Sect. A 96 (984), no. -, [3] F. Benaych-Georges, S. Péché: Poisson statistics for matrix ensembles at large temperature, J. Stat. Phys. 6 (5), no. 3, [4] I. Dumitriu and A. Edelman: Matrix models for beta ensembles, J. Math. Phys. 43 (), no., [5] I. Dumitriu and A. Edelman: Global spectrum fluctuations for the β-hermite and β-laguerre ensembles via matrix models, J. Math. Phys. 47 (6), 633. [6] T.K. Duy and T. Shirai: The mean spectral measures of random Jacobi matrices related to Gaussian beta ensembles, Electron. Commun. Probab. (5), no. 68, 3. [7] K. Johansson: On fluctuations of eigenvalues of random Hermitian matrices, Duke Math. J. 9 (998), no., 5 4.
16 CONCENTRATION FOR FIRST PASSAGE PERCOLATION. First passage percolation [] First passage percolation. ({η(j, x)} (j,x) N Z d, P); P(η(j, x) = ) = p (, ]. ω p (.) ω p = (k,x) N Z d η(k, x)δ (k,p /d x). (P, ω ) intensity counting measure Lebesgue measure N R d Poisson point process (.) p /d ω p p ω ω p, n { n } (.) T n (ω p ) = min x k x k α : x = and {(k, x k )} n k= ω p k= subadditive ergodic theorem (.3) µ p = lim n n T n(ω p ) P-a.s. µ p time constant 3. Theorem (concentration around the mean). For any χ < /, there exist C, C, λ > independent of p such that for any n N, P( n T n (ω p ) n E[T n (ω p )] > n χ ) < C exp{ C n λ }. Theorem (rate of convergence). For any γ (, ), there exists C > independent of p such that for any n N, µ p n E[T n (ω p )] < Cn γ. Theorem 3 (continuity). µ p is continuous in p [, ]. Theorem Concentration First Passage Percolation T n α > jump References [] F. Comets, R. Fukushima, S. Nakajima, N. Yoshida. Limiting results for the free energy of directed polymers in random environment with unbounded jumps Journal of Statistical Physics.
17 , [, ] Dirichlet Varadhan,. ([, Theorem.7]). (E, B, µ) σ-, (E, F) L (E, µ) Dirichlet, (T t ) t> L (E, µ) Markov. A, B, P t (A, B) = T t B dµ. lim t log P t(a, B) = d(a, B). t d(a, B) A, B, (E, F). A ( ). (E, B, µ) Dirichlet (E, F), Hilbert (H, (, ) H ) : F L (E H; µ) E (f, g) = ( f, g) H dµ E. E H- b, c, ϵ >, A ϵ > { } max b Hf dµ, c Hf dµ ϵe (f, f) + A ϵ f dµ, f F., E(f, g) = ( f, g) H dµ + E E E E (b, f) H g dµ + (c, g) H f dµ, E E f, g F (T t ) t>, (E, F) d Varadhan., A, B B, P t (A, B) = T t B dµ, A ( d(a, B) = sup essinf f F x B lim t log P t(a, B) = d(a, B) t {., F = f F b : h F b, ( f, f)h dµ. E f(x) esssup E x A ) f(x) } h dµ, essinf, esssup µ [] Ariyoshi, T. and Hino, M., Small-time asymptotic estimate in local Dirichlet spaces, Electron. J. Probab. (5), [] Hino, M. and Ramírez, J. A., Small-time Gaussian behavior of symmetric diffusion semigroups, Ann. Probab. 4 (3),
18 Parametrix method for simulation of stochastic differential equations Tomooki Yuasa Abstract In this talk, we introduce an exact simulation method for stochastic differential equations (SDE s) via Parametrix method with localization. This method can be expressed as two methods: the Forward method and the Backward method. In particular, the Backward method does not require the regularity of the coefficients of the SDE. Although the method has no bias, the error will be big because the variance is not generally finite. Accordingly, we consider an important sampling method together with and a new idea called the localization technique. Preparation Here, we simply explain the -dimensional SDE with the Forward method only. The general case can be obtained similarly. Let X t X t (x) be the weak solution of -dimensional SDE { dxt = σ(x t )dw t + b(x t )dt, t [, T ], () X = x, where W t is a -dimensional Wiener process. Let ξ i, i Z > be a sequence of i.i.d. random variables with common density f ξ which satisfies f ξ (,) (x) = dx a.e. and f ξ [, ) (x) > dx a.e. Also, let τ, τ i = i j= ξ j, i Z > be a point process and let N t = inf{n Z ; τ n t < τ n+ }, t [, T ] be its associated renewal process. Xτ π i, i Z denotes the Euler-Maruyama scheme of () with X π = x and the random partition π = { = τ < τ <... < τ NT T }. We define q n, n Z, for every n Z > and s = < s < < s n < T, q n (s, s,..., s n ) = T s n dxf ξ (x) n i= f ξ(s i+ s i ) and q (s ). Basic form of Parametrix method Theorem. (P. Andersson, V. Bally and A. Kohatsu-Higa (5)). Suppose that the coefficients of () satisfy σ Cb (R), b C b (R) and a := σ is uniformly elliptic. Then, for every f Cc (R), [ ] f(xt π E[f(X T )] = E ) N T θ τi+ τ q NT (τ, τ,..., τ NT ) i (Xτ π i, Xτ π i+ ). i= Dept. of Mathematics, Ritsumeikan University, -- Nojihigashi, Kusatsu, Shiga , Japan, ra9ff@ed.ritsumei.ac.jp
19 where for every t [, T ] and x, y R, θ t (x, y) = ( ) ( ) a (y) + (a(y) a(x))h () t (x, y) + a (y)h () t (x, y) b (y) + (b(y) b(x))h () t (x, y), h () t (x, y) = y x b(x)t a(x)t and h () t (x, y) = ( y x b(x)t a(x)t ) a(x)t. Also, if there exists ω such that N T (ω) =, then we understand N T (ω) i= θ τi+ τ i (Xτ π i (ω), Xτ π i+ (ω)) =. 3 Samplings of the sequence ξ i Once one defines f ξ, we can obtain the concrete Parametrix method. 3. Exponential sampling We arbitrarily choose λ >. Let ξ i, i Z > be a sequence of i.i.d. random variables with common density f ξ (x) = λe λx [, ) (x), x R. Namely, the common distribution of ξ i is the Exponential distribution with parameter λ and N t is a Poisson process. Also, for every n Z > and = s < s < < s n < T, q n (s, s,..., s n ) = e λt λ N T. 3. Beta sampling We arbitrarily choose τ > T and γ (, ). Let ξ i, i Z > be a sequence of i.i.d. random variables with common density f ξ (x) = γ x γ τ γ (,τ) (x), x R. Then, the common distribution F ξ (x) = ( x τ ) γ (,τ) (x) + [τ, ) (x), x R. Also, for every n Z > and = s < s < < s n < T, q n (s, s,..., s n ) = ( ( T sn τ ) τ )( γ ) n n τ γ i= (s i+ s i). γ 3.3 Gamma sampling We arbitrarily choose γ ( α, ) and ϑ >. Now, we need to choose α (, ) in Forward method, too. Let ξ i, i Z > be a sequence of i.i.d. random variables with common density f ξ (x) = Γ( γ)ϑ γ x e x γ ϑ (, ) (x), x R. Namely, the common distribution of ξ i is the Gamma distribution with parameter γ and ϑ. Also, for every n Z > and = s < s < < s n < T, q n (s, s,..., s n ) = Γ( γ, T s n ϑ ) Γ( γ) ( Γ( γ)ϑ ) n e sn n γ ϑ i= (s i+ s i ). γ Here, Beta sampling and Gamma sampling are two example of importance samplings which lead to methods with finite variance. Note that we must choose parameters for smaller variance and shorter trial time. However, the variance is still big even if we use importance sampling and its application although successful has some limitations. Thus, we suggest a new localization technique and give their theoretical properties in this talk. This technique can theoretically improve the importance sampling methods proposed above. In some sense, it becomes a space importance sampling method while the ones described above are time importance sampling methods.
untitled
3 3. (stochastic differential equations) { dx(t) =f(t, X)dt + G(t, X)dW (t), t [,T], (3.) X( )=X X(t) : [,T] R d (d ) f(t, X) : [,T] R d R d (drift term) G(t, X) : [,T] R d R d m (diffusion term) W (t)
More information43433 8 3 . Stochastic exponentials...................................... 3. Girsanov s theorem......................................... 4 On the martingale property of stochastic exponentials 5. Gronwall
More information2 (March 13, 2010) N Λ a = i,j=1 x i ( d (a) i,j x j ), Λ h = N i,j=1 x i ( d (h) i,j x j ) B a B h B a = N i,j=1 ν i d (a) i,j, B h = x j N i,j=1 ν i
1. A. M. Turing [18] 60 Turing A. Gierer H. Meinhardt [1] : (GM) ) a t = D a a xx µa + ρ (c a2 h + ρ 0 (0 < x < l, t > 0) h t = D h h xx νh + c ρ a 2 (0 < x < l, t > 0) a x = h x = 0 (x = 0, l) a = a(x,
More information第5章 偏微分方程式の境界値問題
October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ
More informationuntitled
2 : n =1, 2,, 10000 0.5125 0.51 0.5075 0.505 0.5025 0.5 0.4975 0.495 0 2000 4000 6000 8000 10000 2 weak law of large numbers 1. X 1,X 2,,X n 2. µ = E(X i ),i=1, 2,,n 3. σi 2 = V (X i ) σ 2,i=1, 2,,n ɛ>0
More information25 II :30 16:00 (1),. Do not open this problem booklet until the start of the examination is announced. (2) 3.. Answer the following 3 proble
25 II 25 2 6 13:30 16:00 (1),. Do not open this problem boolet until the start of the examination is announced. (2) 3.. Answer the following 3 problems. Use the designated answer sheet for each problem.
More informationDesign of highly accurate formulas for numerical integration in weighted Hardy spaces with the aid of potential theory 1 Ken ichiro Tanaka 1 Ω R m F I = F (t) dt (1.1) Ω m m 1 m = 1 1 Newton-Cotes Gauss
More information浜松医科大学紀要
On the Statistical Bias Found in the Horse Racing Data (1) Akio NODA Mathematics Abstract: The purpose of the present paper is to report what type of statistical bias the author has found in the horse
More informationADM-Hamiltonian Cheeger-Gromov 3. Penrose
ADM-Hamiltonian 1. 2. Cheeger-Gromov 3. Penrose 0. ADM-Hamiltonian (M 4, h) Einstein-Hilbert M 4 R h hdx L h = R h h δl h = 0 (Ric h ) αβ 1 2 R hg αβ = 0 (Σ 3, g ij ) (M 4, h ij ) g ij, k ij Σ π ij = g(k
More information1. Introduction Palatini formalism vierbein e a µ spin connection ω ab µ Lgrav = e (R + Λ). 16πG R µνab µ ω νab ν ω µab ω µac ω νcb + ω νac ω µcb, e =
Chiral Fermion in AdS(dS) Gravity Fermions in (Anti) de Sitter Gravity in Four Dimensions, N.I, Takeshi Fukuyama, arxiv:0904.1936. Prog. Theor. Phys. 122 (2009) 339-353. 1. Introduction Palatini formalism
More information80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = i=1 i=1 n λ x i e λ i=1 x i! = λ n i=1 x i e nλ n i=1 x
80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = n λ x i e λ x i! = λ n x i e nλ n x i! n n log l(λ) = log(λ) x i nλ log( x i!) log l(λ) λ = 1 λ n x i n =
More information,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,.
9 α ν β Ξ ξ Γ γ o δ Π π ε ρ ζ Σ σ η τ Θ θ Υ υ ι Φ φ κ χ Λ λ Ψ ψ µ Ω ω Def, Prop, Th, Lem, Note, Remark, Ex,, Proof, R, N, Q, C [a, b {x R : a x b} : a, b {x R : a < x < b} : [a, b {x R : a x < b} : a,
More informationT rank A max{rank Q[R Q, J] t-rank T [R T, C \ J] J C} 2 ([1, p.138, Theorem 4.2.5]) A = ( ) Q rank A = min{ρ(j) γ(j) J J C} C, (5) ρ(j) = rank Q[R Q,
(ver. 4:. 2005-07-27) 1 1.1 (mixed matrix) (layered mixed matrix, LM-matrix) m n A = Q T (2m) (m n) ( ) ( ) Q I m Q à = = (1) T diag [t 1,, t m ] T rank à = m rank A (2) 1.2 [ ] B rank [B C] rank B rank
More informationJOURNAL OF THE JAPANESE ASSOCIATION FOR PETROLEUM TECHNOLOGY VOL. 66, NO. 6 (Nov., 2001) (Received August 10, 2001; accepted November 9, 2001) Alterna
JOURNAL OF THE JAPANESE ASSOCIATION FOR PETROLEUM TECHNOLOGY VOL. 66, NO. 6 (Nov., 2001) (Received August 10, 2001; accepted November 9, 2001) Alternative approach using the Monte Carlo simulation to evaluate
More informationII Brown Brown
II 16 12 5 1 Brown 3 1.1..................................... 3 1.2 Brown............................... 5 1.3................................... 8 1.4 Markov.................................... 1 1.5
More informationtokei01.dvi
2. :,,,. :.... Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 3. (probability),, 1. : : n, α A, A a/n. :, p, p Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN
More informationTrapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x
University of Hyogo 8 8 1 d x(t) =f(t, x(t)), dt (1) x(t 0 ) =x 0 () t n = t 0 + n t x x n n x n x 0 x i i = 0,..., n 1 x n x(t) 1 1.1 1 1 1 0 θ 1 θ x n x n 1 t = θf(t n 1, x n 1 ) + (1 θ)f(t n, x n )
More information149 (Newell [5]) Newell [5], [1], [1], [11] Li,Ryu, and Song [2], [11] Li,Ryu, and Song [2], [1] 1) 2) ( ) ( ) 3) T : 2 a : 3 a 1 :
Transactions of the Operations Research Society of Japan Vol. 58, 215, pp. 148 165 c ( 215 1 2 ; 215 9 3 ) 1) 2) :,,,,, 1. [9] 3 12 Darroch,Newell, and Morris [1] Mcneil [3] Miller [4] Newell [5, 6], [1]
More informationI
I 6 4 10 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............
More information数学の基礎訓練I
I 9 6 13 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 3 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............
More informationh23w1.dvi
24 I 24 2 8 10:00 12:30 1),. Do not open this problem booklet until the start of the examination is announced. 2) 3.. Answer the following 3 problems. Use the designated answer sheet for each problem.
More informationC 1 -path x t x 1 (f(x u), dx u ) rough path analyi p-variation (1 < p < 2) rough path 2 Introduction f(x) = (fj i(x)) 1 i n,1 j d (x R d ) (n, d) Cb
Rough path analyi 1 x t ( t 1) R d path f(x) = t (f 1 (x),, f d (x)) R d R d - C x t 1 (f(x u), dx u ) Stieltje path x t p-variation norm (1 < p < 2) x p := { up D } 1/p N x ti x ti 1 p (D = { = t <
More informationxia2.dvi
Journal of Differential Equations 96 (992), 70-84 Melnikov method and transversal homoclinic points in the restricted three-body problem Zhihong Xia Department of Mathematics, Harvard University Cambridge,
More informationJFE.dvi
,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho
More informationNumerical Analysis II, Exam End Term Spring 2017
H. Ammari W. Wu S. Yu Spring Term 2017 Numerical Analysis II ETH Zürich D-MATH End Term Spring 2017 Problem 1 Consider dx = f(t, x), t [0, T ] dt x(0) = x 0 R [28 Marks] with f C subject to the Lipschitz
More information,,,17,,, ( ),, E Q [S T F t ] < S t, t [, T ],,,,,,,,
14 5 1 ,,,17,,,194 1 4 ( ),, E Q [S T F t ] < S t, t [, T ],,,,,,,, 1 4 1.1........................................ 4 5.1........................................ 5.........................................
More information1. R n Ω ε G ε 0 Ω ε B n 2 Ωε = with Bu = 0 on Ω ε i=1 x 2 i ε +0 B Bu = u (Dirichlet, D Ω ε ), Bu = u ν (Neumann, N Ω ε ), Ω ε G ( ) / 25
.. IV 2012 10 4 ( ) 2012 10 4 1 / 25 1. R n Ω ε G ε 0 Ω ε B n 2 Ωε = with Bu = 0 on Ω ε i=1 x 2 i ε +0 B Bu = u (Dirichlet, D Ω ε ), Bu = u ν (Neumann, N Ω ε ), Ω ε G ( ) 2012 10 4 2 / 25 1. Ω ε B ε t
More informationI (Analysis I) Lebesgue (Lebesgue Integral Theory) 1 (Seiji HIRABA) 1 ( ),,, ( )
I (Analysis I) Lebesgue (Lebesgue Integral Theory) 1 (Seiji HIRABA) 1 ( ),,, ( ) 1 (Introduction) 1 1.1... 1 1.2 Riemann Lebesgue... 2 2 (Measurable sets and Measures) 4 2.1 σ-... 4 2.2 Borel... 5 2.3...
More information1 # include < stdio.h> 2 # include < string.h> 3 4 int main (){ 5 char str [222]; 6 scanf ("%s", str ); 7 int n= strlen ( str ); 8 for ( int i=n -2; i
ABC066 / ARC077 writer: nuip 2017 7 1 For International Readers: English editorial starts from page 8. A : ringring a + b b + c a + c a, b, c a + b + c 1 # include < stdio.h> 2 3 int main (){ 4 int a,
More informationL. S. Abstract. Date: last revised on 9 Feb translated to Japanese by Kazumoto Iguchi. Original papers: Received May 13, L. Onsager and S.
L. S. Abstract. Date: last revised on 9 Feb 01. translated to Japanese by Kazumoto Iguchi. Original papers: Received May 13, 1953. L. Onsager and S. Machlup, Fluctuations and Irreversibel Processes, Physical
More informationPage 1 of 6 B (The World of Mathematics) November 20, 2006 Final Exam 2006 Division: ID#: Name: 1. p, q, r (Let p, q, r are propositions. ) (10pts) (a
Page 1 of 6 B (The World of Mathematics) November 0, 006 Final Exam 006 Division: ID#: Name: 1. p, q, r (Let p, q, r are propositions. ) (a) (Decide whether the following holds by completing the truth
More information空間多次元 Navier-Stokes 方程式に対する無反射境界条件
81 Navier-Stokes Poinsot Lele Poinsot Lele Thompson Euler Navier-Stokes A Characteristic Nonreflecting Boundary Condition for the Multidimensional Navier-Stokes Equations Takaharu YAGUCHI, Kokichi SUGIHARA
More informationTest IV, March 22, 2016 6. Suppose that 2 n a n converges. Prove or disprove that a n converges. Proof. Method I: Let a n x n be a power series, which converges at x = 2 by the assumption. Applying Theorem
More information2
2011 8 6 2011 5 7 [1] 1 2 i ii iii i 3 [2] 4 5 ii 6 7 iii 8 [3] 9 10 11 cf. Abstracts in English In terms of democracy, the patience and the kindness Tohoku people have shown will be dealt with as an exception.
More information3. ( 1 ) Linear Congruential Generator:LCG 6) (Mersenne Twister:MT ), L 1 ( 2 ) 4 4 G (i,j) < G > < G 2 > < G > 2 g (ij) i= L j= N
RMT 1 1 1 N L Q=L/N (RMT), RMT,,,., Box-Muller, 3.,. Testing Randomness by Means of RMT Formula Xin Yang, 1 Ryota Itoi 1 and Mieko Tanaka-Yamawaki 1 Random matrix theory derives, at the limit of both dimension
More informationBasic Math. 1 0 [ N Z Q Q c R C] 1, 2, 3,... natural numbers, N Def.(Definition) N (1) 1 N, (2) n N = n +1 N, (3) N (1), (2), n N n N (element). n/ N.
Basic Mathematics 16 4 16 3-4 (10:40-12:10) 0 1 1 2 2 2 3 (mapping) 5 4 ε-δ (ε-δ Logic) 6 5 (Potency) 9 6 (Equivalence Relation and Order) 13 7 Zorn (Axiom of Choice, Zorn s Lemma) 14 8 (Set and Topology)
More informationIsogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206,
H28. (TMU) 206 8 29 / 34 2 3 4 5 6 Isogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206, http://link.springer.com/article/0.007/s409-06-0008-x
More informationver.1 / c /(13)
1 -- 11 1 c 2010 1/(13) 1 -- 11 -- 1 1--1 1--1--1 2009 3 t R x R n 1 ẋ = f(t, x) f = ( f 1,, f n ) f x(t) = ϕ(x 0, t) x(0) = x 0 n f f t 1--1--2 2009 3 q = (q 1,..., q m ), p = (p 1,..., p m ) x = (q,
More information24 Depth scaling of binocular stereopsis by observer s own movements
24 Depth scaling of binocular stereopsis by observer s own movements 1130313 2013 3 1 3D 3D 3D 2 2 i Abstract Depth scaling of binocular stereopsis by observer s own movements It will become more usual
More information(check matrices and minimum distances) H : a check matrix of C the minimum distance d = (the minimum # of column vectors of H which are linearly depen
Hamming (Hamming codes) c 1 # of the lines in F q c through the origin n = qc 1 q 1 Choose a direction vector h i for each line. No two vectors are colinear. A linearly dependent system of h i s consists
More information(Basics of Proability Theory). (Probability Spacees ad Radom Variables,, (Ω, F, P ),, X,. (Ω, F, P ) (probability space) Ω ( ω Ω ) F ( 2 Ω ) Ω σ (σ-fi
I (Basics of Probability Theory ad Radom Walks) 25 4 5 ( 4 ) (Preface),.,,,.,,,...,,.,.,,.,,. (,.) (Basics of Proability Theory). (Probability Spacees ad Radom Variables...............2, (Expectatios,
More information1 Fourier Fourier Fourier Fourier Fourier Fourier Fourier Fourier Fourier analog digital Fourier Fourier Fourier Fourier Fourier Fourier Green Fourier
Fourier Fourier Fourier etc * 1 Fourier Fourier Fourier (DFT Fourier (FFT Heat Equation, Fourier Series, Fourier Transform, Discrete Fourier Transform, etc Yoshifumi TAKEDA 1 Abstract Suppose that u is
More information4/15 No.
4/15 No. 1 4/15 No. 4/15 No. 3 Particle of mass m moving in a potential V(r) V(r) m i ψ t = m ψ(r,t)+v(r)ψ(r,t) ψ(r,t) = ϕ(r)e iωt ψ(r,t) Wave function steady state m ϕ(r)+v(r)ϕ(r) = εϕ(r) Eigenvalue problem
More informationLebesgue Fubini L p Banach, Hilbert Höld
II (Analysis II) Lebesgue (Applications of Lebesgue Integral Theory) 1 (Seiji HIABA) 1 ( ),,, ( ) 1 1 1.1 1 Lebesgue........................ 1 1.2 2 Fubini...................... 2 2 L p 5 2.1 Banach, Hilbert..............................
More information1 4 1 ( ) ( ) ( ) ( ) () 1 4 2
7 1995, 2017 7 21 1 2 2 3 3 4 4 6 (1).................................... 6 (2)..................................... 6 (3) t................. 9 5 11 (1)......................................... 11 (2)
More informationAuthor Workshop 20111124 Henry Cavendish 1731-1810 Biot-Savart 26 (1) (2) (3) (4) (5) (6) Priority Proceeding Impact factor Full paper impact factor Peter Drucker 1890-1971 1903-1989 Title) Abstract
More informationFeynman Encounter with Mathematics 52, [1] N. Kumano-go, Feynman path integrals as analysis on path space by time slicing approximation. Bull
Feynman Encounter with Mathematics 52, 200 9 [] N. Kumano-go, Feynman path integrals as analysis on path space by time slicing approximation. Bull. Sci. Math. vol. 28 (2004) 97 25. [2] D. Fujiwara and
More information²ÄÀÑʬΥ»¶ÈóÀþ·¿¥·¥å¥ì¡¼¥Ç¥£¥ó¥¬¡¼ÊýÄø¼°¤ÎÁ²¶á²òÀÏ Asymptotic analysis for the integrable discrete nonlinear Schrödinger equation
Asymptotic analysis for the integrable discrete nonlinear Schrödinger equation ( ) ( ) 2016 12 17 1. Schrödinger focusing NLS iu t + u xx +2 u 2 u = 0 u(x, t) =2ηe 2iξx 4i(ξ2 η 2 )t+i(ψ 0 +π/2) sech(2ηx
More information(2 X Poisso P (λ ϕ X (t = E[e itx ] = k= itk λk e k! e λ = (e it λ k e λ = e eitλ e λ = e λ(eit 1. k! k= 6.7 X N(, 1 ϕ X (t = e 1 2 t2 : Cauchy ϕ X (t
6 6.1 6.1 (1 Z ( X = e Z, Y = Im Z ( Z = X + iy, i = 1 (2 Z E[ e Z ] < E[ Im Z ] < Z E[Z] = E[e Z] + ie[im Z] 6.2 Z E[Z] E[ Z ] : E[ Z ] < e Z Z, Im Z Z E[Z] α = E[Z], Z = Z Z 1 {Z } E[Z] = α = α [ α ]
More informationI, II 1, A = A 4 : 6 = max{ A, } A A 10 10%
1 2006.4.17. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 1. 1. 2. 3. 4. 5. 2. ɛ-δ 1. ɛ-n
More information確率論と統計学の資料
5 June 015 ii........................ 1 1 1.1...................... 1 1........................... 3 1.3... 4 6.1........................... 6................... 7 ii ii.3.................. 8.4..........................
More informationcompact compact Hermann compact Hermite ( - ) Hermann Hermann ( ) compact Hermite Lagrange compact Hermite ( ) a, Σ a {0} a 3 1
014 5 4 compact compact Hermann compact Hermite ( - ) Hermann Hermann ( ) compact Hermite Lagrange compact Hermite ( ) 1 1.1. a, Σ a {0} a 3 1 (1) a = span(σ). () α, β Σ s α β := β α,β α α Σ. (3) α, β
More informationuntitled
1 25/5/3-6/3 1 1 1.1.................................. 1 1.2.................................. 4 2 5 2.1.............................. 5 2.2.............................. 6 3 Black Scholes 7 3.1 BS............................
More information2
p1 i 2 = 1 i 2 x, y x + iy 2 (x + iy) + (γ + iδ) = (x + γ) + i(y + δ) (x + iy)(γ + iδ) = (xγ yδ) + i(xδ + yγ) i 2 = 1 γ + iδ 0 x + iy γ + iδ xγ + yδ xδ = γ 2 + iyγ + δ2 γ 2 + δ 2 p7 = x 2 +y 2 z z p13
More information1 M = (M, g) m Riemann N = (N, h) n Riemann M N C f : M N f df : T M T N M T M f N T N M f 1 T N T M f 1 T N C X, Y Γ(T M) M C T M f 1 T N M Levi-Civi
1 Surveys in Geometry 1980 2 6, 7 Harmonic Map Plateau Eells-Sampson [5] Siu [19, 20] Kähler 6 Reports on Global Analysis [15] Sacks- Uhlenbeck [18] Siu-Yau [21] Frankel Siu Yau Frankel [13] 1 Surveys
More information微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.
微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. ttp://www.morikita.co.jp/books/mid/00571 このサンプルページの内容は, 初版 1 刷発行時のものです. i ii 014 10 iii [note] 1 3 iv 4 5 3 6 4 x 0 sin x x 1 5 6 z = f(x, y) 1 y = f(x)
More information201711grade1ouyou.pdf
2017 11 26 1 2 52 3 12 13 22 23 32 33 42 3 5 3 4 90 5 6 A 1 2 Web Web 3 4 1 2... 5 6 7 7 44 8 9 1 2 3 1 p p >2 2 A 1 2 0.6 0.4 0.52... (a) 0.6 0.4...... B 1 2 0.8-0.2 0.52..... (b) 0.6 0.52.... 1 A B 2
More informationDevelopement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Sum
Developement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Summary Previously, the authors developed the plastic
More information磁性物理学 - 遷移金属化合物磁性のスピンゆらぎ理論
email: takahash@sci.u-hyogo.ac.jp May 14, 2009 Outline 1. 2. 3. 4. 5. 6. 2 / 262 Today s Lecture: Mode-mode Coupling Theory 100 / 262 Part I Effects of Non-linear Mode-Mode Coupling Effects of Non-linear
More information03.Œk’ì
HRS KG NG-HRS NG-KG AIC Fama 1965 Mandelbrot Blattberg Gonedes t t Kariya, et. al. Nagahara ARCH EngleGARCH Bollerslev EGARCH Nelson GARCH Heynen, et. al. r n r n =σ n w n logσ n =α +βlogσ n 1 + v n w
More information.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T
NHK 204 2 0 203 2 24 ( ) 7 00 7 50 203 2 25 ( ) 7 00 7 50 203 2 26 ( ) 7 00 7 50 203 2 27 ( ) 7 00 7 50 I. ( ν R n 2 ) m 2 n m, R = e 2 8πε 0 hca B =.09737 0 7 m ( ν = ) λ a B = 4πε 0ħ 2 m e e 2 = 5.2977
More information基礎数学I
I & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............
More information318 T. SICE Vol.52 No.6 June 2016 (a) (b) (c) (a) (c) ) 11) (1) (2) 1 5) 6) 7), 8) 5) 20 11) (1
Vol.52, No.6, 317/329 2016 Mathematical Model and Analysis on Human s Handwork and Their Application to Brush Motion and Junkatsu of Japanese Calligraphy Tadamichi Mawatari and Kensuke Tsuchiya There is
More information(Basic of Proability Theory). (Probability Spacees ad Radom Variables , (Expectatios, Meas) (Weak Law
I (Radom Walks ad Percolatios) 3 4 7 ( -2 ) (Preface),.,,,...,,.,,,,.,.,,.,,. (,.) (Basic of Proability Theory). (Probability Spacees ad Radom Variables...............2, (Expectatios, Meas).............................
More informationII (Percolation) ( 3-4 ) 1. [ ],,,,,,,. 2. [ ],.. 3. [ ],. 4. [ ] [ ] G. Grimmett Percolation Springer-Verlag New-York [ ] 3
II (Percolation) 12 9 27 ( 3-4 ) 1 [ ] 2 [ ] 3 [ ] 4 [ ] 1992 5 [ ] G Grimmett Percolation Springer-Verlag New-York 1989 6 [ ] 3 1 3 p H 2 3 2 FKG BK Russo 2 p H = p T (=: p c ) 3 2 Kesten p c =1/2 ( )
More informationDecember 28, 2018
e-mail : kigami@i.kyoto-u.ac.jp December 28, 28 Contents 2............................. 3.2......................... 7.3..................... 9.4................ 4.5............. 2.6.... 22 2 36 2..........................
More information28 Horizontal angle correction using straight line detection in an equirectangular image
28 Horizontal angle correction using straight line detection in an equirectangular image 1170283 2017 3 1 2 i Abstract Horizontal angle correction using straight line detection in an equirectangular image
More information1 : ( ) ( ) ( ) ( ) ( ) etc (SCA)
START: 17th Symp. Auto. Decentr. Sys., Jan. 28, 2005 Symplectic cellular automata as a test-bed for research on the emergence of natural systems 1 : ( ) ( ) ( ) ( ) ( ) etc (SCA) 2 SCA 2.0 CA ( ) E.g.
More information通信容量制約を考慮したフィードバック制御 - 電子情報通信学会 情報理論研究会(IT) 若手研究者のための講演会
IT 1 2 1 2 27 11 24 15:20 16:05 ( ) 27 11 24 1 / 49 1 1940 Witsenhausen 2 3 ( ) 27 11 24 2 / 49 1940 2 gun director Warren Weaver, NDRC (National Defence Research Committee) Final report D-2 project #2,
More information医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.
医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009192 このサンプルページの内容は, 第 2 版 1 刷発行時のものです. i 2 t 1. 2. 3 2 3. 6 4. 7 5. n 2 ν 6. 2 7. 2003 ii 2 2013 10 iii 1987
More information鉄鋼協会プレゼン
NN :~:, 8 Nov., Adaptive H Control for Linear Slider with Friction Compensation positioning mechanism moving table stand manipulator Point to Point Control [G] Continuous Path Control ground Fig. Positoining
More information1..FEM FEM 3. 4.
008 stress behavior at the joint of stringer to cross beam of the steel railway bridge 1115117 1..FEM FEM 3. 4. ABSTRACT 1. BackgroundPurpose The occurrence of fatigue crack is reported in the joint of
More informationOn the Wireless Beam of Short Electric Waves. (VII) (A New Electric Wave Projector.) By S. UDA, Member (Tohoku Imperial University.) Abstract. A new e
On the Wireless Beam of Short Electric Waves. (VII) (A New Electric Wave Projector.) By S. UDA, Member (Tohoku Imperial University.) Abstract. A new electric wave projector is proposed in this paper. The
More information02-量子力学の復習
4/17 No. 1 4/17 No. 2 4/17 No. 3 Particle of mass m moving in a potential V(r) V(r) m i ψ t = 2 2m 2 ψ(r,t)+v(r)ψ(r,t) ψ(r,t) Wave function ψ(r,t) = ϕ(r)e iωt steady state 2 2m 2 ϕ(r)+v(r)ϕ(r) = εϕ(r)
More informationseminar0220a.dvi
1 Hi-Stat 2 16 2 20 16:30-18:00 2 2 217 1 COE 4 COE RA E-MAIL: ged0104@srv.cc.hit-u.ac.jp 2004 2 25 S-PLUS S-PLUS S-PLUS S-code 2 [8] [8] [8] 1 2 ARFIMA(p, d, q) FI(d) φ(l)(1 L) d x t = θ(l)ε t ({ε t }
More informationBlack-Scholes [1] Nelson [2] Schrödinger 1 Black Scholes [1] Black-Scholes Nelson [2][3][4] Schrödinger Nelson Parisi Wu [5] Nelson Parisi-W
003 7 14 Black-Scholes [1] Nelson [] Schrödinger 1 Black Scholes [1] Black-Scholes Nelson [][3][4] Schrödinger Nelson Parisi Wu [5] Nelson Parisi-Wu Nelson e-mail: takatoshi-tasaki@nifty.com kabutaro@mocha.freemail.ne.jp
More informationV(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H
199 1 1 199 1 1. Vx) m e V cos x π x π Vx) = x < π, x > π V i) x = Vx) V 1 x /)) n n d f dξ ξ d f dξ + n f = H n ξ) ii) H n ξ) = 1) n expξ ) dn dξ n exp ξ )) H n ξ)h m ξ) exp ξ )dξ = π n n!δ n,m x = Vx)
More information28
y i = Z i δ i +ε i ε i δ X y i = X Z i δ i + X ε i [ ] 1 δ ˆ i = Z i X( X X) 1 X Z i [ ] 1 σ ˆ 2 Z i X( X X) 1 X Z i Z i X( X X) 1 X y i σ ˆ 2 ˆ σ 2 = [ ] y i Z ˆ [ i δ i ] 1 y N p i Z i δ ˆ i i RSTAT
More informationAnderson ( ) Anderson / 14
Anderson 2008 12 ( ) Anderson 2008 12 1 / 14 Anderson ( ) Anderson 2008 12 2 / 14 Anderson P.W.Anderson 1958 ( ) Anderson 2008 12 3 / 14 Anderson tight binding Anderson tight binding Z d u (x) = V i u
More informationIII III 2010 PART I 1 Definition 1.1 (, σ-),,,, Borel( ),, (σ-) (M, F, µ), (R, B(R)), (C, B(C)) Borel Definition 1.2 (µ-a.e.), (in µ), (in L 1 (µ)). T
III III 2010 PART I 1 Definition 1.1 (, σ-),,,, Borel( ),, (σ-) (M, F, µ), (R, B(R)), (C, B(C)) Borel Definition 1.2 (µ-a.e.), (in µ), (in L 1 (µ)). Theorem 1.3 (Lebesgue ) lim n f n = f µ-a.e. g L 1 (µ)
More informationI A A441 : April 21, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) Google
I4 - : April, 4 Version :. Kwhir, Tomoki TA (Kondo, Hirotk) Google http://www.mth.ngoy-u.c.jp/~kwhir/courses/4s-biseki.html pdf 4 4 4 4 8 e 5 5 9 etc. 5 6 6 6 9 n etc. 6 6 6 3 6 3 7 7 etc 7 4 7 7 8 5 59
More informationGelfand 3 L 2 () ix M : ϕ(x) ixϕ(x) M : σ(m) = i (λ M) λ (L 2 () ) ( 0 ) L 2 () ϕ, ψ L 2 () ((λ M) ϕ, ψ) ((λ M) ϕ, ψ) = λ ix ϕ(x)ψ(x)dx. λ /(λ ix) ϕ,
A spectral theory of linear operators on Gelfand triplets MI (Institute of Mathematics for Industry, Kyushu University) (Hayato CHIBA) chiba@imi.kyushu-u.ac.jp Dec 2, 20 du dt = Tu. (.) u X T X X T 0 X
More information(1) (2) (3) (4) 1
8 3 4 3.................................... 3........................ 6.3 B [, ].......................... 8.4........................... 9........................................... 9.................................
More informationmf.dvi
21 9 29 1 2 3....................................... 3 :......................... 3....................................... 4................................ 4..................................... 5................................
More informationn ξ n,i, i = 1,, n S n ξ n,i n 0 R 1,.. σ 1 σ i .10.14.15 0 1 0 1 1 3.14 3.18 3.19 3.14 3.14,. ii 1 1 1.1..................................... 1 1............................... 3 1.3.........................
More informationFig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels).
Fig. 1 The scheme of glottal area as a function of time Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels). Fig, 4 Parametric representation
More informationCOE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho
COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate School of Economics and Institute of Economic Research
More information( ) Loewner SLE 13 February
( ) Loewner SLE 3 February 00 G. F. Lawler, Conformally Invariant Processes in the Plane, (American Mathematical Society, 005)., Summer School 009 (009 8 7-9 ) . d- (BES d ) d B t = (Bt, B t,, Bd t ) (d
More informationばらつき抑制のための確率最適制御
( ) http://wwwhayanuemnagoya-uacjp/ fujimoto/ 2011 3 9 11 ( ) 2011/03/09-11 1 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 2 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 3 / 46 (1/2) r + Controller - u Plant y
More informationGreen
28 6/-3 2 ax, bx R, L = d2 ax 2 dx 2 + bx d dx 2 L X = {Xt, t }. σ y = inf{t > ; Xt = y} t Xt L y = sup{t > ; Xt = y} 2. Gauss Gamma 3. tree 4. t Black-Scholes Madan-Roynette-Yor [6] σ y E x [exp λσ y
More informationIA 2013 : :10722 : 2 : :2 :761 :1 (23-27) : : ( / ) (1 /, ) / e.g. (Taylar ) e x = 1 + x + x xn n! +... sin x = x x3 6 + x5 x2n+1 + (
IA 2013 : :10722 : 2 : :2 :761 :1 23-27) : : 1 1.1 / ) 1 /, ) / e.g. Taylar ) e x = 1 + x + x2 2 +... + xn n! +... sin x = x x3 6 + x5 x2n+1 + 1)n 5! 2n + 1)! 2 2.1 = 1 e.g. 0 = 0.00..., π = 3.14..., 1
More information1 Fig. 1 Extraction of motion,.,,, 4,,, 3., 1, 2. 2.,. CHLAC,. 2.1,. (256 ).,., CHLAC. CHLAC, HLAC. 2.3 (HLAC ) r,.,. HLAC. N. 2 HLAC Fig. 2
CHLAC 1 2 3 3,. (CHLAC), 1).,.,, CHLAC,.,. Suspicious Behavior Detection based on CHLAC Method Hideaki Imanishi, 1 Toyohiro Hayashi, 2 Shuichi Enokida 3 and Toshiaki Ejima 3 We have proposed a method for
More informationudc-2.dvi
13 0.5 2 0.5 2 1 15 2001 16 2009 12 18 14 No.39, 2010 8 2009b 2009a Web Web Q&A 2006 2007a20082009 2007b200720082009 20072008 2009 2009 15 1 2 2 2.1 18 21 1 4 2 3 1(a) 1(b) 1(c) 1(d) 1) 18 16 17 21 10
More information2012専門分科会_new_4.pptx
d dt L L = 0 q i q i d dt L L = 0 r i i r i r r + Δr Δr δl = 0 dl dt = d dt i L L q i q i + q i i q i = q d L L i + q i i dt q i i q i = i L L q i L = 0, H = q q i L = E i q i i d dt L q q i i L = L(q
More information¼§À�ÍýÏÀ – Ê×ÎòÅŻҼ§À�¤È¥¹¥Ô¥ó¤æ¤é¤® - No.7, No.8, No.9
No.7, No.8, No.9 email: takahash@sci.u-hyogo.ac.jp Spring semester, 2012 Introduction (Critical Behavior) SCR ( b > 0) Arrott 2 Total Amplitude Conservation (TAC) Global Consistency (GC) TAC 2 / 25 Experimental
More information1 = = = (set) (element) a A a A a A a A a A {2, 5, (0, 1)}, [ 1, 1] = {x; 1 x 1}. (proposition) A = {x; P (x)} P (x) x x a A a A Remark. (i) {2, 0, 0,
2005 4 1 1 2 2 6 3 8 4 11 5 14 6 18 7 20 8 22 9 24 10 26 11 27 http://matcmadison.edu/alehnen/weblogic/logset.htm 1 1 = = = (set) (element) a A a A a A a A a A {2, 5, (0, 1)}, [ 1, 1] = {x; 1 x 1}. (proposition)
More information_念3)医療2009_夏.indd
Evaluation of the Social Benefits of the Regional Medical System Based on Land Price Information -A Hedonic Valuation of the Sense of Relief Provided by Health Care Facilities- Takuma Sugahara Ph.D. Abstract
More informationMacdonald, ,,, Macdonald. Macdonald,,,,,.,, Gauss,,.,, Lauricella A, B, C, D, Gelfand, A,., Heckman Opdam.,,,.,,., intersection,. Macdona
Macdonald, 2015.9.1 9.2.,,, Macdonald. Macdonald,,,,,.,, Gauss,,.,, Lauricella A, B, C, D, Gelfand, A,., Heckman Opdam.,,,.,,., intersection,. Macdonald,, q., Heckman Opdam q,, Macdonald., 1 ,,. Macdonald,
More information1 I
1 I 3 1 1.1 R x, y R x + y R x y R x, y, z, a, b R (1.1) (x + y) + z = x + (y + z) (1.2) x + y = y + x (1.3) 0 R : 0 + x = x x R (1.4) x R, 1 ( x) R : x + ( x) = 0 (1.5) (x y) z = x (y z) (1.6) x y =
More information16_.....E...._.I.v2006
55 1 18 Bull. Nara Univ. Educ., Vol. 55, No.1 (Cult. & Soc.), 2006 165 2002 * 18 Collaboration Between a School Athletic Club and a Community Sports Club A Case Study of SOLESTRELLA NARA 2002 Rie TAKAMURA
More information