- PDF Free Download

Size: px

Start display at page:

Download ""

원태 범
5 years ago
Views:

2 예제 1.1 ( 행벡터만들기 ) >> a=[ ] % a 는 1 에서 5 까지 a = >> b=1:2:9 % b 는 1 에서 2 씩증가시켜 9 까지 b = >> c=[b a] c = >> d=[1 0 1 b(3) a(1:2:5)] d = 예제 1.2 ( 복소수사용하기 ) >> [1-2i i] % contains 5 elements i i >> [(1-2i) i] % contains 4 elements i i >> [1-2i i] % contains 4 elements i i 예제 1.3 ( 열벡터만들기 ) >> c=[1;2;3;4;5] c = 1 2

3 3 4 5 예제 1.4 ( 행렬만들기 ) >> g=[ ; ] % g는 2 by 4 matrix g = >> g=[g, g; 9:12, 9:12] % g는 6행 4열로확장됨 g = >> h=[1 2 3; ] % 모든 row 가같은수의 column 을가져야만한다.??? All rows in the bracketed expression must have the same number of column π π πππ

4 π π π π π π π π π π 예제 1.5 (Sampling) >> x=[0.1*pi.2*pi.3*pi.4*pi.5*pi.6*pi.7*pi.8*pi.9*pi pi] x = Column 1 through Column 8 through >> y=sin(x) y = Column 1 through Column 8 through >> plot(x,y) % y = sin(x) 의그래프 예제 1.6 (linspace 혹은 logspace 이용하기 ) >> x=(0:0.1:1)*pi % 0 에서 0.1*pi 씩증가하여 1*pi 까지 11 개의원소 x = Column 1 through

5 Column 8 through >> x=linspace(0,pi,11) % 0 에서 pi 까지 11 개의원소를표현 x = Column 1 through Column 8 through >> logspace(0,2,11) % 10^0 부터 10^2 까지 11 개의점이포함되도록만듬 Column 1 through Column 8 through

6 예제 2.1 ( 저장된변수명과변수의크기와길이 ) >> whos Name Size Elements Bytes Density Complex a 1 by Full No b 5 by Full No c 1 by Full No d 1 by Full Yes e 5 by Full Yes f 5 by Full Yes g 3 by Full No Grand total is 42 elements using 456 bytes leaving bytes of memory free. >> A=[ ; ] A = >> s=size(a) % size 함수는행의개수와열의개수를알려줌 s = 2 4 >> [r,c]=size(a) r = 2 c = 4 >> r=size(a,1), c=size(a,2) % 두개이상의명령어를쓰려면, 혹은 ; 를사용 r = 2 c = 4 >> length(a) % length 함수는행의개수와열의개수중더큰값을준다 4 >> size([]) 0 0

7 예제 2.2 ( 특수행렬만들기 ) >> zeros(2) % a 2-by-2 matrix of zeros >> ones(2,4)*pi % a 2-by-4 matrix of ones * pi >> A=[1 2 3; 4 5 6]; % 명령문마지막에 ; 를사용하면결과가출력안됨 >> ones(size(a)) >> eye(3) >> eye(3,2)

8 >> rand(3,1) % 0 과 1 사이의균일하게분포된임의의수들의 3 by 1 행렬 >> randn(2) % 평균 (mean) 이 0 이고표준편차 (variance) 가 1 인 % 정규분포를가진수들의 2 by 2 행렬 >> a=[1 2 3;4 5 6]; >> b=find(a>10) b = [ ]

9 예제 3.1 ( 열벡터만들기 ) >> a=1:5 a = >> b=a' b = >> c=b' c = 예제 3.2 ( 전치행렬만들기 ) >> c=a.' % 실수데이터에서는 와. 는서로동일하다. c =

10 >> d=a+i*a % 여기서 i=sqrt(-1) d = Columns 1 through i i i i Columns i >> e=d' % complex conjugate transpose of d e = i i i i i >> f=d.' % noncomplex conjugate transpose of d f = i i i i i 예제 3.3 ( 대칭행렬과반대칭행렬 ) >> a = rand(3), S = (a+a')/2, K = (a-a')/2, sum = S+K a = S = K =

