IV (2)

Size: px
Start display at page:

Download "IV (2)"

Transcription

1 COMPUTATIONAL FLUID DYNAMICS (CFD) IV (2) The Analysis of Numerical Schemes (2) 11. Iterative methods for algebraic systems Reima Iwatsu, iwatsu@cck.dendai.ac.jp Winter Semester 2007, Graduate School, Tokyo Denki University 19, F,

2 IV (2)

3 11. Iterative methods for algebraic systems CFD 3D RANS (7 100 = 700 ),,, :, Varga (1962), Dahlquist & Bjork (1974), Saad (2003) Dongarra 3

4 u = f 0 x, y L (1) u = g on the bounradies (2) (u i+1,j 2u ij + u i 1,j ) + (u i,j+1 2u ij + u i 1,j ) = f ij x 2 (3) M j 1 0 y i M 6 x 4

5 U = (u 11, u 21,, u M,1, u 1,2, u 2,2,, u 1,j, u 2,j,, u M,j, ) (4) (i) S SU = F x 2 + G Q (5) (ii) S (, ) N = M 2, I = i + (j 1)M N J=1 u 11. u i1. = x2 f 11. f i1. g 10 + g 01. g i0. s IJ u J = Q I I = 1,, N (7) (6) 5

6 / ( ) u n ij u n+1 un+1 ij (n ) ij = 1 4 (un i+1,j + un i 1,j + un i,j+1 + un i,j 1 ) qn ij (8) u n+1 i = 1 N q n s i s ij u n j (9) ij j=1,j i : Unknown at iteration n + 1 : Known from iteration n 6

7 D:, E:, F : S = D + E + F (10) DU n+1 = Q n (E + F )U n (11) (D + E + F )U = Q (12) - R n R n R(U n ) = SU n + Q n = (D + E + F )U n + Q n (13) D U n = R n (14) U n U n+1 U n (15) 4(u n+1 ij u n ij ) = (un i+1,j + un i 1,j + un i,j+1 + un i,j 1 4un ij ) + qn ij (16) 7

8 u n+1 ij u n+1 I = 1 s II = 1 4 (un+1 i+1,j + un i 1,j + un i,j+1 + un+1 i,j 1 ) qn ij I 1 qi n N s IJ u n+1 J s IJ u n J J=1 J=I+1 (17a) (17b) (D + E)U n+1 = Q n F U n (18) (D + E) U n = R n (19) 6 6 : Unknown at iteration n + 1 : Known from iteration n + 1 : Known from iteration n 8

9 n e n 0, e n = U n Ū Ū : (20) S S = P + A (21) P U n+1 = Q n AU n (22) P U n = R n (23) P : ( (22) ) SU n+1 = Q n or S U n = R n (24) (25) (22) SŪ = (P + A)Ū = Q (25) 9

10 P e n+1 = Ae n (26) e n+1 = (P 1 A)e n = (P 1 A) n e 1 = (1 P 1 S) n e 1 (27) n e n+1 0 P 1 A ( ) G G = 1 P 1 S (28) ρ(g) 1 (29) λ J (G) 1 for all J (30) S P ˆP 1 O(N) ˆP S 10

11 R e R n = Se n (31) 0 0 ( ) R R U n Ū SU 2 x = 10 2 U 10 4 U n+1 = (1 P 1 S)U n P 1 Q = GU n P 1 Q (32) e n+1 = Ge n = G n+1 e 0 (33) n 11

12 G = D 1 (E + F ) = 1 D 1 S G = (D + E) 1 F = 1 (D + E) 1 S, : S S : S S, (2D S) (S 2D) : S S, (2D S) 12

13 n (33) Varga (1962) e n ( e n ) 1/n e 0 (34) e n G n 1/n (35) log s n 1 (10 ) s lim n Gn 1/n = ρ(g) (36) ρ(g) < 1 13 ( ) 1 n s = ρ(g) (37) 10 n 1/ log ρ(g) (38)

14 G (P 1 S ) U n = J=1 (λ J ) n U 0 V (J) J=1 Q J Ω J [1 (λ J ) n ]V (J) (39) V (J) : S, U 0 :, Ω J : S, S, G, (λ J : G ) lim n U n = J=1 U = S 1 Q Q J Ω J V (J) (40) Ū = J=1 U J V (J) (41) SŪ = J=1 Ω J U J V (J) = J=1 Q J V (J) (42) U J = Q J Ω J (43) 14

15 n (e 0J : ) e n = J=1(λ J ) n e 0J V (J) (44) V (J) ( ) R n = J=1 (λ J ) n e 0J SV (J) = J=1(λ J ) n Ω J e 0J V (J) (45) L 2 R n L2 = J=1[(λ J ) n Ω J e 0J ] 2 (46) n (λ J ) (λ J ) 1 R b a n 15

16 G S λ J = λ(ω J ) D G J = 1 D 1 S (47) λ J = 1 1 d J Ω J (S) (48) 16

17 S, G S Ω = 4(sin 2 φ x /2 + sin 2 φ y /2) φ x = lπ/m φ y = mπ/m J = l+(m 1)M l, m = 1,, M (49) G J λ(g J ) = 1 (sin 2 φ x /2 + sin 2 φ y /2) = 1 2 (cos φ x + cos φ y ) (50) (φ x, φ y 0) (λ 1) (Ω- ) x = y ρ(g J ) = cos πm (51) ρ(g J ) 1 π2 2M 2 = 1 π2 2N = 1 O( x2 ) (52) 17

18 (4 e iφ x e iφ y )λ GS = e iφ x + e iφ y (53) 2 (Young, 1971) λ(g GS ) = 1 4 (cos φ x + cos φ y ) 2 = λ(g J ) 2 (54) ρ(g GS ) = ρ(g J ) 2 = cos 2 π/m (55) M ρ(g GS ) 1 π2 M 2 = 1 π2 N 2 1 N/(0.43π 2 ) Varga (1962) (56) 18

19 10.2 U n = U n+1 U n (Frankel & Young, 1950) U n+1 ω U n+1 = ωū n+1 + (1 ω)u n (57) U n = U n+1 U n with U n = ω U n (58) ω 19

