IV (2)
|
|
- ほだか こいたばし
- 5 years ago
- Views:
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
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 informationall.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 information24 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 informationII ( ) (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) (
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 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 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 informationnumb.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 informationx 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 informationIII 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 informationx () 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 informationm 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 informationChap10.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 informationUntitled
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)
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 information20 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
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
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 informationIA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................
More informationGauss 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章 微分方程式と近似解法
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 informationII 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 information2011de.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
(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 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日本統計学会誌, 第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 information2016 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 informationmeiji_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 informationD = [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 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 informationPart () () Γ 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 informationhttp://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 からご覧いただけます. 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
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 ( ) ( ) () () ( ) ( ) ( ) ( )
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 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 informationSO(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 informationchap9.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章 偏微分方程式の境界値問題
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 information18 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 informationuntitled
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 information4. ϵ(ν, 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 informationI 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.
2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................
More information66 σ σ (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 informationkawa (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 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 informationK E N Z OU
K E N Z OU 11 1 1 1.1..................................... 1.1.1............................ 1.1..................................................................................... 4 1.........................................
More information006 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 informationv 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 informationTOP 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,,.
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 information20 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›ð’Í
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 information1 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 informationma22-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 informationmain.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
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 information4 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 からご覧いただけます. 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 informationGmech08.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: 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 informationmain.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 informationx,, 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 informationkeisoku01.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 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 information= 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 information2 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 informationgr09.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 information5 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 informationTOP 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 information5. [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>
単純適応制御 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 informationmain.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 informationHanbury-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)
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 informationII 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 informationuntitled
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 information1 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 informationohpr.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 informationSO(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 information1. 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 information1 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 information44 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 informationOHP.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 informationi 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 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 informationH 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 & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............
More information