11 sum = 예제 4.1 ( 뎃셈과상수연산 ) >> g, h % recall previous array g = h = >> g+ones(3,3)??? Error using ==> + Matrix dimensions must agree. >> g+h % add h to g on an element-by-element basis >> ans-h % subtract h from the previous answer to get g back

12 >> g-2 % g의모든원소에서 2를뺀경우 >> 2*g >> 2*g-h % multiplies g by 2 and subtracts h from the result 예제 4.2 ( 원소끼리의곱셈과나눗셈 ) >> g.*h >> g./h >> h.\g % 슬래쉬 (slash) 나백슬래쉬 (back slash) 모두나눗셈을정의할때, % 슬래쉬아래쪽이슬래쉬위쪽을나눈다.

13 >> g*h % 행렬의곱셈은매우특별한방식으로정의되어있다.??? Error using ==> * Inner matrix dimensions must agree. >> g/h Warning: Rank deficient, rank = 1 tol = e-15 >> h\g % 행렬의두나눗셈이서로다른결과를가지고있다. Warning: Rank deficient, rank = 1 tol = e 예제 4.3 ( 원소끼리의지수승 ) >> g,h g = h = >> g.^ >> g.^-1 % MATLAB doesn't like this syntax

14 ??? g.^- Missing variable or function >> g.^(-1) % However, with parentheses it works! >> 2.^g >> g.^h >> g.^(h-1)

15 예제 5.1 ( 행렬의곱와내적을이용한사이각구하기 ) >> A=[2 1 5; 1 3 2], B=[3 4; -1 2; 2 1] A = B = >> A*B

16 >> x = 1:4, y=2:2:8 x = y = >> x * y??? Error using ==> * Inner matrix dimensions must agree. >> x * y' 60 >> a=[-1 2 4]; b=[ ]; % angle between vectors a and b >> theta = acos( a*b' / sqrt(a*a'*b*b') ) 5.2 ( 행렬곱연산의교환법칙 ) >> C=[1 3; -1-1] C = >> C*ones(2,2) % 행렬의곱셈에대한교환법칙은성립하지않는다 >> ones(2,2)*c 0 2

17 0 2 >> AB=A*B, AB*C % 행렬의곱셈에대한결합법칙은성립한다. AB = >> BC=B*C, A*BC BC = 예제 5.3 ( 행렬의나눗셈 ) >> A=[2 1 5; 1 3 2], B=[3 4; -1 2; 2 1] % 정방행렬이아닌경우의나눗셈 A = B = >> (A*B)/B >> A\(A*B)

18 >> A=rand(3,3), B=rand(3,3) % 일반적정방행렬의나눗셈 A = B = >> (A*B)/B >> A\(A*B) 예제 6.1 (1 차원배열의주소할당 ) >> x=(0:0.1:1)*pi; % 0 에서 0.1*pi 씩증가하여 1*pi 까지 11 개의원소 >> y=sin(x); >> x(3) % The third element of x >> y(5) % The fifth element of y >> x(1:5) % 첫번째에서다섯번째원소까지나타냄

19 >> y(3:-1:1) >> x(2:2:7) % 두번째, 네번째, 여섯번째원소. % 즉, 두번째원소부터 7보다작은원소까지 2씩증가시킴 >> y([ ]) % [ ] 자체가 y의원소들의주소의배열이다 예제 6.2 (2차원배열의주소할당 ) >> A=[1 2 3; 4 5 6; 7 8 9] A = >> A(2,3)=0 % A(2,3) 의원소 [2nd row, 3rd column] 가 0으로바뀜 A = >> A(2,6)=1 % 원래 A가 6개의열을가지지않음으로 0을채워서늘려놓았다. A = >> A=[1 2 3; 4 5 6; 7 8 9]; % restore original data >> B=A(3:-1:1,1:3) % row의순서를거꾸로쓰는법 B = >> B=A(3:-1:1,:) % 마지막의싱글콜론 (single colon) 은모든열을가지는 % 것을의미한다. 즉, : 은 1:3 을짧게줄여쓴것이다.

