이전 포스트들에서 설명한 plot(a,b) 명령어에서는 a벡터를 x축에, b벡터를 y축에 선형 스케일로 2차원 그래프를 그린다. 수학이나 공학에서는 넓은 영역의 값을 한 그래프에 나타내기 위해 때때로 로그 단위의 그래프를 그릴 필요가 있다. 특히 주파수를 기본 축으로 하는 함수에서는 그 값이 일반적으로 log 단위로 표시되는데 이때 사용되는 명령어로는 아래와 같이 semilogx(), semilogy(), loglog() 등이 있다. 그리고 극좌표 스케일의 그래프를 그릴 때는 polar() 함수를 사용한다.

[표 1] 로그 스케일 그래프를 생성하는 함수

loglog(a,b)

x축과 y축 모두 로그 스케일인 그래프 생성

semilogx(a,b)

x축만 로그 스케일인 그래프 생성

semilogy(a,b)

y축만 로그 스케일인 그래프 생성

polar(a,b)

극좌표 그래프를 그린다

 예를 들어 다음의 식을 고려해보자.

먼저 x축의 값인 주파수를 로그단위로 나타내 보자. 함수 logspace()를 사용하면 일정한 상용 로그값 차이를 갖는 벡터를 만들어 준다.


>> w = logspace(-1,2,100)
>> H = 1 ./(1+w.^2).^0.5
>> semilogx(w,H)
>> grid on


벡터 w에는 10-1(rad/sec)에서 102 (rad/sec)까지의 값을 100등분하여 로그값을 저장한다. 즉, logspace() 함수는 다음과 완전히 동일하다.


>> w = 10.^linspace(-1,2,100)


x축을 log 스케일로 그리는 그래프는 semilogx() 함수를 이용하면 된다.

[그림 1] 실행 결과 (gnu octave 실행 결과를 캡쳐)


 polar()함수는 극좌표계의 그래프를 그려준다. 입력 인수로 거리와 각도를 각각 받는다. 즉 polar(r, a) 에서 r은 거리의 벡터, a는 각도(단위는 라디안)의 벡터로 주어져야 한다.


>> polar(0:10, 0:10)


[그림 2] 실행 결과 (gnu octave 실행 결과를 캡쳐)


polar()함수가 그림을 그리는 원리는 plot()함수와 동일하다. (거리, 각도) 쌍의 점을 찍은 후 그것들을 선으로 이어주는 것이다. 선의 색/모양을 지정하는 방법은 polt()함수와 동일하다.


>> t=linspace(0,pi,100)
>> polar(sin(t), cos(t))


[그림 2] 실행 결과 (gnu octave 실행 결과를 캡쳐)


이렇게 MATLAB에서 로그스케일 그래프와 극좌표 그래프를 쉽게 생서할 수 있다.



Posted by 살레시오

 여기에서는 이전 포스트들 (part 1, part 2, part 3)에 이어서 plot()함수에 대해서 더 알아보도록 하겠다.  plot()함수에서 그림의 제목과 x축 라벨 그리고 y축 라벨을 써주는 함수는 각각 title(), xlabel(), ylabel()함수이며 이들은 문자열을 입력 파라메터로 받아들여서 그래프에 도시하는 함수이다. 예로서 함수 sin(x)/x 그래프를 그려보자.


>> x=0:0.1:10
>> y=sin(x)./x
>> plot(x,y)


여기서 두 번째 줄에서 y값을 계산하는데 (./)연산자를 사용했음에 주의해야 한다. sin(x)도 벡터이고 x도 벡터이므로 각각의 요소간 나눗셈을 수행해야 하므로 (./)연산자를 사용한다. 이제 이렇게 생성된 그래프에 제목과 x축 라벨 그리고 y축 라벨을 추가해 보자.


>> title('graph of sin(x)/x')
>> xlabel('x')
>> ylabel('y')


결과 그림을 확인해 보면 그래프 상단에 타이틀이 하단에 x축 라벨이 그리고 오른편에 y축 라벨이 표기되어 있음을 확인할 수 있다.

[그림 1] 실행 결과 (gnu octave 를 캡처)


각 축의 화면에 보여지는 범위를 조절하는 함수로서 axis()함수가 있으며 이 함수의 문법은 다음과 같다.


>> axis([XMIN XMAX YMIN YMAX])


여기서 XMIN와 XMAX는 각각 x축의 최소값과 최대값을 그리고 YMIN와 YMAX는 각각 y축의 최소값과 최대값을 표기한다. 예를 들면 다음과 같다.