20 u n+1 ij = ω 4 (un i+1,j + un i 1,j + un i,j+1 + un i,j 1 + qn ij ) + (1 ω)un ij (59) ( ) DU n+1 = ωq n ω(e+f )U n +(1 ω)du n = ω(su n +Q n )+DU n (60) D U n = ωr n (61) G J (ω) = (1 ω)i + ωg J = 1 ωd 1 S (62) ρ(g J (ω)) 1 ω + ωρ(g J ) < 1 (63) 20

21 0 < ω < ρ(g J ) (64) λ(g J (ω)) = (1 ω) + ωλ(g J ) (65) ω ω = 1/(1 λ(g J )) λ(g J ) = 0 λ 1, ω < 1 λ 1, ω ω max 21

22 ω ω opt 1 + ω opt (1 λ min ) = 1 ω opt (1 λ max ) (66) λ min λ max ω opt = (50) ω opt = (λ min + λ max ) (67) λ max = λ min = cos π/m (68) λ(g J (ω)) Unstable λ min ω opt 1 λ max ω 22

23 : SOR (Successive Overrelaxation) ( ) ū n+1 ij = ω 4 (un i+1,j + un+1 i 1,j + un i,j+1 + un+1 i,j 1 + qn ij ) + (1 ω)un ij u n+1 ij = ωū n+1 ij + (1 ω)u n ij (69) D U n = R n E U n (70) - ( (D + ωe)) (D + ωe) U n = ωr n (71) G SOR (ω) = 1 ω(d + ωe) 1 S = (D + ωe) 1 [(1 ω)d ωf ] (72) 23

24 det G SOR (ω) = det(i + ωd 1 E) 1 det[(1 ω)i ωd 1 F ] = I det[(1 ω)i ωd 1 F ] = (1 ω) N (73) ρ(g SOR (ω)) N det G SOR (ω) = (1 ω) N (74) 1 (1 ω) ρ(g SOR (ω)) < 1 (75) SOR (Young, 1971) 0 < ω < 0 < ω < 2 (76) ρ(g GS ) (77) 24

25 ( ) D S = 1 F F T (78) D 2 D 1, D 2 :, λ(g SOR (ω)) = λ(ω) i) SOR λ(ω) = 1 ω + ωλ 1/2 (ω)λ(g J ) (79) ii) ω = 1 λ(ω = 1) λ(g GS ) = λ 2 (G J ) (80) iii) ω opt = ρ 2 (G J ) (81) ρ(g SOR (ω opt )) = ω opt 1 (82) 25

26 ω opt = ρ(g SOR (ω opt )) = ρ(g J ) = cos π/m (83) sin π/m 2(1 π M + π2 M 2) (84) 1 sin π/m 1 + sin π/m 1 2π M + O( 1 M 2) (85) 1 2.3M/π N SOR k kn N ρ(g) = lim n e n+1 e n n (79) ρ(g J ) (81) ω (86) 26

27 SSOR (Symmetric Successive Overrelaxation) SOR U n+1 SOR (D + ωe) U n+1/2 = ωr n (87) (D + ωf ) U n = ωr(ū n+1/2 ) (88) U n+1 = Ū n+1/2 + Ū n+1/2 = U n + U n = U n + U n+1/2 U n+1/2 + U n (90) SSOR SOR ω opt (1 ρ 2 (G J )) (91) SOR 2 ;2 27

28 SLOR (Successive Line Overrelaxation) U n SLOR (VLOR) ū n+1 ij = ω 4 (un i+1,j + un+1 i 1,j + ūn i,j+1 + ūn i,j qn ij ) u n+1 ij = ωū n+1 ij + (1 ω)u n ij (92) 4 Uij n U i,j 1 n U i,j+1 n ω U i 1,j n = ωrn ij (93) SLOR ω opt 2 6 : Unknown at iteration n + 1 : Known from iteration n : Known from iteration n

29 - : I = i + (j 1)M, : I : Red point : Black point 29

30 10.3 ( -form) n P U n τ = ω(su n + Q n ) = ωr n (94) τ P du dt = ω(su + Q) (95) ( ) P 1 G = 1 ωτp 1 S (96) ρ(g) 1 ( ωτp 1 )S ( ) P = (ωτs), G = 0 1 P/ωτ S P 30

31 P U n τ = ω(su n + Q n ) = ωr n P/τ = 1 ( S ) U n+1 = U n + ω(su n + Q n ) = (1 + ωs)u n + ωq n G R U n + ωq n (97) ( ) G R = 1 + ωs G R S Ω J λ J = 1 + ωω J (98) 0 < ω < 2 = 2 Ω J max ρ(s) (99) 31

32 λ 1 Ω ω 1 max opt ω opt = Unstable Ω min ( ) S ρ(g R ) = 1 2 Ω J max + Ω J min (100) 2 Ω J max = κ(s) 1 Ω J max + Ω J min κ(s) + 1 ω (101) κ(s) = Ω J max Ω J min (102) Ω J max = 8 Ω J min = 8 sin 2 π 2M 2π2 N 32 (103)

33 κ(s) = 1 sin 2 π 4N 2M π 2 (104) M κ(s) 1000 ( ) ρ(g R ) = 1 2π2 N = 1 2 κ(s) (105) ω 1/ Ω J max 1/ Ω J min non-stationary cf. stationary 33

34 ADI (Alternating Direction Implicit) P ADI (1 τs x )(1 τs y )(1 τs z ) = τω(su n + Q n ) (106) (1 τs x ) U = τω(su n + Q n ) (1 τs y ) U = U (1 τs z ) U = U (107) ω, τ ADI N log N ( ) ω opt 2 τ opt = τ 1 ΩJ min Ω J max (108) λ(g ADI ) = (1 4τ 1 sin 2 φ x /2)(1 4τ 1 sin 2 φ y /2) (1 + 4τ 1 sin 2 φ x /2)(1 + 4τ 1 sin 2 φ y /2) (109) 34

35 P P U n τ = ω(su n + Q n ) = ωr n U n+1 = U n τωp 1 (SU n + Q n ) B = τp 1 S G = 1 τp 1 S, G ω opt = λ(g) = 1 ωλ(τp 1 S) (110) 2 λ min (τp 1 S) + λ max (τp 1 S) ρ(g) = κ(τp 1 s) 1 κ(τp 1 s) + 1 (111) = : τp 1 S P 35

