ばらつき抑制のための確率最適制御

Size: px
Start display at page:

Download "ばらつき抑制のための確率最適制御"

Transcription

1 ( ) fujimoto/ ( ) 2011/03/ / 46

2 Outline ( ) 2011/03/ / 46

3 Outline ( ) 2011/03/ / 46

4 (1/2) r + Controller - u Plant y Plant Controller Plant: dx = f(x(t), u(t)) dt y(t) = h(x(t), u(t)) lim y(t) r(t) = 0 t x(t + 1) = A x(t) + B u(t) y(t) = C x(t) + D u(t) ( ) 2011/03/ / 46

5 (2/2) r + - Controller u Plant y x(t + 1) = A x(t) + B u(t) y(t) = C x(t) + D u(t) ( x(t) ) u(t) = K x(t) ( ) 2011/03/ / 46

6 vs ( ) vs ( ) y y o x o x ( ) 2011/03/ / 46

7 vs ( ) ( ) ( ) ( ) ( ) 2011/03/ / 46

8 vs ( ) ( ) ( ) EM ( + + ) ( + ) ( + ) ( + + ) ( ) 2011/03/ / 46

9 ?! ( ) 2011/03/ / 46

10 Outline ( ) 2011/03/ / 46

11 X p(y X) Y X Y Y X p(x Y) X, Y p(x, Y) = p(x Y)p(Y) = p(y X)p(X) p(x Y) { }} { p(x Y) = ( ) { }} { p(y X) p(y) }{{} = p(y X)p(X)dX {}}{ p(x) ( ) 2011/03/ / 46

12 { }} { p(x Y) = ( ) { }} { p(y X) {}}{ p(x) p(y X) p(x) dx : p(x Y)! p(x) ( ) p(y X) p(x) ( ) p(y X) ( ) 2011/03/ / 46

13 y = θu + w w N(0, σ w ) ( ) p(w) = p(y θu) = p(u, y θ) 1 (y θu)2 exp 2πσw 2σ w p(θ)! p(θ) = 1 exp (θ µ θ) 2 2πσθ 2σ θ p(u, y θ)p(θ) = ( ) exp(θ 2 )! ( ) 2011/03/ / 46

14 { xt+1 = Ax t + Bu t + w t, w t N(0, Q), x 0 N(µ 0, Σ 0 ) y t = Cx t + Du t + v t, v t N(0, R) Θ := {A, B, C, D, Q, R} X := {x 0, x 1,, x t, µ 0, Σ 0 } ( ) : Θ X! p(θ, X) p(θ)p(x) p(θ) p(x) ( ) 2011/03/ / 46

15 : p(θ, X Y) q(θ, X Y) p q KL KL(q, p) := q log q p dxdθ q F(q) q log p(y) = F(q) + KL(q, p) } {{ } ( ) q(θ, X Y) = q(θ Y)q(X Y) F(q(Θ, X Y)) Θ, X ( ) ( ) 2011/03/ / 46

16 VB-E VB-M q(x) (new) p(x) exp E q(θ) (old) [log p(x, Y Θ)] q(θ) (new) p(θ) exp E q(x) (old) [log p(θ, Y X)] E( ) p( ), q( ) ( ) ( ) 2011/03/ / 46

17 ( ) ( ) 2011/03/ / 46

18 Outline ( ) 2011/03/ / 46

19 ( ) { xt+1 = Ax t + Bu t + w t, w t N(0, Q), x 0 N(µ 0, Σ 0 ) y t = Cx t + Du t + v t, v t N(0, R) Θ I := {A, B, Q} Θ O := {C, D, R} X := {x 0, x 1,, x t, µ 0, Σ 0 } p(θ I, Θ O, X) p(θ I )p(θ O )p(x) ( ) 2011/03/ / 46

20 : p(y X) p(x) p(x Y) = p(y X)p(X)dX [Beal 03, Barber & Chiappa 07] x t 1 Q, R ( ) p(vec(a, B) µ AB, Q, G) = N(vec(A, B) µ AB, G Q) p(vec(c, D) µ CD, R, H) = N(vec(C, D) µ CD, H R) { p(q ν, S Q ) Q 1 (ν n 1)/2 exp 1 } 2 tr(q 1 S 1 Q ) { p(r η, S R ) R 1 (η l 1)/2 exp 1 } 2 tr(r 1 S 1 R ) ( ) 2011/03/ / 46

21 { ẋ = Ax + Bu y = Cx + Du x = Tx { x = T AT 1 x + TBu y = CT 1 x + D? ( ) 2011/03/ / 46

22 p(a, B, C, D Q 1, R 1 ) q(a, B, C, D Q 1, R 1 ) f T p(ā, B, C, D Q 1, R 1 )? f T q(ā, B, C, D Q 1, R 1 ) ( ) 2011/03/ / 46

23 p(a, B, C, D Q 1, R 1 ) f T p(ā, B, C, D Q 1, R 1 ) q(a, B, C, D Q 1, R 1 ) f T q(ā, B, C, D Q 1, R 1 ) ( ) 2011/03/ / 46

24 p(θ) x = Tx p(a, B, C, D, Q, R) = p(t AT 1, TB, CT 1, D, TQT T, R)! d(a, B, C, D, Q, R) = d(t AT 1, TB, CT 1, D, TQT T, R)! ( ) 2011/03/ / 46

25 X ( ) Θ? ( ( ) x = Tx?) Yes! ( ) 2011/03/ / 46

26 ( A = Ā = ( ) (, B = ) ( 1, B = ), C = ( 1, 0 ), D = 0 ), C = ( 1, 1 ), D = 1 ( ) 2011/03/ / 46

27 (l 2 ) l 2 : : ± 2σ 5 x ( ) 2011/03/ / 46

28 (50 ) (50 ) : : : ± 2σ ( ) 2011/03/ / 46

29 (200 ) (200 ) : : : ± 2σ ( ) 2011/03/ / 46

30 Outline ( ) 2011/03/ / 46

31 ?! : : ( ) 2011/03/ / 46

32 x t+1 = Ax t + Bu t T 1 ( ): J T (u) = ( ): J (u) = t=0 t=0 J T, J! ( x T t Qx t + u T t Ru ) t + x T T Fx T ( x T t Qx t + u T t Ru ) t u t = u t (x t, t) ( ) 2011/03/ / 46

33 V(x 0 ) ( ) ( V(x t ) := min J(u x t ) = min V(xt+1 ) + x T u t,,u T 1 u t t Qx ) t + u T t Ru t 2 V(x t ) = x T t Π t x t Π t x T t Π ( t x t = min x T Π u t t+1 t+1x t+1 + x T t Qx ) t + u T t Ru t ( = min (Axt + Bu t ) T Π t+1 (Ax t + Bu t ) + x T u t t Qx ) t + u T t Ru t u t 2 ( ) u t = (R + B T Π t+1 B) 1 (B T Π t+1 A)x t Π t =Q + A T Π t+1 A (A T Π t+1 B)(R + B T Π t+1 B) 1 (B T Π t+1 A) : Π N = F : Π t+1 = Π t = Π Π ( ) 2011/03/ / 46

34 LQG LQG(Linear Quadratic Gaussian) ϵ t T 1 ( ): J T (u) = E x t+1 = Ax t + Bu t + Gϵ t t=0 ( ): J (u) = E t=0 ( x T t Qx t + u T t Ru ) t + x T ( x T t Qx t + u T t Ru t) x 0 T Fx T x 0 ( ) 2011/03/ / 46

35 x t+1 = Ax t + Bu t + Gϵ t MCV(Minimum Cost variance) [Sain 66] J(u) =E[Ĵ] + λ var[ĵ] T 1 Ĵ = (x T t Qx t + u T t Ru t) + x T T Fx T t=0 RS(Risk Sensitive) [Whittle 81] J(u) = 2θ 1 log E[exp( θ 2 Ĵ)] =E[Ĵ] θ 4 var[ĵ] + O(θ2 ) A, B ( ) 2011/03/ / 46

36 [De Koning 82, F 10] x t+1 = A t x t + B t u t + G t ϵ t A t, B t, G t, ϵ t T 1 ( ( ): J T (u) = E x T t Qx t + u T t Ru t + tr (S cov[x t+1 x t ]) ) + x T T Fx T t=0 ( ( ): J (u) = E x T t Qx t + u T t Ru t + tr (S cov[x t+1 x t ]) ) t=0 tr (S cov[x t+1 x t ]) = E[x T Sx t+1 t+1 x t ] E[x t+1 ] T SE[x t+1 x t ] ( ) 2011/03/ / 46

37 S = 0 V ( ) V(x t ) := min J(u x t ) = min u t,,u T 1 u t ( E[V(xt+1 ) x t ] + x T t Qx t + u T t Ru ) t 2 V Π t V(x t ) = x T t Π t x t + β t ( ) 2011/03/ / 46

38 u t = (R + Σ BB + E[B T Π t+1 B]) 1 (Σ BA + E[B T Π t+1 A])x t Π t =Q + Σ AA + E[A T Π t+1 A] (Σ AB + E[A T Π t+1 B]) (R + Σ BB + E[B T Π t+1 B]) 1 (Σ BA + E[B T Π t+1 A]) Π T =F Σ XY := E[X T SY] E[X] T SE[Y] ( ) 2011/03/ / 46

39 u t = (R + B T ΠB) 1 B T ΠAx t Π =Q + A T ΠA A T ΠB(R + B T ΠB) 1 B T ΠA (S = 0, T = ) u t = (R + E[B T ΠB]) 1 E[B T ΠA]x t Π =Q + E[A T ΠA] E[A T ΠB](R + E[B T ΠB]) 1 E[B T ΠA] 2 ( ) 2011/03/ / 46

40 (1/2) ( ) 1 01 E[A] = ( ) 0 E[B] = cov[vec( A, B)] = Q = ( ) ( ) , R = 1, F = ( ) 2011/03/ / 46

41 数値例 (2/2) 従来法と提案法 (S = 100) の状態 x state X1 state X time time 従来法と提案法 (S = 100) の入力 u 藤本 健治 (名古屋大学) ばらつき抑制のための確率最適制御 /03/09-11 力学系研究会 40 / 46

42 ( ) x k+1 = A t x t lim E[x t x 0 ] = 0 t lim E[ x t 2 x 0 ] = 0 t [De Koning 82] ( ) 2011/03/ / 46

43 x t+1 = A t x t + B t u t u t = (R + Σ BB + E[B T ΠB]) 1 (Σ BA + E[B T ΠA])x t Π =Q + Σ AA + E[A T ΠA] (Σ AB + E[A T ΠB]) (R + Σ BB + E[B T ΠB]) 1 (Σ BA + E[B T ΠA])! ( ) 2011/03/ / 46

44 : T 1 ( J T (u) =E x T t Qx t + u T t Ru t + tr (S cov[x t+1 x t ]) ) + x T T Fx T t=0 B L B L Π :=Q + L T RL + E[(A BL) T (Π + S)(A BL)] E[(A BL)] T SE[(A BL)] B L u t = Lx t J t (u) := x T 0 B L t Fx 0 B L J J J ( ) 2011/03/ / 46

45 Outline ( ) 2011/03/ / 46

46 ? or vs? (2 ) ( ) 2011/03/ / 46

47 [1] D Barber and S Chiappa Unified inference for variational Bayesian linear Gaussian state-space models In Advances in Neural Information Processing Systems 19 (NIPS 20), pp The MIT Press, 2007 [2] M J Beal Variational Algorithms for Approximate Bayesian inference PhD thesis, University of Londong, London, UK, 2003 [3] W L De Koning Infinite horizon optimal control of linear discrete time systems with stochastic parameters Automatica, Vol 18, No 4, pp , 1982 [4] K Fujimoto, S Ogawa, Y Ota, and M Nakayama Optimal control of linear systems with stochastic parameters for variance suppression: The finite time horizon case In Proceedings of the 18th IFAC World Congress, 2011 [5] K Fujimoto, A Satoh, and S Fukunaga System identification based on variational bayes method and the inavriance under coordinate transformations Submitted, 2011 [6] M K Sain Control of linear systems according to the minimal variance criterion: A new approach to the disturbance problem IEEE Trans Autom Contr, Vol 11, No 1, pp , 1966 [7] P Whittle Risk-sensitive Linear/Quadratic/Gaussian control Advances in Applied Probability, Vol 13, pp , 1981 [8],, Vol 20, No 10, pp , 2007 [9], H D, Vol J91-D, No 6, pp , 2008 ( ) 2011/03/ / 46

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi) 0. A A = 4 IC () det A () A () x + y + z = x y z X Y Z = A x y z ( 5) ( s5590) 0. a + b + c b c () a a + b + c c a b a + b + c 0 a b c () a 0 c b b c 0 a c b a 0 0. A A = 7 5 4 5 0 ( 5) ( s5590) () A ()

More information

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P 9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)

More information

通信容量制約を考慮したフィードバック制御 - 電子情報通信学会 情報理論研究会(IT) 若手研究者のための講演会

通信容量制約を考慮したフィードバック制御 -  電子情報通信学会 情報理論研究会(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

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,. 24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)

More information

Part () () Γ Part ,

Part () () Γ Part , Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35

More information

CVaR

CVaR CVaR 20 4 24 3 24 1 31 ,.,.,. Markowitz,., (Value-at-Risk, VaR) (Conditional Value-at-Risk, CVaR). VaR, CVaR VaR. CVaR, CVaR. CVaR,,.,.,,,.,,. 1 5 2 VaR CVaR 6 2.1................................................

More information

数値計算:有限要素法

数値計算:有限要素法 ( ) 1 / 61 1 2 3 4 ( ) 2 / 61 ( ) 3 / 61 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61 u(x) 0 L x ( ) 6 / 61 P(0) P(L) f d dx ( EA du dx

More information

( )/2 hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1

( )/2   hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1 ( )/2 http://www2.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html 1 2011 ( )/2 2 2011 4 1 2 1.1 1 2 1 2 3 4 5 1.1.1 sample space S S = {H, T } H T T H S = {(H, H), (H, T ), (T, H), (T, T )} (T, H) S

More information

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11

More information

renshumondai-kaito.dvi

renshumondai-kaito.dvi 3 1 13 14 1.1 1 44.5 39.5 49.5 2 0.10 2 0.10 54.5 49.5 59.5 5 0.25 7 0.35 64.5 59.5 69.5 8 0.40 15 0.75 74.5 69.5 79.5 3 0.15 18 0.90 84.5 79.5 89.5 2 0.10 20 1.00 20 1.00 2 1.2 1 16.5 20.5 12.5 2 0.10

More information

K E N Z OU

K E N Z OU K E N Z OU 11 1 1 1.1..................................... 1.1.1............................ 1.1..................................................................................... 4 1.........................................

More information

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0 1 1 1.1 1.) T D = T = D = kn 1. 1.4) F W = F = W/ = kn/ = 15 kn 1. 1.9) R = W 1 + W = 6 + 5 = 11 N. 1.9) W b W 1 a = a = W /W 1 )b = 5/6) = 5 cm 1.4 AB AC P 1, P x, y x, y y x 1.4.) P sin 6 + P 1 sin 45