>> axis([2 8 -0.2 0.6])


이렇게 하면 그래프의 x축 범위가 2에서 8로 그리고 y축 범위가 -0.2에서 0.6으로 조정된다. 또한 axis()함수가 입력파라메터 없이 사용되면 현재 그래프의 xy축의 범위를 반환한다.


>> axis
2.0000 8.0000 -0.2000 0.6000

 MATLAB이 자동으로 축의 범위를 설정하게 하려면 axis auto 명령을 사용한다.


 이제 하나의 그림창에 여러 함수의 그래프를 도시하는 다른 방법을 살펴보자. 먼저 다음과 같은 두함수를 x∈[0,10]의 구간에 대해서 도시하고자 한다.

도시할 데이터를 다음과 같이 생성한다.


>> x = linspace(0,10,50)
>> y1= -0.5*x.^2+1;
>> y2= 2*sqrt(x)-10;


그리고 다음과 같이 그래프를 생성한다.


>> plot(x,y1,'k')
>> hold on
>> plot(x,y2,'k--')
>> grid on
>> legend('-0,5*x^2+1','2*sqrt(x)-10')


여기서 ‘hold on’ 이라는 명령어가 사용되었다. 이 명령어는 지금 활성화된 그래프 창의 내용을 그대로 유지하라는 명령어이다. 그렇기 때문에 두 번째 plot()함수에 의해서 처음에 그려진 그림이 지워지지 않고 그 위에 덧그려진다. 만약 hold on 명령어가 없었다면 첫 번째 그림은 지워지고 두 번째 그림만 남았을 것이다. hold off명령어를 내리지 않는 이상 이후에 호출되는 plot()함수에 의한 결과는 계속 그 위에 덧그려진다. 그리고 grid on 이라는 명령어는 그림에 격자 표시를 하라는 명령어이다. 격자를 지우고 싶으면 grid off라고 하면 된다. 마지막으로 legend()함수는 문자열들을 입력으로 받아서 그래프들에 대한 설명을 생성하는 명령어이다.

[그림 2] 실행 결과 (gnu octave 를 캡처)

 현재 그래프 창의 내용을 지우는 명령은 clf이고 새로운 그래프 창을 생성하는 함수는 figure()함수이다. 위의 예로 다시 설명하면 두 함수를 각각 다른 창에 도시하고 싶을 때는 다음과 같이 하면 된다.


>> plot(x,y1)
>> figure(2), plot(x,y2,'k')


첫 번째 plot()함수로 첫 번째 그림창이 생성되므로 이 그래프 창이 1번의 번호를 가지게 된다. figure(2)라는 명령어로 두 번째 새로운 그래프 창이 생성이 되며 이후의 plot()함수에 의한 그림은 두 번째 창에 그려지게 된다. (캡션바에 그래프 창 번호가 표시됨)


[그림 3] 실행 결과 (gnu octave 를 캡처)

 현재 열려있는 특정한 그림창을 없애는 명령어는 close(n)명령어이다. 이 명령어로 n번 그래프 창이 소멸된다. 현재 생성되어 있는 모든 그림창들을 닫고 싶으면 close all 명령을 내리면 된다.



Posted by 살레시오

 이전 포스트들(part 1, part 2)에 이어서 plot()함수데 대해서 더 알아보도록 하겠다. plot()함수의  자세한 온라인 설명은 help plot 명령으로 볼 수 있으며 다음 표에 옵션을 정리하였다.


[표 1] plot() 함수의 그래프 모양을 정하는 옵션

심벌

색상

심벌

표시자

심벌

선의종류

y

노란색

.

.

-

실선

m

진홍색

o

:

점선

c

청록색

x

x

-.

실선과 점선

g

녹색

+

+

--

실선과실선

b

청색

*

*

k

검은색

^

r

빨간색

v

w

흰색

<

>

s(quare)

d(iamond)

p(entagram)

h(exagram)


위의 표에 소개된 옵션은 다음과 같이 지정할 수 있다.


>> plot(x,y,'option')


만약 option을 지정하지 않으면 기본 설정으로 ‘파란색 실선’의 그래프가 그려지게 되며 위의 option자리에 [표 1]에 소개된 색과 표시자 그리고 선의 종류를 지정하면 된다. 다음의 예를 보자.


>> a=1:0.1:10;
>> b=cos(a);
>> c=sin(a);
>> plot(a,b,'*',a,c,'k')


