倒立振子で学ぶ制御工学 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

Size: px
Start display at page:

Download "倒立振子で学ぶ制御工学 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです."

Transcription

1

2 倒立振子で学ぶ制御工学 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

3

4 i (PID ) MATLAB/Simulink book_ip_page.html ()

5 ii I () II () //

6 iii MATLAB/Simulink MATLAB/Simulink ( mfiles.zip ) mfiles p1c2 2.. p1c8 8 p2c1 1 p2c2 2 p2c3 3 (3.2 ) AutoGenU_InvPend 3 (3.3 ) page.html (Windows R2013a R2016a ) ip_toolbox_1.0.2.zip ip_toolbox_1.0.2 iptools odqlab_2.1.3 cdip_sample adip_sample / M Simulink ODQ Toolbox/Lab ( 2.3 ) C hoge ip_toolbox_1.0.2 iptools odqlab_2.1.3 >> addpath( C:Y=hogeY=ip_toolbox_1.0.2Y=iptools ) >> addpath( C:Y=hogeY=ip_toolbox_1.0.2Y=odqlab_2.1.3 ) iptools 3.4 (p. 56) ODQ Toolbox/Lab 2.3 (p. 173) 1 (p. 135) 2.3 (p. 173) SeDuMi YALMIP

7 iv f(s) =Lˆf(t) f(t) =L 1ˆf(s) f(z) =Zˆf[k] f(t) Laplace Laplace t s f[k] :=f(kt s) (k =0, 1,...) Z R R n n R m n m n C C n n C m n m n I (I n) (n n ) 0(0 m n) (m n ) M M M 1 M M + (m n ) M M 1 2 M M =(M 1 2 ) M 1 2 M M rank(m) M () tr(m) M diag{a 1,..., a n} M 0 M (M ) M 0 M (M ) M 0 M (M ) M 0 M (M ) M N M N Kronecker x x Euclid x x M F M Frobenius M M j (j 2 = 1) Reˆλ λ = α + jβ α Imˆλ λ = α + jβ β min f f () max f f () inf f f (infimum) () sup f f (supremum) () subject to

8 v I 1 1 () D/A DC PID () PID P PD PI PID PID P P D I PD PID ()

9 vi MATLAB/Simulink () () () ()

10 vii 8 () Z H II LMI () LMI LMI () ()

11 viii

12 I I PID

13

14 3 1 () ( 1.1) 1) 2) 3) (Segway) 4) (1.3 (a)) ibot 5) 6) (1.3 (b)) 1.2

15 4 1 (a) 1.3 (b) () LEGO MINDSTORMS NXT/EV3 7) 11) 12) 17) Lagrange () PID

16 ) 15) 13) 12), 13) 14), 15) 13) 20) 1.4 ( ) 1.4 Furuta Pendulum 12) 14), 16) 1.5 LEGO MINDSTORMS NXT/EV3 LEGO PowerFunctions XL 18) mindsensors.com GlideWheel-M 19) 10), 11) 1.5 LEGO MINDSTORMS 21) 1.6 Pendulum Robot ( ) Pendubot ), 14)

17

18 37 3 () () () () 1) 4) 5) 8) 2) 4),9),10) (p. 8) Newton-Euler Lagrange 5) 7) () (a) Newton-Euler Newton-Euler 5) 7) Newton-Euler 3.1 F (t) =M z(t) (3.1) T (t) =J θ(t) (3.2)

