. (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u(

Size: px
Start display at page:

Download ". (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u("

Transcription

1 3 8. (.8.) Nermark β (.8.) (3.4.6) (.8.3) (.8.8) (Case ) (3.4.7) (..5)

2 . (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u(t) + ( t) 3 ü(t) + ( t)3 3! ü(t) + ( t) 6... u (t) + () ü(t + t) (3) t t + t ü u u t 3... ü(t + t) ü(t) u (t) = t (4) () u(t + t) = u(t) + t+ t t ü(t)dt = u(t) + t (ü(t + t) + ü(t)) (5) (3) ( f(t + t) c u(t) + t ) ( ) ü(t) k u(t) + t u(t) + ( t) ü(t) 3 ü(t + t) = m + c t + k ( t) 6 (6) (4) ü(t + t) = f(t + t) c u(t + t) = u(t) + t ( ü(t) + t u(t + t) =u(t) + t u(t) + ( t) 3 u(t) + t ) ( ) ü(t) k u(t) + t u(t) + ( t) ü(t) 3 m + c t + k ( t) 6 ü(t + t) (8) ü(t) + ( t) 6 (7) ü(t + t) (9)

3 .3 () t + t u Taylor t u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u(t) + ( t) 4 ü(t) + () ü(t) + ( t) 4 t ü(t + t) () ü(t) = ü(t + t) + ü(t) () () u(t + t) = u(t) + t+ t t ü(t)dt = u(t) + t (ü(t + t) + ü(t)) (3) (3) ( f(t + t) c u(t) + t ) ( ) ü(t) k u(t) + t u(t) + ( t) ü(t) 4 ü(t + t) = m + c t + k ( t) 4 (4) (4) ü(t + t) = f(t + t) c u(t + t) = u(t) + t ( ü(t) + t u(t + t) =u(t) + t u(t) + ( t) 4 u(t) + t ) ( ) ü(t) k u(t) + t u(t) + ( t) ü(t) 4 m + c t + k ( t) 4 (5) ü(t + t) (6) ü(t) + ( t) 4 ü(t + t) (7)

4 .4 Nermark β () t + t u u Taylor u(t + t) = u(t) + t! u(t + t) = u(t) + t! u(t) + ( t)! ü(t) + ( t)! ü(t) + ( t)3 3!... u (t) + (8)... u (t) + (9) 4 /3! = β 3 /! = γ u(t + t) = u(t) + t u(t) + ( t) ü(t) + β ( t) 3... u (t) () u(t + t) = u(t) + t ü(t) + γ ( t)... u (t) ()... u... ü(t + t) ü(t) u (t) = t () ( ) u(t + t) = u(t) + t u(t) + β ( t) ü(t) + β( t) ü(t + t) (3) u(t + t) = u(t) + ( γ) t ü(t) + γ t ü(t + t) (4) γ = / Newmark β Newmark β ü(t + t) = f(t + t) c ( u(t) + t ) ( ( ) ) ü(t) k u(t) + t u(t) + β ( t) ü(t) m + c t + k β ( t) u(t + t) = u(t) + t t ü(t) + ü(t + t) (6) ( ) u(t + t) =u(t) + t u(t) + β ( t) ü(t) + β ( t) ü(t + t) (7) (5) β = 4 β = 6 3

5 () Newmark β γ β t T min π γ/ β (8) T min β γ (9) γ / γ = / γ = / β = /4 γ = / β = /6 t T min π γ/ β =.553 (3) t =.(sec) T min.8(sec).(sec).(sec) FEM.8(sec) 4

6 .. (.8.) φ u m ( φ + ü) + c u + k u = (3) φ ü u u φ + ü u u m ü + h ω u + ω u = φ (3) h ω m c k T c = h k k m ω = m T = π m = π ω k φ h T m = (33) m ü + c u + k u = m φ (34) k = 4π m T c = h k m (35) t T = t (sec) h =.5 ü + φ ü φ subroutine IACC subroutine CRAC FORTRAN Fortran9 ( ) KiK-net sabroutine ERES Fortran9 5

7 Acceleration (gal) Acceleration (gal) Acceleration (gal) 96.7 TCGH6_EW time (sec) TCGH6_NS time (sec) 87.7 TCGH6_UD time (sec) Damping factor h=.5 Response acceleration (gal) TCGH6_EW TCGH6_NS TCGH6_UD.. Period (sec) 6

8 . 3 (3.4.6) 3 S a π T S v S d T π S v (36) S a : S v S d : : 5 5 Damping factor h=.5 Response acceleration (gal) 5 5 Response velocity (kine) Response displacement (cm) Period (sec) TCGH6_EW TCGH6_NS TCGH6_UD

9 3. (.8.3) 3. () x(t) N t x m m =,,,, N N N = L L x m k = N/ N/ { ( ) ( )} πk πk x m = A k cos N t t + B k sin N t t k= = A + N/ k= ( A k cos πkm N + B k sin πkm ) + A n/ cos π(n/)m N N m =,,,, N t = m t x(t) A k, B k (Fourier transform) (Fourier inverse transform) A k = N B k = N N m= N m= x m cos πkm N x m sin πkm N k =,,,, N/, N/ (37) k =,,, N/ (38) () A = N x m ( ) (39) N m= e ±iθ = cos θ ± i sin θ C k = N N m= x m e i(πkm/n) k =,,,, N (4) N x m = C k e +i(πkm/n) m =,,,, N (4) k= (3) C k = A k ib k A k + ib k = A N k ib N k A k = Re(C k ) B k = Im(C k ) k =,,,, N/ 8

10 (4) f k = k N t F k = N t φ k = arctan N/ k =,,,, N/ Ak + B k = N t C k = N t Re(C k ) + Im(C k ) ( B ) k = arctan A k ˆx(t) = cos(π f k t + φ k ) k= ( ) Im(Ck ) Re(C k ) ( π < φ k < π) ( ) 3. Fourier data Fortran9 random number 6 Re(Ck) Im(Ck) FFT abs(ck) fp(ck) Re(xm) Im(xm) i data Re(Ck) Im(Ck) abs(ck) fp(ck) Re(xm) Im(xm) Ave..478 (i=) n n/+(i=8) (i=) n n/+(i=8) 9

11 BASIC Fortran9 (4) N m= x m e i(πkm/n) N (4) (4) module d e f p i implicit none real ( 8 ), parameter : : p i = D end module d e f p i program f9 FFT!! FFT and I n v e r s e FFT! use d e f p i implicit none integer : : i! integer : : ndata! Number o f o r i g i n a l data integer : : nn! Number o f c a l c u l a t e d data ( powers o f ) real ( 8 ) : : dt! Time increment real ( 8 ), allocatable : : accinp ( : )! I n p u t t e d o r i g i n a l data real ( 8 ), allocatable : : xr ( : )! Real p a r t o f c a l c u l a t e d data real ( 8 ), allocatable : : x i ( : )! Imaginary p a r t o f c a l c u l a t e d data real ( 8 ), allocatable : : cr ( : )! Memory o f r e a l p a r t real ( 8 ), allocatable : : c i ( : )! Remory o f imaginary p a r t real ( 8 ), allocatable : : ck ( : )! A bsolute o f complex Fourier c o e f f i c i e n t real ( 8 ), allocatable : : fp ( : )! Fourier phase spectrum character : : fnamer 5, fnamew 5 character : : strcom,dummy 5, fmt character : : l i n e b u f f c a l l g e t a r g (, fnamer )! Input data f i l e name c a l l g e t a r g (,fnamew)! Input output f i l e name open (, f i l e=fnamer, status= old ) read (, ( a ) ) strcom read (, ) dummy, dt read (, ) dummy, ndata! To s e t number o f data! (nn i s e q u a l to to t h e power o f p o s i t i v e i n t e g e r number! and nn i s g r e a t e r than or e q u a l to number o f data ) nn= do nn=nn i f ( ndata <=nn ) exit allocate ( accinp ( : nn ) ) allocate ( xr ( : nn ) ) allocate ( x i ( : nn ) ) allocate ( cr ( : nn ) ) allocate ( c i ( : nn ) ) allocate ( ck ( : nn ) ) allocate ( fp ( : nn ) )

12 stop do do i =,nn accinp ( i )=.D xr ( i )=.D x i ( i )=.D i =, ndata read (, ) accinp ( i ) xr ( i )= accinp ( i ) close ( ) c a l l FFT( nn, xr, x i )! Fourier transform do i =,nn! To d i v i d e r eturned v a l u e s by number o f data xr ( i )=xr ( i )/ dble ( nn ) x i ( i )= x i ( i )/ dble ( nn ) c r ( i )=xr ( i )! To memory r e a l p a r t o f transformed v a l u e s c i ( i )= x i ( i )! To memory imaginal p a r t o f transformed v a l u e s ck ( i )= s q r t ( cr ( i ). D+c i ( i ). D) fp ( i )=.D i f (. D 3<abs ( cr ( i ) ) ) then fp ( i )=atan ( c i ( i )/ cr ( i ) ) / pi 8.D i f ( cr ( i ) <.D. and. c i ( i ) <.D) fp ( i )= fp ( i ) 8.D i f ( cr ( i ) <.D. and.. D<c i ( i ) ) fp ( i )=8.D+fp ( i ) else i f (. D<c i ( i ). and.. D<cr ( i ) ) fp ( i )=9.D i f ( c i ( i ) <.D. and. cr ( i ) <.D) fp ( i )= 9.D i f (. D<c i ( i ). and. cr ( i ) <.D) fp ( i )=9.D i f ( c i ( i ) <.D. and.. D<cr ( i ) ) fp ( i )= 9.D end i f do i =,nn! To change t h e s i g n o f imaginary p a r t s f o r i n v e r s e FFT xr ( i )=xr ( i ) x i ( i )= x i ( i ) Call FFT( nn, xr, x i )! I n v e r s e Fourier transform fmt= ( i5,,, e5. 7, 4 (,, e5. 7 ), (,, e5. 7 ) ) open (, f i l e=fnamew, status= r e p l a c e ) write (, ( a ) ) i, data, Re(Ck), Im(Ck), abs (Ck), fp (Ck), Re(xm), Im(xm) do i =,nn write ( l i n e b u f f,fmt=fmt ) i, accinp ( i ), cr ( i ), c i ( i ), ck ( i ), fp ( i ), xr ( i ), x i ( i ) c a l l d e l s p a c e s ( l i n e b u f f ) write (, ( a ) ) trim ( l i n e b u f f ) close ( ) contains subroutine FFT( nn, xr, x i )!! Fast Fourier Transform, I n v e r s e transform!! nn : Number o f data ( powers o f )! xr ( ) : Real p a r t o f i /o data! x i ( ) : Imaginary p a r t o f i /o data integer, intent ( in ) : : nn real ( 8 ), intent ( inout ) : : xr ( : nn ) real ( 8 ), intent ( inout ) : : x i ( : nn ) integer : : g, h, i, j, k integer : : l,m, n, p, q real ( 8 ) : : a, b, xd real ( 8 ), allocatable : : s ( : ) real ( 8 ), allocatable : : c ( : ) n=nn allocate ( s ( : n/+)) allocate ( c ( : n/+)) i =; j =;k=; l =;p=;h=;g=;q= m=i n t ( l o g ( dble ( n ) ) / l o g (. D)+.D) a=.d

13 b=p i. D/ dble ( n ) do i =,n/ s ( i +)= s i n ( a ) c ( i +)=cos ( a ) a=a+b l=n h= do g=,m l=l / k= do q=,h p= do i=k, l+k j=i+l a=xr ( i +) xr ( j +) b=x i ( i +) x i ( j +) xr ( i +)=xr ( i +)+xr ( j +) x i ( i +)=xi ( i +)+xi ( j +) i f ( p==)then xr ( j +)=a x i ( j +)=b else xr ( j +)=a c ( p+)+b s ( p+) x i ( j +)=b c ( p+) a s ( p+) end p=p+h k=k+l+l h=h+h j=n/ do i f i =,n k=n i f ( j<i ) then xd=xr ( i +) xr ( i +)=xr ( j +) xr ( j +)=xd xd=x i ( i +) x i ( i +)=xi ( j +) x i ( j +)=xd end i f k=k/ do while ( j>=k ) j=j k k=k/ i f ( k==)exit j=j+k end subroutine FFT subroutine d e l s p a c e s ( s ) character ( ), intent ( inout ) : : s character ( len=len ( s ) ) tmp integer i, j j = do i =, len ( s ) i f ( s ( i : i )== ) cycle tmp( j : j ) = s ( i : i ) j = j + s = tmp ( : j ) end subroutine d e l s p a c e s end program f9 FFT

14 3.3 f k = k (k = N/) N t F k = N t Re(C k ) + Im(C k ) (k = N/) 4 4 Parzen window Parzen winndow subroutine SWIN FORTRAN Fortran9 3

15 Band=.Hz TCGH6_EW Fourier spectrum (gal*sec). 5 5 Frequency (Hz) Band=.Hz TCGH6_NS Fourier spectrum (gal*sec). 5 5 Frequency (Hz) Band=.Hz TCGH6_UD Fourier spectrum (gal*sec). 5 5 Frequency (Hz) 4 4

16 4. (.8.8) % Fourier S a ( ) S a r k = S a /S a e r (r k ) ɛ Yes n ( r k ) k= e r (r k ) = n ɛ e rror S a r k = S a S a ɛ No Fourier S a n Fourier C k Fourier C k = C k r k Fourier 5 5

17 4. ( ) ( ) 4 3 Case a(sec) b(sec) c(sec) e(t) = (t/a) e(t) =. e(t) = exp{ α (t b)} e(t) =. at t = c α = ln(.) c b Envelop function of acceleration wave due to earthquake e(t) a=5 b=5 b=5 b=35 c=3 c=6 c= e= time t (sec) 6 Fortran9 random number t =.(sec) gal 3gal 6

18 Acceleration (gal) Acceleration (gal) 4 Original time (sec) 4 Simulated time (sec) Damping factor h=.5 Response acceleration (gal). Target Simulated Original... Period (sec) 7 Case ( C=3 sec) 7

19 Acceleration (gal) Acceleration (gal) 4 Original time (sec) 4 Simulated time (sec) Damping factor h=.5 Response acceleration (gal). Target Simulated Original... Period (sec) 8 Case ( C=6 sec) 8

20 Acceleration (gal) Acceleration (gal) 4 Original time (sec) 4 Simulated time (sec) Damping factor h=.5 Response acceleration (gal). Target Simulated Original... Period (sec) 9 Case 3 ( C= sec) 9

21 4.3 (Case ) (3.4.7) Phase angle (degree) 8 Fourier Phase Spectrum (simulated wave) Frequency (Hz) Acceleration (gal) Acceleration data (input wave) Time (sec) Acceleration (gal) Phase wave Acceleration data (simulated wave) Time (sec) (Case )

22 5. (..5) ü(t) + h ω u(t) + ω u(t) = a(t) (4) u a h ω a(t) = π u(t) = π u(t) = π ü(t) = π (43)(44)(45)(46) (4) A(ω) U(ω) = (ω ω ) hω ω i u(t) = π u(t) = π ü(t) = π H(ω) = a A(ω) H(ω) A(ω)e iωt dω (43) U(ω)e iωt dω (44) iω U(ω)e iωt dω (45) ω U(ω)e iωt dω (46) (47) A(ω) (ω ω ) hω ω i eiωt dω (48) iωa(ω) (ω ω ) hω ω i eiωt dω (49) ω A(ω) (ω ω ) hω ω i eiωt dω (5) (ω ω ) hω ω i 3 A(ω) H(ω) u(t) 4 iω A(ω) H(ω) u(t) 5 ω A(ω) H(ω) ü(t) N N t f k ω N N/+ N/+ N/+ N f k = k N t (5) ω k = πf k k = N/ (5) N/ +

23 a(t) Input acceleration (gal) Time (sec) 5 4 A(ω) A Frequency (Hz) H(ω) H Frequency (Hz) 5 4 ω A(ω) H(ω) ω *A*H Frequency (Hz) a(t) + ü(t) Responce acceleration (gal) Time (sec) ( h=.5 Hz)

24 Damping factor h=.5 Response acceleration (gal) TCGH6_EW TCGH6_NS TCGH6_UD.. Period (sec) ( ) KiK-net sabroutine ERES Fortran9 3

1. ( ) 1.1 t + t [m]{ü(t + t)} + [c]{ u(t + t)} + [k]{u(t + t)} = {f(t + t)} (1) m ü f c u k u 1.2 Newmark β (1) (2) ( [m] + t ) 2 [c] + β( t)2

1. ( ) 1.1 t + t [m]{ü(t + t)} + [c]{ u(t + t)} + [k]{u(t + t)} = {f(t + t)} (1) m ü f c u k u 1.2 Newmark β (1) (2) ( [m] + t ) 2 [c] + β( t)2 212 1 6 1. (212.8.14) 1 1.1............................................. 1 1.2 Newmark β....................... 1 1.3.................................... 2 1.4 (212.8.19)..................................

More information

(5 B m e i 2π T mt m m B m e i 2π T mt m m B m e i 2π T mt B m (m < 0 C m m (6 (7 (5 g(t C 0 + m C m e i 2π T mt (7 C m e i 2π T mt + m m C m e i 2π T

(5 B m e i 2π T mt m m B m e i 2π T mt m m B m e i 2π T mt B m (m < 0 C m m (6 (7 (5 g(t C 0 + m C m e i 2π T mt (7 C m e i 2π T mt + m m C m e i 2π T 2.6 FFT(Fast Fourier Transform 2.6. T g(t g(t 2 a 0 + { a m b m 2 T T 0 2 T T 0 (a m cos( 2π T mt + b m sin( 2π mt ( T m 2π g(t cos( T mtdt m 0,, 2,... 2π g(t sin( T mtdt m, 2, 3... (2 g(t T 0 < t < T

More information

Microsoft Word - 資料 (テイラー級数と数値積分).docx

Microsoft Word - 資料 (テイラー級数と数値積分).docx δx δx n x=0 sin x = x x3 3 + x5 5 x7 7 +... x ak = (-mod(k,2))**(k/2) / fact_k ( ) = a n δ x n f x 0 + δ x a n = f ( n) ( x 0 ) n f ( x) = sin x n=0 58 I = b a ( ) f x dx ΔS = f ( x)h I = f a h h I = h

More information

Microsoft Word - 03-数値計算の基礎.docx

Microsoft Word - 03-数値計算の基礎.docx δx f x 0 + δ x n=0 a n = f ( n) ( x 0 ) n δx n f x x=0 sin x = x x3 3 + x5 5 x7 7 +... x ( ) = a n δ x n ( ) = sin x ak = (-mod(k,2))**(k/2) / fact_k 10 11 I = f x dx a ΔS = f ( x)h I = f a h I = h b (

More information

untitled

untitled Fortran90 ( ) 17 12 29 1 Fortran90 Fortran90 FORTRAN77 Fortran90 1 Fortran90 module 1.1 Windows Windows UNIX Cygwin (http://www.cygwin.com) C\: Install Cygwin f77 emacs latex ps2eps dvips Fortran90 Intel

More information

power.tex

power.tex Contents ii 1... 1... 1... 7... 7 3 (DFFT).................................... 8 4 (CIFT) DFFT................................ 10 5... 13 6... 16 3... 0 4... 0 5... 0 6... 0 i 1987 SN1987A 0.5 X SN1987A

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

1 u t = au (finite difference) u t = au Von Neumann

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

Microsoft Word - 信号処理3.doc

Microsoft Word - 信号処理3.doc Junji OHTSUBO 2012 FFT FFT SN sin cos x v ψ(x,t) = f (x vt) (1.1) t=0 (1.1) ψ(x,t) = A 0 cos{k(x vt) + φ} = A 0 cos(kx ωt + φ) (1.2) A 0 v=ω/k φ ω k 1.3 (1.2) (1.2) (1.2) (1.1) 1.1 c c = a + ib, a = Re[c],

More information

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

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 + ( IA 2013 : :10722 : 2 : :2 :761 :1 23-27) : : 1 1.1 / ) 1 /, ) / e.g. Taylar ) e x = 1 + x + x2 2 +... + xn n! +... sin x = x x3 6 + x5 x2n+1 + 1)n 5! 2n + 1)! 2 2.1 = 1 e.g. 0 = 0.00..., π = 3.14..., 1

More information

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

() 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 ( 3 n nc k+ k + 3 () n C r n C n r nc r C r + C r ( r n ) () 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 + (4) n C n n C + n C + n C + + n C n (5) k k n C k n C k (6) n C + nc

More information

[1] 1.1 x(t) t x(t + n ) = x(t) (n = 1,, 3, ) { x(t) : : 1 [ /, /] 1 x(t) = a + a 1 cos πt + a cos 4πt + + a n cos nπt + + b 1 sin πt + b sin 4πt = a

[1] 1.1 x(t) t x(t + n ) = x(t) (n = 1,, 3, ) { x(t) : : 1 [ /, /] 1 x(t) = a + a 1 cos πt + a cos 4πt + + a n cos nπt + + b 1 sin πt + b sin 4πt = a 13/7/1 II ( / A: ) (1) 1 [] (, ) ( ) ( ) ( ) etc. etc. 1. 1 [1] 1.1 x(t) t x(t + n ) = x(t) (n = 1,, 3, ) { x(t) : : 1 [ /, /] 1 x(t) = a + a 1 cos πt + a cos 4πt + + a n cos nπt + + b 1 sin πt + b sin

More information

コンピュータ概論

コンピュータ概論 4.1 For Check Point 1. For 2. 4.1.1 For (For) For = To Step (Next) 4.1.1 Next 4.1.1 4.1.2 1 i 10 For Next Cells(i,1) Cells(1, 1) Cells(2, 1) Cells(10, 1) 4.1.2 50 1. 2 1 10 3. 0 360 10 sin() 4.1.2 For

More information

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

3. :, 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 information

数値計算:フーリエ変換

数値計算:フーリエ変換 ( ) 1 / 72 1 8 2 3 4 ( ) 2 / 72 ( ) 3 / 72 ( ) 4 / 72 ( ) 5 / 72 sample.m Fs = 1000; T = 1/Fs; L = 1000; t = (0:L-1)*T; % Sampling frequency % Sample time % Length of signal % Time vector y=1+0.7*sin(2*pi*50*t)+sin(2*pi*120*t)+2*randn(size(t));

More information

SFGÇÃÉXÉyÉNÉgÉãå`.pdf

SFGÇÃÉXÉyÉNÉgÉãå`.pdf SFG 1 SFG SFG I SFG (ω) χ SFG (ω). SFG χ χ SFG (ω) = χ NR e iϕ +. ω ω + iγ SFG φ = ±π/, χ φ = ±π 3 χ SFG χ SFG = χ NR + χ (ω ω ) + Γ + χ NR χ (ω ω ) (ω ω ) + Γ cosϕ χ NR χ Γ (ω ω ) + Γ sinϕ. 3 (θ) 180

More information

Microsoft Word - 資料 docx

Microsoft Word - 資料 docx y = Asin 2πt T t t = t i i 1 n+1 i i+1 Δt t t i = Δt i 1 ( ) y i = Asin 2πt i T 29 (x, y) t ( ) x = Asin 2πmt y = Asin( 2πnt + δ ) m, n δ (x, y) m, n 30 L A x y A L x 31 plot sine curve program sine implicit

More information

Part () () Γ Part ,

Part () () Γ Part , Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35

More information

LLG-R8.Nisus.pdf

LLG-R8.Nisus.pdf d M d t = γ M H + α M d M d t M γ [ 1/ ( Oe sec) ] α γ γ = gµ B h g g µ B h / π γ g = γ = 1.76 10 [ 7 1/ ( Oe sec) ] α α = λ γ λ λ λ α γ α α H α = γ H ω ω H α α H K K H K / M 1 1 > 0 α 1 M > 0 γ α γ =

More information

chap1.dvi

chap1.dvi 1 1 007 1 e iθ = cos θ + isin θ 1) θ = π e iπ + 1 = 0 1 ) 3 11 f 0 r 1 1 ) k f k = 1 + r) k f 0 f k k = 01) f k+1 = 1 + r)f k ) f k+1 f k = rf k 3) 1 ) ) ) 1+r/)f 0 1 1 + r/) f 0 = 1 + r + r /4)f 0 1 f

More information

1 7 ω ω ω 7.1 0, ( ) Q, 7.2 ( Q ) 7.1 ω Z = R +jx Z 1/ Z 7.2 ω 7.2 Abs. admittance (x10-3 S) RLC Series Circuit Y R = 20 Ω L = 100

1 7 ω ω ω 7.1 0, ( ) Q, 7.2 ( Q ) 7.1 ω Z = R +jx Z 1/ Z 7.2 ω 7.2 Abs. admittance (x10-3 S) RLC Series Circuit Y R = 20 Ω L = 100 7 7., ) Q, 7. Q ) 7. Z = R +jx Z / Z 7. 7. Abs. admittance x -3 S) 5 4 3 R Series ircuit Y R = Ω = mh = uf Q = 5 5 5 V) Z = R + jx 7. Z 7. ) R = Ω = mh = µf ) 7 V) R Z s = R + j ) 7.3 R =. 7.4) ) f = π.

More information

3. :, 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 +

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

) ] [ 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

) ] [ 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 information

1. (8) (1) (x + y) + (x + y) = 0 () (x + y ) 5xy = 0 (3) (x y + 3y 3 ) (x 3 + xy ) = 0 (4) x tan y x y + x = 0 (5) x = y + x + y (6) = x + y 1 x y 3 (

1. (8) (1) (x + y) + (x + y) = 0 () (x + y ) 5xy = 0 (3) (x y + 3y 3 ) (x 3 + xy ) = 0 (4) x tan y x y + x = 0 (5) x = y + x + y (6) = x + y 1 x y 3 ( 1 1.1 (1) (1 + x) + (1 + y) = 0 () x + y = 0 (3) xy = x (4) x(y + 3) + y(y + 3) = 0 (5) (a + y ) = x ax a (6) x y 1 + y x 1 = 0 (7) cos x + sin x cos y = 0 (8) = tan y tan x (9) = (y 1) tan x (10) (1 +

More information

all.dvi

all.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 information

A = A x x + A y y + A, B = B x x + B y y + B, C = C x x + C y y + C..6 x y A B C = A x x + A y y + A B x B y B C x C y C { B = A x x + A y y + A y B B

A = A x x + A y y + A, B = B x x + B y y + B, C = C x x + C y y + C..6 x y A B C = A x x + A y y + A B x B y B C x C y C { B = A x x + A y y + A y B B 9 7 A = A x x + A y y + A, B = B x x + B y y + B, C = C x x + C y y + C..6 x y A B C = A x x + A y y + A B x B y B C x C y C { B = A x x + A y y + A y B B x x B } B C y C y + x B y C x C C x C y B = A

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

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

(Compton Scattering) Beaming 1 exp [i (k x ωt)] k λ k = 2π/λ ω = 2πν k = ω/c k x ωt ( ω ) k α c, k k x ωt η αβ k α x β diag( + ++) x β = (ct, x) O O x

(Compton Scattering) Beaming 1 exp [i (k x ωt)] k λ k = 2π/λ ω = 2πν k = ω/c k x ωt ( ω ) k α c, k k x ωt η αβ k α x β diag( + ++) x β = (ct, x) O O x Compton Scattering Beaming exp [i k x ωt] k λ k π/λ ω πν k ω/c k x ωt ω k α c, k k x ωt η αβ k α x β diag + ++ x β ct, x O O x O O v k α k α β, γ k γ k βk, k γ k + βk k γ k k, k γ k + βk 3 k k 4 k 3 k

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

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A .. Laplace ). A... i),. ω i i ). {ω,..., ω } Ω,. ii) Ω. Ω. A ) r, A P A) P A) r... ).. Ω {,, 3, 4, 5, 6}. i i 6). A {, 4, 6} P A) P A) 3 6. ).. i, j i, j) ) Ω {i, j) i 6, j 6}., 36. A. A {i, j) i j }.

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

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

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