이 예는 코사인함수는 별표로 그리고 사인함수는 검은색 실선으로 그리는 예제이다. 여기에서 심벌로 지정한 경우는 점의 위치에 그 기호를 표시한다. (선의 경우는 점과 점사이를 잇는다.)

[그림 1] 실행 결과 (gnu octave에서 캡처)

지정자는 색과 표시자 그리고 선의 종류를 동시에 지정할 수도 있다.


>> plot(a,b,'b:',a,c,'m+') 󰎠


이 예에서 ‘b:'지정자는 ‘파란색 점선’으로 그리라는 것이고 ’m+'는 ‘진홍색 십자’로 그리라는 지정이다.


[그림 2] 실행 결과 (gnu octave에서 캡처)

 다음에 몇가지 예제를 더 들었다. 각각을 직접 실행시켜서 짐작대로 그림이 그려지는지 확인해 보자.


plot(x,y,'-.c')

plot(x,y,':r')

plot(x,y,'-sg')

plot(x,y1,'-r',x,y2,'--c')

plot(x,y1,'b',x,y2,'oy')

plot(x,sin(x),'-kp')


위에서 마지막 예를 이용해 보면 다음과 같다.


>> x=linspace(0,2*pi,30);
>> plot(x,sin(x),'-kp')


[그림 3] 실행 결과 (gnu octave에서 캡처)


여기에서 ‘실선’과 ‘펜타그렘’이 동시에 지정되어 있으므로 점의 위치에는 펜타그램이 그려지고 그것들을 실선으로 잇게 되는 것이다. (‘k’지정자에 의해서 색은 검은색)




'MATLAB' 카테고리의 다른 글

MATLAB의 로그 스케일 그래프와 극좌표 그래프  (0) 2015.07.02
MATLAB의 plot()함수 (part 4)  (0) 2015.06.30
MATLAB의 plot()함수 (part 3)  (0) 2015.06.30
MATLAB의 plot()함수 (part 2)  (0) 2015.06.30
MATLAB의 plot()함수 (part 1)  (0) 2015.06.30
MATLAB의 ezplot()함수  (1) 2015.06.27
Posted by 살레시오

 여기에서는 이전 포스트에 이어서 계속 plot()함수에 대해서 설명하도록 하겠다. 만약 sin(x)cos(x)함수를 그리고자 한다면 다음과 같이 하면 될 것이다.


>> x = 0:0.1:2*pi;
>> y = sin(x).*cos(x);
>> plot(x,y)


두 번째 줄에서 sin(x)와 cos(x)두 벡터의 곱에 (.*)연산자를 사용했음에 주의해야 한다. 계산하고자 하는 것이 두 벡터의 요소들 사이의 곱이기 때문이다.


[그림 1] 실행 결과 (gnu octave 에서 캡쳐)


또 다른 예로 x2+xsinx+1 이라는 함수를 도시하고자 한다면 다음과 같이 입력하면 된다.


>> x = -10:0.1:10;
>> y = x.^2 + x .*sin(x) + 1;
>> plot(x,y)


여기에서도 두 번째 줄에 (.^)연산자와 (.*)연산자가 쓰였음을 주의해서 보기 바란다. 요소간 연산을 수행해야 하기 때문이다.


[그림 2] 실행 결과 (gnu octave 에서 캡쳐)


두 번째 입력 파라메터가 행렬인 경우는 각각의 행을 서로 다른 그래프로 보여주게 된다. 다음 예를 보자.


>> x=linspace(0,2*pi,30);
>> y1=sin(x);
>> y2=cos(x);
>> Y=[y1;y2];
>> plot(x,Y)


의 예에서 linspace()함수는 0부터 까지 30개의 등간격 점을 생성하여 반환해주는 함수이다. plot()함수의 두 번째 파라메터 Y는 행렬이다. 아래  결과를 보면 각각의 행에 해당하는 그래프가 따로 도시됨을 알 수 있다.


[그림 3] 실행 결과 (gnu octave 에서 캡쳐)

 

또는 다음과 같이 따로따로 입력해도 같은 결과를 얻을 수 있다.


>> plot(x,y1,x,y2)


다음과 같이 x와 y값을 바꾸어서 입력하면 90도 회전시킨 그림을 도시하게 된다.


>> plot(Y,x)

[그림 4] 실행 결과 (gnu octave 에서 캡쳐)


이제까지 예제에서 그래프의 모양은 기본형이 실선으로 되어 있기 때문에 실선으로 그려졌으나 사용자가 임의로 지정할 수도 있다. 즉 선의 색상, 선의 모양 그리고 선의 종류를 사용자가 직접 지정할 수도 있다. 이것에 대해서는 다음 포스트에서 자세히 알아보도록 하겠다.



