Introduction to SCI-Lab

Introduction to Sci-Lab

• 행렬의 이름은 보통 A, B 이렇게 대문자로 변수를 지정한다.
• 각 행은 띄어쓰기
• 각 열은 세미콜론(;)으로 구분한다.

A = [];
A = [1 2 3; -1 2 4; 5 6 1] 

# A  = 
#
#   1.   2.   3.
#  -1.   2.   4.
#   5.   6.   1.

E = [3 4 -1; 2 7 0; 1 -2 -3]

# E  = 
#
#   3.   4.  -1.
#   2.   7.   0.
#   1.  -2.  -3.

A + E

# ans  =
#
#   4.   6.   2.
#   1.   9.   4.
#   6.   4.  -2.

E - A

# ans  =
#
#   2.   2.  -4.
#   3.   5.  -4.
#  -4.  -8.  -4.

• determinant(행렬식)은 det(A)로 코드를 짤 수 있다.

det(A)
# ans  =
#
#  -28.

• 역행렬을 구하는 법은 inv(A)이런식으로 구할 수 있겠다.

inv(E)
# ans  =
#
#  0.75       -0.5        -0.25     
# -0.2142857   0.2857143   0.0714286
#  0.3928571  -0.3571429  -0.4642857

• Identity Matrix도 정의할 수 있음.

eye(2,2)
# ans  =
#
#  1.   0.
#  0.   1.

• 영행렬도 나타낼 수 있음

zeros(3,3)
# ans  =
#
#  0.   0.   0.
#  0.   0.   0.
#  0.   0.   0.

• 파이(%pi), 복소수(%i), 복소수 모듈러스(abs(z))도 가능하다

z = 3+4*%i
# z  = 
#
#  3. + 4.i

abs(z)
# ans  =
#
#  5.

3 + imult(4)
# ans  =
#
#  3. + 4.i

방정식의 해를 구해보자

• x^2 -5*x + 6 의 해를 구해보면

F = [1;-5; 6]
# F  = 
#
#  1.
# -5.
#  6.

roots(F)
# ans  =
#
#  3. + 0.i
#  2. + 0.i

그래프 그리기

• scf(0) : 그래프창 자체를 띄워준다

---

Fibonacci 수열

Fn+2 = Fn + Fn+1 이런 꼴의 수열을 피보나치 수열이라고 한다.

function y = Fibo(n)
  >     if n == 1 then y =1;
  >     elseif n == 2 then y = 1;
  >     else
  >         y1= 1; // first term
  >         y = 1; // second term
  >         for i = 1:n-2;
  >             y2 = y; // Earlier term of sequence
  >             y = y+y1; // New term of sequence
  >             y1 = y2; // store earlier term in memory
  >         end
  >     end
  > endfunction

Rational Function, Series Function

function y = rat(x)
    y = (x^3-2*x^2+5*x)/(x^4+x^3-2*x-5)
    
endfunction

scf(0)
plot(0:0.01:5,rat)

function y = series(n)
    y = 0;//initialization
    for i =1:n
        y = y+i^4
        
    end
endfunction