More information

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

Korteweg-de Vries

Korteweg-de Vries Korteweg-de Vries 2011 03 29 ,.,.,.,, Korteweg-de Vries,. 1 1 3 1.1 K-dV........................ 3 1.2.............................. 4 2 K-dV 5 2.1............................. 5 2.2..............................

More information

(interferometer) 1 N *3 2 ω λ k = ω/c = 2π/λ ( ) r E = A 1 e iφ1(r) e iωt + A 2 e iφ2(r) e iωt (1) φ 1 (r), φ 2 (r) r λ 2π 2 I = E 2 = A A 2 2 +

(interferometer) 1 N *3 2 ω λ k = ω/c = 2π/λ ( ) r E = A 1 e iφ1(r) e iωt + A 2 e iφ2(r) e iωt (1) φ 1 (r), φ 2 (r) r λ 2π 2 I = E 2 = A A 2 2 + 7 1 (Young) *1 *2 (interference) *1 (1802 1804) *2 2 (2005) (1993) 1 (interferometer) 1 N *3 2 ω λ k = ω/c = 2π/λ ( ) r E = A 1 e iφ1(r) e iωt + A 2 e iφ2(r) e iωt (1) φ 1 (r), φ 2 (r) r λ 2π 2 I = E 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 f f x, y, u, v, r, θ r > = x + iy, f = u + iv C γ D f f D f f, Rm,. = x + iy = re iθ = r cos θ + i sin θ = x iy = re iθ = r cos θ i sin θ x = + = Re, y = = Im i r = = = x + y θ = arg = arctan y x e i =