19 F (t) [N] T (t) [N m] ( 1) M [kg] ( 2) J [kg m 2 ] ( 3) z(t) [m] θ(t) [rad] m c [kg] μ c [kg/s] m p [kg] μ p [kg m 2 /s] J p [kg m 2 ] l p [m] g [m/s 2 ] H(t) V (t) 1), 3), 6) H(t) V (t) f c (t) μ c ż(t) m c z(t) =f c (t) μ c ż(t) H(t) (3.3) ( 4) { { m p Ẍ 2 (t) =H(t) m p Ÿ 2 (t) =V (t) m p g, X 2 (t) =z(t)+l p sin θ(t) Y 2 (t) = l p cos θ(t) (3.4) μ p θ(t) ( 1) ( 2) M ( 3) 2l [m] m [kg] J =(4/3)ml 2 [kg m 2 ] ( 4)

20 J p θ(t) = μp θ(t)+v (t)sinθ(t) lp H(t)cosθ(t) l p (3.5) (3.3) (3.5) H(t), V (t) (m c + m p ) z(t)+m p l p cos θ(t) θ(t) = μ c ż(t)+m p l p θ(t) 2 sin θ(t)+f c (t) (3.6) m p l p cos θ(t) z(t)+(j p + m p lp) θ(t) 2 = μ p θ(t)+mp gl p sin θ(t) (3.7) (b) Lagrange Lagrange 5), 7), 8) 3.3 Maple, Mathematica, MATLAB/Symbolic Math Toolbox 11) W (t) U(t) D(t) W (t) = {}}{ 1 2 m cż(t) 2 + U(t) =m p gy 2 (t) }{{} {}}{ 1 2 m pẋ2(t) 2 +, D(t) = 1 2 μ cż(t) 2 } {{ } {}}{ 1 2 m pẏ2(t) μ p θ(t) 2 }{{} {}}{ 1 2 J p θ(t) 2 (3.8) (3.4) { Ẋ2 (t) =ż(t)+l p θ(t)cosθ(t) Ẏ 2 (t) = l p θ(t)sinθ(t) (3.9) q(t) = [ q 1 (t) q 2 (t) ] = [ z(t) θ(t) ] f(t) = [ f 1 (t) f 2 (t) ] = [ f c (t) 0 ] 3.3 (k z, k θ ) (h(t)) (μ z, μ θ ) Mż(t)2 2 J θ(t) kzz(t)2 2 k θθ(t) 2 Mgh(t) 1 2 μzż(t)2 1 2 μ θ θ(t) 2

21 40 3 Lagrangian L(t) =W (t) U(t) D(t) Lagrange Lagrange ( ) d L(t) L(t) dt q i (t) q i (t) + D(t) q i (t) = f i(t) (i =1, 2) (3.10) (3.6), (3.7) MATLAB/Symbolic Math Toolbox Lagrange (3.10) (3.6), (3.7) M M p1c311_cdip_lagrange.m 1 clear; format compact 2 3 syms m_c mu_c real m c, μ c 4 syms J_p m_p mu_p g l_p real J p, m p, μ p, g, l p 5 syms z th dz dth ddz ddth fc real z, θ, ż, θ, z, θ, f c 6 7 q = [ z th ] ; q =[q 1 q 2 ] =[z θ ] 8 dq = [ dz dth ] ; q =[ q 1 q 2 ] =[ż θ ] 9 ddq = [ ddz ddth ] ; q =[ q 1 q 2 ] =[ z θ ] 10 f = [ fc 0 ] ; f =[f 1 f 2 ] =[f c 0 ] 11 % X2 = q(1) + l_p*sin(q(2)); X 2 = z + l p sin θ 13 Y2 = l_p*cos(q(2)); Y 2 = l p cos θ dx2 = diff(x2,q(1))*dq(1)... Ẋ 2 = X2 z ż + X2 θ θ 16 + diff(x2,q(2))*dq(2); 17 dy2 = diff(y2,q(1))*dq(1)... Ẏ 2 = Y2 z ż + Y2 θ θ 18 + diff(y2,q(2))*dq(2); 19 % W = (1/2)*m_c*dq(1)^2... W = 1 2 mcż mpẋ (1/2)*m_p*dX2^ (1/2)*m_p*dY2^ mpẏ Jp θ (1/2)*J_p*dq(2)^2; 24 U = m_p*g*y2; U = m pgy 2 25 D = (1/2)*mu_c*dq(1)^ (1/2)*mu_p*dq(2)^2; 27 D = 1 2 μcż μp θ 2 28 L = W - U; L = W U 29 % N = length(q); N =2q 31 for i = 1:N 32 dlq(i) = diff(l,dq(i)); L qi = L q i 33 (i =1,,N) 34 temp = 0; 35 for j = 1:N «d L NX Lqi 36 temp = temp + diff(dlq(i),dq(j))*ddq(j)... = q j + Lq i dt q i q 37 + diff(dlq(i),q(j))*dq(j); j=1 j q j 38 end 39 ddlq(i) = temp; 40 «d L 41 eq(i) = ddlq(i) - diff(l,q(i))... L + D f i dt q i q i q i 42 + diff(d,dq(i)) - f(i); 43 end q j «

22 eq = simplify(eq ) eq M p1c311_cdip_lagrange.m eq = - l_p*m_p*sin(th)*dth^2 - fc + dz*mu_c + ddz*(m_c + m_p) + ddth*l_p*m_p*cos(th) (3.6) J_p*ddth + dth*mu_p + ddth*l_p^2*m_p + ddz*l_p*m_p*cos(th) - g*l_p*m_p*sin(th) (3.7) f c (t) v(t) (a) DC 3.2 DC di v a (t) =R a i a (t)+l a (t) a + e b (t), dt e b (t) =k b θm (t) (3.11) J m θm (t) =τ m (t) τ L (t) μ m θm (t), τ m (t) =k t i a (t) (3.12) 8), 12) DC 3.4 L a di a (t)/dt 0 (3.11), (3.12) DC J m θm (t) = μ m θm (t)+ k t R a v a (t) τ L (t), μ m = μ m + k tk b R a (3.13) 3.2 DC 3.4 DC v a(t) i a(t) R a L a e b (t) k b θ m(t) τ m(t) J m μ m τ L (t) k t (b) DC 1.12 (p. 8) DC DC

23

24 II I LMI II LMI 1 2 3

25

26 135 1 LMI (LMI: linear matrix inequality) LMI LMI LMI ( ) 1), 2) y(s) =P(s)u(s), P(s) = ω 2 n s 2 +2ζω n s + ω 2 n (1.1) ω n > 0 0 <ζ<1 2 u(t) =1(t 0) y(t) % 3/(ζω n ) T s 3/(ζω n ) <T s

27 136 1 LMI (0 0 <ζ<1) ζω n > 3 T s (1.2) ζ> 1 2 (1.3) <ζ<1 λ = ζω n ± jω n 1 ζ 2 λ λ A(α) := { λ C Re[λ] < α } (1.4) A(α) α- ( 1.3) α (1.2) 5%α =3/T s θ ζ >cos θ tan θ> Im[λ] Re[λ] = 1 ζ 2 λ S(k) := { λ C Im[λ] <k Re[λ], Re[λ] < 0, k =tanθ>0 } (1.5) S(k) k ( 1.4) (1.3) ζ >cos θ =1/ 2 k =1 ζ

28 A(α)α- 1.4 S(k) (k k =tanθ) 1.5 B(β) ( 1.5) λ B(β) := { λ C Re[λ] > β } (1.6) ẋ(t) =Ax(t)+Bu(t), x(0) = x 0 (1.7) y(t) =Cx(t) (1.8) x(t) R n u(t) R m y(t) R p (A, B, C) A λ (1.7) (1.4) (1.5) (LMI: linear matrix inequality) ( 1) (1.4) LMI 3), 4) 1.1 A(α) (α-) LMI A R n n α R (i) (ii) (i) (ii) A λ λ A(α) ( 1.3) LMI n P 0 A P + PA+2αP 0 (1.9) A λ α- LMI P 0 (1.9) P ( 1) LMI ( ) LMI

