C 2. /

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

sakigake1.dvi

ẋ = ax + y ẏ = x x by 2 Griffith a b Saddle Node Saddle-Node (phase plane) Griffith mrna(y) Protein(x) (nullcline) 0 (nullcline) (


Part () () Γ Part ,

phs.dvi


Shunsuke Kobayashi 1 [6] [11] [7] u t = D 2 u 1 x 2 + f(u, v) + s L u(t, x)dx, L x (0.L), t > 0, Neumann 0 v t = D 2 v 2 + g(u, v), x (0, L), t > 0. x

r 1 m A r/m i) t ii) m i) t B(t; m) ( B(t; m) = A 1 + r ) mt m ii) B(t; m) ( B(t; m) = A 1 + r ) mt m { ( = A 1 + r ) m } rt r m n = m r m n B

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

1 a b = max{a, b}, a b = mi{a, b} a 1 a 2 a a 1 a = max{a 1,... a }, a 1 a = mi{a 1,... a }. A sup A, if A A A A A sup A sup A = + A if A = ± y = arct

(Basic of Proability Theory). (Probability Spacees ad Radom Variables , (Expectatios, Meas) (Weak Law

A11 (1993,1994) 29 A12 (1994) 29 A13 Trefethen and Bau Numerical Linear Algebra (1997) 29 A14 (1999) 30 A15 (2003) 30 A16 (2004) 30 A17 (2007) 30 A18

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


x (x, ) x y (, y) iy x y z = x + iy (x, y) (r, θ) r = x + y, θ = tan ( y ), π < θ π x r = z, θ = arg z z = x + iy = r cos θ + ir sin θ = r(cos θ + i s

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

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


DVIOUT

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

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

i 18 2H 2 + O 2 2H 2 + ( ) 3K

86 6 r (6) y y d y = y 3 (64) y r y r y r ϕ(x, y, y,, y r ) n dy = f(x, y) (6) 6 Lipschitz 6 dy = y x c R y(x) y(x) = c exp(x) x x = x y(x ) = y (init

() n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (5) (6 ) n C + nc + 3 nc n nc n (7 ) n C + nc + 3 nc n nc n (


1 Abstract 2 3 n a ax 2 + bx + c = 0 (a 0) (1) ( x + b ) 2 = b2 4ac 2a 4a 2 D = b 2 4ac > 0 (1) 2 D = 0 D < 0 x + b 2a = ± b2 4ac 2a b ± b 2

D 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

tokei01.dvi

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

(2 X Poisso P (λ ϕ X (t = E[e itx ] = k= itk λk e k! e λ = (e it λ k e λ = e eitλ e λ = e λ(eit 1. k! k= 6.7 X N(, 1 ϕ X (t = e 1 2 t2 : Cauchy ϕ X (t

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

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

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

6. Euler x

II (Percolation) ( 3-4 ) 1. [ ],,,,,,,. 2. [ ],.. 3. [ ],. 4. [ ] [ ] G. Grimmett Percolation Springer-Verlag New-York [ ] 3

(Basics of Proability Theory). (Probability Spacees ad Radom Variables,, (Ω, F, P ),, X,. (Ω, F, P ) (probability space) Ω ( ω Ω ) F ( 2 Ω ) Ω σ (σ-fi

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

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

TOP URL 1

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

B [ 0.1 ] x > 0 x 6= 1 f(x) µ 1 1 xn 1 + sin sin x 1 x 1 f(x) := lim. n x n (1) lim inf f(x) (2) lim sup f(x) x 1 0 x 1 0 (


i

Gmech08.dvi

1

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

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

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


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

1 Introduction 1 (1) (2) (3) () {f n (x)} n=1 [a, b] K > 0 n, x f n (x) K < ( ) x [a, b] lim f n (x) f(x) (1) f(x)? (2) () f(x)? b lim a f n (x)dx = b

IA 2013 : :10722 : 2 : :2 :761 :1 (23-27) : : ( / ) (1 /, ) / e.g. (Taylar ) e x = 1 + x + x xn n! +... sin x = x x3 6 + x5 x2n+1 + (

数学Ⅱ演習(足助・09夏)

1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2

数学の基礎訓練I

B 38 1 (x, y), (x, y, z) (x 1, x 2 ) (x 1, x 2, x 3 ) 2 : x 2 + y 2 = 1. (parameter) x = cos t, y = sin t. y = f(x) r(t) = (x(t), y(t), z(t)), a t b.

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

入試の軌跡

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

LINEAR ALGEBRA I Hiroshi SUZUKI Department of Mathematics International Christian University

1 1 n 0, 1, 2,, n n 2 a, b a n b n a, b n a b (mod n) 1 1. n = (mod 10) 2. n = (mod 9) n II Z n := {0, 1, 2,, n 1} 1.

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

III III 2010 PART I 1 Definition 1.1 (, σ-),,,, Borel( ),, (σ-) (M, F, µ), (R, B(R)), (C, B(C)) Borel Definition 1.2 (µ-a.e.), (in µ), (in L 1 (µ)). T

03.Œk’ì

Transcription:

C 2. / 24 8 5-6 9-2

i.................................................................................2.................................. 3..3................................. 5..4..................................... 6..5........................................ 7..6.................................. 7.2......................................... 9.2.................................... 9.2.2 DYNASCAL....................................... 2 2 22 2............................................ 22 2......................................... 22 2..2.................................. 23 2..3............................. 24 2..4................................. 26 2..5..................... 26 2..6....................................... 27 2.2............................................ 29 2.2.................. 29 2.2.2................................... 29 2.2.3......................... 29 A 3

... (differential equation (solution curve dx = ax, (. dt (a solution k x = f(t =e ta k, (.2 (. k (. (an initial value problem x = ax, x( = k, (.3 (. (a unique solution. a =. (.2 k 2,,, 2 Maple 9 (uniqueness of solutions t (local existence t (existence in the large Hartman (982 Hartman [ ] y = y 2, y( =. (.4 c t c = t c =(t c 2 /4 [.] (. (.2

2 4 2 - -5 t 5-2 -4.: The simplest differential equations with initial values (.3 (.2. a> k > x k < x 2. a = x = k(constant 3. a< x Fig. a =. k =2,, 2 (. (linear homogeneous first-order differential equations dx/dt = A(t x (homogeneous equation (. x (nonhomogeneous equation (. x B(t dx/dt = A(t x + B(t (higher-order differential equations d 2 x dt 2 + a dx + b(t = dt

3..2 (state space (a system of linear first-order differential equations dx dt = Ax,,x( = k Rn. (.5 (.5 x = f(t =e ta k = e ta x. (.6 e ta = k= k! Ak (Hirsch & Smale, 974, p.86. dx dt = g(x, x R n g(x = x = x (an equilibrium point x = φ t (x t φ t (x =x. (.7 x (a stationary point (a fixed point (a zero point (a singular point, or singularity 2 - -.5.5 - -2.2: The vector field and some solution curves near a saddle with λ < <λ 2 (.5 A dx/dt = Ax x = Py P (.5 dy dt = By, B = P AP, y( = P x(. (.8

4 P x( = k [.2] dx/dt = Ax (.8 P A A A (real canonical form.3.2. -2 - -. 2 -.2 -.3.3: The vector field and some solution curves near a saddle with λ 2 < <λ.5 - -.5.5 -.5 -.4: The vector field and some solution curves near an inward node with λ <λ 2 <

5 n A J O J 2... P AP = F F 2. O.., (.9 J J 2 (elementary Jordan matrices λ (elementary λ blocks λ O... J l =, l =, 2,...,p. (....... O λ F F 2 D O. I.. 2 F m =...... O I 2 D, D, m =, 2,...,q. (. ( D = a b b a (, I 2 = (.2 D λ = a ± bi.5 - -.5.5 -.5 -.5: The vector field and some solution curves near an inward node with λ 2 <λ < A A

6 (a ( (b ( ( λ λ 2 λ λ λ λ,,, (c ( a b b a, b > ( a b b a, b > 2 -.4 -.2.2.4 - -2.6: The vector field and some solution curves near an outward node with <λ <λ 2 (.8 =e tb k, (.3 ( λ Case (a B =, λ 2 (.8 ( λ y = λ 2 y 2 ( dy dt dy 2 dt, (.4

7 (.2 ( ( y (t e λt k =. (.5 y 2 (t e λ2t k 2 (.5 y 2 (t y (t = e(λ2 λ t k 2. (.6 k (.5 (.6 λ λ 2 λ λ 2 (saddle (.2.3 (inward node.4.5 (outward node.6.7 λ = λ 2 (.6 y 2 (t y (t = k 2 k, (.7 t y (t y 2 (t k k 2.8.9.4.2-2 - 2 -.2 -.4.7: The vector field and some solution curves near an outward node with <λ 2 <λ ( ( λ λ Case (b B =, B =, λ λ ( λ B =, =e tb k λ ( ( y (t = y 2 (t e λt k e λt tk + e λt k 2. (.8 (.8

8.5 - -.5.5 -.5 -.8: The vector field and some solution curves near an inward focus with λ = λ 2 <.4.2 -.4 -.2.2.4 -.2 -.4.9: The vector field and some solution curves near an outward focus with λ = λ 2 > ( tλ Case (b B tb = t tλ ( tb =(tλi + N, N =. t. e tb = e tλi e N. e X = r= X r r!,

9 e tλi = e tλ I N 2 = N 3 =...= O e N = I + N ( e tb = e tλ I (I + N =e tλ (I + N =e tλ. t ( e tb = e tλ te tλ e tλ..5 - -.5.5 -.5 -.: The vector field and some solution curves near an inward improper node with eqs.(.2,λ <.5 - -.5.5 -.5 -.: The vector field and some solution curves near an inward improper node with eqs.(.8,λ < (.8 y 2 (t y (t = t + k 2 k. (.9

( λ B = λ ( λ λ ( e tb = e tλ te tλ e tλ. (.2 ( y (t y 2 (t ( e λt k + e λt tk 2 =. (.2 e λt k 2 y (t y 2 (t = t + k k 2. (.22 (.9 (.22 t y (t y 2 (t, y 2 (t. (.23 y (t.. λ >.2.3..5 -. -.5.5. -.5 -..2: The vector field and some solution curves near an outward improper node with eqs.(.2,λ > ( ( a b a b Case (c B =, B =, b > b a b a ( a b B =, b a ( ta tb tb = = tai + tbl. (.24 tb ta L = (. (.25

(.24 e tb = e tai e tbl = e ta e tbl. (.26 e tbl L ( e tbl cos tb sin tb =. (.27 sin tb cos tb e tbl = I + tb! L + (tb2 L 2 +, 2! L 2 = I, L 3 = L, L 4 = I, L 5 = L, (.28 L cos(tb sin(tb.5 - -.5.5 -.5 -.3: The vector field and some solution curves near an outward improper node with eqs.(.8,λ > (.26 (.27 e tb = e ta ( cos tb sin tb sin tb cos tb. (.29 ( y (t y 2 (t ( e ta (k cos tb k 2 sin tb = e ta (k sin tb + k 2 cos tb. (.3 y 2 (t y (t = k sin tb + k 2 cos tb k cos tb k 2 sin tb. (.3 (.3 cos tb sin tb a < lim t =

2 y ( a b = b a ( dy dt dy 2 dt ( y ( ay = by (.32 y 2 dy 2 dt = by. (.33 b > y > dy 2 /dt > y < dy 2 /dt < spiral.4 b < spiral.5 y 2 = dy 2 /dt = by a > lim t =.6.7.8.4 -.8 -.6 -.4 -.2.2.4 -.4.4: The vector field and some solution curves near a (counterclockwise spiral sink.8.4 -.4 -.2.2.4.6.8 -.4.5: The vector field and some solution curves near a (clockwise spiral sink

3 ( ( b b Case (c B =, B =, b > b b Case (c a = (.3.8.4 -.4 -.2.2.4.6.8 -.4.6: The vector field and some solution curves near a (counterclockwise spiral source.8.4 -.8 -.6 -.4 -.2.2.4 -.4.7: The vector field and some solution curves near a (clockwise spiral source ( y (t y 2 (t ( = (k cos tb k 2 sin tb (k sin tb + k 2 cos tb. (.34 (.34 y 2(t+y2 2 (t y 2 (t+y2 2 (t =k2 + k2 2 = r2 (k,k 2. (.35 (k,k 2 Case (b.8.9

4 A (negative real parts sink inward node inward improper node inward focus spiral sink sink A (positive real parts source outward node outward improper node outward focus spiral source source.5.5 -.5 - -.5.5.5 -.5 - -.5.8: The vector field and some solution curves near a (counterclockwise center.5.5 -.5 - -.5.5 -.5.5 - -.5.9: The vector field and some solution curves near a (clockwise center 99, Chino, 987, 99.2.3 saddle

5.4.5 spiral sink.8 focus.6.7.9 focus focus.. improper node.2.3 improper node.4.5 spiral sinks.6.7 spiral source.8.9 centers center..3..2 (.5 A (a system of nonlinear first-order differential equations (linearize Guckenheimer & Holmes, 983, pp.2-3 dx dt = f(x, x Rn, x( = k. (.36 (.36 x f(x = (.36 f(x x (the Jacobian matrix Jx dx dt = J x x, (.37 Jx { } = fi = x j x=x f f x f 2 f 2 x x 2 f x n f 2 x n x 2... f n x f n x 2 f n x n. (.38 x=x (non-zero real part

6 center.8.9..4 (limit cycles spiral sink spiral source (closed orbit (α limit cycle (ω limit cycle spiral sink spiral source.2.5.5 -.5.5 -.5 -.2: An ω-limit cycle Hirsch & Smale (974, pp.227-228 µ =.5 Maple ( ( µx x 3 + y =. (.39 x source µ

7 8 Heath, 2 (99 Chino (987, 99.2 spiral sink spiral source..5 (Hartman & Grobman s theorem (non-zero real part (structurally stable (structurally unstable.8.9 center..6 (bifurcation (bifurcation parameters c dx dt = f(x, c. (.4 (.39 c = µ (local bifurcation Guckenheimer & Holmes, 983 (saddle-node bifurcation (pitchfork bifurcation

8.4.2 -.2 -...2.3.4.5 -.2.2: A spiral sink before an ω-limit cycle appears (transcritical bifurcation (Hopf bifurcation saddle node spiral sink spiral source center center spiral sink spiral source spiral source spiral sink.2.2 spiral sink Maple (.39 µ =.5.2 µ =.5 µ =. (saddle loop (saddle connection spiral sink saddle spiral sink center center source spiral sink saddle saddles saddle saddle saddles (99 Chino (987, 99

9.2.2 spiral sink sink source spiral source spiral sink.2 PsycINFO dynamical system Tobler (976-977 (the wind model (flow (curl free part (divergence free part (gradient vector field (scalar potential 8 (bifurcation theory DYNASCAL (Dynamical System Scaling Chino & Nakagawa, 983, 99 DYNASCAL Yadohisa and Niki (999 DYNASCAL.2. (997 t jk = d jk r + c jk. (.4 t jk j k d jk r c jk

2 (.4 c jk = c kj t kj = d kj /(r + c kj c jk = r t kj t jk t jk + t kj. (.42 t jk (j,k (997.2.2 DYNASCAL DYNASCAL Guttman (968 SSA-II (nonautonomous system (autonomous system 4 2-8 -4 4 8-2.22: Illustration of the estimated vector field at a time by DYNASCAL.22 spiral sinks saddles DYNASCAL DYNASCAL (.4 x t dx dt = f(x,t, (.43

2 dx /dt f (x,x 2,x 3 dx 2 /dt = f 2 (x,x 2,x 3, (.44 dx 3 /dt DYNASCAL (.4 c DYNASCAL 23. 2. DYNASCAL 3. (limit cycles DYNASCAL Newcomb (96 Chino and Nakagawa (99 DYNASCAL...3 DYNASCAL. saddle.2 saddle 4 spiral source.3

22 2 2. 2.. (difference equation x(n +=a(n x(n, x(n =x, n n, (2. (a linear homogeneous first-order difference equation ( Elaydi (999 ( n x(n = a(i x, (2.2 i=n x(n +=a(n x(n+g(n, x(n =x, n n, (2.3 (a linear nonhomogeneous first-order difference equation ( n ( n n x(n = a(i x + a(i g(r, (2.4 i=n i=r+ k- Elaydi, 999, p.54 r=n y(n + k+p (n y(n + k + + p k (n y(n =g(n. (2.3 a(n =a g(n =b x(n +=ax(n+b, x(n =x, n n, (2.5 { ( a n a x + b n a, a x(n = x + bn, a = (2.6 (2.6 n. a> = x.

2 23 2. a = = b> x, b = x = x, b< x. 3. <a< = x b a. 4. a = = x x b x 5. a< = x ±. 2..2 k (a system of linear first-order difference equations x(n +=Ax(n, x(n =x. (2.7 A = a a 2 a k a 2 a 22 a 2k.... (2.8 a k a k2 a kk (real nonsingular matrix (2.7 x(n, n, x =A n n x. (2.9 y(n x =x(n y(n +=Ay(n y( = x(n (2.9 y(n =A n y(, (2. (2.7 A A n (the Putzer algorithm Elaydi, 999, p.6 (2.7 x(n =(x (n,x 2 (n,...,x k (n t k (states s,s 2,...,s k x(n =p(n =(p (n,p 2 (n,...,p k (n t A = P = {p ij = p(s i s j } p ij s j s i (the Markov chain p(n +=Sp(n, (2. S = {p ij }, (2.2 (the transition matrix A = {a ij } A a ij k k A j =, 2,...,k k i= a ij = (Markov (stochastic λ λ Elaydi, 999, p.42 λ = λ ρ(a ρ(a k k (the spectral radius ρ(a =max i k λ i (the Perron s theorem (Elaydi, 999, p.42

2 24 2..3 x(n +=f(x(n, x( = x, (2.3 x f (the fixed point f(x =x, (2.4 (2.3 (the equilibrium point (2.4 (.7 x(n x(n + x(n (2.4 x(n + x(n = (a x ɛ > x x <δ(δ > n> f n (x x <ɛ δ (stable (b x x x <η(η > lim n x n = x η (attracting x η = (a global attractor (globally attracting (c x (asymptotically stable equilibrium η = (globally asymptotically stable x x x δ> n> x n x ɛ> x n x x x x η> x n n x x x x n n x x x x η> n> x n x ɛ> n x Elaydi, 999 x x(n +=f(x(n, (2.5 f x (i f (x < x

2 25 (ii f (x > x x f (x (hyperbolic (2.5 x f (x = (i f (x x (ii f (x = f (x > x (iii f (x = f (x < x (2.5 x f (x = (i Sf(x < x (ii Sf(x > x Sf(x (the Schwarzian derivative of a function f Sf(x = f (x f (x 3 ( f 2 (x 2 f. (2.6 (x k k (a periodic point Elaydi, 999 b f (i b k f k (b =b f (2.5 x(n +=g(x(n, g = f k, (2.7 f k k- (k-periodic b O(b = { b, f(b,...,f k (b } k- (k-cycle (ii b m f m (b k- k- (eventually k-periodic (995 k- b (988, p.54 k- b f k- b (i f k (ii f k (iii f k

2 26 b k- { x( = b, x( = f(b,...,x(k = f k (b } k- O(b ={b = x(,x(,...,x(k } f k- k- (i f (x(f (x( f (x(k < (ii f (x(f (x( f (x(k > 2..4 dx/dt = f(x,t,b dx/dt = f(x, b b x(n +=f(x(n,n,b b Peitgen and Richter (986 (Verhulst dynamics x( n x(n R R = r( x(n x(n + = ( + Rx(n x(n + = ( + rx(n rx 2 (n. r r r r n r (chaos 2..5 x(n +=f(n, x(n, x(n =x, (2.8 x(n R k f : z + R k R k f(n, x(n x (2.8 n (nonaoutonomous (2.8 n x(n +=f(x(n (autonomous (2.8 n Z N f(n + N,x =f(n, x, (2.9

2 27 (periodic R k x f(n, x =x n n (2.8 (an equilibrium point x y(n =x(n x (2.8 y(n +=f(n, y(n+x x = g(n, y(n, y = x = x (2.8 Elaydi (999, p.57 (2.8 x (i ɛ > n n n x x < δ x(n, n, x x <ɛ δ = δ(ɛ, n (stable δ n (uniformly stable (unstable (ii x x <µ lim n x(n, n, x =x µ = µ(n (attractive µ n (uniformly attractive ɛ n x x <µ n n + N x(n, n, x x <ɛ n N = N(ɛ µ> (iii x (asymptotically stable ( (uniformly asymptotically stable (iv x x <δ x(n, n, x x M x x η n n δ> M> η (, (exponentially stable (v M n n x(n, n, x M x(n, n, x (bounded (ii (iii µ = (iv δ = (global 2..6 2..3 k- 2 x x(n +=Ax(n, x(n =x. (2.2 A

2 28 Ax = x (A Ix = (2.2 A I x = A I = c x = c y = x(n x (2.2 y(n +=Ay(n x x = (2.2 A x(n =Py(n P (2.2 y(n +=By(n, B = P AP, y( = P x(. (2.2 y( = k =(k,k 2 t (..2 A B (..2 (a (b (c ( λ Case (a B =, λ 2 (2.2 ( ( y (n + λ y (n =, (2.22 y 2 (n + λ 2 y 2 (n (2.6 b = ( ( y (n λ n = k y 2 (n λ n 2 k. (2.23 2 (2.23 ( n y 2 (n y (n = λ2 k 2. (2.24 λ k (2.23 (2.24 λ λ 2 λ λ 2 (saddle (asymptotically stable node (unstable node

2 29 2.2 (the Julia set (the filled Julia set (the Fatou set 2.2. f(z p 2 f(z =a n z p + a p z p + + a z + a, a n, (2.25. K f z z C 2. J K (the boundary K 3. F f n : S S n f n (the domain of normality Maple 2.2.2 (2.25 p =2 z(n +=z 2 (n+c, (2.26 c (the Mandelbrot set 2.2.3 Chino (22, 23 MDS HFM

2 3 References Chino, N. (987. A bifurcation model of changes in interdependence structure among objects. Bulletin of The Faculty of Humanities of Aichigakuin University, 7, 85-9. (99. II (pp.385-43 Chino, N. (22. Complex space models for the analysis of asymmetry. In S. Nishisato, Y. Baba, H. Bozdogan, & K. Kanefuji (Eds. Measurement and Multivariate Analysis, Tokyo: Springer. pp.7-4. Chino, N. (23. Complex difference system models for the analysis of asymetry. In H. Yanai, A. Okada, K. Shigemasu, Y. Kano, & J. J. Meulman (Eds. New Developments in Psychometrics, Tokyo: Springer. pp.479-486. (23. Chino, N., & Nakagawa, M. (983. A vector field model for sociometric data. Proceedings of the th annual meeting of the Behaviormetric Society of Japan (pp.9-, Kyoto, September. Chino, N., & Nakagawa, M. (99. A bifurcation model of change in group structure. The Japanese Journal of Experimental Social Psychology, 29, No.3, 25-38. Chino, N., & Shiraiwa, K. (993. Geometrical structures of some non-distance models for asymmetric MDS. Behaviormetrika, 2, 35-47. Elaydi, S. N. (999. An introduction to differential equations. 2nd Ed., New York: Springer-Verlag. Guckenheimer, J. & Holmes, P. (983. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Berlin: Springer-Verlag. Guttman, L. (968. A general nonmetric technique for finding the smallest coordinate space for a configuration of points. Psychometrika, 33, 469-56. Hartman, P. (982. Ordinary Differential Equations, Second Edition, Boston: Birkhäuser. Hirsch, M. W., & Smale, S. (974. Differential equations, dynamical systems, and linear algebra. New York: Academic Press. Milnor, J. (2. Dynamics in One Complex Variable. 2nd edition. Wiesbaden:Vieweg & Sohn. Newcomb, T. M. (96. The acquaintance process. New York: Holt. Rinehart and Winston. Peitgen, H.-O., & Richter, P. H. (986. The Beauty of Fractals. New York: Springer-Verlag. Tobler, W. (976-77. Spatial interaction patterns. Journal of Environmental Systems, 6, 27-3. (995. Yadohisa, H., & Niki, N. (999. Vector field representation of asymmetric proximity data. Communications in Statistics, 28, 35-48.

3 A. g(t d dt (g(t e at =g (t e at + g(t( ae at =ag(t e at ag(t e at =. (A. g(t e at k g(t =e at k.2 dx/dt = Ax x = Py P dx/dt = d(py/dt = P dy/dt dx/dt = P dy/dt = A(Py dy/dt = P AP y