More information

#A A A F, F d F P + F P = d P F, F y P F F x A.1 ( α, 0), (α, 0) α > 0) (x, y) (x + α) 2 + y 2, (x α) 2 + y 2 d (x + α)2 + y 2 + (x α) 2 + y 2 =

#A A A F, F d F P + F P = d P F, F y P F F x A.1 ( α, 0), (α, 0) α > 0) (x, y) (x + α) 2 + y 2, (x α) 2 + y 2 d (x + α)2 + y 2 + (x α) 2 + y 2 = #A A A. F, F d F P + F P = d P F, F P F F A. α, 0, α, 0 α > 0, + α +, α + d + α + + α + = d d F, F 0 < α < d + α + = d α + + α + = d d α + + α + d α + = d 4 4d α + = d 4 8d + 6 http://mth.cs.kitmi-it.c.jp/

More information

main.dvi

main.dvi 4 DFT DFT Fast Fourier Transform: FFT 4.1 DFT IDFT X(k) = 1 n=0 x(n)e j2πkn (4.1) 1 x(n) = 1 X(k)e j2πkn (4.2) k=0 x(n) X(k) DFT 2 ( 1) 2 4 2 2(2 1) 2 O( 2 ) 4.2 FFT 4.2.1 radix2 FFT 1 (4.1) 86 4. X(0)

