数値計算:有限要素法
|
|
- ありあ ほがり
- 5 years ago
- Views:
Transcription
1 ( ) 1 / 61
2 ( ) 2 / 61
3 ( ) 3 / 61
4 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61
5 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61
6 u(x) 0 L x ( ) 6 / 61
7 P(0) P(L) f d dx ( EA du dx ) = 0 u(0) = 0 E(L)A(L) du dx (L) = f (boundary value problem) ( ) 7 / 61
8 (finite element method; FEM) ( ) ( ) 8 / 61
9 U U = f W L 0 ( ) 2 1 du 2 EA dx, dx W = f u(l) min I = U W, subject to u(0) = 0 ( ) 9 / 61
10 [0, L] 6 h = L/6 Pi Pj 0 xi xj L x (nodal point) x 0 = 0, x 1 = h, x 2 = 2h,, x 6 = L u i = u(xi ) u(x) 7 u 0, u 1, u 2,, u 6 ( ) 10 / 61
11 [x i, x j ] u(x) u(x) = u i N i,j (x) + u j N j,i (x), (x i x x j ) u i, u j x i, x j xi xj L N i,j (x) = x j x { h 1 (x = xi ) = 0 (x = x j ) x 0 0 xi xj L 1 N j,i (x) = x x i { h 1 (x = xj ) = 0 (x = x i ) x ( ) 11 / 61
12 u(x) = u(x) (0 x L) u 0 N 0,1 (x) + u 1 N 1,0 (x) (x 0 x x 1 ) u 1 N 1,2 (x) + u 2 N 2,1 (x) (x 1 x x 2 ) u 2 N 2,3 (x) + u 3 N 3,2 (x) (x 2 x x 3 ). u 5 N 5,6 (x) + u 6 N 6,5 (x) (x 5 x x 6 ) u(x) 7 u 0, u 1,, u 6 ( ) 12 / 61
13 u(x) 0 h 2h 3h 4h 5h L P0 P1 P2 P3 P4 P5 P6 x ( ) 13 / 61
14 u(x) u5 u6 u2 u4 u0 u1 u3 0 h 2h 3h 4h 5h L P0 P1 P2 P3 P4 P5 P6 x ( ) 13 / 61
15 [x i, x j ] N i,j (x) = x j x, N j,i (x) = x x i h h N i,j(x) = 1 h, N j,i(x) = 1 h u(x) = u i N i,j (x) + u j N j,i (x) du = u i N dx i,j(x) + u j N j,i(x) 1 = u i h + u 1 j h = u i + u j h ( ) 14 / 61
16 = 7 u 0, u 1,, u 6 L 0 = u N = u 0 u 1 u 2. u 6 x1 x 0 x2 x 1 x3 + + x x6 x 5 ( ) 15 / 61
17 E xj ( ) 2 1 du xj ( ) 2 x i 2 EA 1 dx = dx x i 2 EA ui + u j dx h = 1 E xj 2 h ( u 2 i + u j ) 2 A dx x i = 1 [ [ ] E V ui u i,j j 2 h 2 V i,j V i,j = xj x i A dx = [x i, x j ] V i,j V i,j ] [ ui u j ], ( ) 16 / 61
18 L ( ) 2 1 du U = 0 2 EA dx dx = 1 [ ] [ [ ] E V u0 u 0,1 V 0,1 u0 1 2 h 2 V 0,1 V 0,1 u 1 [ ] [ 1 [ ] E V u1 u 1,2 V 1,2 u1 2 2 h 2 V 1,2 V 1,2 u 2 [ ] [ 1 [ ] E V u2 u 2,3 V 2,3 u2 3 2 h 2 V 2,3 V 2,3 u 3 [ ] [ 1 [ ] E V u5 u 5,6 V 5,6 u5 6 2 h 2 V 5,6 V 5,6 u 6 ] + ] + ] + + ] ( ) 17 / 61
19 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61
20 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61
21 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61
22 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61
23 A V i,j = Ah K = EA h ( ) 19 / 61
24 A V i,j = Ah K = EA h ( ) 20 / 61
25 A V i,j = Ah K = EA h ( ) 20 / 61
26 f W = f u(l) = f u 6 W = f T u N f = f ( ) 21 / 61
27 u(0) = u 0 = 0 a T u N = 0 a = ( ) 22 / 61
28 min I = U W, subject to u(0) = 0 min I (u N ) = 1 2 ut N K u N f T u N, subject to a T u N = 0 ( ) 23 / 61
29 ( ) min I (u N ) = 1 2 ut N K u N f T u N, subject to a T u N = 0 min J(u N, λ) = I (u N ) λa T u N = 1 2 ut N K u N f T u N λa T u N ( ) 24 / 61
30 ( ) J u N = Ku N f λa = 0 J λ = at u N = 0 [ K a a T 0 ] [ un λ ] = [ f 0 ] ( ) 25 / 61
31 ( ) J u N = Ku N f λa = 0 J λ = at u N = 0 [ K a a T 0 ] [ un λ ] = [ f 0 ] ( ) 25 / 61
32 ( ) J u N = Ku N f λa = 0 J λ = at u N = 0 [ K a a T 0 ] [ un λ ] = [ f 0 ] = u N λ ( ) 25 / 61
33 A(x) = a 2bx (a, b a 2bL > 0 ) E V i,j = xj x i (a 2bx) dx = [ ax bx 2] xj x i = {a b(x i + x j )} h V 0,1 = (a bh)h, V 1,2 = (a 3bh)h, V 2,3 = (a 5bh)h, [0, L] 4 a bh K = E h a + bh a + bh 2a 4bh a + 3bh a + 3bh 2a 8bh a + 5bh a + 5bh 2a 12bh a + 7bh a + 7bh a 7bh ( ) 26 / 61
34 P(x) r(x) = 3(a bx) (a, b a bl > 0 ) A(x) = πr 2 = 3π(a bx) 2 xj ] x=xj V i,j = 3π(a bx) 2 (a bx)3 dx = [π x i b x=x i = π { 3a 2 3ab(x i + x j ) + b 2 (xi 2 + x i x j + xj 2 ) } h [0, L] 4 (x i = ih h = L/4) V 0,1 = π(3a 2 3abh + b 2 h 2 )h V 1,2 = π(3a 2 9abh + 7b 2 h 2 )h V 2,3 = π(3a 2 15abh + 19b 2 h 2 )h V 3,4 = π(3a 2 21abh + 37b 2 h 2 )h ( ) 27 / 61
35 ( ) 28 / 61
36 P(0) P(x) u(x,t) 0 P(L) f(t) P(0) P(x) t P(L) ( ) 29 / 61
37 T = L 0 ( ) 2 1 u L 2 ρa 1 dx = t 0 2 ρa u2 dx 7 u 0, u 1,, u 6 [x i, x j ] u(x, t) u(x, t) = u i (t) N i,j (x) + u j (t) N j,i (x), (x i x x j ) u(x, t) = u i (t) N i,j (x) + u j (t) N j,i (x), (x i x x j ) ( ) 30 / 61
38 xj x i 1 2 ρa u2 dx = 1 2 m i,j; i,j = xj = 1 2 m i,j; j,i = m j,i; i,j = xj ρa { u i (t) N i,j (x) + u j (t) N j,i (x)} 2 dx x i [ ] [ ] [ ] m ui u i,j; i,j m i,j; j,i ui j m j,i; i,j m j,i; j,i u j x i ρa (N i,j ) 2 dx, m j,i; j,i = xj x i ρa N i,j N j,i dx xj x i ρa (N j,i ) 2 dx, ( ) 31 / 61
39 L T = 1 ρa u 2 dx 2 0 = 1 [ ] [ ] [ m u0 u 0,1; 0,1 m 0,1; 1,0 u0 1 2 m 1,0; 0,1 m 1,0; 1,0 u 1 1 [ ] [ ] [ m u1 u 1,2; 1,2 m 1,2; 2,1 u1 2 2 m 2,1; 1,2 m 2,1; 2,1 u 2 1 [ ] [ ] [ m u2 u 2,3; 2,3 m 2,3; 3,2 u2 3 2 m 3,2; 2,3 m 3,2; 3,2 u 3 1 [ ] [ ] [ m u5 u 5,6; 5,6 m 5,6; 6,5 u5 6 2 m 6,5; 5,6 m 6,5; 6,5 u 6 ] + ] + ] + + ] = 1 2 u N T M u N ( ) 32 / 61
40 M = m 0,1; 0,1 m 0,1; 1,0 m 1,0; 0,1 m 1,0; 1,0 + m 1,2; 1,2 m 1,2; 2,1. m 2,1; 1,2.. m4,5; 5,4 m 5,4; 4,5 m 5,4; 5,4 + m 5,6; 5,6 m 5,6; 6,5 m 6,5; 5,6 m 6,5; 6,5 ( ) 33 / 61
41 A ρ xj x i (N i,j ) 2 dx = xj x i N i,j N j,i dx = M = ρah 6 xj x i (N j,i ) 2 dx = h 3 xj x i N j,i N i,j dx = h ( ) 34 / 61
42 A ρ M = ρah ( ) 35 / 61
43 A ρ M = ρah ( ) 36 / 61
44 A ρ M = ρah ( ) 36 / 61
45 A(x) = a 2bx (a, b a 2bL > 0 ) ρ xj x (N i,j ) 2 dx = h x i 12 (3x i + x j ), xj x (N j,i ) 2 dx = h x i 12 (x i + 3x j ), xj x i x N i,j N j,i dx = xj x i x N j,i N i,j dx = h 12 (x i + x j ) m i,j; i,j = ρh 6 {2a b(3x i + x j )}, m j,i; j,i = ρh 6 {2a b(x i + 3x j )}, m i,j; j,i = m j,i; i,j = ρh 6 {a b(x i + x j )} ( ) 37 / 61
46 [0, L] 4 2a bh M = ρh 6 a bh a bh 4a 8bh a 3bh a 3bh 4a 16bh a 5bh a 5bh 4a 24bh a 7bh a 7bh 2a 15bh ( ) 38 / 61
47 L(u N, u N ) = T U + W + λr = 1 2 u N T M u N 1 2 u N T Ku N + f T u N + λa T u N L d L = Ku N + f + λa Mü N = 0 u N dt u N ( ) 39 / 61
48 R = a T u N Ṙ = a T u N R = a T ü N a T ü N + 2νa T u N + ν 2 a T u N = 0 ( ) 40 / 61
49 Mü N λa = Ku N + f a T ü N = a T (2ν u N + ν 2 u N ) u N = v N [ M a a T ] [ vn λ ] [ = Ku N + f a T (2νv N + ν 2 u N ) ] ( ) 41 / 61
50 λ, µ ( ) 42 / 61
51 P i P i P i Pk Pi Pj ( ) 43 / 61
52 P i P i P i Pk Pi Pj ( ) 43 / 61
53 P i P i P i Pk Pi uk ui Pj uj P i P i P i u i = [ ui v i ] [ uj, u j = v j ] [ uk, u k = v k ] ( ) 43 / 61
54 U i,j,k = 1 [ u T 2 i uj T uk T K i,j,k (6 6) ] Ki,j,k u i u j u k K i,j,k = λj i,j,k λ J i,j,k λ, Jµ i,j,k (6 6) + µj i,j,k µ ( ) 44 / 61
55 1 P2 P2 1 1 P3 P0 1 P1 P1 ( ) 45 / 61
56 J 0,1,2 λ = J 3,2,1 λ = 1 2 J 0,1,2 µ = J 3,2,1 µ = ( ) 46 / 61
57 P2 P3 P2 P2 P3 1 = P0 1 P1 P0 1 P1 P1 ( ) 47 / 61
58 U = U 0,1,2 + U 3,2,1 = 1 2 [ u T 0 u T 1 u T 2 u T 3 ] K u 0 u 1 u 2 u 3 K (8 8) J λ, J µ (8 8) K = λj λ + µj µ ( ) 48 / 61
59 J λ (0, 0) J 0,1,2 λ = ( ) 49 / 61
60 J λ (0, 1) J 0,1,2 λ = ( ) 49 / 61
61 J λ (0, 2) J 0,1,2 λ = ( ) 49 / 61
62 J λ (1, 0) J 0,1,2 λ = ( ) 49 / 61
63 J λ (1, 1) J 0,1,2 λ = ( ) 49 / 61
64 J λ (1, 2) J 0,1,2 λ = ( ) 49 / 61
65 J λ (2, 0) J 0,1,2 λ = ( ) 49 / 61
66 J λ (2, 1) J 0,1,2 λ = ( ) 49 / 61
67 J λ (2, 2) J 0,1,2 λ = ( ) 49 / 61
68 J 0,1,2 λ J λ J λ = ( ) 50 / 61
69 J λ (3, 3) J 3,2,1 λ = ( ) 51 / 61
70 J λ (3, 2) J 3,2,1 λ = ( ) 51 / 61
71 J λ (3, 1) J 3,2,1 λ = ( ) 51 / 61
72 J λ (2, 3) J 3,2,1 λ = ( ) 51 / 61
73 J λ (2, 2) J 3,2,1 λ = ( ) 51 / 61
74 J λ (2, 1) J 3,2,1 λ = ( ) 51 / 61
75 J λ (1, 3) J 3,2,1 λ = ( ) 51 / 61
76 J λ (1, 2) J 3,2,1 λ = ( ) 51 / 61
77 J λ (1, 1) J 3,2,1 λ = ( ) 51 / 61
78 J 3,2,1 λ J λ J λ = ( ) 52 / 61
79 J 0,1,2 λ J 3,2,1 λ J λ = ( ) 53 / 61
80 J 0,1,2 µ J µ J µ = ( ) 54 / 61
81 J 3,2,1 µ J µ J µ = ( ) 55 / 61
82 J 0,1,2 µ J 1,3,2 µ J µ = ( ) 56 / 61
83 U = 1 2 u N T Ku N, T = 1 2 u N T M u N = = ( ) 57 / 61
84 x = x 1 x 2 x 3 y x y x y 1 = y x x 2 y x 3 ( ) 58 / 61
85 b = b 1 b 2 b 3 y = b T x = b 1 x 1 + b 2 x 2 + b 3 x 3 y x 1 = b 1, y x 2 = b 2, y x 3 = b 3 ( b T x ) = b x ( ) 59 / 61
86 A = a 11 a 12 a 13 a 12 a 22 a 23 a 13 a 23 a 33 y = 1 2 x T Ax = 1 2 { a11 x a 22 x a 33 x a 12 x 1 x 2 + 2a 13 x 1 x 3 + 2a 23 x 2 x 3 } y x 1 = 1 2 {a 11 2x 1 + 2a 12 x 2 + 2a 13 x 3 } = a 11 x 1 + a 12 x 2 + a 13 x 3 ( ) 60 / 61
87 y x 1 = a 11 x 1 + a 12 x 2 + a 13 x 3 y x 2 = a 12 x 1 + a 22 x 2 + a 23 x 3 y x 3 = a 13 x 1 + a 23 x 2 + a 33 x 3 ( ) 1 x 2 x T Ax = Ax ( ) 61 / 61
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 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 information73
73 74 ( u w + bw) d = Ɣ t tw dɣ u = N u + N u + N 3 u 3 + N 4 u 4 + [K ] {u = {F 75 u δu L σ (L) σ dx σ + dσ x δu b δu + d(δu) ALW W = L b δu dv + Aσ (L)δu(L) δu = (= ) W = A L b δu dx + Aσ (L)δu(L) Aσ
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 informationall.dvi
5,, Euclid.,..,... Euclid,.,.,, e i (i =,, ). 6 x a x e e e x.:,,. a,,. a a = a e + a e + a e = {e, e, e } a (.) = a i e i = a i e i (.) i= {a,a,a } T ( T ),.,,,,. (.),.,...,,. a 0 0 a = a 0 + a + a 0
More informationA 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 informationW 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 informationAHPを用いた大相撲の新しい番付編成
5304050 2008/2/15 1 2008/2/15 2 42 2008/2/15 3 2008/2/15 4 195 2008/2/15 5 2008/2/15 6 i j ij >1 ij ij1/>1 i j i 1 ji 1/ j ij 2008/2/15 7 1 =2.01/=0.5 =1.51/=0.67 2008/2/15 8 1 2008/2/15 9 () u ) i i i
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 information1 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 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 information1 8, : 8.1 1, 2 z = ax + by + c ax by + z c = a b +1 x y z c = 0, (0, 0, c), n = ( a, b, 1). f = n i=1 a ii x 2 i + i<j 2a ij x i x j = ( x, A x), f =
1 8, : 8.1 1, z = ax + by + c ax by + z c = a b +1 x y z c = 0, (0, 0, c), n = ( a, b, 1). f = a ii x i + i
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 informationA
A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................
More information6.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 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 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 informationct_root.dvi
7 1 1 3 ( ) (i) ( ) (ii) ( ) (iii) 1.1 2 1.2 1.3 1.4 1.3 1.4 1.1 1.1.1 ( ) ( ) 1) 1 x x X R X =[0, x] x b(x) b( ) 2 b(0) = 0 x 0 b (0) = + b (x) =0 x 0 b (x) < 0 x 2 ( ) Θ Θ ={θ 0,θ 1 } 0
More informationall.dvi
29 4 Green-Lagrange,,.,,,,,,.,,,,,,,,,, E, σ, ε σ = Eε,,.. 4.1? l, l 1 (l 1 l) ε ε = l 1 l l (4.1) F l l 1 F 30 4 Green-Lagrange Δz Δδ γ = Δδ (4.2) Δz π/2 φ γ = π 2 φ (4.3) γ tan γ γ,sin γ γ ( π ) γ tan
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 informationnakata/nakata.html p.1/20
http://www.me.titech.ac.jp/ nakata/nakata.html p.1/20 1-(a). Faybusovich(1997) Linear systems in Jordan algebras and primal-dual interior-point algorithms,, Euclid Jordan p.2/20 Euclid Jordan V Euclid
More informationばらつき抑制のための確率最適制御
( ) http://wwwhayanuemnagoya-uacjp/ fujimoto/ 2011 3 9 11 ( ) 2011/03/09-11 1 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 2 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 3 / 46 (1/2) r + Controller - u Plant y
More 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 informationkou05.dvi
2 C () 25 1 3 1.1........................................ 3 1.2..................................... 4 1.3..................................... 7 1.3.1................................ 7 1.3.2.................................
More information変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy,
変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, z + dz) Q! (x + d x + u + du, y + dy + v + dv, z +
More information直交座標系の回転
b T.Koama x l x, Lx i ij j j xi i i i, x L T L L, L ± x L T xax axx, ( a a ) i, j ij i j ij ji λ λ + λ + + λ i i i x L T T T x ( L) L T xax T ( T L T ) A( L) T ( LAL T ) T ( L AL) λ ii L AL Λ λi i axx
More information6.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 information1 1 u m (t) u m () exp [ (cπm + (πm κ)t (5). u m (), U(x, ) f(x) m,, (4) U(x, t) Re u k () u m () [ u k () exp(πkx), u k () exp(πkx). f(x) exp[ πmxdx
1 1 1 1 1. U(x, t) U(x, t) + c t x c, κ. (1). κ U(x, t) x. (1) 1, f(x).. U(x, t) U(x, t) + c κ U(x, t), t x x : U(, t) U(1, t) ( x 1), () : U(x, ) f(x). (3) U(x, t). [ U(x, t) Re u k (t) exp(πkx). (4)
More information106 4 4.1 1 25.1 25.4 20.4 17.9 21.2 23.1 26.2 1 24 12 14 18 36 42 24 10 5 15 120 30 15 20 10 25 35 20 18 30 12 4.1 7 min. z = 602.5x 1 + 305.0x 2 + 2
105 4 0 1? 1 LP 0 1 4.1 4.1.1 (intger programming problem) 1 0.5 x 1 = 447.7 448 / / 2 1.1.2 1. 2. 1000 3. 40 4. 20 106 4 4.1 1 25.1 25.4 20.4 17.9 21.2 23.1 26.2 1 24 12 14 18 36 42 24 10 5 15 120 30
More information3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t +
B: 2016 12 2, 9, 16, 2017 1 6 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : t = ν 2 u x 2, (1), c. t + c x = 0, (2). e-mail: iwayama@kobe-u.ac.jp,. 1 3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t +
More information6 2 2 x y x y t P P = P t P = I P P P ( ) ( ) ,, ( ) ( ) cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ y x θ x θ P
6 x x 6.1 t P P = P t P = I P P P 1 0 1 0,, 0 1 0 1 cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ x θ x θ P x P x, P ) = t P x)p ) = t x t P P ) = t x = x, ) 6.1) x = Figure 6.1 Px = x, P=, θ = θ P
More information< 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3)
< 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3) 6 y = g(x) x = 1 g( 1) = 2 ( 1) 3 = 2 ; g 0 ( 1) =
More information1 Ricci V, V i, W f : V W f f(v ) = Imf W ( ) f : V 1 V k W 1
1 Ricci V, V i, W f : V W f f(v = Imf W ( f : V 1 V k W 1 {f(v 1,, v k v i V i } W < Imf > < > f W V, V i, W f : U V L(U; V f : V 1 V r W L(V 1,, V r ; W L(V 1,, V r ; W (f + g(v 1,, v r = f(v 1,, v r
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 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 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 informationsimx simxdx, cosxdx, sixdx 6.3 px m m + pxfxdx = pxf x p xf xdx = pxf x p xf x + p xf xdx 7.4 a m.5 fx simxdx 8 fx fx simxdx = πb m 9 a fxdx = πa a =
II 6 ishimori@phys.titech.ac.jp 6.. 5.4.. f Rx = f Lx = fx fx + lim = lim x x + x x f c = f x + x < c < x x x + lim x x fx fx x x = lim x x f c = f x x < c < x cosmx cosxdx = {cosm x + cosm + x} dx = [
More information1 1.1 / Fik Γ= D n x / Newton Γ= µ vx y / Fouie Q = κ T x 1. fx, tdx t x x + dx f t = D f x 1 fx, t = 1 exp x 4πDt 4Dt lim fx, t =δx 3 t + dxfx, t = 1
1 1.1......... 1............. 1.3... 1.4......... 1.5.............. 1.6................ Bownian Motion.1.......... Einstein.............. 3.3 Einstein........ 3.4..... 3.5 Langevin Eq.... 3.6................
More informationランダムウォークの境界条件・偏微分方程式の数値計算
B L06(2018-05-22 Tue) : Time-stamp: 2018-05-22 Tue 21:53 JST hig,, 2, multiply transf http://hig3.net L06 B(2018) 1 / 38 L05-Q1 Quiz : 1 M λ 1 = 1 u 1 ( ). M u 1 = u 1, u 1 = ( 3 4 ) s (s 0)., u 1 = 1
More information18 ( ) 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,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), Mellor and Yamada 1974) 4 2 3, 2 4,
Mellor and Yamada1974) The Turbulence Closure Model of Mellor and Yamada 1974) Kitamori Taichi 2004/01/30 ,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), 4 1 4 Mellor and Yamada 1974) 4 2 3, 2
More informationf(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第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 informationSturm-Liouville Green KEN ZOU Hermite Legendre Laguerre L L [p(x) d2 dx 2 + q(x) d ] dx + r(x) u(x) = Lu(x) = 0 (1) L = p(x) d2 dx
Sturm-Liouville Green KEN ZOU 2006 4 23 1 Hermite Legendre Lguerre 1 1.1 2 L L p(x) d2 2 + q(x) d + r(x) u(x) = Lu(x) = 0 (1) L = p(x) d2 2 + q(x) d + r(x) (2) L = d2 2 p(x) d q(x) + r(x) (3) 2 (Self-Adjoint
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
[ ] 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 informationD v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco
post glacial rebound 3.1 Viscosity and Newtonian fluid f i = kx i σ ij e kl ideal fluid (1.9) irreversible process e ij u k strain rate tensor (3.1) v i u i / t e ij v F 23 D v D F v/d F v D F η v D (3.2)
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 information211 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 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 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 informationuntitled
- k k k = y. k = ky. y du dx = ε ux ( ) ux ( ) = ax+ b x u() = ; u( ) = AE u() = b= u () = a= ; a= d x du ε x = = = dx dx N = σ da = E ε da = EA ε A x A x x - σ x σ x = Eε x N = EAε x = EA = N = EA k =
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 informationMicrosoft Word - 表紙.docx
黒住英司 [ 著 ] サピエンティア 計量経済学 訂正および練習問題解答 (206/2/2 版 ) 訂正 練習問題解答 3 .69, 3.8 4 (X i X)U i i i (X i μ x )U i ( X μx ) U i. i E [ ] (X i μ x )U i i E[(X i μ x )]E[U i ]0. i V [ ] (X i μ x )U i i 2 i j E [(X i
More information1. 1 BASIC PC BASIC BASIC BASIC Fortran WS PC (1.3) 1 + x 1 x = x = (1.1) 1 + x = (1.2) 1 + x 1 = (1.
Section Title Pages Id 1 3 7239 2 4 7239 3 10 7239 4 8 7244 5 13 7276 6 14 7338 7 8 7338 8 7 7445 9 11 7580 10 10 7590 11 8 7580 12 6 7395 13 z 11 7746 14 13 7753 15 7 7859 16 8 7942 17 8 Id URL http://km.int.oyo.co.jp/showdocumentdetailspage.jsp?documentid=
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 information構造と連続体の力学基礎
II 37 Wabash Avenue Bridge, Illinois 州 Winnipeg にある歩道橋 Esplanade Riel 橋6 6 斜張橋である必要は多分無いと思われる すぐ横に道路用桁橋有り しかも塔基部のレストランは 8 年には営業していなかった 9 9. 9.. () 97 [3] [5] k 9. m w(t) f (t) = f (t) + mg k w(t) Newton
More information名称未設定
! = ( u v w = u i u = u 1 u u 3 u = ( u 1 u = ( u v = u i! 11! 1! 13 % T = $! 1!! 3 ' =! ij #! 31! 3! 33 & 1 0 0%! ij = $ 0 1 0 ' # 0 0 1& # % 1! ijm = $ 1 & % 0 (i, j,m = (1,,3, (, 3,1, (3,1, (i, j,m
More information: 1g99p038-8
16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................
More information23 7 28 i i 1 1 1.1................................... 2 1.2............................... 3 1.2.1.................................... 3 1.2.2............................... 4 1.2.3 SI..............................
More information微分方程式の解を見る
( ) norikazu[at]ms.u-tokyo.ac.jp http://www.infsup.jp/saito/materials/829gairon.pdf 28 2 9 NS ( ) 28 2 9 / 5 Newton 2 3 Navier-Stokes 4 5 6 NS ( ) 28 2 9 2 / 5 2 x 2 + bx + c = x = b ± b 2 4c 2 3 x 3 +
More informationlecture
5 3 3. 9. 4. x, x. 4, f(x, ) :=x x + =4,x,.. 4 (, 3) (, 5) (3, 5), (4, 9) 95 9 (g) 4 6 8 (cm).9 3.8 6. 8. 9.9 Phsics 85 8 75 7 65 7 75 8 85 9 95 Mathematics = ax + b 6 3 (, 3) 3 ( a + b). f(a, b) ={3 (a
More information, = = 7 6 = 42, =
http://www.ss.u-tokai.ac.jp/~mahoro/2016autumn/alg_intro/ 1 1 2016.9.26, http://www.ss.u-tokai.ac.jp/~mahoro/2016autumn/alg_intro/ 1.1 1 214 132 = 28258 2 + 1 + 4 1 + 3 + 2 = 7 6 = 42, 4 + 2 = 6 2 + 8
More information2013 25 9 i 1 1 1.1................................... 1 1.2........................... 2 1.3..................................... 3 1.4..................................... 4 2 6 2.1.................................
More informationSample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1
Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) µ = lim xk( k = autocorrelation function R( t, t + τ) = lim ( ) ( + τ) xk t xk t k = V p o o R p o, o V S M R realization
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 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 informationn ( (
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 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 informationB 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 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 informationIV (2)
COMPUTATIONAL FLUID DYNAMICS (CFD) IV (2) The Analysis of Numerical Schemes (2) 11. Iterative methods for algebraic systems Reima Iwatsu, e-mail : iwatsu@cck.dendai.ac.jp Winter Semester 2007, Graduate
More information3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,.,
B:,, 2017 12 1, 8, 15, 22 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : u t = ν 2 u x 2, (1), c. u t + c u x = 0, (2), ( ). 1 3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2,
More informationii 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 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 informationp = mv p x > h/4π λ = h p m v Ψ 2 Ψ
II p = mv p x > h/4π λ = h p m v Ψ 2 Ψ Ψ Ψ 2 0 x P'(x) m d 2 x = mω 2 x = kx = F(x) dt 2 x = cos(ωt + φ) mω 2 = k ω = m k v = dx = -ωsin(ωt + φ) dt = d 2 x dt 2 0 y v θ P(x,y) θ = ωt + φ ν = ω [Hz] 2π
More informationn Y 1 (x),..., Y n (x) 1 W (Y 1 (x),..., Y n (x)) 0 W (Y 1 (x),..., Y n (x)) = Y 1 (x)... Y n (x) Y 1(x)... Y n(x) (x)... Y n (n 1) (x) Y (n 1)
D d dx 1 1.1 n d n y a 0 dx n + a d n 1 y 1 dx n 1 +... + a dy n 1 dx + a ny = f(x)...(1) dk y dx k = y (k) a 0 y (n) + a 1 y (n 1) +... + a n 1 y + a n y = f(x)...(2) (2) (2) f(x) 0 a 0 y (n) + a 1 y
More informationuntitled
40 4 4.3 I (1) I f (x) C [x 0 x 1 ] f (x) f 1 (x) (3.18) f (x) ; f 1 (x) = 1! f 00 ((x))(x ; x 0 )(x ; x 1 ) (x 0 (x) x 1 ) (4.8) (3.18) x (x) x 0 x 1 Z x1 f (x) dx ; Z x1 f 1 (x) dx = Z x1 = 1 1 Z x1
More informationHolton semigeostrophic semigeostrophic,.., Φ(x, y, z, t) = (p p 0 )/ρ 0, Θ = θ θ 0,,., p 0 (z), θ 0 (z).,,,, Du Dt fv + Φ x Dv Φ + fu +
Holton 9.2.2 semigeostrophic 1 9.2.2 semigeostrophic,.., Φ(x, y, z, t) = (p p 0 )/ρ 0, Θ = θ θ 0,,., p 0 (z), θ 0 (z).,,,, Du Dt fv + Φ x Dv Φ + fu + Dt DΘ Dt + w dθ 0 dz = 0, (9.2) = 0, (9.3) = 0, (9.4)
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 informationy π π 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( ) 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 information05-5.dvi
131 71 71 71 71 71 71 71 7 1 71 71 71 71 71 7 1 71 71 71 71 71 71 7 1 71 71 71 71 71 71 71 71 7 1 1 71 71 71 71 71 71 71 71 71 7 1 75(468) 1 71 71 7 517.95 1 7.1 7. 1 71 71 71 71 71 71 71 7, 1 7.1 7. 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 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 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 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 informationtokei01.dvi
2. :,,,. :.... Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 3. (probability),, 1. : : n, α A, A a/n. :, p, p Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN
More information) ] [ h m x + y + + V x) φ = Eφ 1) z E = i h t 13) x << 1) N n n= = N N + 1) 14) N n n= = N N + 1)N + 1) 6 15) N n 3 n= = 1 4 N N + 1) 16) N n 4
1. k λ ν ω T v p v g k = π λ ω = πν = π T v p = λν = ω k v g = dω dk 1) ) 3) 4). p = hk = h λ 5) E = hν = hω 6) h = h π 7) h =6.6618 1 34 J sec) hc=197.3 MeV fm = 197.3 kev pm= 197.3 ev nm = 1.97 1 3 ev
More informationIA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................
More informationPart y mx + n mt + n m 1 mt n + n t m 2 t + mn 0 t m 0 n 18 y n n a 7 3 ; x α α 1 7α +t t 3 4α + 3t t x α x α y mx + n
Part2 47 Example 161 93 1 T a a 2 M 1 a 1 T a 2 a Point 1 T L L L T T L L T L L L T T L L T detm a 1 aa 2 a 1 2 + 1 > 0 11 T T x x M λ 12 y y x y λ 2 a + 1λ + a 2 2a + 2 0 13 D D a + 1 2 4a 2 2a + 2 a
More information() Remrk I = [0, ] [x i, x i ]. (x : ) f(x) = 0 (x : ) ξ i, (f) = f(ξ i )(x i x i ) = (x i x i ) = ξ i, (f) = f(ξ i )(x i x i ) = 0 (f) 0.
() 6 f(x) [, b] 6. Riemnn [, b] f(x) S f(x) [, b] (Riemnn) = x 0 < x < x < < x n = b. I = [, b] = {x,, x n } mx(x i x i ) =. i [x i, x i ] ξ i n (f) = f(ξ i )(x i x i ) i=. (ξ i ) (f) 0( ), ξ i, S, ε >
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 informationMicrosoft 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 information2 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a 2 + b 2 α (norm) N(α) = a 2 + b 2 = αα = α 2 α (spure) (trace) 1 1. a R aα =
1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a + b α (norm) N(α) = a + b = αα = α α (spure) (trace) 1 1. a R aα = aα. α = α 3. α + β = α + β 4. αβ = αβ 5. β 0 6. α = α ( ) α = α
More informationdifgeo1.dvi
1 http://matlab0.hwe.oita-u.ac.jp/ matsuo/difgeo.pdf ver.1 8//001 1 1.1 a A. O 1 e 1 ; e ; e e 1 ; e ; e x 1 ;x ;x e 1 ; e ; e X x x x 1 ;x ;x X (x 1 ;x ;x ) 1 1 x x X e e 1 O e x x 1 x x = x 1 e 1 + x
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 informationf : R R f(x, y) = x + y axy f = 0, x + y axy = 0 y 直線 x+y+a=0 に漸近し 原点で交叉する美しい形をしている x +y axy=0 X+Y+a=0 o x t x = at 1 + t, y = at (a > 0) 1 + t f(x, y
017 8 10 f : R R f(x) = x n + x n 1 + 1, f(x) = sin 1, log x x n m :f : R n R m z = f(x, y) R R R R, R R R n R m R n R m R n R m f : R R f (x) = lim h 0 f(x + h) f(x) h f : R n R m m n M Jacobi( ) m n
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 informationIntroduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math))
Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math)) 2001 1 e-mail:s00x0427@ip.media.kyoto-u.ac.jp 1 1 Van der Pol 1 1 2 2 Bergers 2 KdV 2 1 5 1.1........................................
More informationA
A04-164 2008 2 13 1 4 1.1.......................................... 4 1.2..................................... 4 1.3..................................... 4 1.4..................................... 5 2
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 information