29 138 1 LMI LMI P A α- 1.1 LMI A P, α LMI A A(α) (α-) [ ] [ ] 0 1 ξ 1 ξ 2 A =, P =, α =1 (T s = 3) (1.10) ξ 2 1 A λ =( 9 ± 3 7j)/8 Re[λ] = 9/8 α = 1 (A A(α) ) A A(α) LMI (P 0 (1.9) ) ( 2) [ ] P 0 = ξ 1 ξ 2 0 ξ 2 1 [ ] (1.11) A P + PA+2αP 0 = 8ξ 1 +18ξ 2 9 4ξ 1 + ξ 2 9 4ξ 1 + ξ ξ 1 0 (1.12) LMI (ξ 1,ξ 2 ) LMI LMI LMI» ( 2) ξ1 ξ 2 ξ 3 P = ξ 3 =1 ξ 2 ξ 3 α-

30 (ξ 1,ξ 2 ) λ α- LMI SeDuMi (Ver ( )) LMI YALMIP (Release ) ( 3) P 0 (1.9) ((1.11), (1.12) ) LMI M M p2c112_ex1_alpha_stability.m 1 clear; format compact 2 % A = [ ]; 5 n = length(a); 6 alpha = 1; 7 ep = 1e-6; 8 % xi1 = sdpvar(1); xi2 = sdpvar(1); 10 P = [ xi1 xi2 ξ 1, ξ 2 P = 11 xi2 1 ]; 12 % LMI = [ P >= ep*eye(n) ]; 14 LMI = LMI + [ A *P + P*A + 2*alpha*P <= - ep*eye(n) ]; 15 % solvesdp(lmi) 17 % P = double(p)» A = nx(t) α =1 ε =10 6» ξ1 ξ 2 ξ 2 1 P εi ( 0) ( 4) A P + PA+2αP εi ( 0) LMI P YALMIP LMI ε>0 F (ξ) 0 F (ξ) εi ( 0) (1.13) M p2c112_ex1_alpha_stability.m SeDuMi 1.32 by AdvOL, and Jos F. Sturm, Alg = 2: xz-corrector, theta = 0.250, beta = eqs m = 2, order n = 5, dim = 9, blocks = 3 nnz(a) = 7 + 0, nnz(ada) = 4, nnz(l) = 3 it : b*y gap delta rate t/tp* t/td* feas cg cg prec 0 : 9.31E : 0.00E E E-10 ( 3) SeDuMi YALMIP johanl/yalmip/ SeDuMi YALMIP http: // ( 4) ξ1 ξ 2 P = 9 11 P = sdpvar(n,n, sy ); ( sy ξ 2 ξ 3 )

31

32 223 A Ackermann 87 A/D 169, 173 AutoGenU 220 B Butterworth 28 C Cayley-Hamilton 81 Cholesky 73 D D/A 7, 9, 11, 46, 120, 169, 173 DC 12, 41 DC 41 Dirac 65 DMC 216 D 21, 23, 25, 26 E Euclid 51, 68 F Frobenius 144 F/V 13 G GPC 216 H Hankel 181 Hankel 63 Hurwitz 71 H 13 H 195 I inf 143, 145 ITAE 28 I 21, 25, 26 I PD 32 I PD J Jordan () 88 K Kronecker 160 L Lagrange 39 Lagrangian 40 Laplace 64 Lie 84 LMI 137, 141 LPV , , PDLMI LPV 158 LMI 161 LTI Luenberger 109 Lyapunov 69 () , 89, 203, , , 72 M Möbius 127 MPC M 52 N Newton-Euler 37 Nyquist 69 O ODQ Toolbox 186 P PDLMI PD 23 PD PFC 216 PID 21, 22 PID 20, PI 25, 96 PI , 99 PWM () 12 P 22, 42, 45, 95 P 22, , 98 P 21, 25 P D 31 P D 30, R Receding Horizon Riccati 89, 101, 143, 205 Routh-Hurwitz 69 S Schur 170 Schur 144 SeDuMi 139, 186 SOS 161 sup 143, 145 T Tustin 126, 127

33 224 Y YALMIP 139, 161 Youla 113 Z Z α- 136 LMI , () Lyapunov 69 Nyquist 69 Routh-Hurwitz , , , 54 39, , 199 Lagrange 40, 54 Newton-Euler Luenberger , , 143, , 105 () , 25, 26 PI 99 P , 9, 10, 169, 173 inf 63, 89, 100, 106, , () , 146, 163 Z , 118 () , 76, 85, 100, 112, , 77 81, 82, , 80 () () PD 23 PID 25 PI 25 P 22 P D , , , 72, , 85, 108, 137 Ackermann LMI 146 PDLMI LMI 137, , 47, LMI , , 47, , , 45, , 143, 205, 214, 217 LMI 145

34 225 16, 88, 101, 143, 147, 152, 164, () 169, , , 11, 119, , sup , () , () , 71, , 71, , , 147, I 136 LMI , 120, , , 83 69, 70, 71 () , 43, 55, 92 LMI 44, , 161 LPV 126, , , , , , I PD 32 PD 23 PID 25 PI 25, 96, 99 P 22, 29, 45, 65, 95, 98 P D 30 27, , 146 Dirac , 51 27, I PD 32 PD 23 PID 25 PI 25 P 23 P D 31 21, 22, 29, 32, 93, 95, 99, 102 PI 96, 99

35 226 P 95, 98 9, 11 Dirac () , 175, 181 5, 14, 78, 87, 90, 146, 151, 156, 162, , 7, 103, 113, 124, 130, 182, 193, 207 7, , 98, , , LMI PDLMI 16, , 71 69, , , 71 70, H 201, , 66, 76, PD P D D P PID I PD , 85, 94 33, 84, , 80 24, , 71, , 71, , 68, 92, , , 12, , , , , 169, 187 Riccati , , , , , 10, 24, 50, 107, , 24, PDLMI , 201, 222

36 SVBL Carnegie Mellon University 2007/102009/

37 i

38 c FAX PrintedinJapan ISBN

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

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

More information

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

<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