More information

Kroneher Levi-Civita 1 i = j δ i j = i j 1 if i jk is an even permutation of 1,2,3. ε i jk = 1 if i jk is an odd permutation of 1,2,3. otherwise. 3 4

Kroneher Levi-Civita 1 i = j δ i j = i j 1 if i jk is an even permutation of 1,2,3. ε i jk = 1 if i jk is an odd permutation of 1,2,3. otherwise. 3 4 [2642 ] Yuji Chinone 1 1-1 ρ t + j = 1 1-1 V S ds ds Eq.1 ρ t + j dv = ρ t dv = t V V V ρdv = Q t Q V jdv = j ds V ds V I Q t + j ds = ; S S [ Q t ] + I = Eq.1 2 2 Kroneher Levi-Civita 1 i = j δ i j =

More information

[ ] 0.1 lim x 0 e 3x 1 x IC ( 11) ( s114901) 0.2 (1) y = e 2x (x 2 + 1) (2) y = x/(x 2 + 1) 0.3 dx (1) 1 4x 2 (2) e x sin 2xdx (3) sin 2 xdx ( 11) ( s

[ ] 0.1 lim x 0 e 3x 1 x IC ( 11) ( s114901) 0.2 (1) y = e 2x (x 2 + 1) (2) y = x/(x 2 + 1) 0.3 dx (1) 1 4x 2 (2) e x sin 2xdx (3) sin 2 xdx ( 11) ( s [ ]. lim e 3 IC ) s49). y = e + ) ) y = / + ).3 d 4 ) e sin d 3) sin d ) s49) s493).4 z = y z z y s494).5 + y = 4 =.6 s495) dy = 3e ) d dy d = y s496).7 lim ) lim e s49).8 y = e sin ) y = sin e 3) y =

