Transcription of シミュレーション基礎(1) - cml-seikei.com
1 Mail : Tel : 0422-37-3722.. Mail : Tel : 0422-37-3706. 1. C OS .. MATLAB. 2. C . Y sin( X ) 0 X .. MATLAB . X=[ 0: : 2*pi ];. Y=sin(X);. plot(X,Y);. 3. C .. A=[1 2 3 ; 2 3 1; 3 1 2] ;. B =[-1 0 1]' ;. X=inv(A)*B ;. 4. Matlab .. MATLAB .. 2000 . ISBN978 4 434 15591 8.. 5.. 6.. MATLAB . 1+2 . >> 1+2.. ans =. 3. 7.. M .. >>ans = 3. ex1102 . >> ex1102. ans = 3. 8.. M . X=0:pi/100:2*pi;. Y=sin(X);. plot(X, Y);.. MATLAB . Sin SIN .. 9.. a+b a-b a .*b a ./ b .. ( .^). 2^10. 5! . prod(1:5). 10. #2.. i j 5. - 2*i 5. 2i .. X=10+6i; Y=5*2i Z1=X+Y, Z2=X.*Y, Z3= .. X=10+6i; Y=X'. 11.. yugekohei2341_vab . yuge yuge . Matlab .. 12.. ax 2 bx c 0.. b b 4ac 2. x . 2a . A= ; B= ; C= ;. X1=(-B+(B.^2-4.*A.*C).^ )./(2.*A). X2=(-B-(B.^2-4.*A.*C).^ )./(2.*A). EX10301 .. Matlab . 13.. A= ; B= ; C= ;. X1=(-B+(B.^2-4.*A.*C).^ )./(2.*A). X2=(-B-(B.^2-4.*A.*C).^ )./(2.*A). EX10301 .. Matlab.
2 14.. roots . P=[ ];. R=roots(P). ( ex1302 .m) . [ ] [ ] Matlab .. Roots . 15. [ ] . [ ] [ ] . A=[1 2 3], B=[-1 2 -3], c= . c+A=[3 4 5] A c . c.*A=[2 4 6] A c . [2 1 ] c A ( [ c/A(1) c/A(2) c/A(3) ] ). [ 1 ] A c ( [ A(1)/c A(2)/c A(3)/c ] ). A.^c=[1 4 9] A c . A+B=[0 4 0] . A.*B=[a1*b1 a2*b2 a3*b3]=[-1 4 -9] . A' A . C=dot(A,B) A B . D=cross(A,B) A B . 17.. V=(a b) c a=[2 1 0];. (1,1,2) b=[1 2 -1];. (1,2,-1) c=[1 1 2];.. b v=dot(cross(a,b),c). (0,0,0) a (2,1,0). 18.. 100. S 1 2 . 2 2. 100 i 2. 2. i 1. S=0;. for I=1:100. S=S+I.^2;. end S. ( : ) .. S=0; for I=1:100; S=S+I.^2; end S. 19. #2. sum . X=1:100; X=X.^2; S=sum(X).. X=1:100; -> 1 100 X . X=X.^2; -> X . 20.. A=100 R= N=10 .. N=0:10; A=1000000; R= ; F = A.*(1+R).^N , plot(N,F). 2 .. title( FUKURI KEISAN');xlabel( Kikan');ylabel( Ganri Goukei'). 22.. 1. MATLAB a = [ 9 ] . 2. a sin . 3. a 3 . 4. MATLAB b=[ 2 ] . 5.
3 A b . 6. a b . 7. a 2 . 8. 0 10 1 c . 9. 0 10 2 d . 10.. 20. S (k 4 3k 3 2). k 1. 23.. Ax = B. 1 2 3 1 .. A = 2 3 1 b = 0 .. 3 1 2 1 .. A=[1 2 3 ; 2 3 1; 3 1 2] ;. B =[-1 0 1]' ;. X=inv(A)*B. ; .. inv(A) . 24. #2.. A=[1 2 3 ; 2 3 1; 3 1 2] ;. B =[-1 0 1]' ;. X=A B.. 25. 1 .. mx kx 0. m, k . A, A . j nt j nt n k / m . x Ae Ae t 0 x x0 , x v0 . x0 v x v A j 0 , A 0 j 0 k 2 2 n 2 2 n m x v0. x x0 cos nt sin nt n 1 . 1 .. 0 5 . t=0 : : ; t . m= ; k=320;.. x0= ; v0= ;. wn=sqrt(k/m); v x x0 cos nt 0 sin nt x=x0*cos(wn*t)+v0/wn*sin(wn*t); n plot(t,x),xlabel( 'Time (s)' ),ylabel( 'Displacement (m)' ),grid plot(,) .. xlabel() . 34. 1 .. 35.