레이블이 CAS인 게시물을 표시합니다. 모든 게시물 표시
레이블이 CAS인 게시물을 표시합니다. 모든 게시물 표시

2015년 7월 29일 수요일

맥시마(maxima)에서 wxMaxima의 대화창으로 그래프 그리는 법

 wxMaxima에서는 그래프를 그리는 명령으로 wxplot2d 함수와 wxplot3d 함수가 있는데 각각 2차원 그래프와 3차원 그래프를 그려준다. 이번 절에서는 2차원 그래프를 그리는 함수의 사용법에 대해서 알아보자. 일단 메뉴창에서 이 함수를 선택하면 (Plot > Plot 2d...) 아래와 같은 대화창이 화면에 나타난다.
[그림 1] wxMaxima에서 plot>Plot 2d 메뉴 선택

[그림 2] Plot 2d 메뉴의 정보 입력 창
위에서 부터 차례대로

①도시할 함수(expression)
② x축의 범위
③ y축의 범위
④ 그래프의 정밀도
⑤ 표시형식(format)
⑥ 선택사항
⑦ 그래프를 저장할 화일명 (eps)

들을 입력하게 되어 있다. 일단 처음이니까 [그림 2]와 같이 첫 번째 칸만 채우고 [OK]를 클릭(엔터키를 쳐도 된다)하면 다음과 같은 그래프가 입력창에 나타나게 된다. 결과를 보면 알겠지만 대화창은 단순히 wxplot2d 함수의 파라메터를 결정해서 함수를 호출하는 역할을 한다. 따라서 사용자가 직접 이 함수를 입력해도 동일한 결과가 나타나므로 앞으로는 wxplot2d함수를 직접 입력해 보겠다.


맥시마(maxima)에서 수학 함수 정의

 맥시마에서 수학 함수는 := 기호를 써서 정의한다. 변수에 어떤 값을 대입하는 연산자는 콜론(:)이고 함수를 정의하는 기호는 (:=) 이라는 것에 유의해야 한다.


다음 예제를 고려하자.

.

이것을 맥시마로 풀면 다음과 같다.


 조건명령 if를 이용하여 함수를 정의할 때 구간 별로 다른 규칙을 정의할 수 있다. 예를 들어 다음과 같은 불연속 함수를 고려하자.
이 함수는 다음과 같이 입력할 수 있다.


다음 함수도 유사하게 입력할 수 있다.

이와 같이 맥시마에서 함수를 정의할 경우에는 := 연산자를 사용하며 if 조건문을 이용하여 불연속 함수도 자유롭게 생성할 수 있다.




맥시마(maxima)의 if 조건문

 맥시마에서 조건 검사를 하는데 if 명령어가 사용된다. 문법은 다음과 같다.

if <조건식> then <expr1> else <expr2>

만약 조건식이 참이면 <expr1>이 실행되고 거짓이면 <expr2>가 실행된다. 조건식을 구성하는데 사용되는 조건 연산자는 다음과 같다.

[표 1] 조건식에 사용되는 연산자들
연산자
의미
비고
a = b, equal(a,b)
같다
관계 연산자
a # b, notequal(a,b)
다르다(≠)
a < b
크다
a > b
작다
a <= b
크거나 같다
a >= b
작거나 같다
and
and
논리 연산자
or
or
not
not

예를 들면 다음과 같다.
변수 a에 %pi 값이 들어있으니 조건은 거짓이 되므로 b:0 이 실행된다. 위의 예제는 다음과 같이 조금 더 간단하게 바꿀 수 있다.

위에서 (if a<%e else 0)가 단일 표현식이며 그 값은 조건에 따라서 1이거나 0이 된다. 따라서 조건에 따라 그 값을 변수 b에 저장시키는 것이다. and연산자를 이용한 예를 들어보면 다음과 같다.

또한 elseif 명령을 이용하여 if문을 중첩하여 쓸 수 있다.

if <조건식1> then <expr1>
elseif <조건식2> then <expr2>
elseif <조건식3> then <expr3>
...
else <exprN>

이 문법에서 만약 <조건식1> 이 참이라면 <expr1>을 수행하고 빠져나간다. <조건식1>이 거짓이고 <조건식2>가 참이라면 <expr2>를 수행한다. <exprN>은 모든 조건이 거짓일 때 수행된다. 다음 예를 보자.
이 예제에서 변수 b에는 2가 대입된다.


2015년 7월 28일 화요일

맥시마(maxima)의 for 반복문

 맥시마는 반복 수행을 위한 명령어 for 가 있다. 문법은 다음 표에 정리한 바와 같다.