More information

(Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1

(Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1 (Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1 17 Fortran Formular Tranlator Lapack Fortran FORTRAN, FORTRAN66, FORTRAN77, FORTRAN90, FORTRAN95 17.1 A Z ( ) 0 9, _, =, +, -, *,

More information

情報活用資料

情報活用資料 y = Asin 2πt T t t = t i i 1 n+1 i i+1 Δt t t i = Δt i 1 ( ) y i = Asin 2πt i T 21 (x, y) t ( ) x = Asin 2πmt y = Asin( 2πnt + δ ) m, n δ (x, y) m, n 22 L A x y A L x 23 ls -l gnuplot gnuplot> plot "sine.dat"

More information

29

29 9 .,,, 3 () C k k C k C + C + C + + C 8 + C 9 + C k C + C + C + C 3 + C 4 + C 5 + + 45 + + + 5 + + 9 + 4 + 4 + 5 4 C k k k ( + ) 4 C k k ( k) 3 n( ) n n n ( ) n ( ) n 3 ( ) 3 3 3 n 4 ( ) 4 4 4 ( ) n n

More information

. ev=,604k m 3 Debye ɛ 0 kt e λ D = n e n e Ze 4 ln Λ ν ei = 5.6π / ɛ 0 m/ e kt e /3 ν ei v e H + +e H ev Saha x x = 3/ πme kt g i g e n

. ev=,604k m 3 Debye ɛ 0 kt e λ D = n e n e Ze 4 ln Λ ν ei = 5.6π / ɛ 0 m/ e kt e /3 ν ei v e H + +e H ev Saha x x = 3/ πme kt g i g e n 003...............................3 Debye................. 3.4................ 3 3 3 3. Larmor Cyclotron... 3 3................ 4 3.3.......... 4 3.3............ 4 3.3...... 4 3.3.3............ 5 3.4.........

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

1. 4cm 16 cm 4cm 20cm 18 cm L λ(x)=ax [kg/m] A x 4cm A 4cm 12 cm h h Y 0 a G 0.38h a b x r(x) x y = 1 h 0.38h G b h X x r(x) 1 S(x) = πr(x) 2 a,b, h,π

1. 4cm 16 cm 4cm 20cm 18 cm L λ(x)=ax [kg/m] A x 4cm A 4cm 12 cm h h Y 0 a G 0.38h a b x r(x) x y = 1 h 0.38h G b h X x r(x) 1 S(x) = πr(x) 2 a,b, h,π . 4cm 6 cm 4cm cm 8 cm λ()=a [kg/m] A 4cm A 4cm cm h h Y a G.38h a b () y = h.38h G b h X () S() = π() a,b, h,π V = ρ M = ρv G = M h S() 3 d a,b, h 4 G = 5 h a b a b = 6 ω() s v m θ() m v () θ() ω() dθ()

More information

SO(3) 49 u = Ru (6.9), i u iv i = i u iv i (C ) π π : G Hom(V, V ) : g D(g). π : R 3 V : i 1. : u u = u 1 u 2 u 3 (6.10) 6.2 i R α (1) = 0 cos α

SO(3) 49 u = Ru (6.9), i u iv i = i u iv i (C ) π π : G Hom(V, V ) : g D(g). π : R 3 V : i 1. : u u = u 1 u 2 u 3 (6.10) 6.2 i R α (1) = 0 cos α SO(3) 48 6 SO(3) t 6.1 u, v u = u 1 1 + u 2 2 + u 3 3 = u 1 e 1 + u 2 e 2 + u 3 e 3, v = v 1 1 + v 2 2 + v 3 3 = v 1 e 1 + v 2 e 2 + v 3 e 3 (6.1) i (e i ) e i e j = i j = δ ij (6.2) ( u, v ) = u v = ij

More information

85 4

85 4 85 4 86 Copright c 005 Kumanekosha 4.1 ( ) ( t ) t, t 4.1.1 t Step! (Step 1) (, 0) (Step ) ±V t (, t) I Check! P P V t π 54 t = 0 + V (, t) π θ : = θ : π ) θ = π ± sin ± cos t = 0 (, 0) = sin π V + t +V

More information

2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h)

2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h) 1 16 10 5 1 2 2.1 a a a 1 1 1 2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h) 4 2 3 4 2 5 2.4 x y (x,y) l a x = l cot h cos a, (3) y = l cot h sin a (4) h a