'MATLAB' 카테고리의 다른 글

MATLAB의 plot()함수 (part 4)  (0) 2015.06.30
MATLAB의 plot()함수 (part 3)  (0) 2015.06.30
MATLAB의 plot()함수 (part 2)  (0) 2015.06.30
MATLAB의 plot()함수 (part 1)  (0) 2015.06.30
MATLAB의 ezplot()함수  (1) 2015.06.27
MATLAB의 기본적인 수학 함수들  (0) 2015.06.27
Posted by 살레시오

 여기에서는 2차원 그래프를 작성하는데 많이 쓰이는 plot()함수의 사용법에 대해서 자세히 알아보도록 하겠다.(이후 예제는 gnu octave로 모두 실행 가능함) 먼저 다음 예제를 실행해 보자.


>> y=(-3:3).^3; 󰎠
>> plot(y); 󰎠

[그림 1] 실행 결과 (gnu octave 에서 캡쳐)

plot()함수는 입력 파라메터(일단은 벡터라고 가정한다)가 하나일 때는 x축에는 인덱스 값을, y축에는 그 입력 파라메터를 도시한다. 즉 위 예제의 경우

y = [ -27 -8 -1 0 1 8 27]
인덱스  1  2  3  4 5 6 7

이므로 이를 (x,y)쌍으로 나타내면 (1,-27), (2,-8), (3,-1), (4,0), (5,1), (6,8), (7,27) 이고 plot()함수는 이 위치에 점을 찍은 후 이들을 선으로 연결한다. [그림 1]을 자세히 관찰해 보면 이해가 갈 것이다. 만약 점들의 간격이 충분히 좁다면 마치 곡선처럼 보일 것이다.


>> y=(-3:0.1:3).^3;
>> plot(y) 󰎠


[그림 2] 실행 결과 (gnu octave 에서 캡쳐)


점과 점의 간격이 작기 때문에 마치 곡선처럼 보이지만 기본적인 원리는 앞에서 설명한 바와 같다.

 또 다른 예로서 한 주기 동안 시간 t와 시간함수인 사인 함수 y=sin(t)를 그래프로 그려보자. 먼저 기준 축인 시간함수의 벡터를 정해야 한다.


>> t = 0:0.1:2*pi ; %t축 데이터 생성
>> y = sin(t); %y축 데이터 생성
>> plot(t,y)

[그림 3] 실행 결과 (gnu octave 에서 캡쳐)


입력 파라메터가 두개인 경우 두 입력은 크기가 같은 벡터여야 하며 첫 번째 파라메터를 x축에 그리고 두 번째 파라메터를 y축에 대응시켜서 그래프를 그리게 된다.


>> x=-2:5;
>> y=exp(x);
>> plot(x,y)


위의 예에서 다음과 같은 행벡터들을 얻을 수 있다.

x=[ -2 -1 0 1 2 3 4 5]

y=[ 0.14 0.36 1 2.72 7.39 20.09 54.60 148.41]

따라서 plot()함수는 (x,y)순서쌍 (-2,0.14),  (-1,0.36), (0,1), (1,2.72), (2,7.39), (3,20.09), (4,54.60), (5,148.41)의 위치에 점을 찍은 후 이들을 서로 실선으로 연결하게 된다. [그림 4]에 결과가 도시되어 있다.


[그림 4] 실행 결과 (gnu octave 에서 캡쳐)

좀 더 부드러운 곡선을 얻기 위해서는 행벡터 x의 간격만 줄여주면 된다.


>> x=-2:0.01:5;
>> y=exp(x);
>> plot(x,y)

[그림 5] 실행 결과 (gnu octave 에서 캡쳐)



'MATLAB' 카테고리의 다른 글

MATLAB의 plot()함수 (part 3)  (0) 2015.06.30
MATLAB의 plot()함수 (part 2)  (0) 2015.06.30
MATLAB의 plot()함수 (part 1)  (0) 2015.06.30
MATLAB의 ezplot()함수  (1) 2015.06.27
MATLAB의 기본적인 수학 함수들  (0) 2015.06.27
MATLAB의 다항식(polynimial) 다루기 (part 2)  (0) 2015.06.27
Posted by 살레시오

MATLAB의 ezplot()함수

