시스템이 식(1)과 같이 표현된다.
dtndny(t) + an−1dtn−1dn−1y(t) +… + a1dtdy(t) + a0y(t) = bn−1dtn−1dn−1u(t) + bn−2dtn−2dn−2u(t) +… + b1dtdu(t) + b0u(t)(1)
이를 전달함수로 표현하면 식(2)와 같아진다.
G(s) = U(s)Y(s)=sn+an−1sn−1+…+a1s+a0bn−1sn−1+bn−2sn−2+…+b1s+b0(2)
식(2)를 관측 가능 표준형으로 표현하기 위해 식(3)과 같이 식(2)를 크로스로 곱한다.
(sn+an−1sn−1+…+a1s+a0)Y(s)=(bn−1 sn−1+bn−2sn−2+…+b1s+b0)U(s)(3)
식(3)에서 s가 없는 항을 왼쪽에 나머지를 오른쪽으로 옮기고 상태변수를 식(4)와 같이 정의한다.
b0U(s)−a0Y(s)==(sn+an−1sn−1+…+a1s)Y(s)−(bn−1 sn−1+bn−2sn−2+…+b1s)U(s)=sX1(s)(4)
식(4)에서 첫번째 식과 마지막 식을 모아서 식(5)와 같이 만든다.
b0U(s)−a0Y(s)=sX1(s)(5)
식(4)에서 두번째 식과 마지막 식으로 하나의 방정식을 만들고 이 식을 s로 나누면 식(6)과 같아진다.
X1=(sn−1+an−1sn−2+…+a1)Y(s)−(bn−1 sn−2+bn−2sn−3+…+b1)U(s)(6)
식(6)에서 s가 없는 항을 왼쪽으로 이항하고 식(7)과 같이 상태변수를 정의한다.
X1(s)+b1U(s)−a1Y(s)=(sn−1+an−1sn−2+…+a2s)Y(s)−(bn−1 sn−2+bn−2sn−3+…+b2s)U(s)=sX2(s)(7)
식(7)의 첫번째 식과 마지막 식으로 식(8)과 같이 정리한다.
X1(s)+b1U(s)−a1Y(s)=sX2(s)(8)
식(7)의 두번째 식과 마지막 식을 s로 나줘주고 정리하면 식(9)와 같아진다.
X2=(sn−2+an−1sn−3+…+a2)Y(s)−(bn−1 sn−3+bn−2sn−4+…+b2)U(s)(9)
식(9)에서 s가 없는 항을 왼쪽으로 이항하여 식(10)을 만든다.
X2(s)+b2U(s)−a2Y(s)=(sn−2+an−1sn−3+…+a3s)Y(s)−(bn−1 sn−2+bn−2sn−3+…+b3s)U(s)=sX3(s)(10)
식(10)의 첫번째 식과 마지막 식을 정리하면 식(11)과 같아진다. 위의 과정을 반복한다.
X2(s)+b2U(s)−a2Y(s)=sX3(s)(11)
⋮
Xn−1(s) + bn−1 U(s) − an−1Y(s) = sY(s) = sXn(s)(12)
Xn−1(s)+bn−1U(s)−an−1Y(s)=sXn(s)=sY(s)(13)
식(5), (8), (11), ... 을 라플라스 역변환하고 식(13)을 적용하면 식(14)와 같이 쓸 수 있다.
x1(t)˙=−a0y(t)+b0u(t)⇒x1(t)˙=−a0xn(t)+b0u(t)
x2(t)˙=x1(t)−a1y(t)+b1u(t)⇒x2(t)˙=x1(t)−a1xn(t)+b1u(t)x3(t)˙=x2(t)−a2y(t)+b2u(t)⇒x3(t)˙=x2(t)−a2xn(t)+b2u(t)⋮xn(t)˙=xn−1(t)−an−1y(t)+bn−1u(t)⇒xn(t)˙=xn−1(t)−an−1xn(t)+bn−1u(t)(14)
위 식을 행렬식으로 표현하면 식(15)와 같다.
⎣⎡x1(t)˙x2(t)˙x3(t)˙⋮xn(t)˙⎦⎤=⎣⎡01000010………⋮…0001−a0−a1−a2−an−1⎦⎤⎣⎡x1(t)x2(t)x3(t)⋮xn(t)⎦⎤+⎣⎡b0b1b2⋮bn−1⎦⎤u(t)(15)
y(t)=[000…1]⎣⎡x1(t)x2(t)x3(t)⋮xn(t)⎦⎤+[0]u(t)
여기서,
A =⎣⎡01000010………⋮…0001−a0−a1−a2−an−1⎦⎤ , B =⎣⎡b0b1⋮bn−1⎦⎤, C=[0000…1], D=[0]
Login to write a comment.