More information

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

I A A441 : April 15, 2013 Version : 1.1 I   Kawahira, Tomoki TA (Shigehiro, Yoshida ) I013 00-1 : April 15, 013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida) http://www.math.nagoya-u.ac.jp/~kawahira/courses/13s-tenbou.html pdf * 4 15 4 5 13 e πi = 1 5 0 5 7 3 4 6 3 6 10 6 17

More information

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

ω 0 m(ẍ + γẋ + ω0x) 2 = ee (2.118) e iωt x = e 1 m ω0 2 E(ω). (2.119) ω2 iωγ Z N P(ω) = χ(ω)e = exzn (2.120) ϵ = ϵ 0 (1 + χ) ϵ(ω) ϵ 0 = 1 +

ω 0 m(ẍ + γẋ + ω0x) 2 = ee (2.118) e iωt x = e 1 m ω0 2 E(ω). (2.119) ω2 iωγ Z N P(ω) = χ(ω)e = exzn (2.120) ϵ = ϵ 0 (1 + χ) ϵ(ω) ϵ 0 = 1 + 2.6 2.6.1 ω 0 m(ẍ + γẋ + ω0x) 2 = ee (2.118) e iωt x = e 1 m ω0 2 E(ω). (2.119) ω2 iωγ Z N P(ω) = χ(ω)e = exzn (2.120) ϵ = ϵ 0 (1 + χ) ϵ(ω) ϵ 0 = 1 + Ne2 m j f j ω 2 j ω2 iωγ j (2.121) Z ω ω j γ j f j

More information

1 3 1.1.......................... 3 1............................... 3 1.3....................... 5 1.4.......................... 6 1.5........................ 7 8.1......................... 8..............................

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

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

2014 3 10 5 1 5 1.1..................................... 5 2 6 2.1.................................... 6 2.2 Z........................................ 6 2.3.................................. 6 2.3.1..................

More information

ii

ii ii iii 1 1 1.1..................................... 1 1.2................................... 3 1.3........................... 4 2 9 2.1.................................. 9 2.2...............................

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

[1.1] r 1 =10e j(ωt+π/4), r 2 =5e j(ωt+π/3), r 3 =3e j(ωt+π/6) ~r = ~r 1 + ~r 2 + ~r 3 = re j(ωt+φ) =(10e π 4 j +5e π 3 j +3e π 6 j )e jωt