analog-control-mod : 2007/2/4(8:44) 2 E8 P M () r e K P ( ) T I u K M T M K D E8.: DC PID K D E8. (E8.) P M () E8.2 K P D () ( T ) (E8.2) K M T M K, T

analog-control-mod : 2007/2/4(8:44) 2 E8 P M () r e K P ( ) T I u K M T M K D E8.: DC PID K D E8. (E8.) P M () E8.2 K P D () ( T ) (E8.2) K M T M K, T analog-control-mod : 2007/2/4(8:44) E8 E8. PID DC. PID 2. DC PID 3. E8.2 DC PID C8 E8. DC PID E6 DC P M () K M ( T M ) (E8.) DC PID C8 E8. r e u E8.2 PID E8. PID analog-control-mod : 2007/2/4(8:44) 2 E8

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

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

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

More information

<4D F736F F D B B BB2D834A836F815B82D082C88C602E646F63>

<4D F736F F D B B BB2D834A836F815B82D082C88C602E646F63> 入門モーター工学 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/074351 このサンプルページの内容は, 初版 1 刷発行当時のものです. 10 kw 21 20 50 2 20 IGBT IGBT IGBT 21 (1) 1 2 (2) (3) ii 20 2013 2 iii iv...

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

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

More information

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

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

More information