MATLAB 2015.06.27 12:13

 여기에서는 MATLAB에서 그래픽을 그려주는 함수들중 ezplot() 함수에 대해서 설명하겠다. MATLAB에서 자주 사용되는  그래픽 함수는 다음과 같다. (아래 함수들은 gnu octave 4.x 이상 버전에서 모두 지원한다.)


[표 1] 그래프 함수들

함수명

비고

ezplot(), fplot()

plot()

ezpolar(), polar()

loglog(), semilogx(), semilogy()

2차원 그래프

ezplot3(), ezmesh(), ezsurf()

plot3(), mesh(), surf()

3차원 그래프


 이것들 중 ezplot()함수는 사용자가 쉽게 (easy) 수학 함수의 그래프를 그릴 수 있도록 해주는 함수이고 fplot()함수는 함수(function)의 그래프를 그려주는 함수이다. 두 함수 모두 어떤 함수의 2차원 그래프를 그려주는 함수로서 명시적 함수( y=f(x) 의 형식) 과 암시적 함수 ( f(x,y)=0 의 형식)을 도시하는데 있어서 plot()함수보다는 보다 편하게 사용할 수 있다.

 가장 간단한 사용법은 도시하고자 하는 함수를 문자열로 입력하면 된다. 예를 들어 y=cos(x) 함수를 그리고 싶다면


>> ezplot('cos(x)')


[그림 1] gnu octave 4.x에서의 실행 결과


명시적 함수 y=f(x) 에서 f(x) 부분만 문자열로 입력하면 된다. y=sqrt(1-x2) 을 그리고 싶다면 다음과 같이 입력한다.


>> ezplot('sqrt(1-x^2)')

[그림 2] gnu octave 4.x 에서의 실행 결과


 만약 x축의 범위를 명시하고 싶다면 두 번째 파라메터로 x축의 최소값, 최대값을 가지는 행벡터를 입력하면 된다.


>> ezplot('sin(x)*cos(x)', [0 pi])

[그림 3] gnu octave 4.x 에서의 실행 결과


위의 명령은  sin(x)cos(x) 그래프를 그리되 x축의 범위를 0과 π사이로 하라는 것이다.


 암시적 함수도 그대로 문자열로 입력하면 된다. 예를 들어서 원의 방정식 x2+y2=16 을 그리고 싶다면 우변을 0으로 만든 다음 x2+y2-16=0 그 좌변을 입력하면 된다.


>> ezplot('x^2+y^2-16')


[그림 4] gnu octave 4.x 에서의 실행 결과


 이 경우도 두 번째 파라메터로 x축의 범위를 입력할 수 있다. 예를 들어 x3+y3-5xy+0.2=0 의 그래프를 그리되 x축의 범위를 [-3 3]로 하고 싶다면 다음과 같이 하면 된다.


>> ezplot('x^3 + y^3 - 5*x*y + 0.2',[-3,3])

[그림 5] gnu octave 4.x 에서의 실행 결과


 ezplot()함수는 x와 y의 관계식이 직접적으로 나타내지지 않고 매개변수에 의해서 관계 지어지는 경우도 도시할 수 있다. 예를 들어서 매개변수 t 에 대해서 x=sin(t) , y=cos(t)의 관계식이 주어진 경우 이들을 순서대로 입력하면 된다.


>> ezplot('sin(t)','cos(t)')

[그림 6] gnu octave 4.x 에서의 실행 결과


또한 다음과 같은 관계식을 가지는 x-y그래프를 t가 [0 π]의 범위를 가질 경우에 대해서 도시하려면

다음과 같이 하면 된다.


>> ezplot('sin(3*t)*cos(t)', 'sin(3*t)*sin(t)', [0,pi])

[그림 7] gnu octave 4.x 에서의 실행 결과


즉, 매개변수를 가지는 xy그래프를 그릴 때 세 번째 입력벡터는 매개변수의 범위가 된다.



Posted by 살레시오

 여기에서는 MATLAB에 내장된 수학 함수들을 소개한다.


  MATLAB의 모든 삼각함수는 입력과 출력의 단위로 라디안(radian)을 사용한다. 라디안을 도(degree)로 혹은 도를 라디안으로 변환하려면 다음과 같이 하면 된다.


>> angle_deg = angle_rad*(180/pi)
>> angle_rad = angle_deg*(pi/180)


[표 1] 삼각함수 관련 함수들

함수

기능

sin(x)

x의 sine값 계산

cos(x)

x의 cosine값 계산

tan(x)

x의 tangent값 계산

asin(x)

x의 arcsine값 계산. 단 x는 -1과 1사이의 값이어야 하고 결과값은 라디안으로 -π/2 와 π/2 사이의 값이다.