[1.1] r 1 =10e j(ωt+π/4), r 2 =5e j(ωt+π/3), r 3 =3e j(ωt+π/6) ~r = ~r 1 + ~r 2 + ~r 3 = re j(ωt+φ) =(10e π 4 j +5e π 3 j +3e π 6 j )e jωt 3.4.7 [.] =e j(t+/4), =5e j(t+/3), 3 =3e j(t+/6) ~ = ~ + ~ + ~ 3 = e j(t+φ) =(e 4 j +5e 3 j +3e 6 j )e jt = e jφ e jt cos φ =cos 4 +5cos 3 +3cos 6 =.69 sin φ =sin 4 +5sin 3 +3sin 6 =.9 =.69 +.9 =7.74 [.]

More information

1F90/kouhou_hf90.dvi

1F90/kouhou_hf90.dvi Fortran90 3 33 1 2 Fortran90 FORTRAN 1956 IBM IBM704 FORTRAN(FORmula TRANslation ) 1965 FORTRAN66 1978 FORTRAN77 1991 Fortran90 Fortran90 Fortran Fortran90 6 Fortran90 77 90 90 Fortran90 [ ] Fortran90

More information

all.dvi

all.dvi fortran 1996 4 18 2007 6 11 2012 11 12 1 3 1.1..................................... 3 1.2.............................. 3 2 fortran I 5 2.1 write................................ 5 2.2.................................

More information

N88 BASIC 0.3 C: My Documents 0.6: 0.3: (R) (G) : enterreturn : (F) BA- SIC.bas 0.8: (V) 0.9: 0.5:

N88 BASIC 0.3 C: My Documents 0.6: 0.3: (R) (G) : enterreturn : (F) BA- SIC.bas 0.8: (V) 0.9: 0.5: BASIC 20 4 10 0 N88 Basic 1 0.0 N88 Basic..................................... 1 0.1............................................... 3 1 4 2 5 3 6 4 7 5 10 6 13 7 14 0 N88 Basic 0.0 N88 Basic 0.1: N88Basic

More information

I

I I 6 4 10 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

untitled

untitled SPring-8 RFgun JASRI/SPring-8 6..7 Contents.. 3.. 5. 6. 7. 8. . 3 cavity γ E A = er 3 πε γ vb r B = v E c r c A B A ( ) F = e E + v B A A A A B dp e( v B+ E) = = m d dt dt ( γ v) dv e ( ) dt v B E v E

More information

TOP URL 1

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

More information

OHP.dvi

OHP.dvi 0 7 4 0000 5.. 3. 4. 5. 0 0 00 Gauss PC 0 Gauss 3 Gauss Gauss 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 u [] u [3] u [4] u [4] P 0 = P 0 (),3,4 (,), (3,), (4,) 0,,,3,4 3 3 3 3 4 4 4 4 0 3 6 6 0 6 3 6 0 6

More information

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

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

More information

2

2 16 1050026 1050042 1 2 1 1.1 3 1.2 3 1.3 3 2 2.1 4 2.2 4 2.2.1 5 2.2.2 5 2.3 7 2.3.1 1Basic 7 2.3.2 2 8 2.3.3 3 9 2.3.4 4window size 10 2.3.5 5 11 3 3.1 12 3.2 CCF 1 13 3.3 14 3.4 2 15 3.5 3 17 20 20 20

More information

Fortran90/95 [9]! (1 ) " " 5 "Hello!"! 3. (line) Fortran Fortran 1 2 * (1 ) 132 ( ) * 2 ( Fortran ) Fortran ,6 (continuation line) 1

Fortran90/95 [9]! (1 )   5 Hello!! 3. (line) Fortran Fortran 1 2 * (1 ) 132 ( ) * 2 ( Fortran ) Fortran ,6 (continuation line) 1 Fortran90/95 2.1 Fortran 2-1 Hello! 1 program example2_01! end program 2! first test program ( ) 3 implicit none! 4 5 write(*,*) "Hello!"! write Hello! 6 7 stop! 8 end program example2_01 1 program 1!

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

untitled

untitled MRR Physical Basis( 1.8.4) METEK MRR 1 MRR 1.1 MRR 24GHz FM-CW(frequency module continuous wave) 30 r+ r f+ f 1.2 1 4 MRR 24GHz 1.3 50mW 1 rf- (waveguide) (horn) 60cm ( monostatic radar) (continuous wave)

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

2000年度『数学展望 I』講義録