[표 1] for반복문의 문법
생략 가능
(기본값은 1)
다음 중 택일
(반드시 하나는 있어야 함>
for <변수>:<초기값>
step <증감값>
thru <끝값>
do <수행문>
while <조건식>
unless <조건식>
for <변수> in <list>
do <수행문>

먼저 반복문에서 사용할 변수를 초기화 시켜야 한다. 그 다음에 증감값을 step명령 다음에 지정해 준다. 이 증감값은 생략가능하며 이 경우 기본값인 1로 지정된다.

가장 기본적인 사용법의 예는 다음과 같다. (display함수는 단순히 변수의 값을 화면에 출력하는 함수임)
만약 1부터 100까지의 합을 구하고 싶다면 다음과 같다.
증감값을 정하고 싶다면 step 명령어를 다음과 같이 이용한다.
for문을 이용하여 x5+x4+x3+x2+x+1을 생성하면 다음과 같다.
같은 명령을 thru 대신 while을 이용하면 다음과 같다.

또 다른 예로서 리스트 [0, 1, 2, 3, … , 29, 30]를 생성하려면 다음과 같이 하면 된다. 다음 명령을 입력한 후에 그 결과를 확인해 보자. (/*~*/사이의 주석문은 입력할 필요가 없다.)

똑같은 명령을 unless를 이용하여 작성하면 다음과 같다.
이러한 for문을 이용하여 다양한 반복 작업과 많은 요소를 갖는 객체(리스트나 다항식 등)를 생성할 수 있다.


맥시마(maxima)의 괄호를 이용한 표현식 나열

 이전 포스트들에서 설명한 계산식, 변수 대입, 수식 등등이 모두 Maxima의 표현식(expression)이다. 표현식은 그 자체로서 객체를 이루며 또는 정해진 동작을 수행하여 결과 값을 반환해 주기도 한다.
위의 두 예는 동일한 동작을 수행한다. 먼저 x에 2+i 라는 복소수를 대입한 후 그 값을 다시 a변수에 저장한다. 맥시마에서는 표현식을 괄호 ( )로 묶을 수 있으며 단일 표현식에서는 보통 괄호를 생략한다. 따라서 (x:2+%i)는 하나의 표현식이며 그 값은 2+i 라는 값을 가진다. 표현식이 값을 가지는 것이다. 따라서 a변수에 (x:2+%i)의 값인 2+i 가 대입되게 되는 것이다.
표현식이 하나라면 보통 괄호를 생략하지만 여러 개의 표현식을 묶을 때는 괄호를 사용하고 각각의 표현식은 콤마(,)로 구분한다.
위의 예를 보면 첫 번째 줄에서 네 개의 표현식이 괄호로 묶여있다. 이 명령을 내리면 각각의 표현식이 하나씩 실행이 된다. 즉 a변수에 cos2(x)+sin(x)수식이 입력되고 b변수에는 3/sqrt(2)값이, 변수 c에는 3값이 저장된다. 그런데 결과를 보면 (%o1) 3이 나온다. 즉 (...) 안의 표현식은 몇 개가 되었든지 하나의 값을 가지는데 제일 마지막 표현식의 값을 가지며 그 결과가 찍혔다는 것을 알 수 있다. 따라서 다음과 같은 대입이 가능하다.
위의 예를 보면 x:( ... ) 명령어에서 ( ) 안에 두개의 표현식이 있다. 각각은 실행되며 전체 표현식의 값은 마지막 표현식의 값이 되므로 그 값이 변수 x에 저장된다. 따라서 이 예제를 실행시키면 변수 a, b, x 세 개에 각각 값이 저장된다.
 이제 아래 예들의 동작을 구분해 보자.


첫 번째 예는 세 개의 표현식 각각을 독립적으로 실행시킨 것이다. 따라서 그 결과 값도 각각 독립적으로 저장된다. (%o1, %o2, %o3) 두 번째 예는 첫 번째 것과 동작은 똑같지만 결과값은 하나만 생성된다.(%o4) 괄호로 묶여서 전체가 하나의 표현식이고 그 값은 마지막 표현식의 그것이기 때문이다.


맥시마(maxima)의 리스트(list) 자료형

 리스트(list)는 맥시마의 기본 자료형 중 하나로서 다양한 표현식을 하나의 집합으로 묶어놓은 것이다. 기본적인 정의는 대괄호 [ ] 를 이용하며 각각의 요소는 콤마(,)로 구분한다. (공백 문자가 아님에 주의하자.)
위의 예는 lsA라는 변수에 리스트 [1, 2, 7, x+y]를 대입한 것이다. 리스트명은 임의의 식별자가 될 수 있다. 만약 x나 y가 어떤 값을 갖는 변수라면 그 값이 대신 들어간다. 그렇지 않다면 그냥 문자 상수로 인식된다.

 리스트 내부의 하나의 요소를 참조할 때에는 리스트명 뒤에 역시 대괄호를 이용하여 인덱스를 지정한다.


리스트의 요소로는 표현식(상수, 변수, 수학 함수 등등)과 리스트 자체도 리스트의 요소가 될 수 있다. 다음 예를 보자.
lsA는 네 개의 요소를 갖는 리스트인데 네 번째 요소는 리스트이다. 따라서 lsA[4]는 리스트 [a,b,c]인 것이다. 리스트 안의 리스트를 인덱싱하기 위해서는 위의 마지막 예와 같이 하면 된다.

 만약 값을 가진 변수를 이용하여 리스트를 정의한 후에 사용된 변수의 값을 바꾸면 어떻게 될까? 변경된 값이 리스트에 영향을 미칠까? 다음 예를 살펴보자.


만약 위와 같이 변수 a, b값을 이용하여 lsA를 정의했다면 정의된 그 시점의 a, b값이 lsA에 사용되는 것이지 변수 a, b 자체를 사용한 것은 아니다. 따라서 이후에 변수 a, b값을 바꾸는 것은 lsA에 아무런 영향을 미치지 않는다는 사실에 유의해야 한다.


그리고 리스트와 스칼라와의 산술 연산은 리스트 내의 모든 요소에 적용된다.


그리고 리스트끼리의 산술연산은 크기가 똑같아야 성립된다.

 리스트를 다루는 여러 함수 중에 많이 사용되는 몇 개만 소개하면 다음과 같다. 먼저 append() 함수는 두개의 리스트를 입력받아서 그 두개를 붙인 리스트를 생성한다.
리스트의 크기를 구하는 함수로서 length() 함수가 있다.

그리고 첫 번째 요소와 마지막 요소를 반환하는 함수는 각각 first(), last() 함수이다.


[표 1] 리스트 관련 함수 정리
기호
의미
append(lsA, lsB)
lsA와 lsB를 붙여서 새로운 리스트 생성
length(lsA)
lsA의 요소의 개수를 반환
first(lsA)
lsA의 첫 번째 요소
last(lsA)
lsA의 마지막 요소