Skip to main content

Controllable and observable forms (denominator 0.36z⁻²)

Short answer

Controllable form: G=[[0,1],[-0.36,1.368]], H=[[0],[1]], C=[0.264,0.368]. Observable form: G=[[0,-0.36],[1,1.368]], H=[[0.264],[0.368]], C=[0,1]. Both give Y/U=(0.368z+0.264)/(z^2-1.368z+0.36). Poles at 0.3556, 1.0124 (unstable).

Obtain the state space representation of following pulse transfer function in (i) controllable canonical form (ii) observable canonical form: Y(z)/U(z)=(0.368z⁻¹+0.264z⁻²)/(1-1.368z⁻¹+0.36z⁻²).

Here the constant term of the denominator is 0.36, not 0.368. The method is the same; only the entry −a2-a_2 changes.

Coefficients

Y(z)U(z)=0.368z−1+0.264z−21−1.368z−1+0.36z−2=0.368z+0.264z2−1.368z+0.36\frac{Y(z)}{U(z)} = \frac{0.368z^{-1}+0.264z^{-2}}{1-1.368z^{-1}+0.36z^{-2}} = \frac{0.368z+0.264}{z^2-1.368z+0.36}
a1a_1a2a_2b0b_0b1b_1b2b_2
−1.3680.3600.3680.264

(i) Controllable canonical form

Let Q(z)U(z)=1z2−1.368z+0.36\dfrac{Q(z)}{U(z)} = \dfrac{1}{z^2-1.368z+0.36} with x1(k)=q(k)x_1(k) = q(k), x2(k)=q(k+1)x_2(k) = q(k+1):

x1(k+1)=x2(k)x2(k+1)=−0.36x1(k)+1.368x2(k)+u(k)y(k)=0.264x1(k)+0.368x2(k)\begin{aligned} x_1(k+1) &= x_2(k)\\ x_2(k+1) &= -0.36x_1(k) + 1.368x_2(k) + u(k)\\ y(k) &= 0.264x_1(k) + 0.368x_2(k) \end{aligned} G=[01−0.361.368],H=[01],C=[0.2640.368],D=0\mathbf{G} = \begin{bmatrix} 0 & 1 \\ -0.36 & 1.368 \end{bmatrix},\quad \mathbf{H} = \begin{bmatrix} 0 \\ 1 \end{bmatrix},\quad \mathbf{C} = \begin{bmatrix} 0.264 & 0.368 \end{bmatrix},\quad D = 0

(ii) Observable canonical form

From Y=z−1[0.368U+1.368Y+z−1(0.264U−0.36Y)]Y = z^{-1}\left[0.368U + 1.368Y + z^{-1}(0.264U - 0.36Y)\right] with x2=yx_2 = y:

x1(k+1)=−0.36x2(k)+0.264u(k)x2(k+1)=x1(k)+1.368x2(k)+0.368u(k)y(k)=x2(k)\begin{aligned} x_1(k+1) &= -0.36x_2(k) + 0.264u(k)\\ x_2(k+1) &= x_1(k) + 1.368x_2(k) + 0.368u(k)\\ y(k) &= x_2(k) \end{aligned} G=[0−0.3611.368],H=[0.2640.368],C=[01],D=0\mathbf{G} = \begin{bmatrix} 0 & -0.36 \\ 1 & 1.368 \end{bmatrix},\quad \mathbf{H} = \begin{bmatrix} 0.264 \\ 0.368 \end{bmatrix},\quad \mathbf{C} = \begin{bmatrix} 0 & 1 \end{bmatrix},\quad D = 0

Check and remark

  • Both give C(zI−G)−1H=0.368z+0.264z2−1.368z+0.36\mathbf{C}(z\mathbf{I}-\mathbf{G})^{-1}\mathbf{H} = \dfrac{0.368z+0.264}{z^2-1.368z+0.36}.
  • The poles are z=0.3556z = 0.3556 and z=1.0124z = 1.0124. The second lies just outside the unit circle, so this system is (slightly) unstable, unlike the 0.368 version, which has a pole at exactly z=1z = 1.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 5 questions →