acos(x)

x의 arccosine값 계산. 단 x는 -1과 1사이의 값이어야 하고 결과값은 라디안으로 -π/2 와 π/2 사이의 값이다.

atan(x)

x의 arctangent값 계산. 결과값은 라디안으로 -π/2 와 π/2 사이의 값이다.

atan2(y,x)

y/x 의 arctangent값 계산. 결과값은 라디안으로 -π/2 와 π/2 사이의 값이다.


[표 2] 쌍곡선 함수들

함수

기능

sinh(x)

hyperbolic sin 함수

cosh(x)

hyperbolic cos 함수

tanh(x)

hyperbolic tan 함수

asinh(x)

각각의 역함수들

acosh(x)

atanh(x)


[표 3] 공통 함수

함수

기능

abs(x)

절대값(=|x|)

sqrt(x)

sqrtm(A)

제곱근( square root) 요소간 연산

제곱근. 행렬 연산

round(x)

반올림

fix(x)

x를 0의 방향의 가장 가까운 정수로

floor(x)

x를 -inf 의 방향의 가장 가까운 정수로

ceil(x)

x를 +inf 의 방향의 가장 가까운 정수로

sign(x)

x<0이면 -1, x=0이면 0, x>0이면 1을 반환

rem(x,y)

x를 y로 나눈 나머지

exp(x)

expm(A)

지수 함수, 요소간 연산

지수 함수, 행렬 연산

log(x)

logm(A)

밑이 e인 로그(자연 로그) 요소 간 연산

밑이 e인 로그(자연 로그) 행렬 연산

log10(x)

밑이 10인 로그

MATLAB에서 밑이 e인 로그함수는 log()이다. (ln()이 아님에 유의할 것) 밑이 10일 로그 함수는 log10 이다. 보통 수학책에서는 log 가 밑이 10인 로그 함수를 표기하는데 MATLAB에서는 밑이 e인 함수명이라는 것에 유의해야 한다.


[표 4] 복소수 관련 함수

함수

기능

abs(s)

복소수 s의 크기를 계산 (=|s|)

real(s)

복소수 s의 실수부 계산

imag(s)

복소수 s의 허수부 계산

conj(s)

켤레복소수

angle(s)

복소수의 각도(radian)를 계산


여기에서 소개한 함수들은 MATLAB의 기본적인 내장 함수이므로 별다른 조치 없이 바로 사용할 수 있다.



Posted by 살레시오

 여기서는 이전 포스트에 이어서 다항식 관련 함수들에 대해서 알아보겠다. 두 다항식의 곱셈은 conv()함수에 의해서 구해진다. 두 함수 a(x)=x3+2x2+3x+4 와 b(x)=x3+4x2+9x+16의 곱을 구해보자.


>> a=[1 2 3 4], b=[1 4 9 16]
a =
  1   2   3   4
b =
   1    4    9   16

>> c=conv(a,b)
c =
   1    6   20   50   75   84   64


이 결과로부터

c(x) = x6+6x5+20x4+50x3+75x2+84x+64

임을 알 수 있다.


 방정식간의 나눗셈은 deconv()함수를 이용한다. 앞에서 구한 방정식 c(x)와 b(x)를 이용하여 를 구하면 다음과 같다.


>> deconv(c,b) 󰎠
1.0000 2.0000 3.0000 4.0000


결과가 a(x)와 같음을 알 수 있다.


다항식의 미분을 계산하는 함수로 polyder()함수가 있다. 예를 들어 p(x)=x3-2x-5 의 미분함수를 계산하려고 한다면


>> p=[1 0 -2 -5]
p =
1 0 -2 -5
>> q=polyder(p)
q =
3.0000 0 -2.0000


여기서 p’(x)=3x2-2임을 알 수 있다. 만약 함수 polyder()로의 입력벡터가 두 개면 두 다항식의 곱의 미분을 구해준다.


>> a=[1 3 5], b=[2 4 6]
a =
1 3 5
b =
2 4 6

>> polyder(a,b)
8.0000 30.0000 56.0000 38.0000


여에서 a(x)=x2+3x+5와 b(x)=2x2+4x+6의 곱함수 a(x)b(x)을 미분한 결과는 8x3+30x2+56x+38 이라는 것을 알 수 있다. 마지막으로 입력도 두개이고 출력도 두개라면, 예를 들어서


>> [q,d]=polyder(a,b) 󰎠
q =
-2 -8 -2
d =
4 16 40 48 36