1 [ 1] (1) MKS? (2) MKS? [ 2] (1) (42.195k) k 2 (2) (3) k/hr [ 3] t = 0 10 ( 1 velocity [/s] 8 4 O

1 [ 1] (1) MKS? (2) MKS? [ 2] (1) (42.195k) k 2 (2) (3) k/hr [ 3] t = 0 10 ( 1 velocity [/s] 8 4 O : 2014 4 10 1 2 2 3 2.1...................................... 3 2.2....................................... 4 2.3....................................... 4 2.4................................ 5 2.5 Free-Body

More information

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10%

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10% 1 2006.4.17. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 1. 1. 2. 3. 4. 5. 2. ɛ-δ 1. ɛ-n

More information

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 3 版 1 刷発行時のものです. 最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/052093 このサンプルページの内容は, 第 3 版 1 刷発行時のものです. i 3 10 3 2000 2007 26 8 2 SI SI 20 1996 2000 SI 15 3 ii 1 56 6

More information

(iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y = 0., y x, y = x. (v) 1x = x. (vii) (α + β)x = αx + βx. (viii) (αβ)x = α(βx)., V, C.,,., (1)

(iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y = 0., y x, y = x. (v) 1x = x. (vii) (α + β)x = αx + βx. (viii) (αβ)x = α(βx)., V, C.,,., (1) 1. 1.1...,. 1.1.1 V, V x, y, x y x + y x + y V,, V x α, αx αx V,, (i) (viii) : x, y, z V, α, β C, (i) x + y = y + x. (ii) (x + y) + z = x + (y + z). 1 (iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y

More information

1 1.1 [ 1] velocity [/s] 8 4 (1) MKS? (2) MKS? 1.2 [ 2] (1) (42.195k) k 2 (2) (3) k/hr [ 3] t = 0

1 1.1 [ 1] velocity [/s] 8 4 (1) MKS? (2) MKS? 1.2 [ 2] (1) (42.195k) k 2 (2) (3) k/hr [ 3] t = 0 : 2016 4 1 1 2 1.1......................................... 2 1.2................................... 2 2 2 2.1........................................ 2 2.2......................................... 3 2.3.........................................

More information

基礎から学ぶトラヒック理論 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

基礎から学ぶトラヒック理論 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. 基礎から学ぶトラヒック理論 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/085221 このサンプルページの内容は, 初版 1 刷発行時のものです. i +α 3 1 2 4 5 1 2 ii 3 4 5 6 7 8 9 9.3 2014 6 iii 1 1 2 5 2.1 5 2.2 7

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

08-Note2-web

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

More information

* 1 2014 7 8 *1 iii 1. Newton 1 1.1 Newton........................... 1 1.2............................. 4 1.3................................. 5 2. 9 2.1......................... 9 2.2........................

More information

t = h x z z = h z = t (x, z) (v x (x, z, t), v z (x, z, t)) ρ v x x + v z z = 0 (1) 2-2. (v x, v z ) φ(x, z, t) v x = φ x, v z

t = h x z z = h z = t (x, z) (v x (x, z, t), v z (x, z, t)) ρ v x x + v z z = 0 (1) 2-2. (v x, v z ) φ(x, z, t) v x = φ x, v z I 1 m 2 l k 2 x = 0 x 1 x 1 2 x 2 g x x 2 x 1 m k m 1-1. L x 1, x 2, ẋ 1, ẋ 2 ẋ 1 x = 0 1-2. 2 Q = x 1 + x 2 2 q = x 2 x 1 l L Q, q, Q, q M = 2m µ = m 2 1-3. Q q 1-4. 2 x 2 = h 1 x 1 t = 0 2 1 t x 1 (t)

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C602E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C602E646F63> スピントロニクスの基礎 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/077461 このサンプルページの内容は, 初版 1 刷発行時のものです. i 1 2 ii 3 5 4 AMR (anisotropic magnetoresistance effect) GMR (giant magnetoresistance

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

振動と波動

振動と波動 Report JS0.5 J Simplicity February 4, 2012 1 J Simplicity HOME http://www.jsimplicity.com/ Preface 2 Report 2 Contents I 5 1 6 1.1..................................... 6 1.2 1 1:................ 7 1.3

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

Note.tex 2008/09/19( )

Note.tex 2008/09/19( ) 1 20 9 19 2 1 5 1.1........................ 5 1.2............................. 8 2 9 2.1............................. 9 2.2.............................. 10 3 13 3.1.............................. 13 3.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

pdf

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

More information

基礎数学I

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

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

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

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

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

More information

01_.g.r..

01_.g.r.. I II III IV V VI VII VIII IX X XI I II III IV V I I I II II II I I YS-1 I YS-2 I YS-3 I YS-4 I YS-5 I YS-6 I YS-7 II II YS-1 II YS-2 II YS-3 II YS-4 II YS-5 II YS-6 II YS-7 III III YS-1 III YS-2

More information

数値計算:常微分方程式

数値計算:常微分方程式 ( ) 1 / 82 1 2 3 4 5 6 ( ) 2 / 82 ( ) 3 / 82 C θ l y m O x mg λ ( ) 4 / 82 θ t C J = ml 2 C mgl sin θ θ C J θ = mgl sin θ = θ ( ) 5 / 82 ω = θ J ω = mgl sin θ ω J = ml 2 θ = ω, ω = g l sin θ = θ ω ( )

More information

1 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω 1 ω α V T m T m 1 100Hz m 2 36km 500Hz. 36km 1

1 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω 1 ω α V T m T m 1 100Hz m 2 36km 500Hz. 36km 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω ω α 3 3 2 2V 3 33+.6T m T 5 34m Hz. 34 3.4m 2 36km 5Hz. 36km m 34 m 5 34 + m 5 33 5 =.66m 34m 34 x =.66 55Hz, 35 5 =.7 485.7Hz 2 V 5Hz.5V.5V V

More information

1 variation 1.1 imension unit L m M kg T s Q C QT 1 A = C s 1 MKSA F = ma N N = kg m s 1.1 J E = 1 mv W = F x J = kg m s 1 = N m 1.

1 variation 1.1 imension unit L m M kg T s Q C QT 1 A = C s 1 MKSA F = ma N N = kg m s 1.1 J E = 1 mv W = F x J = kg m s 1 = N m 1. 1.1 1. 1.3.1..3.4 3.1 3. 3.3 4.1 4. 4.3 5.1 5. 5.3 6.1 6. 6.3 7.1 7. 7.3 1 1 variation 1.1 imension unit L m M kg T s Q C QT 1 A = C s 1 MKSA F = ma N N = kg m s 1.1 J E = 1 mv W = F x J = kg m s 1 = N

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> 新 Excel コンピュータシミュレーション サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/084871 このサンプルページの内容は, 初版 1 刷発行当時のものです. Microsoft Excel Excel Visual Basic Visual Basic 2007 Excel Excel

More information

4................................. 4................................. 4 6................................. 6................................. 9.................................................... 3..3..........................

More information

u = u(t, x 1,..., x d ) : R R d C λ i = 1 := x 2 1 x 2 d d Euclid Laplace Schrödinger N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3

u = u(t, x 1,..., x d ) : R R d C λ i = 1 := x 2 1 x 2 d d Euclid Laplace Schrödinger N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3 2 2 1 5 5 Schrödinger i u t + u = λ u 2 u. u = u(t, x 1,..., x d ) : R R d C λ i = 1 := 2 + + 2 x 2 1 x 2 d d Euclid Laplace Schrödinger 3 1 1.1 N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3,... } Q

More information

<4D F736F F D B B BB2D834A836F815B82D082C88C602E646F63>

<4D F736F F D B B BB2D834A836F815B82D082C88C602E646F63> 信号処理の基礎 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/081051 このサンプルページの内容は, 初版 1 刷発行時のものです. i AI ii z / 2 3 4 5 6 7 7 z 8 8 iii 2013 3 iv 1 1 1.1... 1 1.2... 2 2 4 2.1...

More information

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

ルベーグ積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. ルベーグ積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/005431 このサンプルページの内容は, 初版 1 刷発行時のものです. Lebesgue 1 2 4 4 1 2 5 6 λ a

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

December 28, 2018

December 28, 2018 e-mail : kigami@i.kyoto-u.ac.jp December 28, 28 Contents 2............................. 3.2......................... 7.3..................... 9.4................ 4.5............. 2.6.... 22 2 36 2..........................

More information

??

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

More information

X G P G (X) G BG [X, BG] S 2 2 2 S 2 2 S 2 = { (x 1, x 2, x 3 ) R 3 x 2 1 + x 2 2 + x 2 3 = 1 } R 3 S 2 S 2 v x S 2 x x v(x) T x S 2 T x S 2 S 2 x T x S 2 = { ξ R 3 x ξ } R 3 T x S 2 S 2 x x T x S 2

More information

( ) ( )

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

More information

II 1 II 2012 II Gauss-Bonnet II

II 1 II 2012 II Gauss-Bonnet II II 1 II 212 II Gauss-Bonnet II 1 1 1.1......................................... 1 1.2............................................ 2 1.3.................................. 3 1.4.............................................

More information

DE-resume

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

More information

C : q i (t) C : q i (t) q i (t) q i(t) q i(t) q i (t)+δq i (t) (2) δq i (t) δq i (t) C, C δq i (t 0 )0, δq i (t 1 ) 0 (3) δs S[C ] S[C] t1 t 0 t1 t 0

C : q i (t) C : q i (t) q i (t) q i(t) q i(t) q i (t)+δq i (t) (2) δq i (t) δq i (t) C, C δq i (t 0 )0, δq i (t 1 ) 0 (3) δs S[C ] S[C] t1 t 0 t1 t 0 1 2003 4 24 ( ) 1 1.1 q i (i 1,,N) N [ ] t t 0 q i (t 0 )q 0 i t 1 q i (t 1 )q 1 i t 0 t t 1 t t 0 q 0 i t 1 q 1 i S[q(t)] t1 t 0 L(q(t), q(t),t)dt (1) S[q(t)] L(q(t), q(t),t) q 1.,q N q 1,, q N t C :

More information

c 2006 Yoneda norimasa All rights reserved

c 2006 Yoneda norimasa All rights reserved LEGO 011242 2006 2 Copyright c 2006 by Yoneda norimasa c 2006 Yoneda norimasa All rights reserved , LEGO,., Matlab Simlink,,., LEGO., 2, 2.,, LEGOMINDSTORMS, LOBOLAB.,,., Matlab,.,,,,,.,,,,, 2.,. i 1

More information

Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009631 このサンプルページの内容は, 初版 1 刷発行時のものです. Excel URL http://www.morikita.co.jp/books/mid/009631 i Microsoft Windows

More information

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

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

More information

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =,

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =, [ ] IC. r, θ r, θ π, y y = 3 3 = r cos θ r sin θ D D = {, y ; y }, y D r, θ ep y yddy D D 9 s96. d y dt + 3dy + y = cos t dt t = y = e π + e π +. t = π y =.9 s6.3 d y d + dy d + y = y =, dy d = 3 a, b

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

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

(u(x)v(x)) = u (x)v(x) + u(x)v (x) ( ) u(x) = u (x)v(x) u(x)v (x) v(x) v(x) 2 y = g(t), t = f(x) y = g(f(x)) dy dx dy dx = dy dt dt dx., y, f, g y = f (g(x))g (x). ( (f(g(x)). ). [ ] y = e ax+b (a, b )

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

量子力学 問題

量子力学 問題 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

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

More information

Untitled

Untitled http://www.ike-dyn.ritsumei.ac.jp/ hyoo/dynamics.html 1 (i) (ii) 2 (i) (ii) (*) [1] 2 1 3 1.1................................ 3 1.2..................................... 3 2 4 2.1....................................

More information

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

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

More information

214 March 31, 214, Rev.2.1 4........................ 4........................ 5............................. 7............................... 7 1 8 1.1............................... 8 1.2.......................

More information

z f(z) f(z) x, y, u, v, r, θ r > 0 z = x + iy, f = u + iv C γ D f(z) f(z) D f(z) f(z) z, Rm z, z 1.1 z = x + iy = re iθ = r (cos θ + i sin θ) z = x iy

z f(z) f(z) x, y, u, v, r, θ r > 0 z = x + iy, f = u + iv C γ D f(z) f(z) D f(z) f(z) z, Rm z, z 1.1 z = x + iy = re iθ = r (cos θ + i sin θ) z = x iy z fz fz x, y, u, v, r, θ r > z = x + iy, f = u + iv γ D fz fz D fz fz z, Rm z, z. z = x + iy = re iθ = r cos θ + i sin θ z = x iy = re iθ = r cos θ i sin θ x = z + z = Re z, y = z z = Im z i r = z = z

More information

2012 IA 8 I p.3, 2 p.19, 3 p.19, 4 p.22, 5 p.27, 6 p.27, 7 p

2012 IA 8 I p.3, 2 p.19, 3 p.19, 4 p.22, 5 p.27, 6 p.27, 7 p 2012 IA 8 I 1 10 10 29 1. [0, 1] n x = 1 (n = 1, 2, 3,...) 2 f(x) = n 0 [0, 1] 2. 1 x = 1 (n = 1, 2, 3,...) 2 f(x) = n 0 [0, 1] 1 0 f(x)dx 3. < b < c [, c] b [, c] 4. [, b] f(x) 1 f(x) 1 f(x) [, b] 5.

More information

213 March 25, 213, Rev.1.5 4........................ 4........................ 6 1 8 1.1............................... 8 1.2....................... 9 2 14 2.1..................... 14 2.2............................

More information

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

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

More information

B ver B

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

More information

II Karel Švadlenka * [1] 1.1* 5 23 m d2 x dt 2 = cdx kx + mg dt. c, g, k, m 1.2* u = au + bv v = cu + dv v u a, b, c, d R

II Karel Švadlenka * [1] 1.1* 5 23 m d2 x dt 2 = cdx kx + mg dt. c, g, k, m 1.2* u = au + bv v = cu + dv v u a, b, c, d R II Karel Švadlenka 2018 5 26 * [1] 1.1* 5 23 m d2 x dt 2 = cdx kx + mg dt. c, g, k, m 1.2* 5 23 1 u = au + bv v = cu + dv v u a, b, c, d R 1.3 14 14 60% 1.4 5 23 a, b R a 2 4b < 0 λ 2 + aλ + b = 0 λ =

More information

Z: Q: R: C: sin 6 5 ζ a, b

Z: Q: R: C: sin 6 5 ζ a, b Z: Q: R: C: 3 3 7 4 sin 6 5 ζ 9 6 6............................... 6............................... 6.3......................... 4 7 6 8 8 9 3 33 a, b a bc c b a a b 5 3 5 3 5 5 3 a a a a p > p p p, 3,

More information

e a b a b b a a a 1 a a 1 = a 1 a = e G G G : x ( x =, 8, 1 ) x 1,, 60 θ, ϕ ψ θ G G H H G x. n n 1 n 1 n σ = (σ 1, σ,..., σ N ) i σ i i n S n n = 1,,

e a b a b b a a a 1 a a 1 = a 1 a = e G G G : x ( x =, 8, 1 ) x 1,, 60 θ, ϕ ψ θ G G H H G x. n n 1 n 1 n σ = (σ 1, σ,..., σ N ) i σ i i n S n n = 1,, 01 10 18 ( ) 1 6 6 1 8 8 1 6 1 0 0 0 0 1 Table 1: 10 0 8 180 1 1 1. ( : 60 60 ) : 1. 1 e a b a b b a a a 1 a a 1 = a 1 a = e G G G : x ( x =, 8, 1 ) x 1,, 60 θ, ϕ ψ θ G G H H G x. n n 1 n 1 n σ = (σ 1,

More information

sec13.dvi

sec13.dvi 13 13.1 O r F R = m d 2 r dt 2 m r m = F = m r M M d2 R dt 2 = m d 2 r dt 2 = F = F (13.1) F O L = r p = m r ṙ dl dt = m ṙ ṙ + m r r = r (m r ) = r F N. (13.2) N N = R F 13.2 O ˆn ω L O r u u = ω r 1 1:

More information

n=1 1 n 2 = π = π f(z) f(z) 2 f(z) = u(z) + iv(z) *1 f (z) u(x, y), v(x, y) f(z) f (z) = f/ x u x = v y, u y = v x

n=1 1 n 2 = π = π f(z) f(z) 2 f(z) = u(z) + iv(z) *1 f (z) u(x, y), v(x, y) f(z) f (z) = f/ x u x = v y, u y = v x n= n 2 = π2 6 3 2 28 + 4 + 9 + = π2 6 2 f(z) f(z) 2 f(z) = u(z) + iv(z) * f (z) u(x, y), v(x, y) f(z) f (z) = f/ x u x = v y, u y = v x f x = i f y * u, v 3 3. 3 f(t) = u(t) + v(t) [, b] f(t)dt = u(t)dt

More information

( ) ( 40 )+( 60 ) Schrödinger 3. (a) (b) (c) yoshioka/education-09.html pdf 1

( ) ( 40 )+( 60 ) Schrödinger 3. (a) (b) (c)   yoshioka/education-09.html pdf 1 2009 1 ( ) ( 40 )+( 60 ) 1 1. 2. Schrödinger 3. (a) (b) (c) http://goofy.phys.nara-wu.ac.jp/ yoshioka/education-09.html pdf 1 1. ( photon) ν λ = c ν (c = 3.0 108 /m : ) ɛ = hν (1) p = hν/c = h/λ (2) h

More information

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

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

More information

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

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> 高速流体力学 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/067361 このサンプルページの内容は, 第 1 版発行時のものです. i 20 1999 3 2 2010 5 ii 1 1 1.1 1 1.2 4 9 2 10 2.1 10 2.2 12 2.3 13 2.4 13 2.5

More information

I, II 1, 2 ɛ-δ 100 A = A 4 : 6 = max{ A, } A A 10

I, II 1, 2 ɛ-δ 100 A = A 4 : 6 = max{ A, } A A 10 1 2007.4.13. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 0. 1. 1. 2. 3. 2. ɛ-δ 1. ɛ-n

More information

4.2.................... 20 4.3.................. 21 4.4 ( )............... 22 4.5 ( )...... 24 4.6 ( )........ 25 4.7 ( )..... 26 5 28 5.1 PID........

4.2.................... 20 4.3.................. 21 4.4 ( )............... 22 4.5 ( )...... 24 4.6 ( )........ 25 4.7 ( )..... 26 5 28 5.1 PID........ version 0.01 : 2004/04/16 1 2 1.1................. 2 1.2.......................... 3 1.3................. 5 1.4............... 6 1.5.............. 7 2 9 2.1........................ 9 2.2......................

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

<4D F736F F D B B BB2D834A836F815B82D082C88C60202D B2E646F63>

<4D F736F F D B B BB2D834A836F815B82D082C88C60202D B2E646F63> 例題で学ぶはじめての塑性力学 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/066721 このサンプルページの内容は, 初版 1 刷発行当時のものです. http://www.morikita.co.jp/support/ 03 3817 5670 FAX 03 3815 8199 i 1

More information

untitled

untitled i ii iii iv v 43 43 vi 43 vii T+1 T+2 1 viii 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 a) ( ) b) ( ) 51

More information

2

2 1 2 3 4 5 6 7 8 9 10 I II III 11 IV 12 V 13 VI VII 14 VIII. 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 _ 33 _ 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 VII 51 52 53 54 55 56 57 58 59

More information

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2 II No.1 [n/] [1]H n x) H n x) = 1) r n! r!n r)! x)n r r= []H n x) n,, H n x) = 1) n H n x) [3] H n x) = 1) n dn x e dx n e x [4] H n+1 x) = xh n x) nh n 1 x) ) d dx x H n x) = H n+1 x) d dx H nx) = nh