36 (Choleski) Meijerink & Van der Vorst U n+1 = U n τωp 1 (SU n + Q n ) (112) Strongly Implicit Procedure (SIP) Stone, Schneider & Zedan λ(g) = 1 ωλ(τp 1 S) (113) CG (Conjugate Gradient) Reid, Concus, Kershaw ω opt = 2 λ min (τp 1 S) + λ max (τp 1 S) GMRES (Generalized Minimum Residual) (114) Saad & Schultz ρ(g) = κ(τp 1 s) 1 κ(τp 1 s) + 1 (115) 36

37 10.4 S(U) = Q (116) S(U n+1 ) = S(U n + U) = S(U n ) + (Jacobian) ( ) J(U) Ū ( ) S U ( ) S U = Q (115) U (116) J(U) U n = R n (117) Ū = U n J(U) 1 R n (118) 37

38 e n = U n Ū J(U)e n = R n (119) (117) J(U) U n = R n P/τ P τ U n = R n (120) e n G e n+1 = U n+1 Ū = en + U = e n τp 1 R n = (1 τp 1 J)e n Ge n P 1 (121) G = 1 τp 1 J (122) 38

39 10.5 Brandt (1972,1977,1982), Thomas et al. (2003), Briggs (2000) = = G 1 = S Ω J min φ x = φ y = π/m S (λ J ) n (G ) λ J = 0.5 n 2 n λ J = Error after smoothing x 39

40 (smoothing) µ = max π/2 φ π : φ x = φ y = π 1 λ(g J ) = 1, +1/2 λ(g) (123) µ(g J (ω)) = max[ 1 2ω, 1 ω/2 ] (124) ω = 4/5 µ =

41 CGC (Coarse Grid Correction) h:, H:, ex. H = 2h (, residual form) S h U h = Q h (125) P U h = R h (126) 1. (restriction) 2. R H = I H h R h (127) S H U H = Q H (128) S H U H = R H (129) 3. U H (prolongation) U H = I h H U H (130) 41

42 Two-grid Iteration Method h, H 1. Uh n S 1 n 1 2. U n+1 h = Uh n + U h 3. U n+1 h S 2 n 2 S n 1 S 2 n 2 6 coarse grid correstion h ( N) N ( N) 42

43 FAS (Full Approximation Scheme) Ih H, Îh H S h (U h ) = Q h (131) S h (U h + U h ) S h (U h ) = R h (132) R H = Ih H R h (133) U H = ÎH h U H (134) S H (U H + U H ) S H (U H ) = R H (135) 43

44 ( ) P G = 1 P 1 S G S O(N) 44

45 n ˆ 1 ˆ 2 n ˆ ( 6 ), I, Toro ˆ Chorin (1967,1968), Harlow & Welch (1965) Patankar & Spalding (1972) 45

untitled

untitled 1 1 Ax = b A R m m A b R m x R m A shift-and invert Lanczos - LU CG A = LU LU Ly = b Ux = y A LU A A = LL T 1 LU b,, Vol. 11, No. 4, pp. 14 18 (2006). x * x (0), x (1), x (2), A Ap A # x (n+1) = Cx (n)

More information

all.dvi

all.dvi 72 9 Hooke,,,. Hooke. 9.1 Hooke 1 Hooke. 1, 1 Hooke. σ, ε, Young. σ ε (9.1), Young. τ γ G τ Gγ (9.2) X 1, X 2. Poisson, Poisson ν. ν ε 22 (9.) ε 11 F F X 2 X 1 9.1: Poisson 9.1. Hooke 7 Young Poisson G

More information

24 I ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x

24 I ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x 24 I 1.1.. ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x 1 (t), x 2 (t),, x n (t)) ( ) ( ), γ : (i) x 1 (t),

More information

量子力学 問題