와 같이 사용한다면 a(x)/b(x)를 미분한 결과를 q(x)/d(x)로 반환해 준다. 즉 다음과 같은 다항식

을 미분한 결과는 다음과 같다는 것을 알 수 있다.


 이제 역라플라스 변환이나 역z변환에서 자주 쓰이는 부분분수전개 (partial fraction expansion) 에 대해서 살펴보자. 부분분수 전개를 수행하는 함수는 residue()함수이다. 어떤 다항식 a(x)와 b(x)에 대해서 a(x)가 다중근을 갖지 않는다면 분수식은 다음과 같이 전개될 수 있다.

예를 들어서

를 부분분수전개한 결과를 보려면 다음과 같이 하면 된다.


>> a=[1 6 8], b=[-4 8]
a =
1 6 8
b =
-4 8

>> [r1,p1,k1]=residue(b,a)
r1 =
-12
8
p1 =
-4
-2
k1 =
Empty matrix: 0-by-0


이것으로 부터 다음을 알 수 있다.

다음과 같이 다중근을 가지는 경우를 살펴보자.

이 식을 부분분수전개를 하려면 예를 들어 다음과 같이 하면 된다.


>> [r2, p2, k2] = residue( 1, poly([1 0.5 0.5]) )
r2 =
4.0000
-4.0000
-2.0000
p2 =
1.0000
0.5000
0.5000
k2 =
Empty matrix: 0-by-0


이 예에서 분모항 을 입력하기 위해서 poly()함수를 사용다. 반환된 벡터 r2,p2 그리고 k2로부터 다음과 같은 결과로 해독한다.

여기서 p2벡터의 세 번째 요소인 두 번째 0.5 는 이중근임을 나타냄에 유의해야 한다.


 지금까지 소개한 다항식에 관련된 함수들을 표로 정리하면 다음과 같다.


[표 1] 다항식에 관련된 함수들 정리

함수

기능

roots()

방정식의 해를 구한다

poly()

해로부터 방정식을 구한다

polyval()

특정한 스칼라값에서의 방정식의 값을 구한다

polyvalm()

특정한 행렬을 대입한 방정식의 값을 구한다

residue()

부분분수전개를 한다.

polyder()

방정식의 미분을 구한다.

conv()

방정식의 곱을 구한다.

deconv()

방정식의 나눗셈을 수행한다.



Posted by 살레시오

 MATLAB은 다항식을 다루는 몇 개의 함수를 제공한다. 예를 들어 방정식의 근을 구한다던가, 어떤 다항식의 미분 혹은 적분을 구하는 함수등이 제공된다.

 예를 들어서 x4-12x3+25x+116 이라는 다항식은 다음과 같이 계수들을 순서대로 가지는 p라는 변수에 저장한다.


>> p=[1 -12 0 25 116]
p =
    1   -12     0    25   116


여기서 계수들을 입력할 때 다항식이 사실 x4-12x3+0x2+25x+116 이므로  x2항의 계수인 0도 반드시 써야 한다는 것에 주의하자. 이제 이 방정식의 근은 roots()함수를 이용하여 구할 수 있다.


>> r=roots(p)
r =
11.747282870823
2.702820743841
-1.22505180733187 + 1.46720800872890i
-1.22505180733187 - 1.46720800872890i


위에서 보듯이 roots()함수는 해들은 열벡터에 저장하여 반환한다. 반대로 해로부터 방정식을 구할 수도 있다. 예를 들어 어떤 방정식의 해가 1,2,3,4,5 일 때 이들을 해로 가지는 방정식은 (x-1)(x-2)(x-3)(x-4)(x-5) 일 것이다. 이 방정식의 계수를 구하는 함수가 poly()함수이다.


>> poly([1 2 3 4 5])
ans =
    1   -15    85  -225   274  -120


위에서 1, 2, 3, 4, 5를 해로 가지는 방정식은 x5-15x4+85x3-225x2+274x-120 이라는 것을 알 수 있다. 처음 나왔던 예제에서 roots()함수를 이용하여 방정식의 근을 r에 저장하였다. 이것을 이용하여 원래의 방정식을 구할 수 있다.


>> pp=poly(r)
pp =
1.0000 -12.0000 -0.0000 25.0000 116.0000


입력변수가 nxn 크기의 행렬인 경우 poly()함수는 이 행렬의 특성방정식(characteristic equation)의 계수를 반환한다.


>> A=[1.2 3 -0.9; 5 1.75 6; 9 0 1] 󰎠
A =
1.2000 3.0000 -0.9000
5.0000 1.7500 6.0000
9.0000 0 1.0000
>> poly(A) 󰎠
1.0000 -3.9500 -1.8500 -163.2750