20 B = >> C=[A B(:,[1 3])] C = >> B=A(1:2,2:3) B = >> C=[1 3] C = 1 3 >> B=A(C,C) % 여기서는배열 C 가 A 의인덱스 (index) 로사용되어진다. B = 예제 6.3 (2 차원배열의조작 ) >> A = [1 2 3; 4 5 6; 7 8 9] >> sum(a) % 합을구하는함수들은열단위로연산한다.

21 >> sum(sum(a)) 45 >> A(1:9) % A 를행벡터로바꿈 A = >> B=A(:) % A 의열들을하나의열로펼쳐놓음 B = >> B=B.' % dot-transpose operation 을이용하여 row 로바꿈 B = >> B=[1 2 3; 4 5 6; 7 8 9]; >> B(:,2)=[] % [](empty matrix) 는앞의 B 의모든행의두번째열을 % 삭제하라는의미이다 B = >> B=B.' B = >> B(2,:)=[] B = >> A(2,:)=B % 두번째열의모든행원소를 B로대체 A =

22 >> B=A(:, [ ]) % 모든행을두번째열에의해네번반복. % 이런것을 Tony's Trick 이라고함 B = 예제 6.4 ( 배열조작시행렬의크기 ) >> A=[1 2 3; 4 5 6; 7 8 9]; >> A(2,2)=[]??? In an assignment A(matrix,matrix) = B, the number of rows in B and the number of elements in the A row index matrix must be the same. >> B B = >> B=A(4,:) % A 가네번째 row 를가지고있지않기때문에 % MATLAB 은어떻게하라는건지알지못한다.??? Index exceeds matrix dimensions >> B(1:2,:)=A % 두행렬의크기 (size) 가다른경우 error 가생긴다??? In an assignment A(matrix,:) = B, the number of columns in A and B must be the same.

23 >> B(3:4,:)=A(2:3,:) % B의크기는 1 by 3이지만, B의 3과 4 행을 % 행렬 A의 2-3 행으로정의함으로서 B의크기가 4 by 3으로바뀐다. % 그러나두번째행에대해서는아무런조건이없으므로 0이대입된다. B = >> G(1:6)=A(:,2:3) % A 의 2 와 3 열을한행에여섯열로늘어뜨림 G = >> D=[ ; ; ], v=[2; 4; 8] D = v = >> E=[D(:,1)-v D(:,2)-v D(:,3)-v D(:,4)-v] E = >> E=D-[v v v v] % 바로위의방법보다더빠르게구할수있다. E = >> % 그러나 D의열의수를모르는경우 Tony's Trick 을사용하는것이좋다. >> c=size(d,2); % 여기서 c는 D의 column 수를나타낸다. >> E=D-v(:,ones(1,c)) % ones(1,c) 는원소가모두 1인 1-by-4 vector 를만든다. E =

24 예제 7.1 ( 논리배열 ) >> x=-3:3 % Create data x = >> abs(x)>1 % x 의절대값이 1 보다크면 1( 참 ) 을, 아니면 0( 거짓 ) 을돌려준다 >> y=x(abs(x)>1) % x 의절대값이 1 보다큰경우의 x 값을돌려준다 y =

25 예제 7.2 ( 배열인덱스 ) >> y=x([ ]) % [ ] 은 x 와같은크기가아니므로배열의 % 인덱스로서 0 을인식하지못한다??? Index into matrix is negative or zero. >> x(abs(x)>1)=[] % abs(x)>1 인경우의값들을제외한나머지 x 의원소만나타냄 x = 예제 7.3 ( 행렬의관계연산 ) >> B=[5-3;2-4] B = >> x=abs(b)>2 % 행렬에대한관계연산자의적용 x = >> y=b(abs(b)>2) y = 예제 7.4 (find 함수 ) >> x=-3:3 x = >> k=find(abs(x)>1) % find 함수를이용해서절대값이 1 보다큰경우의 % 배열 x 의인덱스를찾음 k = >> A=[1 2 3; 4 5 6; 7 8 9]