More information

Kalman ( ) 1) (Kalman filter) ( ) t y 0,, y t x ˆx 3) 10) t x Y [y 0,, y ] ) x ( > ) ˆx (prediction) ) x ( ) ˆx (filtering) )

Kalman ( ) 1) (Kalman filter) ( ) t y 0,, y t x ˆx 3) 10) t x Y [y 0,, y ] ) x ( > ) ˆx (prediction) ) x ( ) ˆx (filtering) ) 1 -- 5 6 2009 3 R.E. Kalman ( ) H 6-1 6-2 6-3 H Rudolf Emil Kalman IBM IEEE Medal of Honor(1974) (1985) c 2011 1/(23) 1 -- 5 -- 6 6--1 2009 3 Kalman ( ) 1) (Kalman filter) ( ) t y 0,, y t x ˆx 3) 10) t

More information

4 4 4 a b c d a b A c d A a da ad bce O E O n A n O ad bc a d n A n O 5 {a n } S n a k n a n + k S n a a n+ S n n S n n log x x {xy } x, y x + y 7 fx

4 4 4 a b c d a b A c d A a da ad bce O E O n A n O ad bc a d n A n O 5 {a n } S n a k n a n + k S n a a n+ S n n S n n log x x {xy } x, y x + y 7 fx 4 4 5 4 I II III A B C, 5 7 I II A B,, 8, 9 I II A B O A,, Bb, b, Cc, c, c b c b b c c c OA BC P BC OP BC P AP BC n f n x xn e x! e n! n f n x f n x f n x f k x k 4 e > f n x dx k k! fx sin x cos x tan

