본문 바로가기

programming language484

MATLAB rot90(), triu(), transpose, hermitian transpose, upper(), lower() MATLAB 의 rot90() 함수는 matrix 를 시계 반대방향으로 돌리는 기능을 하는 함수입니다. x=magic(4) x = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 위 x 에 대하여 다음과 같이 하면 시계 방대방향으로 90 도 돌리게 됩니다. x1=rot90(x) x1 = 13 8 12 1 3 10 6 15 2 11 7 14 16 5 9 4 위 x 에 90 도씩 3번 시계 방대방향으로 돌리고 싶다면 다음과 같이 하면 됩니다. x2=rot90(x,3) x2 = 4 9 5 16 14 7 11 2 15 6 10 3 1 12 8 13 다음으로 설명한 명령어는 triu(), tril() 함수 입니다. triu() 함수는 upper triangular 부분을 구하는 함수이며 tr.. 2011. 4. 3.
MATLAB diff(), sortrows(), 오늘의 포스팅에서는 MATLAB 기본 함수 몇 가지에 대해 설명해 보려 합니다. diff(X, n) 함수는 인자간의 차이를 구하는 함수 입니다. n 값은 차이를 구하는 횟수를 의미 합니다. 아래 예를 보면 쉽게 이해하실 겁니다. diff(x, 2) 는 차이를 두 번 구하는 거죠. diff(x) 의 결과를 보면 [2 2 2 2] 니까 이 벡터에 대해 차이를 구하면 [0 0 0]이 나오게 되는 겁니다. sortrows(X, n) 은 X matrix 를 n column index 기준으로 정렬하는 함수 입니다. default 는 오름차순입니다. 다음 예를 보죠. 1열을 기준으로 오름차순 정렬하는 예입니다. 2열을 기준으로 오름차순 정렬하는 예입니다. 다음 예는 1열을 기준으로 오름차순 정렬한 후에 1열에 같은.. 2011. 4. 3.
MATLAB .* ./ .^ 연산 MATLAB 은 matrix 연산을 기본으로 한다. 즉 다음과 같이 A*B 의 경우 matrix 곱셈을 수행한다는 것이다. 그런데 저 곱하기에 아래와 같이 .(점)만 하나 더 찍어 주면 각 인자간의 곱이 된다. 이처럼 인자간의 연산을 할 때는 당연히 두 벡터 또는 매트릭스의 size 가 같아야 한다. /, ^ 등의 연산자도 마찬가지다. ./, .^ 연산을 하게 되면 각 인자간의 연산을 수행하게 된다. 2011. 4. 1.
MATLAB index 이번 포스팅에서는 MATLAB 인덱스에 대해 설명한다. 다음 matrix A 에 대한 각 인덱스의 접근 방법에 대해 알아보자. : (콜론) 의 사용법에 대해 알아보자. : 은 a:b 와 같은 형태로 쓰이면 a~ b 까지라는 뜻이다. 인덱스에 접근하는 용도로 : 을 쓸 때 a,b 값은 1보다 큰 양의 정수 여야 한다. : 을 이용하여 A 의 3행 1열부터 4열 까지의 값을 B 라는변수에 입력해 보자. 행, 열을 구분하지 않고 다음과 같이 접근할 수도 있다. MATLAB 은 기본적으로 행 벡터 즉 열을 기준으로 연산을 하는 특징이 있다. A 의 5행 모든 열을 D 라는 변수에 입력해 보자. 콜론을 이용하여 인덱스에 접근을 할 때 콜론만 단독으로 쓰일때는 '모든 것' 을 의미 한다. 다음으로 인덱스 접근시의 .. 2011. 4. 1.
MATLAB shortcut, 단축키 설정 오늘은 MATLAB 단축키 설정 방법에 대해 알아 본다. MATLAB desktop 에서 File 메뉴 à 그 다음에 Preferences 메뉴에 들어간다. 그럼 다음과 같은 창이 뜬다. 1번 부분을 선택하여 shortcut 설정 탭으로 들어간다. 2 번 부분에 설정하고자 하는 주제를 입력한다. 위 예에서는 cell 에 대한 단축키를 설정하기 위해서 cell 을 썼다. 3번 부분에 나온 단축키 설정 가능한 리스트들 중에서 설정 하고자 하는 동작을 선택한다. 기존에 단축키가 있었다면 4 번 부분에서 원하는 단축키를 설정한다. 기존에 단축키가 없는 동작에 대해 단축키를 설정하기 위해서는 5번 부분에 있는 + 를 누르고 4번 부분에서 단축키를 설정한다. 마지막으로 다 설정을 했으면 Apply à OK 를 누른다. 2011. 4. 1.
Powershell 특정 확장자 파일의 경로 다 찾기 컴퓨터를 이용하다 보면 특정 확장자의 파일 경로를 다 찾아야 할 경우가 있다. 아래 코드는 powershell 을 이용하여 D 드라이브 및 그 하위 드라이브에 있는 txt 파일들을 다 찾고 이름 및 수정 시간에 따라 정열 한 후에 d:\Temp_file_path.txt 에 그 결과를 저장하는 코드이다. $a = gci D:\* -Include *.txt –Recurse # D 드라이브 및 그 하위 폴더에 있는 txt 확장자의 파일들을 다 찾는다. $h= $a | sort-object -property directoryname -unique | sort-object -property lastwritetime # 이름에 따라 정렬하고 겹치는 이름들은 제외 한다. 그리고 수정 시간에 따라 정렬한다. $n=0 .. 2011. 4. 1.
MATLAB 최대공약수 gcd(), 최소공배수lcm() 오늘의 포스팅에서는 최대 공약수와 최소 공배수 관련 함수에 대해 살펴 본다. MATLAB 내장 함수로 gcd() 함수는 최대 공약수를 계산해 주는 함수이다. scalar 뿐만이 아니라 벡터에 대해서도 최대 공약수를 계산해 준다. 사용법은 다음 예와 같다. 다음으로 최소 공배수를 계산해주는 함수는 lcm() 이다. lcm() 함수 역시 scalar 뿐만이 아니라 벡터에 대해서도 최소 공배수를 계산해 준다. 2011. 3. 30.
MATLAB 소수 관련 함수 primes(), isprime() 이번 포스팅에서는 MATLAB 의 내장 함수 중 primes(), isprime()에 대해 알아보고, 소수인지 판별해주는 함수를 직접 작성해 보겠습니다. 내장 함수 중 소수를 생성해주는 함수로는 primes(n) 라는 함수가 있습니다. 이 함수는 n 보다 작은 소수를 생성해 줍니다. 아래 예는 primes() 함수를 사용하여 30보다 작은 소수들을 생성한 것입니다. 소수인지 판별하는 함수는 isprime(K) 라는 함수가 있습니다. 이 함수는 K 가 소수면 1 아니면 0을 반환하게 됩니다. 사용법은 다음과 같습니다. 그럼 이제 소수인지 판별하는 함수를 직접 코딩 해 보죠. 함수 이름은 My_isprime() 이라 해 보죠. 제가 짠 함수는 아래와 같습니다. 일단 벡터나 매트릭스 입력에 대해서도 동작하게 .. 2011. 3. 30.
계산기 항상 위에 윈도우 계산기 항상 위에~~~ 유용하게 쓰시길~ 위 파일은 윈도우 xp 에서 Autohotkey 로 만들었던 거였는데~ 윈도우 7에서는 현재 안되는군요. 그래서 윈도우 7 에서 다시 만들었습니다. 윈도우 7에서 사용하실 분들은 아래 파일 다운로드 받으세요. 아래 요청에 따라 윈도우 10 에서 만들어 봤습니다. 한글 윈도우 10 에서 정상동작 확인 했습니다. 파일 첨부합니다 2011. 3. 28.
MATLAB diag(), eye(), zeros(), ones() 이번 글에서는 MATLAB 의 기본적인 몇몇 함수들에 대해 알아 보겠습니다. diag() 함수를 이용하여 diagonal matrix 를 만들거나 matrix 의 diagonal 값을 얻을 수 있습니다. eye() 함수를 이용하여 identity matrix 를 만들 수 있습니다. zeros() 함수를 이용하여 zero matrix 를 만들 수 있습니다. ones() 함수를 이용하여 1로 구성된 matrix 를 만들 수 있습니다. blkdiag() 함수를 이용하여 diagonal 값이 matrix 로 구성된 matrix 를 만들 수 있습니다. 2011. 3. 27.
MATLAB persistent MATLAB 함수 내부의 변수는 함수 사용 후에는 메모리 상에서 지워지게 됩니다. 하지만 경우에 따라 함수 내부에서 변수가 계속 유지 되게 할 필요가 있는 경우가 있습니다. 이러한 경우에 persistent keyword 를 사용합니다. 다음 예를 통해 persistent 에 대해 알아보죠. 아래 그림의 y=persistent_func(a,b) 함수는 a, b 값 입력을 받아서 a+b+c 값을 출력하는 함수 입니다. c 값에는 함수 내부에서 메모리상에 계속 남아 있도록 하기 위하여 persistent 로 설정했습니다. 그리고 초기값 설정을 위하여 isempty() 함수를 써서 c 값이 비어 있으면 초기값을 설정하도록 하였습니다. 위 함수에 대하여 다음과 같이 사용해 보도록 하겠습니다. 함수 내부에서 c .. 2011. 3. 27.
MATLAB contour 등고선 그리기 MATLAB contour() 함수를 이용하면 등고선 그래프를 매우 쉽게 그릴 수 있다. 다음과 수식에 따른 값에 대해 등고선을 그려 보자. [X,Y] = meshgrid(-2:2); % X,Y 값 설정 Z = 10*(X/5 - X.^3 - Y.^5).*exp(-X.^2-Y.^2); % Z 값 계산 위 식에 대해 아래 명령으로 등고선을 간단하게 그릴 수 있다. figure, contour(X,Y,Z), colorbar 등고선 사이의 표면에 색을 입히고 싶다면 contourf() 함수를 이용한다. figure, contourf(X,Y,Z), colorbar 다음과 같은 등고선 그래프를 얻을 수 있다. 위 그래프를 보면 간격이 너무 넓어서 등고선이 너무 각진 형태로 나오는 것을 확인 할 수 있다. 이럴 때는.. 2011. 3. 27.
MATLAB, Newton method 를 이용하여 해를 구해 보자 Newton method 는 미분 가능한 연속 함수 f(x) 에 대하여 다음과 같은 방식으로 f(x) =0 의 해를 구하는 방법이다. 위 식에 대하여 Newton method 를 이용하여 f(x) =0 이 되는 실수 해 x 를 구해 보자. 위 식에 대한 미분 결과는 다음과 같다. Newton method 의 최대 반복 횟수는 100 회로 하고 상대 오차는 10-9 으로 설정 한 후 이를 만족하는 해를 찾아 본다. MATLAB code 는 다음과 같다. command 창에서 결과를 확인 하면 Newton method 를 이용한 결과와 roots() 함수를 이용한 결과가 근사적으로 일치 함을 확인 할 수 있다. 2011. 3. 27.
MATLAB QPSK BER simulation in AWGN channel 이번 포스팅에서는 AWGN (Additive white gaussian noise) 채널 환경에서의 QPSK Bits Error Rate(BER)에 성능 검증 simulation 에 대해 설명한다. QPSK 및 잡음의 분산(σ2) 과 No(noise power spectral density) 의 관계에 대한 내용은 다음 포스팅을 참조 하기 바란다. 2011/03/19 - [MATLAB] - MATLAB BPSK BER simulation in AWGN channel 2011/03/08 - [통신] - [디지털 통신] BPSK, QPSK QPSK 에서는 아래 식의 M =4 이다. 따라서 다음 식과 같이 정리 된다. QPSK 시뮬레이션 시의 주의할 점은 신호의 파워를 1 로 정규화 해 줘야 하며 잡음 역시 .. 2011. 3. 27.
MATLAB break, continue 이번 포스팅에서는 루프 문(for, while) 사용시의 제어를 위한 keyword 인 break 와 continue 에 대해 설명 드립니다. 보통 if 문 안에 써서 특정 조건 일 때 루프를 멈추게 하거나 또는 건너뛰게 할 때 이용합니다. break 는 루프 문의 동작을 멈추는 기능을 합니다. 다음 코드를 실행해 보면 for k=1:5 if k==3 break end fprintf('%d\n',k) end 1, 2 만 출력이 되고 k = 3 일 때 break 명령에 의해 for 루프를 빠져 나오게 됩니다. 루프 문 안에 또 루프 문이 있는 nested 루프의 형태에서 break 를 이용할 경우에는 break 를 둘러 싸고 있는 루프문 하나만 빠져 나오게 됩니다. continue 문은 특정 루프를 건너 .. 2011. 3. 26.
MATLAB piecewise 수식 표현 MATLAB 을 이용하여 그림과 같은 piecewise 식을 나타내 보자. 이러한 piecewise 식은 비교 연산자와 논리 연산자를 이용하면 매우 간단하게 표현이 가능하다. 위 그래프에 대한 식은 다음과 같다. 위 식을 y = my_piecewise(t) 라는 함수로 표현해 보자. 위 그림의 오른쪽의 빨간색 네모 칸 친 부분 처럼 각각의 조건을 비교 연산( 2011. 3. 26.
MATLAB global keyword MATLAB 의 함수 내부의 변수 값들은 함수가 사용된 이후에는 자동으로 메모리 상에서 없어진다. 하지만 같은 변수를 다른 함수 또는 파일에서도 이용하려는 경우 함수 내부에서 변화된 값이 계속 남아 있게 해야 하는 경우가 있다. 이러한 경우 global keyword 를 이용한다. global 은 하나의 변수를 여러 파일에서 이용하기 위해서 사용한다. global keyword 의 이용 방법은 변수 사용 전에 선언을 해 주고 사용해야 한다. 다음과 같이 예를 들어 보자. 위 예에서 global_test.m 파일에서 변수 a 값을 my_func() 함수에서도 이용하고 함수 내부에서 a 값을 변경했다. 그리고 global_test.m 파일에서 my_func() 사용 전과 사용 후의 a 값이 다음과 같이 달라.. 2011. 3. 26.
MATLAB Tex character 수학 기호 넣기 MATLAB 이용시 그래프 또는 command 창에 text 출력시 수학 기호등을 넣고 싶을 때가 있다. 다음 예는 sin 그래프를 그리고 x 축과 y 축에 θ 라는 특수 기호를 넣었다. MATLAB 에서 이런 특수 기호를 넣기 위해서는 위 코드에서 보는 것처럼 \character 형태로 입력해 줘야 한다. 아래 표를 참조하면 원하는 기호를 넣을 수 있을 것이다. Reference : MATLAB HELP 2011. 3. 25.
MATLAB Rayleigh random variable 생성 Gaussian r.v.(random variable)을 이용하여 rayleigh r.v. 을 생성하고 그 PDF 및 CDF 를 구해 보자. rayleigh random variable 에 대해서는 다음 포스팅을 참조 2009/12/26 - [수학] - Rayleigh PDF PDF 및 CDF 를 구하는 법에 대해서는 다음 포스팅을 참조 하시길 2011/03/23 - [MATLAB] - MATLAB PDF, CDF, CCDF 서로 독립인 평균이 0 이고 분산(σ2)이 1인 Gaussian r.v. X, Y 를 생성해 보자. rayleigh r.v. Z 는 다음과 같이 발생 시킬 수 있다. 이제 이전 포스팅 에서 설명한 대로 다음과 같이 PDF 를 구해 보자. CDF 는 다음과 같다. 분산이 σ2 인 가우.. 2011. 3. 23.
MATLAB FAST convolution using FFT 이번 포스팅에서는 FFT 를 이용한 고속 convolution 에 대해 알아본다. 일반적으로 conv(x, y) 은 filter() 함수로 구현되는 선형 convolution 이다. 선형 convolution 은 x 또는 y 의 길이가 증가할수록 그 복잡도는 급격히 증가하는 특징이 있다. 이러한 선형 convolution 은 순환(Circular) convolution 을 이용하여 구현이 가능하며, 순환 convolution 은 FFT(Fast Fourier Transform) 와 IFFT(Inverse Fast Fourier Transform)를 이용하여 구현이 가능하다. 일단 선형 convolution 을 순환 convolution 으로 변환하는 과정에 대해 살펴 보자. x=[1 2 3 4] y=[1.. 2011. 3. 23.
MATLAB random number 생성 이번 포스팅에서는 MATLAB 에서 기본으로 제공하는 각종 random number 의 생성 방법에 대해 알아본다. MATLAB 에서도 rand() 함수는 uniform random number 를 발생 시킨다. rand() 함수는 0에서 1 사이 값을 발생 시킨다. 이를 이용하여 다양한 범위의 uniform random number 를 발생 시킬 수 있다. a ~ b 사이의 n행 k열 uniform random number 를 발생시키는 방법은 다음과 같다. Uniform_rand = a + (b-a).*rand(n, k); 예를 들어 -10 ~ 5 사이의 4행 5 열 의 uniform random number 를 발생 시키려면 다음과 같이 하면 된다. 다음으로 정규 분포라 불리는 Gaussian ra.. 2011. 3. 23.
MATLAB PDF, CDF, CCDF 이번 포스팅에서는 각종 랜덤 값들을 이용한 PDF(Probability Density Function), CDF(Cumulative Distribution Function), CCDF(Complementary Cumulative Distribution Function) 를 그리는 방법에 대해 살펴 보려 합니다. PDF 를 나타내는 방법에 대해 모르시는 분들은 다음 포스팅을 참조하시길 ~ 2010/01/03 - [programming language/MATLAB] - MATLAB PDF (Probability Density Function) 평균이 2 이고 분산이 10 인 가우시안 잡음을 생성해 보겠습니다. 평균이 2 이므로 +2 를 하고 분산 값이 10 이므로 sqrt(10) 을 곱해 주게 됩니다. 이제.. 2011. 3. 23.
MATLAB 함수 만들기, inline, function handle 어느 컴퓨터 언어나 반복되는 부분은 함수화 하기 마련 입니다. MATLAB 도 예외는 아니죠. 이번 포스팅에서는 MATLAB 에서 함수 작성 방법에 대해 알아보도록 하겠습니다. MATLAB 에서는 함수를 만들 수 있는 방법이 3 가지 정도 됩니다. 첫 번째로는 함수 파일을 만드는 방법이고, 두 번째로는 inline() 객체를 이용하여 함수화 하는 방법, 세 번째로는 function handle 을 이용하는 방법 입니다. 함수 파일을 만드는 방법부터 설명 드리도록 하겠습니다. 예를 들어 설명해 보죠. a, b 라는 두 값을 입력 받아서 그 합과 차를 출력하는 함수를 만들어 보도록 하죠. 함수는 function 이라는 키워드를 먼저 써줘야 합니다. 그리고 다음 그림과 같이 구성을 합니다. 함수 이름은 fun.. 2011. 3. 21.
MATLAB cross product, 삼각형의 면적, 평행 육면체의 부피 MATLAB 의 cross product함수로 cross() 라는 함수를 제공합니다. 본 함수는 cross(A,B) 형태로 이용을 하며 A,B 는 인자가 3개 여야 합니다. 그렇지 않으면 에러가 납니다. 보통 cross product(외적) 는 x 기호로 나타내고 dot product (내적) 은 점(∙) 기호로 나타냅니다. 일단 그 수학적인 정의에 대해 살펴 보겠습니다. 벡터 A= 이고 벡터 B = 인 경우 cross product 는 A x B 로 표현이 되고 dot product 는 A ∙ B 로 표현이 되며 그 결과는 다음과 같습니다. A x B= = A ∙ B= cross product 는 몇 가지 특징이 있습니다. A x B 는 A, B 에 직교 합니다. 즉 A,B 에 의해 이루어 지는 면적에.. 2011. 3. 21.
MATLAB input() 이번 포스팅에서는 input() 함수에 대해 설명해 보겠습니다. input() 함수는 말 그대로 사용자로부터 입력을 받는 함수 입니다. 사용 방법은 매우 간단합니다. >> A=input('입력해 주세요: ') ; 입력해 주세요: 12 >> A A = 12 다음과 같이 벡터도 입력 가능 합니다. >> A=input('입력해 주세요: ') ; 입력해 주세요: [1 2 3] >> A A = 1 2 3 그럼 다음과 같이 좀 엄한 값을 한번 넣어 보죠. >> A=input('입력해 주세요: ') ; 입력해 주세요: T ??? Error using ==> input Undefined function or variable 'T'. 입력해 주세요: 그럼 위 결과처럼 라는 부분이 다시 나오게 됩니다. 여기서 T 라는 값.. 2011. 3. 21.
MATLAB 기본3, sort, if ~ end, for ~ end, max, min MATLAB 은 공학적으로 굉장히 편리 한 툴입니다. 하지만 편리함 속에는 단점도 있는 게 사실 입니다. 쉽다고 함부로 쓰다가는 정말 예상치 못한 결과에 당황하는 경우가 있곤 합니다. 따라서 MATLAB 함수 등을 쓰실 때는 help 문 이나 예제 등을 통해 정확히 알고 난 이후에 이용하는 것이 좋습니다. 이번 포스팅에서는 그럼 기본 함수 몇 가지를 비롯한 몇몇 주의 할 점 들에 대해 설명 드리겠습니다. 첫 번째로 소개해 드릴 함 수는 sort() 함수 입니다. 함수 이름만 봐도 이게 정렬을 해주는 함수구나 하는 건 아실 겁니다. 예를 통해 설명해 보도록 하죠~ A=magic(3) A = 8 1 6 3 5 7 4 9 2 A 에 대해 sort() 함수를 적용해 보면 다음과 같습니다. A1=sort(A) A.. 2011. 3. 21.
MATLAB erf, erfc, Q- function 통신에 대해 공부하다 보면 erfc(x) 함수나 Q(x), erf(x) 함수에 대한 내용들이 많이 나온다. 각 식들은 다음과 같다. 위 함수들은 각각 MATLAB 에서 erf(), erfc(), qfunc() 라는 명령어로 제공한다. erf(), erfc() 함수는 MATLAB 에서 기본적으로 제공하는 함수들이며 qfunc() 은 communication toolbox 에서 제공해주는 함수이다. 추가적으로 위 함수들과 관련하여 erfcinv() 라는 erfc 의 inverse 를 구해주는 함수와 qfuncinv 라는 qfunc() 의 inverse 를 제공해주는 함수도 있다. 2011. 3. 19.
MATLAB BPSK BER simulation in AWGN channel 오늘은 AWGN (Additive white gaussian noise) 채널 환경에서의 BPSK Bits Error Rate(BER)에 대해 설명한다. 일단 BPSK 가 뭔지 잘 모르는 사람은 다음 페이지를 참조 하기 바란다. 2011/03/08 - [통신] - [디지털 통신] BPSK, QPSK 위 페이지를 참조하면 bit 를 BPSK 심볼로 mapping 하는 방식에 대해 알았을 테고, 이제 잡음을 생성하는 방법을 알아 보자. 잡음을 생성하기 위해서는 몇 지 통신 관련 수식에 대해 알아야 한다. 보통 BER 그래프를 그릴 때는 SNR 축에 따른 그래프가 아닌 Eb/No 에 따른 그래프를 그리곤 한다. 따라서 BER 시뮬레이션을 정확히 하려면 Eb/No 에 따른 잡음을 생성할 수 있어야 하겠다. 잡음.. 2011. 3. 19.
MATLAB 반올림 관련 함수들 round, fix, ceil, floor 이번 포스팅에서는 MATLAB 의 rounding 관련 함수들에 대해 정리한다. 다음과 같은 벡터 값에 대해 위 함수들을 적용해 보고 그 특징을 알아본다. x=[-4.3 -1.8 0.7 3.6 1.2+2.6i ] x = -4.3000 -1.8000 0.7000 3.6000 1.2000 + 2.6000i ceil() : + 무한대 쪽으로 rounding 한다. ceil(x) ans = -4.0000 -1.0000 1.0000 4.0000 2.0000 + 3.0000i floor(): ceil 과는 반대로 - 무한대 쪽으로 rounding 한다. floor(x) ans = -5.0000 -2.0000 0 3.0000 1.0000 + 2.0000i fix(): 0 방향으로 rounding 한다. fix(x)a.. 2011. 3. 19.
MATLAB interpolation, interp1() 실험 데이터가 너무 적을 경우 그 경향성을 알아보기 위하여 실험 데이터의 중간 값들을 알아내야 하는 경우가 있다. 이런 경우 interpolation 을 수행하여 수학적으로 중간의 값들을 추정 하곤 한다. 이러한 수치적인 interpolation 을 해 주는 함수 중 하나가 interp1 함수 이다. 다음과 같은 데이터를 가정해 보자. X=0:10 X = 0 1 2 3 4 5 6 7 8 9 10 Y=sin(X) Y = 0 0.8415 0.9093 0.1411 -0.7568 -0.9589 -0.2794 0.6570 0.9894 0.4121 -0.5440 위 코드의 Y 값은 간단하게 0, 1, 2, … 10 값에 대한 sin() 값을 나타낸다. 위 값들을 이용하여 X 값이 0, 0.01, 0.02, 0.03.. 2011. 3. 19.