量子力学 問題 3 : 203 : 0. H = 0 0 2 6 0 () = 6, 2 = 2, 3 = 3 3 H 6 2 3 ϵ,2,3 (2) ψ = (, 2, 3 ) ψ Hψ H (3) P i = i i P P 2 = P 2 P 3 = P 3 P = O, P 2 i = P i (4) P + P 2 + P 3 = E 3 (5) i ϵ ip i H 0 0 (6) R = 0 0 [H,

More information

「産業上利用することができる発明」の審査の運用指針(案)

「産業上利用することができる発明」の審査の運用指針(案) 1 1.... 2 1.1... 2 2.... 4 2.1... 4 3.... 6 4.... 6 1 1 29 1 29 1 1 1. 2 1 1.1 (1) (2) (3) 1 (4) 2 4 1 2 2 3 4 31 12 5 7 2.2 (5) ( a ) ( b ) 1 3 2 ( c ) (6) 2. 2.1 2.1 (1) 4 ( i ) ( ii ) ( iii ) ( iv)

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

( ) sin 1 x, cos 1 x, tan 1 x sin x, cos x, tan x, arcsin x, arccos x, arctan x. π 2 sin 1 x π 2, 0 cos 1 x π, π 2 < tan 1 x < π 2 1 (1) (

( ) sin 1 x, cos 1 x, tan 1 x sin x, cos x, tan x, arcsin x, arccos x, arctan x. π 2 sin 1 x π 2, 0 cos 1 x π, π 2 < tan 1 x < π 2 1 (1) ( 6 20 ( ) sin, cos, tan sin, cos, tan, arcsin, arccos, arctan. π 2 sin π 2, 0 cos π, π 2 < tan < π 2 () ( 2 2 lim 2 ( 2 ) ) 2 = 3 sin (2) lim 5 0 = 2 2 0 0 2 2 3 3 4 5 5 2 5 6 3 5 7 4 5 8 4 9 3 4 a 3 b

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

JFE.dvi

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

numb.dvi

numb.dvi 11 Poisson kanenko@mbkniftycom alexeikanenko@docomonejp http://wwwkanenkocom/ , u = f, ( u = u+f u t, u = f t ) 1 D R 2 L 2 (D) := {f(x,y) f(x,y) 2 dxdy < )} D D f,g L 2 (D) (f,g) := f(x,y)g(x,y)dxdy (L

More information

x A Aω ẋ ẋ 2 + ω 2 x 2 = ω 2 A 2. (ẋ, ωx) ζ ẋ + iωx ζ ζ dζ = ẍ + iωẋ = ẍ + iω(ζ iωx) dt dζ dt iωζ = ẍ + ω2 x (2.1) ζ ζ = Aωe iωt = Aω cos ωt + iaω sin

x A Aω ẋ ẋ 2 + ω 2 x 2 = ω 2 A 2. (ẋ, ωx) ζ ẋ + iωx ζ ζ dζ = ẍ + iωẋ = ẍ + iω(ζ iωx) dt dζ dt iωζ = ẍ + ω2 x (2.1) ζ ζ = Aωe iωt = Aω cos ωt + iaω sin 2 2.1 F (t) 2.1.1 mẍ + kx = F (t). m ẍ + ω 2 x = F (t)/m ω = k/m. 1 : (ẋ, x) x = A sin ωt, ẋ = Aω cos ωt 1 2-1 x A Aω ẋ ẋ 2 + ω 2 x 2 = ω 2 A 2. (ẋ, ωx) ζ ẋ + iωx ζ ζ dζ = ẍ + iωẋ = ẍ + iω(ζ iωx) dt dζ

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

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

m dv = mg + kv2 dt m dv dt = mg k v v m dv dt = mg + kv2 α = mg k v = α 1 e rt 1 + e rt m dv dt = mg + kv2 dv mg + kv 2 = dt m dv α 2 + v 2 = k m dt d

m dv = mg + kv2 dt m dv dt = mg k v v m dv dt = mg + kv2 α = mg k v = α 1 e rt 1 + e rt m dv dt = mg + kv2 dv mg + kv 2 = dt m dv α 2 + v 2 = k m dt d m v = mg + kv m v = mg k v v m v = mg + kv α = mg k v = α e rt + e rt m v = mg + kv v mg + kv = m v α + v = k m v (v α (v + α = k m ˆ ( v α ˆ αk v = m v + α ln v α v + α = αk m t + C v α v + α = e αk m

More information

Chap10.dvi

Chap10.dvi =0. f = 2 +3 { 2 +3 0 2 f = 1 =0 { sin 0 3 f = 1 =0 2 sin 1 0 4 f = 0 =0 { 1 0 5 f = 0 =0 f 3 2 lim = lim 0 0 0 =0 =0. f 0 = 0. 2 =0. 3 4 f 1 lim 0 0 = lim 0 sin 2 cos 1 = lim 0 2 sin = lim =0 0 2 =0.

More information

Untitled

Untitled II 14 14-7-8 8/4 II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ 6/ ] Navier Stokes 3 [ ] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I 1 balance law t (ρv i )+ j

More information

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

20 9 19 1 3 11 1 3 111 3 112 1 4 12 6 121 6 122 7 13 7 131 8 132 10 133 10 134 12 14 13 141 13 142 13 143 15 144 16 145 17 15 19 151 1 19 152 20 2 21 21 21 211 21 212 1 23 213 1 23 214 25 215 31 22 33

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

.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

(Jacobi Gauss-Seidel SOR ) 1. (Theory of Iteration Method) Jacobi Gauss-Seidel SOR 2. Jacobi (Jacobi s Iteration Method) Jacobi 3. Gauss-Seide

(Jacobi Gauss-Seidel SOR ) 1. (Theory of Iteration Method) Jacobi Gauss-Seidel SOR 2. Jacobi (Jacobi s Iteration Method) Jacobi 3. Gauss-Seide 03 9 (Jacobi Gauss-Seidel SOR (Theory of Iteration Method Jacobi Gauss-Seidel SOR Jacobi (Jacobi s Iteration Method Jacobi 3 Gauss-Seidel (Gauss-Seidel Method Gauss-Seidel 4 SOR (SOR Method SOR 9 Ax =

More information

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

More information

Gauss Gauss ɛ 0 E ds = Q (1) xy σ (x, y, z) (2) a ρ(x, y, z) = x 2 + y 2 (r, θ, φ) (1) xy A Gauss ɛ 0 E ds = ɛ 0 EA Q = ρa ɛ 0 EA = ρea E = (ρ/ɛ 0 )e

Gauss Gauss ɛ 0 E ds = Q (1) xy σ (x, y, z) (2) a ρ(x, y, z) = x 2 + y 2 (r, θ, φ) (1) xy A Gauss ɛ 0 E ds = ɛ 0 EA Q = ρa ɛ 0 EA = ρea E = (ρ/ɛ 0 )e 7 -a 7 -a February 4, 2007 1. 2. 3. 4. 1. 2. 3. 1 Gauss Gauss ɛ 0 E ds = Q (1) xy σ (x, y, z) (2) a ρ(x, y, z) = x 2 + y 2 (r, θ, φ) (1) xy A Gauss ɛ 0 E ds = ɛ 0 EA Q = ρa ɛ 0 EA = ρea E = (ρ/ɛ 0 )e z

More information

第1章 微分方程式と近似解法

第1章 微分方程式と近似解法 April 12, 2018 1 / 52 1.1 ( ) 2 / 52 1.2 1.1 1.1: 3 / 52 1.3 Poisson Poisson Poisson 1 d {2, 3} 4 / 52 1 1.3.1 1 u,b b(t,x) u(t,x) x=0 1.1: 1 a x=l 1.1 1 (0, t T ) (0, l) 1 a b : (0, t T ) (0, l) R, u

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

2011de.dvi

2011de.dvi 211 ( 4 2 1. 3 1.1............................... 3 1.2 1- -......................... 13 1.3 2-1 -................... 19 1.4 3- -......................... 29 2. 37 2.1................................ 37

More information

,, 2. Matlab Simulink 2018 PC Matlab Scilab 2

,, 2. Matlab Simulink 2018 PC Matlab Scilab 2 (2018 ) ( -1) TA Email : ohki@i.kyoto-u.ac.jp, ske.ta@bode.amp.i.kyoto-u.ac.jp : 411 : 10 308 1 1 2 2 2.1............................................ 2 2.2..................................................

More information

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

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

日本統計学会誌, 第44巻, 第2号, 251頁-270頁

日本統計学会誌, 第44巻, 第2号, 251頁-270頁 44, 2, 205 3 25 270 Multiple Comparison Procedures for Checking Differences among Sequence of Normal Means with Ordered Restriction Tsunehisa Imada Lee and Spurrier (995) Lee and Spurrier (995) (204) (2006)

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

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

D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j

D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j 6 6.. [, b] [, d] ij P ij ξ ij, η ij f Sf,, {P ij } Sf,, {P ij } k m i j m fξ ij, η ij i i j j i j i m i j k i i j j m i i j j k i i j j kb d {P ij } lim Sf,, {P ij} kb d f, k [, b] [, d] f, d kb d 6..

More information

1 4 1 ( ) ( ) ( ) ( ) () 1 4 2

1 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 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

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

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

微分積分 サンプルページ この本の定価 判型などは, 以下の 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 information

吸収分光.PDF

吸収分光.PDF 3 Rb 1 1 4 1.1 4 1. 4 5.1 5. 5 3 8 3.1 8 4 1 4.1 External Cavity Laser Diode: ECLD 1 4. 1 4.3 Polarization Beam Splitter: PBS 13 4.4 Photo Diode: PD 13 4.5 13 4.6 13 5 Rb 14 6 15 6.1 ECLD 15 6. 15 6.3

More information

( ) e + e ( ) ( ) e + e () ( ) e e Τ ( ) e e ( ) ( ) () () ( ) ( ) ( ) ( )

( ) e + e ( ) ( ) e + e () ( ) e e Τ ( ) e e ( ) ( ) () () ( ) ( ) ( ) ( ) n n (n) (n) (n) (n) n n ( n) n n n n n en1, en ( n) nen1 + nen nen1, nen ( ) e + e ( ) ( ) e + e () ( ) e e Τ ( ) e e ( ) ( ) () () ( ) ( ) ( ) ( ) ( n) Τ n n n ( n) n + n ( n) (n) n + n n n n n n n n

More information

第3章 非線形計画法の基礎

第3章 非線形計画法の基礎 3 February 25, 2009 1 Armijo Wolfe Newton 2 Newton Lagrange Newton 2 SQP 2 1 2.1 ( ) S R n (n N) f (x) : R n x f R x S f (x ) = min x S R n f (x) (nonlinear programming) x 0 S k = 0, 1, 2, h k R n ɛ k

More information

SO(2)

SO(2) TOP URL http://amonphys.web.fc2.com/ 1 12 3 12.1.................................. 3 12.2.......................... 4 12.3............................. 5 12.4 SO(2).................................. 6

More information

chap9.dvi

chap9.dvi 9 AR (i) (ii) MA (iii) (iv) (v) 9.1 2 1 AR 1 9.1.1 S S y j = (α i + β i j) D ij + η j, η j = ρ S η j S + ε j (j =1,,T) (1) i=1 {ε j } i.i.d(,σ 2 ) η j (j ) D ij j i S 1 S =1 D ij =1 S>1 S =4 (1) y j =

More information

第5章 偏微分方程式の境界値問題

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

18 I ( ) (1) I-1,I-2,I-3 (2) (3) I-1 ( ) (100 ) θ ϕ θ ϕ m m l l θ ϕ θ ϕ 2 g (1) (2) 0 (3) θ ϕ (4) (3) θ(t) = A 1 cos(ω 1 t + α 1 ) + A 2 cos(ω 2 t + α

18 I ( ) (1) I-1,I-2,I-3 (2) (3) I-1 ( ) (100 ) θ ϕ θ ϕ m m l l θ ϕ θ ϕ 2 g (1) (2) 0 (3) θ ϕ (4) (3) θ(t) = A 1 cos(ω 1 t + α 1 ) + A 2 cos(ω 2 t + α 18 I ( ) (1) I-1,I-2,I-3 (2) (3) I-1 ( ) (100 ) θ ϕ θ ϕ m m l l θ ϕ θ ϕ 2 g (1) (2) 0 (3) θ ϕ (4) (3) θ(t) = A 1 cos(ω 1 t + α 1 ) + A 2 cos(ω 2 t + α 2 ), ϕ(t) = B 1 cos(ω 1 t + α 1 ) + B 2 cos(ω 2 t

More information

untitled

untitled 18 1 2,000,000 2,000,000 2007 2 2 2008 3 31 (1) 6 JCOSSAR 2007pp.57-642007.6. LCC (1) (2) 2 10mm 1020 14 12 10 8 6 4 40,50,60 2 0 1998 27.5 1995 1960 40 1) 2) 3) LCC LCC LCC 1 1) Vol.42No.5pp.29-322004.5.

More information

4. ϵ(ν, T ) = c 4 u(ν, T ) ϵ(ν, T ) T ν π4 Planck dx = 0 e x 1 15 U(T ) x 3 U(T ) = σt 4 Stefan-Boltzmann σ 2π5 k 4 15c 2 h 3 = W m 2 K 4 5.

4. ϵ(ν, T ) = c 4 u(ν, T ) ϵ(ν, T ) T ν π4 Planck dx = 0 e x 1 15 U(T ) x 3 U(T ) = σt 4 Stefan-Boltzmann σ 2π5 k 4 15c 2 h 3 = W m 2 K 4 5. A 1. Boltzmann Planck u(ν, T )dν = 8πh ν 3 c 3 kt 1 dν h 6.63 10 34 J s Planck k 1.38 10 23 J K 1 Boltzmann u(ν, T ) T ν e hν c = 3 10 8 m s 1 2. Planck λ = c/ν Rayleigh-Jeans u(ν, T )dν = 8πν2 kt dν c

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

: , 2.0, 3.0, 2.0, (%) ( 2.

: , 2.0, 3.0, 2.0, (%) ( 2. 2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................

More information

66 σ σ (8.1) σ = 0 0 σd = 0 (8.2) (8.2) (8.1) E ρ d = 0... d = 0 (8.3) d 1 NN K K 8.1 d σd σd M = σd = E 2 d (8.4) ρ 2 d = I M = EI ρ 1 ρ = M EI ρ EI

66 σ σ (8.1) σ = 0 0 σd = 0 (8.2) (8.2) (8.1) E ρ d = 0... d = 0 (8.3) d 1 NN K K 8.1 d σd σd M = σd = E 2 d (8.4) ρ 2 d = I M = EI ρ 1 ρ = M EI ρ EI 65 8. K 8 8 7 8 K 6 7 8 K 6 M Q σ (6.4) M O ρ dθ D N d N 1 P Q B C (1 + ε)d M N N h 2 h 1 ( ) B (+) M 8.1: σ = E ρ (E, 1/ρ ) (8.1) 66 σ σ (8.1) σ = 0 0 σd = 0 (8.2) (8.2) (8.1) E ρ d = 0... d = 0 (8.3)

More information

kawa (Spin-Orbit Tomography: Kawahara and Fujii 21,Kawahara and Fujii 211,Fujii & Kawahara submitted) 2 van Cittert-Zernike Appendix A V 2

kawa (Spin-Orbit Tomography: Kawahara and Fujii 21,Kawahara and Fujii 211,Fujii & Kawahara submitted) 2 van Cittert-Zernike Appendix A V 2 Hanbury-Brown Twiss (ver. 1.) 24 2 1 1 1 2 2 2.1 van Cittert - Zernike..................................... 2 2.2 mutual coherence................................. 3 3 Hanbury-Brown Twiss ( ) 4 3.1............................................

More information

春期講座 ~ 極限 1 1, 1 2, 1 3, 1 4,, 1 n, n n {a n } n a n α {a n } α {a n } α lim n an = α n a n α α {a n } {a n } {a n } 1. a n = 2 n {a n } 2, 4, 8, 16,

春期講座 ~ 極限 1 1, 1 2, 1 3, 1 4,, 1 n, n n {a n } n a n α {a n } α {a n } α lim n an = α n a n α α {a n } {a n } {a n } 1. a n = 2 n {a n } 2, 4, 8, 16, 春期講座 ~ 極限 1 1, 1 2, 1 3, 1 4,, 1 n, n n {a n } n a n α {a n } α {a n } α lim an = α n a n α α {a n } {a n } {a n } 1. a n = 2 n {a n } 2, 4, 8, 16, 32, n a n {a n } {a n } 2. a n = 10n + 1 {a n } lim an

More information

多体問題

多体問題 Many Body Problem 997 4, 00 4, 004 4............................................................................. 7...................................... 7.............................................

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

006 11 8 0 3 1 5 1.1..................... 5 1......................... 6 1.3.................... 6 1.4.................. 8 1.5................... 8 1.6................... 10 1.6.1......................

More information

v v = v 1 v 2 v 3 (1) R = (R ij ) (2) R (R 1 ) ij = R ji (3) 3 R ij R ik = δ jk (4) i=1 δ ij Kronecker δ ij = { 1 (i = j) 0 (i

v v = v 1 v 2 v 3 (1) R = (R ij ) (2) R (R 1 ) ij = R ji (3) 3 R ij R ik = δ jk (4) i=1 δ ij Kronecker δ ij = { 1 (i = j) 0 (i 1. 1 1.1 1.1.1 1.1.1.1 v v = v 1 v 2 v 3 (1) R = (R ij ) (2) R (R 1 ) ij = R ji (3) R ij R ik = δ jk (4) δ ij Kronecker δ ij = { 1 (i = j) 0 (i j) (5) 1 1.1. v1.1 2011/04/10 1. 1 2 v i = R ij v j (6) [

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc2.com/ 1 30 3 30.1.............. 3 30.2........................... 4 30.3...................... 5 30.4........................ 6 30.5.................................. 8 30.6...............................

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

20 6 4 1 4 1.1 1.................................... 4 1.1.1.................................... 4 1.1.2 1................................ 5 1.2................................... 7 1.2.1....................................

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

1 n A a 11 a 1n A =.. a m1 a mn Ax = λx (1) x n λ (eigenvalue problem) x = 0 ( x 0 ) λ A ( ) λ Ax = λx x Ax = λx y T A = λy T x Ax = λx cx ( 1) 1.1 Th

1 n A a 11 a 1n A =.. a m1 a mn Ax = λx (1) x n λ (eigenvalue problem) x = 0 ( x 0 ) λ A ( ) λ Ax = λx x Ax = λx y T A = λy T x Ax = λx cx ( 1) 1.1 Th 1 n A a 11 a 1n A = a m1 a mn Ax = λx (1) x n λ (eigenvalue problem) x = ( x ) λ A ( ) λ Ax = λx x Ax = λx y T A = λy T x Ax = λx cx ( 1) 11 Th9-1 Ax = λx λe n A = λ a 11 a 12 a 1n a 21 λ a 22 a n1 a n2

More information

ma22-9 u ( v w) = u v w sin θê = v w sin θ u cos φ = = 2.3 ( a b) ( c d) = ( a c)( b d) ( a d)( b c) ( a b) ( c d) = (a 2 b 3 a 3 b 2 )(c 2 d 3 c 3 d

ma22-9 u ( v w) = u v w sin θê = v w sin θ u cos φ = = 2.3 ( a b) ( c d) = ( a c)( b d) ( a d)( b c) ( a b) ( c d) = (a 2 b 3 a 3 b 2 )(c 2 d 3 c 3 d A 2. x F (t) =f sin ωt x(0) = ẋ(0) = 0 ω θ sin θ θ 3! θ3 v = f mω cos ωt x = f mω (t sin ωt) ω t 0 = f ( cos ωt) mω x ma2-2 t ω x f (t mω ω (ωt ) 6 (ωt)3 = f 6m ωt3 2.2 u ( v w) = v ( w u) = w ( u v) ma22-9

More information

重力方向に基づくコントローラの向き決定方法

重力方向に基づくコントローラの向き決定方法 ( ) 2/Sep 09 1 ( ) ( ) 3 2 X w, Y w, Z w +X w = +Y w = +Z w = 1 X c, Y c, Z c X c, Y c, Z c X w, Y w, Z w Y c Z c X c 1: X c, Y c, Z c Kentaro Yamaguchi@bandainamcogames.co.jp 1 M M v 0, v 1, v 2 v 0 v

More information

main.dvi

main.dvi SGC - 70 2, 3 23 ɛ-δ 2.12.8 3 2.92.13 4 2 3 1 2.1 2.102.12 [8][14] [1],[2] [4][7] 2 [4] 1 2009 8 1 1 1.1... 1 1.2... 4 1.3 1... 8 1.4 2... 9 1.5... 12 1.6 1... 16 1.7... 18 1.8... 21 1.9... 23 2 27 2.1

More information

タイの食品市場(国庫用)訂正.PDF

タイの食品市場(国庫用)訂正.PDF 7 2004 ( ) ... 1... 2... 4... 7...11... 14... 16... 23... 24... 27... 29... 34... 39... 39... 43... 53... 58... 61... 62 ... 70... 78... 79... 83... 83... 86... 92... 93... 95... 101... 102... 104... 105...110...113...116...118...119...

More information

4 2016 3 8 2.,. 2. Arakawa Jacobin., 2 Adams-Bashforth. Re = 80, 90, 100.. h l, h/l, Kármán, h/l 0.28,, h/l.., (2010), 46.2., t = 100 t = 2000 46.2 < Re 46.5. 1 1 4 2 6 2.1............................

More information

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.

医系の統計入門第 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

Gmech08.dvi

Gmech08.dvi 145 13 13.1 13.1.1 0 m mg S 13.1 F 13.1 F /m S F F 13.1 F mg S F F mg 13.1: m d2 r 2 = F + F = 0 (13.1) 146 13 F = F (13.2) S S S S S P r S P r r = r 0 + r (13.3) r 0 S S m d2 r 2 = F (13.4) (13.3) d 2

More information

( ) (, ) arxiv: hgm OpenXM search. d n A = (a ij ). A i a i Z d, Z d. i a ij > 0. β N 0 A = N 0 a N 0 a n Z A (β; p) = Au=β,u N n 0 A

( ) (, ) arxiv: hgm OpenXM search. d n A = (a ij ). A i a i Z d, Z d. i a ij > 0. β N 0 A = N 0 a N 0 a n Z A (β; p) = Au=β,u N n 0 A ( ) (, ) arxiv: 1510.02269 hgm OpenXM search. d n A = (a ij ). A i a i Z d, Z d. i a ij > 0. β N 0 A = N 0 a 1 + + N 0 a n Z A (β; p) = Au=β,u N n 0 A-. u! = n i=1 u i!, p u = n i=1 pu i i. Z = Z A Au

More information

main.dvi

main.dvi 6 FIR FIR FIR FIR 6.1 FIR 6.1.1 H(e jω ) H(e jω )= H(e jω ) e jθ(ω) = H(e jω ) (cos θ(ω)+jsin θ(ω)) (6.1) H(e jω ) θ(ω) θ(ω) = KωT, K > 0 (6.2) 6.1.2 6.1 6.1 FIR 123 6.1 H(e jω 1, ω

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

keisoku01.dvi

keisoku01.dvi 2.,, Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 5 Mon, 2006, 401, SAGA, JAPAN Dept.

More information

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

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

= M + M + M + M M + =.,. f = < ρ, > ρ ρ. ρ f. = ρ = = ± = log 4 = = = ± f = k k ρ. k

= M + M + M + M M + =.,. f = < ρ, > ρ ρ. ρ f. = ρ = = ± = log 4 = = = ± f = k k ρ. k 7 b f n f} d = b f n f d,. 5,. [ ] ɛ >, n ɛ + + n < ɛ. m. n m log + < n m. n lim sin kπ sin kπ } k π sin = n n n. k= 4 f, y = r + s, y = rs f rs = f + r + sf y + rsf yy + f y. f = f =, f = sin. 5 f f =.

More information

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) =

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) = 1 1 1.1 I R 1.1.1 c : I R 2 (i) c C (ii) t I c (t) (0, 0) c (t) c(i) c c(t) 1.1.2 (1) (2) (3) (1) r > 0 c : R R 2 : t (r cos t, r sin t) (2) C f : I R c : I R 2 : t (t, f(t)) (3) y = x c : R R 2 : t (t,

More information

gr09.dvi

gr09.dvi .1, θ, ϕ d = A, t dt + B, t dtd + C, t d + D, t dθ +in θdϕ.1.1 t { = f1,t t = f,t { D, t = B, t =.1. t A, tdt e φ,t dt, C, td e λ,t d.1.3,t, t d = e φ,t dt + e λ,t d + dθ +in θdϕ.1.4 { = f1,t t = f,t {

More information

5 c P 5 kn n t π (.5 P 7 MP π (.5 n t n cos π. MP 6 4 t sin π 6 cos π 6.7 MP 4 P P N i i i i N i j F j ii N i i ii F j i i N ii li i F j i ij li i i i

5 c P 5 kn n t π (.5 P 7 MP π (.5 n t n cos π. MP 6 4 t sin π 6 cos π 6.7 MP 4 P P N i i i i N i j F j ii N i i ii F j i i N ii li i F j i ij li i i i i j ij i j ii,, i j ij ij ij (, P P P P θ N θ P P cosθ N F N P cosθ F Psinθ P P F P P θ N P cos θ cos θ cosθ F P sinθ cosθ sinθ cosθ sinθ 5 c P 5 kn n t π (.5 P 7 MP π (.5 n t n cos π. MP 6 4 t sin π 6

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc2.com/ 1 6 3 6.1................................ 3 6.2.............................. 4 6.3................................ 5 6.4.......................... 6 6.5......................

More information

5. [1 ] 1 [], u(x, t) t c u(x, t) x (5.3) ξ x + ct, η x ct (5.4),u(x, t) ξ, η u(ξ, η), ξ t,, ( u(ξ,η) ξ η u(x, t) t ) u(x, t) { ( u(ξ, η) c t ξ ξ { (

5. [1 ] 1 [], u(x, t) t c u(x, t) x (5.3) ξ x + ct, η x ct (5.4),u(x, t) ξ, η u(ξ, η), ξ t,, ( u(ξ,η) ξ η u(x, t) t ) u(x, t) { ( u(ξ, η) c t ξ ξ { ( 5 5.1 [ ] ) d f(t) + a d f(t) + bf(t) : f(t) 1 dt dt ) u(x, t) c u(x, t) : u(x, t) t x : ( ) ) 1 : y + ay, : y + ay + by : ( ) 1 ) : y + ay, : yy + ay 3 ( ): ( ) ) : y + ay, : y + ay b [],,, [ ] au xx

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> 単純適応制御 SAC サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/091961 このサンプルページの内容は, 初版 1 刷発行当時のものです. 1 2 3 4 5 9 10 12 14 15 A B F 6 8 11 13 E 7 C D URL http://www.morikita.co.jp/support

More information

main.dvi

main.dvi 3 4 2 3 4 2 3 ( ) ( ) 1 1.1 R n R n n Euclid R n R n f =(f 1 ;f 2 ;:::;f n ) T : f(x) =0 (1) x = (x 1 ;x 2 ;:::;x n ) T 2 R n 1 \T " f i (1 i n) R n R (1) f i (x) =0 (1 i n) n 2 n x =cosx f(x) :=x 0 cos

More information

Hanbury-Brown Twiss (ver. 2.0) van Cittert - Zernike mutual coherence

Hanbury-Brown Twiss (ver. 2.0) van Cittert - Zernike mutual coherence Hanbury-Brown Twiss (ver. 2.) 25 4 4 1 2 2 2 2.1 van Cittert - Zernike..................................... 2 2.2 mutual coherence................................. 4 3 Hanbury-Brown Twiss ( ) 5 3.1............................................

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

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K II. () 7 F 7 = { 0,, 2, 3, 4, 5, 6 }., F 7 a, b F 7, a b, F 7,. (a) a, b,,. (b) 7., 4 5 = 20 = 2 7 + 6, 4 5 = 6 F 7., F 7,., 0 a F 7, ab = F 7 b F 7. (2) 7, 6 F 6 = { 0,, 2, 3, 4, 5 },,., F 6., 0 0 a F

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

1 nakayama/print/ Def (Definition ) Thm (Theorem ) Prop (Proposition ) Lem (Lemma ) Cor (Corollary ) 1. (1) A, B (2) ABC

1   nakayama/print/ Def (Definition ) Thm (Theorem ) Prop (Proposition ) Lem (Lemma ) Cor (Corollary ) 1. (1) A, B (2) ABC 1 http://www.gem.aoyama.ac.jp/ nakayama/print/ Def (Definition ) Thm (Theorem ) Prop (Proposition ) Lem (Lemma ) Cor (Corollary ) 1. (1) A, B (2) ABC r 1 A B B C C A (1),(2),, (8) A, B, C A,B,C 2 1 ABC

More information

ohpr.dvi

ohpr.dvi 2003/12/04 TASK PAF A. Fukuyama et al., Comp. Phys. Rep. 4(1986) 137 A. Fukuyama et al., Nucl. Fusion 26(1986) 151 TASK/WM MHD ψ θ ϕ ψ θ e 1 = ψ, e 2 = θ, e 3 = ϕ ϕ E = E 1 e 1 + E 2 e 2 + E 3 e 3 J :

More information

SO(3) 7 = = 1 ( r ) + 1 r r r r ( l ) (5.17) l = 1 ( sin θ ) + sin θ θ θ ϕ (5.18) χ(r)ψ(θ, ϕ) l ψ = αψ (5.19) l 1 = i(sin ϕ θ l = i( cos ϕ θ l 3 = i ϕ

SO(3) 7 = = 1 ( r ) + 1 r r r r ( l ) (5.17) l = 1 ( sin θ ) + sin θ θ θ ϕ (5.18) χ(r)ψ(θ, ϕ) l ψ = αψ (5.19) l 1 = i(sin ϕ θ l = i( cos ϕ θ l 3 = i ϕ SO(3) 71 5.7 5.7.1 1 ħ L k l k l k = iϵ kij x i j (5.117) l k SO(3) l z l ± = l 1 ± il = i(y z z y ) ± (z x x z ) = ( x iy) z ± z( x ± i y ) = X ± z ± z (5.118) l z = i(x y y x ) = 1 [(x + iy)( x i y )

More information

1. 2 P 2 (x, y) 2 x y (0, 0) R 2 = {(x, y) x, y R} x, y R P = (x, y) O = (0, 0) OP ( ) OP x x, y y ( ) x v = y ( ) x 2 1 v = P = (x, y) y ( x y ) 2 (x

1. 2 P 2 (x, y) 2 x y (0, 0) R 2 = {(x, y) x, y R} x, y R P = (x, y) O = (0, 0) OP ( ) OP x x, y y ( ) x v = y ( ) x 2 1 v = P = (x, y) y ( x y ) 2 (x . P (, (0, 0 R {(,, R}, R P (, O (0, 0 OP OP, v v P (, ( (, (, { R, R} v (, (, (,, z 3 w z R 3,, z R z n R n.,..., n R n n w, t w ( z z Ke Words:. A P 3 0 B P 0 a. A P b B P 3. A π/90 B a + b c π/ 3. +

More information

1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0

1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0 1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0 0 < t < τ I II 0 No.2 2 C x y x y > 0 x 0 x > b a dx

More information

44 4 I (1) ( ) (10 15 ) ( 17 ) ( 3 1 ) (2)

44 4 I (1) ( ) (10 15 ) ( 17 ) ( 3 1 ) (2) (1) I 44 II 45 III 47 IV 52 44 4 I (1) ( ) 1945 8 9 (10 15 ) ( 17 ) ( 3 1 ) (2) 45 II 1 (3) 511 ( 451 1 ) ( ) 365 1 2 512 1 2 365 1 2 363 2 ( ) 3 ( ) ( 451 2 ( 314 1 ) ( 339 1 4 ) 337 2 3 ) 363 (4) 46

More information

OHP.dvi

OHP.dvi t 0, X X t x t 0 t u u = x X (1) t t 0 u X x O 1 1 t 0 =0 X X +dx t x(x,t) x(x +dx,t). dx dx = x(x +dx,t) x(x,t) (2) dx, dx = F dx (3). F (deformation gradient tensor) t F t 0 dx dx X x O 2 2 F. (det F

More information

i ii i iii iv 1 3 3 10 14 17 17 18 22 23 28 29 31 36 37 39 40 43 48 59 70 75 75 77 90 95 102 107 109 110 118 125 128 130 132 134 48 43 43 51 52 61 61 64 62 124 70 58 3 10 17 29 78 82 85 102 95 109 iii

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

( ) 2002 1 1 1 1.1....................................... 1 1.1.1................................. 1 1.1.2................................. 1 1.1.3................... 3 1.1.4......................................

More information

H 0 H = H 0 + V (t), V (t) = gµ B S α qb e e iωt i t Ψ(t) = [H 0 + V (t)]ψ(t) Φ(t) Ψ(t) = e ih0t Φ(t) H 0 e ih0t Φ(t) + ie ih0t t Φ(t) = [

H 0 H = H 0 + V (t), V (t) = gµ B S α qb e e iωt i t Ψ(t) = [H 0 + V (t)]ψ(t) Φ(t) Ψ(t) = e ih0t Φ(t) H 0 e ih0t Φ(t) + ie ih0t t Φ(t) = [ 3 3. 3.. H H = H + V (t), V (t) = gµ B α B e e iωt i t Ψ(t) = [H + V (t)]ψ(t) Φ(t) Ψ(t) = e iht Φ(t) H e iht Φ(t) + ie iht t Φ(t) = [H + V (t)]e iht Φ(t) Φ(t) i t Φ(t) = V H(t)Φ(t), V H (t) = e iht V (t)e

More information

基礎数学I

基礎数学I I & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............

More information