More information

( ; ) C. H. Scholz, The Mechanics of Earthquakes and Faulting : - ( ) σ = σ t sin 2π(r a) λ dσ d(r a) =

( ; ) C. H. Scholz, The Mechanics of Earthquakes and Faulting : - ( ) σ = σ t sin 2π(r a) λ dσ d(r a) = 1 9 8 1 1 1 ; 1 11 16 C. H. Scholz, The Mechanics of Earthquakes and Faulting 1. 1.1 1.1.1 : - σ = σ t sin πr a λ dσ dr a = E a = π λ σ πr a t cos λ 1 r a/λ 1 cos 1 E: σ t = Eλ πa a λ E/π γ : λ/ 3 γ =

More information

5 1.2, 2, d a V a = M (1.2.1), M, a,,,,, Ω, V a V, V a = V + Ω r. (1.2.2), r i 1, i 2, i 3, i 1, i 2, i 3, A 2, A = 3 A n i n = n=1 da = 3 = n=1 3 n=1

5 1.2, 2, d a V a = M (1.2.1), M, a,,,,, Ω, V a V, V a = V + Ω r. (1.2.2), r i 1, i 2, i 3, i 1, i 2, i 3, A 2, A = 3 A n i n = n=1 da = 3 = n=1 3 n=1 4 1 1.1 ( ) 5 1.2, 2, d a V a = M (1.2.1), M, a,,,,, Ω, V a V, V a = V + Ω r. (1.2.2), r i 1, i 2, i 3, i 1, i 2, i 3, A 2, A = 3 A n i n = n=1 da = 3 = n=1 3 n=1 da n i n da n i n + 3 A ni n n=1 3 n=1

