next up previous
Next: 連続系の状態変数を離散値系で表示する方法 Up: 状態変数 Previous: 部分分数分解法(Diagonal Form)

繰り返し法(Iterative Programming)

(4.12)式はまた次のように書くことができる。
\begin{displaymath}
\frac{Y(Z)}{U(Z)} = \frac{\alpha(Z+Z_1)(Z+Z_2)\cdots(Z+Z_m)}
{(Z+p_1)(Z+p_2)\cdots(Z+p_n)}
\end{displaymath} (4.44)

これを図4.9のように分けて考える。ただし
\begin{displaymath}
\alpha_1\cdot\alpha_2\cdot\cdots\cdot\alpha_n=\alpha
\end{displaymath} (4.45)

であり、$n > m$であるとする。図より
\begin{displaymath}
X_n(Z) = \frac{\alpha_n}{Z+p_n}U(Z)
\end{displaymath} (4.46)

より
\begin{displaymath}
ZX_n(Z) = -p_nX_n(Z) + \alpha_nU(Z)
\end{displaymath} (4.47)

となる。$X_{m+1}(Z)$のところまでは同様にして
\begin{displaymath}
ZX_{m+1}(Z) = -p_{m+1}X_{m+1}(Z) + \alpha_{m+1}X_{m+2}(Z)
\end{displaymath} (4.48)

となる。$X_m(Z)$のところは
$\displaystyle ZX_m(Z)$ $\textstyle =$ $\displaystyle -p_mX_m(Z) + \alpha_m(Z+Z_m)X_{m+1}(Z)$  
  $\textstyle =$ $\displaystyle -p_mX_m(Z)+\alpha_m(Z_m-p_{m+1})X_{m+1}(Z)+\alpha_m\alpha_{m+1}X_{m+2}(Z)$ (4.49)

となり、以下$X_1(Z)$まで同様の方法で求める。したがって、一般的に状態方程式は
    $\displaystyle \left[ \begin{array}{c}
x_1(k+1) \\
x_2(k+1) \\
\vdots \\
x_m(k+1) \\
\vdots \\
x_n(k+1)
\end{array} \right]
=$  
    $\displaystyle \left[ \begin{array}{ccccccc}
-p_1 & \alpha_1(Z_1-p_2) & \alpha_1...
...}{c}
0 \\
0 \\
\cdot \\
\cdot \\
\cdot \\
\alpha_n
\end{array} \right]u(k)$  

という形になる。

図 4.9:
\begin{figure}\begin{center}
\psbox[scale=0.60]{eps/4-2-7.eps} \end{center} \end{figure}

[例] ${\displaystyle \frac{Y(Z)}{U(Z)} = \frac{3(Z+5/3)}{(Z+2)(Z+3)}}$ の場合

4.10のように状態変数をとると$X_2(Z)$のところは

\begin{displaymath}
ZX_2(Z) = -2X_2(Z) + 3U(Z)
\end{displaymath} (4.50)

となり、$X_1(Z)$のところは
$\displaystyle ZX_1(Z)$ $\textstyle =$ $\displaystyle -3X_1(Z) + (Z +5/3)X_2(Z)$  
  $\textstyle =$ $\displaystyle -3X_1(Z) +\frac{5}{3}X_2(Z) +ZX_2(Z)$  
  $\textstyle =$ $\displaystyle -3X_1(Z) + \frac{5}{3}X_2(Z) -2X_2(Z) + 3U(Z)$  
  $\textstyle =$ $\displaystyle -3X_1(Z) -\frac{1}{3}X_2(Z) +3U(Z)$ (4.51)

となる。したがって状態方程式は
\begin{displaymath}
\begin{array}{rcl}
\left[ \begin{array}{c}
x_1(k+1) \\
...
...
3
\end{array} \right]u(k)\\
y(k) & = & x_1(k)
\end{array}\end{displaymath} (4.52)

となる。
図 4.10:
\begin{figure}\begin{center}
\psbox[scale=0.60]{eps/4-2-8.eps} \end{center} \end{figure}

next up previous
Next: 連続系の状態変数を離散値系で表示する方法 Up: 状態変数 Previous: 部分分数分解法(Diagonal Form)
Yasunari SHIDAMA
平成15年6月30日