More information

December 28, 2018

December 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 information

meiji_resume_1.PDF

meiji_resume_1.PDF β β β (q 1,q,..., q n ; p 1, p,..., p n ) H(q 1,q,..., q n ; p 1, p,..., p n ) Hψ = εψ ε k = k +1/ ε k = k(k 1) (x, y, z; p x, p y, p z ) (r; p r ), (θ; p θ ), (ϕ; p ϕ ) ε k = 1/ k p i dq i E total = E

More information

(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y

(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y [ ] 7 0.1 2 2 + y = t sin t IC ( 9) ( s090101) 0.2 y = d2 y 2, y = x 3 y + y 2 = 0 (2) y + 2y 3y = e 2x 0.3 1 ( y ) = f x C u = y x ( 15) ( s150102) [ ] y/x du x = Cexp f(u) u (2) x y = xey/x ( 16) ( s160101)

More information

( 28 ) ( ) ( ) 0 This note is c 2016, 2017 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercial purp

( 28 ) ( ) ( ) 0 This note is c 2016, 2017 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercial purp ( 28) ( ) ( 28 9 22 ) 0 This ote is c 2016, 2017 by Setsuo Taiguchi. It may be used for persoal or classroom purposes, but ot for commercial purposes. i (http://www.stat.go.jp/teacher/c2epi1.htm ) = statistics

More information

S I. dy fx x fx y fx + C 3 C dy fx 4 x, y dy v C xt y C v e kt k > xt yt gt [ v dt dt v e kt xt v e kt + C k x v + C C k xt v k 3 r r + dr e kt S dt d

S I. dy fx x fx y fx + C 3 C dy fx 4 x, y dy v C xt y C v e kt k > xt yt gt [ v dt dt v e kt xt v e kt + C k x v + C C k xt v k 3 r r + dr e kt S dt d S I.. http://ayapin.film.s.dendai.ac.jp/~matuda /TeX/lecture.html PDF PS.................................... 3.3.................... 9.4................5.............. 3 5. Laplace................. 5....

More information

DE-resume

DE-resume - 2011, http://c-faculty.chuo-u.ac.jp/ nishioka/ 2 11 21131 : 4 1 x y(x, y (x,y (x,,y (n, (1.1 F (x, y, y,y,,y (n =0. (1.1 n. (1.1 y(x. y(x (1.1. 1 1 1 1.1... 2 1.2... 9 1.3 1... 26 2 2 34 2.1,... 35 2.2

More information

DVIOUT

DVIOUT A. A. A-- [ ] f(x) x = f 00 (x) f 0 () =0 f 00 () > 0= f(x) x = f 00 () < 0= f(x) x = A--2 [ ] f(x) D f 00 (x) > 0= y = f(x) f 00 (x) < 0= y = f(x) P (, f()) f 00 () =0 A--3 [ ] y = f(x) [, b] x = f (y)

More information

入試の軌跡

入試の軌跡 4 y O x 4 Typed by L A TEX ε ) ) ) 6 4 ) 4 75 ) http://kumamoto.s.xrea.com/plan/.. PDF) Ctrl +L) Ctrl +) Ctrl + Ctrl + ) ) Alt + ) Alt + ) ESC. http://kumamoto.s.xrea.com/nyusi/kumadai kiseki ri i.pdf

More information