More information

2

2 III 22 7 4 3....................................... 3.2 Kepler ( ).......................... 2 2 4 2.................................. 4 2.2......................................... 8 3 20 3..........................................

More information

Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x

Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x University of Hyogo 8 8 1 d x(t) =f(t, x(t)), dt (1) x(t 0 ) =x 0 () t n = t 0 + n t x x n n x n x 0 x i i = 0,..., n 1 x n x(t) 1 1.1 1 1 1 0 θ 1 θ x n x n 1 t = θf(t n 1, x n 1 ) + (1 θ)f(t n, x n )

More information

genron-3

genron-3 " ( K p( pasals! ( kg / m 3 " ( K! v M V! M / V v V / M! 3 ( kg / m v ( v "! v p v # v v pd v ( J / kg p ( $ 3! % S $ ( pv" 3 ( ( 5 pv" pv R" p R!" R " ( K ( 6 ( 7 " pv pv % p % w ' p% S & $ p% v ( J /

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

() (, y) E(, y) () E(, y) (3) q ( ) () E(, y) = k q q (, y) () E(, y) = k r r (3).3 [.7 ] f y = f y () f(, y) = y () f(, y) = tan y y ( ) () f y = f y

() (, y) E(, y) () E(, y) (3) q ( ) () E(, y) = k q q (, y) () E(, y) = k r r (3).3 [.7 ] f y = f y () f(, y) = y () f(, y) = tan y y ( ) () f y = f y 5. [. ] z = f(, y) () z = 3 4 y + y + 3y () z = y (3) z = sin( y) (4) z = cos y (5) z = 4y (6) z = tan y (7) z = log( + y ) (8) z = tan y + + y ( ) () z = 3 8y + y z y = 4 + + 6y () z = y z y = (3) z =

More information

215 11 13 1 2 1.1....................... 2 1.2.................... 2 1.3..................... 2 1.4...................... 3 1.5............... 3 1.6........................... 4 1.7.................. 4

More information

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

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

More information

.1 z = e x +xy y z y 1 1 x 0 1 z x y α β γ z = αx + βy + γ (.1) ax + by + cz = d (.1') a, b, c, d x-y-z (a, b, c). x-y-z 3 (0,

.1 z = e x +xy y z y 1 1 x 0 1 z x y α β γ z = αx + βy + γ (.1) ax + by + cz = d (.1') a, b, c, d x-y-z (a, b, c). x-y-z 3 (0, .1.1 Y K L Y = K 1 3 L 3 L K K (K + ) 1 1 3 L 3 K 3 L 3 K 0 (K + K) 1 3 L 3 K 1 3 L 3 lim K 0 K = L (K + K) 1 3 K 1 3 3 lim K 0 K = 1 3 K 3 L 3 z = f(x, y) x y z x-y-z.1 z = e x +xy y 3 x-y ( ) z 0 f(x,

More information

2 1 1 (1) 1 (2) (3) Lax : (4) Bäcklund : (5) (6) 1.1 d 2 q n dt 2 = e q n 1 q n e q n q n+1 (1.1) 1 m q n n ( ) r n = q n q n 1 r ϕ(r) ϕ (r)

2 1 1 (1) 1 (2) (3) Lax : (4) Bäcklund : (5) (6) 1.1 d 2 q n dt 2 = e q n 1 q n e q n q n+1 (1.1) 1 m q n n ( ) r n = q n q n 1 r ϕ(r) ϕ (r) ( ( (3 Lax : (4 Bäcklud : (5 (6 d q = e q q e q q + ( m q ( r = q q r ϕ(r ϕ (r 0 5 0 q q q + 5 3 4 5 m d q = ϕ (r + ϕ (r + ( Hooke ϕ(r = κr (κ > 0 ( d q = κ(q q + κ(q + q = κ(q + + q q (3 ϕ(r = a b e br

More information

7-12.dvi

7-12.dvi 26 12 1 23. xyz ϕ f(x, y, z) Φ F (x, y, z) = F (x, y, z) G(x, y, z) rot(grad ϕ) rot(grad f) H(x, y, z) div(rot Φ) div(rot F ) (x, y, z) rot(grad f) = rot f x f y f z = (f z ) y (f y ) z (f x ) z (f z )

More information

r d 2r d l d (a) (b) (c) 1: I(x,t) I(x+ x,t) I(0,t) I(l,t) V in V(x,t) V(x+ x,t) V(0,t) l V(l,t) 2: 0 x x+ x 3: V in 3 V in x V (x, t) I(x, t

r d 2r d l d (a) (b) (c) 1: I(x,t) I(x+ x,t) I(0,t) I(l,t) V in V(x,t) V(x+ x,t) V(0,t) l V(l,t) 2: 0 x x+ x 3: V in 3 V in x V (x, t) I(x, t 1 1 2 2 2r d 2r d l d (a) (b) (c) 1: I(x,t) I(x+ x,t) I(0,t) I(l,t) V in V(x,t) V(x+ x,t) V(0,t) l V(l,t) 2: 0 x x+ x 3: V in 3 V in x V (x, t) I(x, t) V (x, t) I(x, t) V in x t 3 4 1 L R 2 C G L 0 R 0

More information

2 Chapter 4 (f4a). 2. (f4cone) ( θ) () g M. 2. (f4b) T M L P a θ (f4eki) ρ H A a g. v ( ) 2. H(t) ( )

2 Chapter 4 (f4a). 2. (f4cone) ( θ) () g M. 2. (f4b) T M L P a θ (f4eki) ρ H A a g. v ( ) 2. H(t) ( ) http://astr-www.kj.yamagata-u.ac.jp/~shibata f4a f4b 2 f4cone f4eki f4end 4 f5meanfp f6coin () f6a f7a f7b f7d f8a f8b f9a f9b f9c f9kep f0a f0bt version feqmo fvec4 fvec fvec6 fvec2 fvec3 f3a (-D) f3b

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> 通信方式第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/072662 このサンプルページの内容は, 第 2 版発行当時のものです. i 2 2 2 2012 5 ii,.,,,,,,.,.,,,,,.,,.,,..,,,,.,,.,.,,.,,.. 1990 5 iii 1 1

More information

(yx4) 1887-1945 741936 50 1995 1 31 http://kenboushoten.web.fc.com/ OCR TeX 50 yx4 e-mail: yx4.aydx5@gmail.com i Jacobi 1751 1 3 Euler Fagnano 187 9 0 Abel iii 1 1...................................

More information