행렬의 특성방정식이란 det(xI-A) 를 말하며 위의 경우 이 특성방정식이 x3-3.95x2-1.85x-163.275 임을 알 수 있다. 이 특성방정식의 해가 바로 행렬의 고유치(eigenvalue)이다. 행렬의 고유치를 구하는 함수는 따로 마련되어 있으며 eig()함수이다.


>> roots(ans)
7.2826
-1.6663 + 4.4321i
-1.6663 - 4.4321i

>> eig(A)
7.2826
-1.6663 + 4.4321i
-1.6663 - 4.4321i


특정한 값에서의 함수값을 구하는 함수는 polyval()함수이다. 예를 들어 x=5에서 p에 저장된 다항식 값을 구하려면


>> polyval(p,5) 󰎠
-634


정방행렬에 대해서 다항식의 값을 구하는 것도 가능하다. 예를 들어서 q(x)=x3-2x-5 라는 함수가 있을때 행렬 A에 대해서 q(A)=A3-2A-5I 가 되는 것이다. 여기서 I는 항등행렬이다. 다항식에 대해서 행렬 계산을 구하는 함수가 polyvalm()함수이다.


>> A=[2 4 5;-1 0 3; 7 1 5]
A =
2 4 5
-1 0 3
7 1 5

>> q=[1 0 -2 -5]
q =
1 0 -2 -5

>> polyvalm(q,A)
377 179 439
111 81 136
490 253 639


위의 결과에 대해서 검산을 하면 다음과 같다.


>> A^3-2*A-5*eye(3) 󰎠
377 179 439
111 81 136
490 253 639


여기서 eye(3)은 3x3 항등행렬을 생성해주는 함수이다.



Posted by 살레시오

 MATLAB에서 제공하는 논리 연산자는 다음 표와 같다.


[표 1] 논리 연산자

논리연산자

기능

&

AND

|

OR

~

NOT


(부연하자면 C/C++의 비트연산자와 같다. MATLAB은 비트 연산자는 없다.)

 다음 예를 보자


>> A=magic(3)
A =
8 1 6
3 5 7
4 9 2

>> B=pascal(3)
B =
1 1 1
1 2 3
1 3 6

>> ~(A>B)
0 1 0
0 0 0
0 0 1

>> (A>2)&(B<3)
1 0 1
1 1 0
1 0 0


 이전 포스트에서 설명한 기본 논리/관계연산자 외에도 MATLAB은 추가적인 다양한 논리/관계 함수를 제공한다.


[표 2] 관계/논리 연산 관련 함수들 (X, Y는 함수)

함수명

함수의 기능

Z=xor(X,Y)

XOR 연산을 수행한다

Z=any(X)

벡터 X의 요소 중 하나라도 0이 아닌 요소가 있으면 1을 반환, X가 행렬일때는 열단위로 연산된다.

z=all(X)

벡터 X의 요소들 모두 0이 아니면 1을 반환

(X가 행렬일 때는 열단위로 연산된다.)

z=isempty(X)

행렬 X가 빈행렬이면 1을 반환

z=isequal(X,Y)

행렬 X와 Y가 완전히 동일하면 1을 반환

Z=isinf(X)

행렬 X의 요소 중 inf가 있는 곳에 1을 반환

Z=isnan(X)

행렬 X의 요소 중 nan이 있는 곳에 1을 반환

z=isreal(X)

행렬 X의 요소 중 복소수가 하나도 없으면 1을 반환


위의 도표에서 소문자 z가 출력변수이면 그 함수는 스칼라값(0혹은1)을 반환함을 의미하며 대문자 Z가 출력행렬이면 그 함수는 입력 파라메터와 동일한 크기를 갖는 행렬을 반환한다는 것을 의미한다.


 이 중 isequal()함수는 연산자 (==)와 같은 역할을 하는 것 아닌가 생각할수 있으나 (==)연산자는 요소간 비교이고 isequal()함수는 두 행렬 전체를 비교하는 함수이다. 다음 예에서 이를 확인해 보라.


>> A=[1 2;3 4]
A =
1 2
3 4

>> B=[1 2;3 4]
B =
1 2
3 4

>> C=[1 2; 3 5]
C =
1 2
3 5

>> A==B
1 1
1 1

>> A==C
1 1
1 0

>> isequal(A,B)
1

>> isequal(A,C)
0




Posted by 살레시오