26 A = >> [i,j]=find(a>5) % 인덱스 i와 j에각각행과열의인덱스가저장된다 i = j = 예제 8.1 ( 성긴행렬 ) >> As=sparse(1:5,1:5,ones(1,5)) As = (1,1) 1 (2,2) 1 (3,3) 1

27 (4,4) 1 (5,5) 1 >> As=sparse(eye(5)) As = (1,1) 1 (2,2) 1 (3,3) 1 (4,4) 1 (5,5) 1 >> A=full(As) % 성긴행렬의전체적인형태를 full 함수를이용하여볼수있다. A = >> clear % 모든변수들을지움 >> B=eye(200) % FULL 200-by-200 identity matrix >> Bs=sparse(B) % Sparse 200-by-200 identity matrix >> whos % Bs는단지 0.5% 의 0 아닌원소만을포함하고 % 3200 bytes 만을저장하면된다. Name Size Elements Bytes Density Complex B 200 by Full No Bs 200 by No

29 예제 9.1 (Diary 저장과, Mat-File 저장, 그리고그림화일저장 ) >> diary sec19.txt % 아래에나타나는출력을 sec19.txt 에저장한다. >> x = -pi:0.01:pi; >> % My work will be saved on sec19.txt >> y = sin(x); >> save sec19 x y % 변수 x,y 의값을 sec19.mat 파일에저장한다. >> clear % 모든저장변수값을소거한다. >> load sec19 %sec19.mat 에저장된변수값을읽어온다. >> plot(x,y) >> print -dbmp sec19.bmp % 그림창에나타난그래프를 bmp 형식으로저장한다. >> diary off >> A = [2 3-5; 4-1 6]; B = [1-3 4; 3 0 1]; >> A+B, 2*A, 3*B, 3*A-2*B >> X= 3*B - 1/2 * ( 2*A + B) >> A = [ ; ; ]; >> size(a)

30 >> a13=a(1,3), a31=a(3,1), a25=a(2,5), a34=a(3,4) >> [i,j] = find(a==4) >> [i,j] = find(a==7) >> [i,j] = find(a==15) >> A = [2 1 0; 3 2 0; 1 0 1]; B=[ ; ; ]; >> A*B >> % B*A % it's not defined (!) >> A = [1 0 0; 0 1 0; 1 1 0]; B=[1 0 1; 0 1 1; 0 0 1]; >> AB=A*B, BA=B*A >> AB==BA % 정수값을가지는행렬이동치인가를비교하여얻은논리배열 >> all(ab==ba) % 각행의원소가모든같은값을가지면 1 을출력 >> all(all(ab==ba)) % 모든원소가같은값을가지면 1 을출력 >> A = [1 4-1; 3 0 2]; B=[3 0 1; 2 1 3; 4 1 2]; C=[2 1; 0 2; -3 1]; >> (A*B) * C == A * (B*C) >> all(all( (A*B)*C == A*(B*C) )) % integer matrix equality check

31 >> A=rand(2,2); B=rand(2,2); C=rand(2,2); >> A*(B+C) - (A*B+A*C) %#% A*(B+C) == A*B + A*C >> (A+B)*C - (A*C+B*C) %#% (A+B)*C == A*C + B*C >> A = rand(4); A = (A-A')/2; A2 = A^2; >> A2 - A2', max(max(abs(a2-a2'))) >> P = [ ; ; ]; >> di = sum(p(1,:)) >> rc = sum(p(:,2)) >> iu = P(3,4)

34 π π

예제 1.1 ( 관계연산자 ) >> A=1:9, B=9-A A = B = >> tf = A>4 % 4 보다큰 A 의원소들을찾을경우 tf = >> tf = (A==B) % A

예제 1.1 ( 관계연산자 ) >> A=1:9, B=9-A A = 1 2 3 4 5 6 7 8 9 B = 8 7 6 5 4 3 2 1 0 >> tf = A>4 % 4 보다큰 A 의원소들을찾을경우 tf = 0 0 0 0 1 1 1 1 1 >> tf = (A==B) % A 의원소와 B 의원소가똑같은경우를찾을때 tf = 0 0 0 0 0 0 0 0 0 >> tf