6.1 (P (P (P (P (P (P (, P (, P.101

6.1 (P (P (P (P (P (P (, P (, P.101 (008 0 3 7 ( ( ( 00 1 (P.3 1 1.1 (P.3.................. 1 1. (P.4............... 1 (P.15.1 (P.15................. (P.18............3 (P.17......... 3.4 (P................ 4 3 (P.7 4 3.1 ( P.7...........

More information

I A A441 : April 15, 2013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida )

I A A441 : April 15, 2013 Version : 1.1 I   Kawahira, Tomoki TA (Shigehiro, Yoshida ) I013 00-1 : April 15, 013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida) http://www.math.nagoya-u.ac.jp/~kawahira/courses/13s-tenbou.html pdf * 4 15 4 5 13 e πi = 1 5 0 5 7 3 4 6 3 6 10 6 17

More information

B ver B

B ver B B ver. 2017.02.24 B Contents 1 11 1.1....................... 11 1.1.1............. 11 1.1.2.......................... 12 1.2............................. 14 1.2.1................ 14 1.2.2.......................

More information

S I. dy fx x fx y fx + C 3 C vt dy fx 4 x, y dy yt gt + Ct + C dt v e kt xt v e kt + C k x v k + C C xt v k 3 r r + dr e kt S Sr πr dt d v } dt k e kt

S I. dy fx x fx y fx + C 3 C vt dy fx 4 x, y dy yt gt + Ct + C dt v e kt xt v e kt + C k x v k + C C xt v k 3 r r + dr e kt S Sr πr dt d v } dt k e kt S I. x yx y y, y,. F x, y, y, y,, y n http://ayapin.film.s.dendai.ac.jp/~matuda n /TeX/lecture.html PDF PS yx.................................... 3.3.................... 9.4................5..............

More information

6.1 (P (P (P (P (P (P (, P (, P.

6.1 (P (P (P (P (P (P (, P (, P. (011 30 7 0 ( ( 3 ( 010 1 (P.3 1 1.1 (P.4.................. 1 1. (P.4............... 1 (P.15.1 (P.16................. (P.0............3 (P.18 3.4 (P.3............... 4 3 (P.9 4 3.1 (P.30........... 4 3.

More information

??

?? ( ) 2014 2014 1/119 = (ISS) ISS ISS ISS iss-clf iss-clf ISS = (ISS) FB 2014 2/119 = (ISS) ISS ISS ISS iss-clf iss-clf ISS R + : 0 K: γ: R + R + K γ γ(0) = 0 K : γ: R + R + K γ K γ(r) (r ) FB K K K K R

More information

ver.1 / c /(13)

ver.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 information

( )

( ) 7..-8..8.......................................................................... 4.................................... 3...................................... 3..3.................................. 4.3....................................

More information

名古屋工業大の数学 2000 年 ~2015 年 大学入試数学動画解説サイト

名古屋工業大の数学 2000 年 ~2015 年 大学入試数学動画解説サイト 名古屋工業大の数学 年 ~5 年 大学入試数学動画解説サイト http://mathroom.jugem.jp/ 68 i 4 3 III III 3 5 3 ii 5 6 45 99 5 4 3. () r \= S n = r + r + 3r 3 + + nr n () x > f n (x) = e x + e x + 3e 3x + + ne nx f(x) = lim f n(x) lim

More information

A A = a 41 a 42 a 43 a 44 A (7) 1 (3) A = M 12 = = a 41 (8) a 41 a 43 a 44 (3) n n A, B a i AB = A B ii aa

A A = a 41 a 42 a 43 a 44 A (7) 1 (3) A = M 12 = = a 41 (8) a 41 a 43 a 44 (3) n n A, B a i AB = A B ii aa 1 2 21 2 2 [ ] a 11 a 12 A = a 21 a 22 (1) A = a 11 a 22 a 12 a 21 (2) 3 3 n n A A = n ( 1) i+j a ij M ij i =1 n (3) j=1 M ij A i j (n 1) (n 1) 2-1 3 3 A A = a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33

More information

? (EM),, EM? (, 2004/ 2002) von Mises-Fisher ( 2004) HMM (MacKay 1997) LDA (Blei et al. 2001) PCFG ( 2004)... Variational Bayesian methods for Natural

? (EM),, EM? (, 2004/ 2002) von Mises-Fisher ( 2004) HMM (MacKay 1997) LDA (Blei et al. 2001) PCFG ( 2004)... Variational Bayesian methods for Natural SLC Internal tutorial Daichi Mochihashi daichi.mochihashi@atr.jp ATR SLC 2005.6.21 (Tue) 13:15 15:00@Meeting Room 1 Variational Bayesian methods for Natural Language Processing p.1/30 ? (EM),, EM? (, 2004/

More information

08-Note2-web

08-Note2-web r(t) t r(t) O v(t) = dr(t) dt a(t) = dv(t) dt = d2 r(t) dt 2 r(t), v(t), a(t) t dr(t) dt r(t) =(x(t),y(t),z(t)) = d 2 r(t) dt 2 = ( dx(t) dt ( d 2 x(t) dt 2, dy(t), dz(t) dt dt ), d2 y(t) dt 2, d2 z(t)

More information

st.dvi

st.dvi 9 3 5................................... 5............................. 5....................................... 5.................................. 7.........................................................................

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> MATLAB/Simulink による現代制御入門 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/9241 このサンプルページの内容は, 初版 1 刷発行当時のものです. i MATLAB/Simulink MATLAB/Simulink 1. 1 2. 3. MATLAB/Simulink

More information

mugensho.dvi

mugensho.dvi 1 1 f (t) lim t a f (t) = 0 f (t) t a 1.1 (1) lim(t 1) 2 = 0 t 1 (t 1) 2 t 1 (2) lim(t 1) 3 = 0 t 1 (t 1) 3 t 1 2 f (t), g(t) t a lim t a f (t) g(t) g(t) f (t) = o(g(t)) (t a) = 0 f (t) (t 1) 3 1.2 lim

More information

1 y(t)m b k u(t) ẋ = [ 0 1 k m b m x + [ 0 1 m u, x = [ ẏ y (1) y b k m u

1 y(t)m b k u(t) ẋ = [ 0 1 k m b m x + [ 0 1 m u, x = [ ẏ y (1) y b k m u ( ) LPV( ) 1 y(t)m b k u(t) ẋ = [ 0 1 k m b m x + [ 0 1 m u, x = [ ẏ y (1) y b k m u m 1 m m 2, b 1 b b 2, k 1 k k 2 (2) [m b k ( ) k 0 b m ( ) 2 ẋ = Ax, x(0) 0 (3) (x(t) 0) ( ) V (x) V (x) = x T P x >

More information

n ( (

n ( ( 1 2 27 6 1 1 m-mat@mathscihiroshima-uacjp 2 http://wwwmathscihiroshima-uacjp/~m-mat/teach/teachhtml 2 1 3 11 3 111 3 112 4 113 n 4 114 5 115 5 12 7 121 7 122 9 123 11 124 11 125 12 126 2 2 13 127 15 128

More information

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,.

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,. (1 C205) 4 10 (2 C206) 4 11 (2 B202) 4 12 25(2013) http://www.math.is.tohoku.ac.jp/~obata,.,,,..,,. 1. 2. 3. 4. 5. 6. 7. 8. 1., 2007 ( ).,. 2. P. G., 1995. 3. J. C., 1988. 1... 2.,,. ii 3.,. 4. F. ( ),..

More information

QCD 1 QCD GeV 2014 QCD 2015 QCD SU(3) QCD A µ g µν QCD 1

QCD 1 QCD GeV 2014 QCD 2015 QCD SU(3) QCD A µ g µν QCD 1 QCD 1 QCD GeV 2014 QCD 2015 QCD SU(3) QCD A µ g µν QCD 1 (vierbein) QCD QCD 1 1: QCD QCD Γ ρ µν A µ R σ µνρ F µν g µν A µ Lagrangian gr TrFµν F µν No. Yes. Yes. No. No! Yes! [1] Nash & Sen [2] Riemann

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,,.

,. 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 information

Microsoft Word - 11問題表紙(選択).docx

Microsoft Word - 11問題表紙(選択).docx A B A.70g/cm 3 B.74g/cm 3 B C 70at% %A C B at% 80at% %B 350 C γ δ y=00 x-y ρ l S ρ C p k C p ρ C p T ρ l t l S S ξ S t = ( k T ) ξ ( ) S = ( k T) ( ) t y ξ S ξ / t S v T T / t = v T / y 00 x v S dy dx

More information

さくらの個別指導 ( さくら教育研究所 ) A 2 P Q 3 R S T R S T P Q ( ) ( ) m n m n m n n n

さくらの個別指導 ( さくら教育研究所 ) A 2 P Q 3 R S T R S T P Q ( ) ( ) m n m n m n n n 1 1.1 1.1.1 A 2 P Q 3 R S T R S T P 80 50 60 Q 90 40 70 80 50 60 90 40 70 8 5 6 1 1 2 9 4 7 2 1 2 3 1 2 m n m n m n n n n 1.1 8 5 6 9 4 7 2 6 0 8 2 3 2 2 2 1 2 1 1.1 2 4 7 1 1 3 7 5 2 3 5 0 3 4 1 6 9 1

More information

pdf

pdf http://www.ns.kogakuin.ac.jp/~ft13389/lecture/physics1a2b/ pdf I 1 1 1.1 ( ) 1. 30 m µm 2. 20 cm km 3. 10 m 2 cm 2 4. 5 cm 3 km 3 5. 1 6. 1 7. 1 1.2 ( ) 1. 1 m + 10 cm 2. 1 hr + 6400 sec 3. 3.0 10 5 kg

More information

II ( ) (7/31) II ( [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re

II ( ) (7/31) II (  [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re II 29 7 29-7-27 ( ) (7/31) II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I Euler Navier

More information

1W II K =25 A (1) office(a439) (2) A4 etc. 12:00-13:30 Cafe David 1 2 TA appointment Cafe D

1W II K =25 A (1) office(a439) (2) A4 etc. 12:00-13:30 Cafe David 1 2 TA  appointment Cafe D 1W II K200 : October 6, 2004 Version : 1.2, kawahira@math.nagoa-u.ac.jp, http://www.math.nagoa-u.ac.jp/~kawahira/courses.htm TA M1, m0418c@math.nagoa-u.ac.jp TA Talor Jacobian 4 45 25 30 20 K2-1W04-00

More information

tokei01.dvi

tokei01.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 information

1: 3.3 1/8000 1/ m m/s v = 2kT/m = 2RT/M k R 8.31 J/(K mole) M 18 g 1 5 a v t πa 2 vt kg (

1: 3.3 1/8000 1/ m m/s v = 2kT/m = 2RT/M k R 8.31 J/(K mole) M 18 g 1 5 a v t πa 2 vt kg ( 1905 1 1.1 0.05 mm 1 µm 2 1 1 2004 21 2004 7 21 2005 web 2 [1, 2] 1 1: 3.3 1/8000 1/30 3 10 10 m 3 500 m/s 4 1 10 19 5 6 7 1.2 3 4 v = 2kT/m = 2RT/M k R 8.31 J/(K mole) M 18 g 1 5 a v t πa 2 vt 6 6 10

More information

Microsoft Word - 信号処理3.doc

Microsoft Word - 信号処理3.doc Junji OHTSUBO 2012 FFT FFT SN sin cos x v ψ(x,t) = f (x vt) (1.1) t=0 (1.1) ψ(x,t) = A 0 cos{k(x vt) + φ} = A 0 cos(kx ωt + φ) (1.2) A 0 v=ω/k φ ω k 1.3 (1.2) (1.2) (1.2) (1.1) 1.1 c c = a + ib, a = Re[c],

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 3.............................. 3.............................. 4.3 4................... 5.4........................ 6.5........................ 8.6...........................7

More information

‚åŁÎ“·„´Šš‡ðŠp‡¢‡½‹âfi`fiI…A…‰…S…−…Y…•‡ÌMarkovŸA“½fiI›ð’Í

‚åŁÎ“·„´Šš‡ðŠp‡¢‡½‹âfi`fiI…A…‰…S…−…Y…•‡ÌMarkovŸA“½fiI›ð’Í Markov 2009 10 2 Markov 2009 10 2 1 / 25 1 (GA) 2 GA 3 4 Markov 2009 10 2 2 / 25 (GA) (GA) L ( 1) I := {0, 1} L f : I (0, ) M( 2) S := I M GA (GA) f (i) i I Markov 2009 10 2 3 / 25 (GA) ρ(i, j), i, j I

More information

03.Œk’ì

03.Œ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

untitled

untitled II(c) 1 October. 21, 2009 1 CS53 yamamoto@cs.kobe-u.ac.jp 3 1 7 1.1 : : : : : : : : : : : : : : : : : : : : : : 7 1.2 : : : : : : : : : : : : : : : : 8 1.2.1 : : : : : : : : : : : : : : : : : : : 8 1.2.2

More information

n ξ 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 information

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A .. Laplace ). A... i),. ω i i ). {ω,..., ω } Ω,. ii) Ω. Ω. A ) r, A P A) P A) r... ).. Ω {,, 3, 4, 5, 6}. i i 6). A {, 4, 6} P A) P A) 3 6. ).. i, j i, j) ) Ω {i, j) i 6, j 6}., 36. A. A {i, j) i j }.

More information

(5) 75 (a) (b) ( 1 ) v ( 1 ) E E 1 v (a) ( 1 ) x E E (b) (a) (b)

(5) 75 (a) (b) ( 1 ) v ( 1 ) E E 1 v (a) ( 1 ) x E E (b) (a) (b) (5) 74 Re, bondar laer (Prandtl) Re z ω z = x (5) 75 (a) (b) ( 1 ) v ( 1 ) E E 1 v (a) ( 1 ) x E E (b) (a) (b) (5) 76 l V x ) 1/ 1 ( 1 1 1 δ δ = x Re x p V x t V l l (1-1) 1/ 1 δ δ δ δ = x Re p V x t V

More information

17 ( ) II III A B C(100 ) 1, 2, 6, 7 II A B (100 ) 2, 5, 6 II A B (80 ) 8 10 I II III A B C(80 ) 1 a 1 = 1 2 a n+1 = a n + 2n + 1 (n = 1,

17 ( ) II III A B C(100 ) 1, 2, 6, 7 II A B (100 ) 2, 5, 6 II A B (80 ) 8 10 I II III A B C(80 ) 1 a 1 = 1 2 a n+1 = a n + 2n + 1 (n = 1, 17 ( ) 17 5 1 4 II III A B C(1 ) 1,, 6, 7 II A B (1 ), 5, 6 II A B (8 ) 8 1 I II III A B C(8 ) 1 a 1 1 a n+1 a n + n + 1 (n 1,,, ) {a n+1 n } (1) a 4 () a n OA OB AOB 6 OAB AB : 1 P OB Q OP AQ R (1) PQ

More information

QMII_10.dvi

QMII_10.dvi 65 1 1.1 1.1.1 1.1 H H () = E (), (1.1) H ν () = E ν () ν (). (1.) () () = δ, (1.3) μ () ν () = δ(μ ν). (1.4) E E ν () E () H 1.1: H α(t) = c (t) () + dνc ν (t) ν (), (1.5) H () () + dν ν () ν () = 1 (1.6)

More information

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x [ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),

More information

y π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a =

y π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a = [ ] 9 IC. dx = 3x 4y dt dy dt = x y u xt = expλt u yt λ u u t = u u u + u = xt yt 6 3. u = x, y, z = x + y + z u u 9 s9 grad u ux, y, z = c c : grad u = u x i + u y j + u k i, j, k z x, y, z grad u v =

More information

http://www.ike-dyn.ritsumei.ac.jp/ hyoo/wave.html 1 1, 5 3 1.1 1..................................... 3 1.2 5.1................................... 4 1.3.......................... 5 1.4 5.2, 5.3....................

More information

I

I I io@hiroshima-u.ac.jp 27 6 A A. /a δx = lim a + a exp π x2 a 2 = lim a + a = lim a + a exp a 2 π 2 x 2 + a 2 2 x a x = lim a + a Sic a x = lim a + a Rect a Gaussia Loretzia Bilateral expoetial Normalized

More information

& 3 3 ' ' (., (Pixel), (Light Intensity) (Random Variable). (Joint Probability). V., V = {,,, V }. i x i x = (x, x,, x V ) T. x i i (State Variable),

& 3 3 ' ' (., (Pixel), (Light Intensity) (Random Variable). (Joint Probability). V., V = {,,, V }. i x i x = (x, x,, x V ) T. x i i (State Variable), .... Deeping and Expansion of Large-Scale Random Fields and Probabilistic Image Processing Kazuyuki Tanaka The mathematical frameworks of probabilistic image processing are formulated by means of Markov

More information

CALCULUS II (Hiroshi SUZUKI ) f(x, y) A(a, b) 1. P (x, y) A(a, b) A(a, b) f(x, y) c f(x, y) A(a, b) c f(x, y) c f(x, y) c (x a, y b)

CALCULUS II (Hiroshi SUZUKI ) f(x, y) A(a, b) 1. P (x, y) A(a, b) A(a, b) f(x, y) c f(x, y) A(a, b) c f(x, y) c f(x, y) c (x a, y b) CALCULUS II (Hiroshi SUZUKI ) 16 1 1 1.1 1.1 f(x, y) A(a, b) 1. P (x, y) A(a, b) A(a, b) f(x, y) c f(x, y) A(a, b) c f(x, y) c f(x, y) c (x a, y b) lim f(x, y) = lim f(x, y) = lim f(x, y) = c. x a, y b

More information

( ) ( )

( ) ( ) 20 21 2 8 1 2 2 3 21 3 22 3 23 4 24 5 25 5 26 6 27 8 28 ( ) 9 3 10 31 10 32 ( ) 12 4 13 41 0 13 42 14 43 0 15 44 17 5 18 6 18 1 1 2 2 1 2 1 0 2 0 3 0 4 0 2 2 21 t (x(t) y(t)) 2 x(t) y(t) γ(t) (x(t) y(t))

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,. 23(2011) (1 C104) 5 11 (2 C206) 5 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 ( ). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5.. 6.. 7.,,. 8.,. 1. (75%

More information

I A A441 : April 21, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) Google

I 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 information

A 2 3. m S m = {x R m+1 x = 1} U + k = {x S m x k > 0}, U k = {x S m x k < 0}, ϕ ± k (x) = (x 0,..., ˆx k,... x m ) 1. {(U ± k, ϕ± k ) 0 k m} S m 1.2.

A 2 3. m S m = {x R m+1 x = 1} U + k = {x S m x k > 0}, U k = {x S m x k < 0}, ϕ ± k (x) = (x 0,..., ˆx k,... x m ) 1. {(U ± k, ϕ± k ) 0 k m} S m 1.2. A A 1 A 5 A 6 1 2 3 4 5 6 7 1 1.1 1.1 (). Hausdorff M R m M M {U α } U α R m E α ϕ α : U α E α U α U β = ϕ α (ϕ β ϕβ (U α U β )) 1 : ϕ β (U α U β ) ϕ α (U α U β ) C M a m dim M a U α ϕ α {x i, 1 i m} {U,

More information

.3. (x, x = (, u = = 4 (, x x = 4 x, x 0 x = 0 x = 4 x.4. ( z + z = 8 z, z 0 (z, z = (0, 8, (,, (8, 0 3 (0, 8, (,, (8, 0 z = z 4 z (g f(x = g(

.3. (x, x = (, u = = 4 (, x x = 4 x, x 0 x = 0 x = 4 x.4. ( z + z = 8 z, z 0 (z, z = (0, 8, (,, (8, 0 3 (0, 8, (,, (8, 0 z = z 4 z (g f(x = g( 06 5.. ( y = x x y 5 y 5 = (x y = x + ( y = x + y = x y.. ( Y = C + I = 50 + 0.5Y + 50 r r = 00 0.5Y ( L = M Y r = 00 r = 0.5Y 50 (3 00 0.5Y = 0.5Y 50 Y = 50, r = 5 .3. (x, x = (, u = = 4 (, x x = 4 x,

More information

熊本県数学問題正解

熊本県数学問題正解 00 y O x Typed by L A TEX ε ( ) (00 ) 5 4 4 ( ) http://www.ocn.ne.jp/ oboetene/plan/. ( ) (009 ) ( ).. http://www.ocn.ne.jp/ oboetene/plan/eng.html 8 i i..................................... ( )0... (

More information

( ) ) ) ) 5) 1 J = σe 2 6) ) 9) 1955 Statistical-Mechanical Theory of Irreversible Processes )

( ) ) ) ) 5) 1 J = σe 2 6) ) 9) 1955 Statistical-Mechanical Theory of Irreversible Processes ) ( 3 7 4 ) 2 2 ) 8 2 954 2) 955 3) 5) J = σe 2 6) 955 7) 9) 955 Statistical-Mechanical Theory of Irreversible Processes 957 ) 3 4 2 A B H (t) = Ae iωt B(t) = B(ω)e iωt B(ω) = [ Φ R (ω) Φ R () ] iω Φ R (t)

More information

,,,17,,, ( ),, E Q [S T F t ] < S t, t [, T ],,,,,,,,

,,,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 information

ohp_06nov_tohoku.dvi

ohp_06nov_tohoku.dvi 2006 11 28 1. (1) ẋ = ax = x(t) =Ce at C C>0 a0 x(t) 0(t )!! 1 0.8 0.6 0.4 0.2 2 4 6 8 10-0.2 (1) a =2 C =1 1. (1) τ>0 (2) ẋ(t) = ax(t τ) 4 2 2 4 6 8 10-2 -4 (2) a =2 τ =1!! 1. (2) A. (2)

More information

2019 1 5 0 3 1 4 1.1.................... 4 1.1.1......................... 4 1.1.2........................ 5 1.1.3................... 5 1.1.4........................ 6 1.1.5......................... 6 1.2..........................

More information

211 kotaro@math.titech.ac.jp 1 R *1 n n R n *2 R n = {(x 1,..., x n ) x 1,..., x n R}. R R 2 R 3 R n R n R n D D R n *3 ) (x 1,..., x n ) f(x 1,..., x n ) f D *4 n 2 n = 1 ( ) 1 f D R n f : D R 1.1. (x,

More information

( ) Loewner SLE 13 February

( ) 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

K E N Z U 2012 7 16 HP M. 1 1 4 1.1 3.......................... 4 1.2................................... 4 1.2.1..................................... 4 1.2.2.................................... 5................................

More information

) a + b = i + 6 b c = 6i j ) a = 0 b = c = 0 ) â = i + j 0 ˆb = 4) a b = b c = j + ) cos α = cos β = 6) a ˆb = b ĉ = 0 7) a b = 6i j b c = i + 6j + 8)

) a + b = i + 6 b c = 6i j ) a = 0 b = c = 0 ) â = i + j 0 ˆb = 4) a b = b c = j + ) cos α = cos β = 6) a ˆb = b ĉ = 0 7) a b = 6i j b c = i + 6j + 8) 4 4 ) a + b = i + 6 b c = 6i j ) a = 0 b = c = 0 ) â = i + j 0 ˆb = 4) a b = b c = j + ) cos α = cos β = 6) a ˆb = b ĉ = 0 7) a b = 6i j b c = i + 6j + 8) a b a b = 6i j 4 b c b c 9) a b = 4 a b) c = 7

More information

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C 8 ( ) 8 5 4 I II III A B C( ),,, 5 I II A B ( ),, I II A B (8 ) 6 8 I II III A B C(8 ) n ( + x) n () n C + n C + + n C n = 7 n () 7 9 C : y = x x A(, 6) () A C () C P AP Q () () () 4 A(,, ) B(,, ) C(,,

More information

( ) ( ) 1729 (, 2016:17) = = (1) 1 1

( ) ( ) 1729 (, 2016:17) = = (1) 1 1 1729 1 2016 10 28 1 1729 1111 1111 1729 (1887 1920) (1877 1947) 1729 (, 2016:17) 12 3 1728 9 3 729 1729 = 12 3 + 1 3 = 10 3 + 9 3 (1) 1 1 2 1729 1729 19 13 7 = 1729 = 12 3 + 1 3 = 10 3 + 9 3 13 7 = 91

More information

x,, z v = (, b, c) v v 2 + b 2 + c 2 x,, z 1 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) v 1 = ( 1, b 1, c 1 ), v 2 = ( 2, b 2, c 2 ) v

x,, z v = (, b, c) v v 2 + b 2 + c 2 x,, z 1 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) v 1 = ( 1, b 1, c 1 ), v 2 = ( 2, b 2, c 2 ) v 12 -- 1 4 2009 9 4-1 4-2 4-3 4-4 4-5 4-6 4-7 4-8 4-9 4-10 c 2011 1/(13) 4--1 2009 9 3 x,, z v = (, b, c) v v 2 + b 2 + c 2 x,, z 1 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) v 1 = ( 1, b 1, c 1 ), v 2

More information

III 1 (X, d) d U d X (X, d). 1. (X, d).. (i) d(x, y) d(z, y) d(x, z) (ii) d(x, y) d(z, w) d(x, z) + d(y, w) 2. (X, d). F X.. (1), X F, (2) F 1, F 2 F

III 1 (X, d) d U d X (X, d). 1. (X, d).. (i) d(x, y) d(z, y) d(x, z) (ii) d(x, y) d(z, w) d(x, z) + d(y, w) 2. (X, d). F X.. (1), X F, (2) F 1, F 2 F III 1 (X, d) d U d X (X, d). 1. (X, d).. (i) d(x, y) d(z, y) d(x, z) (ii) d(x, y) d(z, w) d(x, z) + d(y, w) 2. (X, d). F X.. (1), X F, (2) F 1, F 2 F F 1 F 2 F, (3) F λ F λ F λ F. 3., A λ λ A λ. B λ λ

More information

W u = u(x, t) u tt = a 2 u xx, a > 0 (1) D := {(x, t) : 0 x l, t 0} u (0, t) = 0, u (l, t) = 0, t 0 (2)

W u = u(x, t) u tt = a 2 u xx, a > 0 (1) D := {(x, t) : 0 x l, t 0} u (0, t) = 0, u (l, t) = 0, t 0 (2) 3 215 4 27 1 1 u u(x, t) u tt a 2 u xx, a > (1) D : {(x, t) : x, t } u (, t), u (, t), t (2) u(x, ) f(x), u(x, ) t 2, x (3) u(x, t) X(x)T (t) u (1) 1 T (t) a 2 T (t) X (x) X(x) α (2) T (t) αa 2 T (t) (4)

More information

A B P (A B) = P (A)P (B) (3) A B A B P (B A) A B A B P (A B) = P (B A)P (A) (4) P (B A) = P (A B) P (A) (5) P (A B) P (B A) P (A B) A B P

A B P (A B) = P (A)P (B) (3) A B A B P (B A) A B A B P (A B) = P (B A)P (A) (4) P (B A) = P (A B) P (A) (5) P (A B) P (B A) P (A B) A B P 1 1.1 (population) (sample) (event) (trial) Ω () 1 1 Ω 1.2 P 1. A A P (A) 0 1 0 P (A) 1 (1) 2. P 1 P 0 1 6 1 1 6 0 3. A B P (A B) = P (A) + P (B) (2) A B A B A 1 B 2 A B 1 2 1 2 1 1 2 2 3 1.3 A B P (A

More information

³ÎΨÏÀ

³ÎΨÏÀ 2017 12 12 Makoto Nakashima 2017 12 12 1 / 22 2.1. C, D π- C, D. A 1, A 2 C A 1 A 2 C A 3, A 4 D A 1 A 2 D Makoto Nakashima 2017 12 12 2 / 22 . (,, L p - ). Makoto Nakashima 2017 12 12 3 / 22 . (,, L p

More information

IA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................

More information

2016 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 1 16 2 1 () X O 3 (O1) X O, O (O2) O O (O3) O O O X (X, O) O X X (O1), (O2), (O3) (O2) (O3) n (O2) U 1,..., U n O U k O k=1 (O3) U λ O( λ Λ) λ Λ U λ O 0 X 0 (O2) n =

More information

6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m f 4

6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m f 4 35-8585 7 8 1 I I 1 1.1 6kg 1m P σ σ P 1 l l λ λ l 1.m 1 6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m

More information

untitled

untitled 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 information

2 1 x 1.1: v mg x (t) = v(t) mv (t) = mg 0 x(0) = x 0 v(0) = v 0 x(t) = x 0 + v 0 t 1 2 gt2 v(t) = v 0 gt t x = x 0 + v2 0 2g v2 2g 1.1 (x, v) θ

2 1 x 1.1: v mg x (t) = v(t) mv (t) = mg 0 x(0) = x 0 v(0) = v 0 x(t) = x 0 + v 0 t 1 2 gt2 v(t) = v 0 gt t x = x 0 + v2 0 2g v2 2g 1.1 (x, v) θ 1 1 1.1 (Isaac Newton, 1642 1727) 1. : 2. ( ) F = ma 3. ; F a 2 t x(t) v(t) = x (t) v (t) = x (t) F 3 3 3 3 3 3 6 1 2 6 12 1 3 1 2 m 2 1 x 1.1: v mg x (t) = v(t) mv (t) = mg 0 x(0) = x 0 v(0) = v 0 x(t)

More information

I ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT

I ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT I (008 4 0 de Broglie (de Broglie p λ k h Planck ( 6.63 0 34 Js p = h λ = k ( h π : Dirac k B Boltzmann (.38 0 3 J/K T U = 3 k BT ( = λ m k B T h m = 0.067m 0 m 0 = 9. 0 3 kg GaAs( a T = 300 K 3 fg 07345

More information

f(x) = f(x ) + α(x)(x x ) α(x) x = x. x = f (y), x = f (y ) y = f f (y) = f f (y ) + α(f (y))(f (y) f (y )) f (y) = f (y ) + α(f (y)) (y y ) ( (2) ) f

f(x) = f(x ) + α(x)(x x ) α(x) x = x. x = f (y), x = f (y ) y = f f (y) = f f (y ) + α(f (y))(f (y) f (y )) f (y) = f (y ) + α(f (y)) (y y ) ( (2) ) f 22 A 3,4 No.3 () (2) (3) (4), (5) (6) (7) (8) () n x = (x,, x n ), = (,, n ), x = ( (x i i ) 2 ) /2 f(x) R n f(x) = f() + i α i (x ) i + o( x ) α,, α n g(x) = o( x )) lim x g(x) x = y = f() + i α i(x )

More information

1 u t = au (finite difference) u t = au Von Neumann

1 u t = au (finite difference) u t = au Von Neumann 1 u t = au 3 1.1 (finite difference)............................. 3 1.2 u t = au.................................. 3 1.3 Von Neumann............... 5 1.4 Von Neumann............... 6 1.5............................

More information

x = a 1 f (a r, a + r) f(a) r a f f(a) 2 2. (a, b) 2 f (a, b) r f(a, b) r (a, b) f f(a, b)

x = a 1 f (a r, a + r) f(a) r a f f(a) 2 2. (a, b) 2 f (a, b) r f(a, b) r (a, b) f f(a, b) 2011 I 2 II III 17, 18, 19 7 7 1 2 2 2 1 2 1 1 1.1.............................. 2 1.2 : 1.................... 4 1.2.1 2............................... 5 1.3 : 2.................... 5 1.3.1 2.....................................

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

(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 information

{ 8. { CHAPTER 8. Å (sampling time) x[k] =x(kå) u(ú) t t + Å (u[k]) x[k + 1] =A d x[k] +B d u[k] (8:) (diãerence equation) A d =e AÅ ; B d = Z Å 0 e A

{ 8. { CHAPTER 8. Å (sampling time) x[k] =x(kå) u(ú) t t + Å (u[k]) x[k + 1] =A d x[k] +B d u[k] (8:) (diãerence equation) A d =e AÅ ; B d = Z Å 0 e A Chapter 8 ( _x = Ax +Bu; y = Cx) 8.1 8.1.1 x(t) =e At x(0) + _x =Ax +Bu (8:1) Z t 0 e A(tÄú) Bu(ú)dú u(t) x(0) t x(t) t =kå (k + 1)Å x(t) x[k + 1] =e AÅ x[k] + Z t+å t { 8.1 { e A(t+ÅÄú) Bu(ú)dú (8:) {

More information

No δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2

No δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2 No.2 1 2 2 δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i δx j (5) δs 2 = δx i δx i + 2 u i δx i δx j = δs 2 + 2s ij δx i δx j

More information

25 7 18 1 1 1.1 v.s............................. 1 1.1.1.................................. 1 1.1.2................................. 1 1.1.3.................................. 3 1.2................... 3

More information

( ) 2.1. C. (1) x 4 dx = 1 5 x5 + C 1 (2) x dx = x 2 dx = x 1 + C = 1 2 x + C xdx (3) = x dx = 3 x C (4) (x + 1) 3 dx = (x 3 + 3x 2 + 3x +

( ) 2.1. C. (1) x 4 dx = 1 5 x5 + C 1 (2) x dx = x 2 dx = x 1 + C = 1 2 x + C xdx (3) = x dx = 3 x C (4) (x + 1) 3 dx = (x 3 + 3x 2 + 3x + (.. C. ( d 5 5 + C ( d d + C + C d ( d + C ( ( + d ( + + + d + + + + C (5 9 + d + d tan + C cos (sin (6 sin d d log sin + C sin + (7 + + d ( + + + + d log( + + + C ( (8 d 7 6 d + 6 + C ( (9 ( d 6 + 8 d

More information