2000年度『数学展望 I』講義録 2000 I I IV I II 2000 I I IV I-IV. i ii 3.10 (http://www.math.nagoya-u.ac.jp/ kanai/) 2000 A....1 B....4 C....10 D....13 E....17 Brouwer A....21 B....26 C....33 D....39 E. Sperner...45 F....48 A....53

More information

2004 2005 2 2 1G01P038-0 1 2 1.1.............................. 2 1.2......................... 2 1.3......................... 3 2 4 2.1............................ 4 2.2....................... 4 2.3.......................

More information

1 1 3 ABCD ABD AC BD E E BD 1 : 2 (1) AB = AD =, AB AD = (2) AE = AB + (3) A F AD AE 2 = AF = AB + AD AF AE = t AC = t AE AC FC = t = (4) ABD ABCD 1 1

1 1 3 ABCD ABD AC BD E E BD 1 : 2 (1) AB = AD =, AB AD = (2) AE = AB + (3) A F AD AE 2 = AF = AB + AD AF AE = t AC = t AE AC FC = t = (4) ABD ABCD 1 1 ABCD ABD AC BD E E BD : () AB = AD =, AB AD = () AE = AB + () A F AD AE = AF = AB + AD AF AE = t AC = t AE AC FC = t = (4) ABD ABCD AB + AD AB + 7 9 AD AB + AD AB + 9 7 4 9 AD () AB sin π = AB = ABD AD

More information

9 5 ( α+ ) = (α + ) α (log ) = α d = α C d = log + C C 5. () d = 4 d = C = C = 3 + C 3 () d = d = C = C = 3 + C 3 =

9 5 ( α+ ) = (α + ) α (log ) = α d = α C d = log + C C 5. () d = 4 d = C = C = 3 + C 3 () d = d = C = C = 3 + C 3 = 5 5. 5.. A II f() f() F () f() F () = f() C (F () + C) = F () = f() F () + C f() F () G() f() G () = F () 39 G() = F () + C C f() F () f() F () + C C f() f() d f() f() C f() f() F () = f() f() f() d =

More information

nosenote3.dvi

nosenote3.dvi i 1 1 2 5 3 Verlet 9 4 18 5 23 6 26 1 1 1 MD t N r 1 (t), r 2 (t),, r N (t) ṙ 1 (t), ṙ 2 (t),, ṙ N (t) MD a 1, a 2, a 3 r i (i =1,,n) 1 2 T =0K r i + m 1 a 1 + m 2 a 2 + m 3 a 3 (m 1,m 2,m 3 =0, ±1, ±2,,

More information

(2-1) x, m, 2 N(m, 2 ) x REAL*8 FUNCTION NRMDST (X, M, V) X,M,V REAL*8 x, m, 2 X X N(0,1) f(x) standard-norm.txt normdist1.f x=0, 0.31, 0.5

(2-1) x, m, 2 N(m, 2 ) x REAL*8 FUNCTION NRMDST (X, M, V) X,M,V REAL*8 x, m, 2 X X N(0,1) f(x) standard-norm.txt normdist1.f x=0, 0.31, 0.5 2007/5/14 II II agata@k.u-tokyo.a.jp 0. 1. x i x i 1 x i x i x i x x+dx f(x)dx f(x) f(x) + 0 f ( x) dx = 1 (Probability Density Funtion 2 ) (normal distribution) 3 1 2 2 ( x m) / 2σ f ( x) = e 2πσ x m

More information

x y x-y σ x + τ xy + X σ y B = + τ xy + Y B = S x = σ x l + τ xy m S y = σ y m + τ xy l σ x σ y τ xy X B Y B S x S y l m δu δv [ ( σx δu + τ )

x y x-y σ x + τ xy + X σ y B = + τ xy + Y B = S x = σ x l + τ xy m S y = σ y m + τ xy l σ x σ y τ xy X B Y B S x S y l m δu δv [ ( σx δu + τ ) 1 8 6 No-tension 1. 1 1.1................................ 1 1............................................ 5.1 - [B].................................. 5................................. 6.3..........................................

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

I-2 (100 ) (1) y(x) y dy dx y d2 y dx 2 (a) y + 2y 3y = 9e 2x (b) x 2 y 6y = 5x 4 (2) Bernoulli B n (n = 0, 1, 2,...) x e x 1 = n=0 B 0 B 1 B 2 (3) co

I-2 (100 ) (1) y(x) y dy dx y d2 y dx 2 (a) y + 2y 3y = 9e 2x (b) x 2 y 6y = 5x 4 (2) Bernoulli B n (n = 0, 1, 2,...) x e x 1 = n=0 B 0 B 1 B 2 (3) co 16 I ( ) (1) I-1 I-2 I-3 (2) I-1 ( ) (100 ) 2l x x = 0 y t y(x, t) y(±l, t) = 0 m T g y(x, t) l y(x, t) c = 2 y(x, t) c 2 2 y(x, t) = g (A) t 2 x 2 T/m (1) y 0 (x) y 0 (x) = g c 2 (l2 x 2 ) (B) (2) (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

untitled

untitled ( ) c a sin b c b c a cos a c b c a tan b a b cos sin a c b c a ccos b csin (4) Ma k Mg a (Gal) g(98gal) (Gal) a max (K-E) kh Zck.85.6. 4 Ma g a k a g k D τ f c + σ tanφ σ 3 3 /A τ f3 S S τ A σ /A σ /A

More information

ohpr.dvi

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

More information

keisoku01.dvi

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

More information

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

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

More information

sp3.dvi

sp3.dvi 3 15 5 22 1 2 1.1... 2 1.2... 4 1.3... 5 1.4... 8 1.5 (Matlab )... 11 2 15 2.1... 15 2.2... 16 2.3... 17 3 19 3.1... 19 3.2... 2 3.3... 21 3.4... 22 3.5... 23 3.6... 24 3.7... 25 3.8 Daubechies... 26 4

More information

数学の基礎訓練I

数学の基礎訓練I I 9 6 13 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 3 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

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

数学Ⅱ演習(足助・09夏) II I 9/4/4 9/4/2 z C z z z z, z 2 z, w C zw z w 3 z, w C z + w z + w 4 t R t C t t t t t z z z 2 z C re z z + z z z, im z 2 2 3 z C e z + z + 2 z2 + 3! z3 + z!, I 4 x R e x cos x + sin x 2 z, w C e z+w

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

xyz,, uvw,, Bernoulli-Euler u c c c v, w θ x c c c dv ( x) dw uxyz (,, ) = u( x) y z + ω( yz, ) φ dx dx c vxyz (,, ) = v( x) zθ x ( x) c wxyz (,, ) =

xyz,, uvw,, Bernoulli-Euler u c c c v, w θ x c c c dv ( x) dw uxyz (,, ) = u( x) y z + ω( yz, ) φ dx dx c vxyz (,, ) = v( x) zθ x ( x) c wxyz (,, ) = ,, uvw,, Bernoull-Euler u v, w θ dv ( ) dw u (,, ) u( ) ω(, ) φ d d v (,, ) v( ) θ ( ) w (,, ) w( ) θ ( ) (11.1) ω φ φ dθ / dφ v v θ u w u w 11.1 θ θ θ 11. vw, (11.1) u du d v d w ε d d d u v ω γ φ w u

More information

B 1 B.1.......................... 1 B.1.1................. 1 B.1.2................. 2 B.2........................... 5 B.2.1.......................... 5 B.2.2.................. 6 B.2.3..................

More information

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

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

More information

m(ẍ + γẋ + ω 0 x) = ee (2.118) e iωt P(ω) = χ(ω)e = ex = e2 E(ω) m ω0 2 ω2 iωγ (2.119) Z N ϵ(ω) ϵ 0 = 1 + Ne2 m j f j ω 2 j ω2 iωγ j (2.120)

m(ẍ + γẋ + ω 0 x) = ee (2.118) e iωt P(ω) = χ(ω)e = ex = e2 E(ω) m ω0 2 ω2 iωγ (2.119) Z N ϵ(ω) ϵ 0 = 1 + Ne2 m j f j ω 2 j ω2 iωγ j (2.120) 2.6 2.6.1 mẍ + γẋ + ω 0 x) = ee 2.118) e iωt Pω) = χω)e = ex = e2 Eω) m ω0 2 ω2 iωγ 2.119) Z N ϵω) ϵ 0 = 1 + Ne2 m j f j ω 2 j ω2 iωγ j 2.120) Z ω ω j γ j f j f j f j sum j f j = Z 2.120 ω ω j, γ ϵω) ϵ

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

25 7 18 1 1 1.1 v.s............................. 1 1.1.1.................................. 1 1.1.2................................. 1 1.1.3.................................. 3 1.2................... 3

More information