next up previous
次へ: 複素PARCOR係数 上へ: 複素自己回帰モデルと係数の高速計算法 戻る: 複素自己回帰モデルと係数の高速計算法

複素自己回帰モデル

平面図形の輪郭線(平面曲線)を追跡して得られる点列を $(x_j,y_j)(j=0,1,...,N-1)$とし、その複素表現を$z_j=x_j+iy_j$とする。

このとき、$m$次の複素自己回帰モデルは、

\begin{displaymath}
\hat{z}_j=\sum_{k=1}^{m}a_kz_{j-k}
\end{displaymath} (298)

のように輪郭点を$m$個前までの輪郭点の線形結合で近似するモデルとして定義され る。ここで、モデルの係数 $\{a_k\}_{k=1}^{m}$ は平均2乗予測誤差 $\varepsilon^2(m)$
\begin{displaymath}
\varepsilon^2(m)=\mbox{E}_j\vert\hat{z}_j-z_j\vert^2=\mbox{...
...boldmath$a$} - \mbox{\boldmath$a$}^T \mbox{\boldmath$r$} + r_0
\end{displaymath} (299)

が最小となるように決められる。ただし、
\begin{displaymath}
R = \left[
\begin{array}{rrrrr}
r_{0}, & \bar{r}_{1}, & ...
...1}, & r_{m-2}, & r_{m-3}, & ..., & r_{0}
\end{array} \right],
\end{displaymath} (300)


\begin{displaymath}
\mbox{\boldmath$a$} = [a_{1}, a_{2}, a_{3},..., a_{m}]^T,
\end{displaymath} (301)


\begin{displaymath}
\mbox{\boldmath$r$} = [r_{1}, r_{2}, r_{3},..., r_{m}]^T,
\end{displaymath} (302)


\begin{displaymath}
r_{k}=\mbox{E}_{j}(z_{j}\bar{z}_{j-k})
\end{displaymath} (303)

である。また、$\mbox{E}_{j}$は平均操作( $(1/N)\sum_{j=0}^{N-1}$)を表わし、 $\bar{ }$ は複素共役、$^T$は転置を表わす。$r_k$は複素相関係数であり、$R$は Hermite 行列 ($R^*=R, \ ^*$は共役転置)となる。

平均2乗予測誤差(式(7.4))を最小とする最適な複素自己回帰係数 $\{a_k\}_{k=1}^{m}$は、

\begin{displaymath}
\mbox{\boldmath$a$}=R^{-1}\mbox{\boldmath$r$}
\end{displaymath} (304)

で与えられる。

また、このとき達成される最小平均2乗予測誤差 $\hat{\varepsilon}^2(m)$

\begin{displaymath}
\hat{\varepsilon}^2(m)=r_0-\mbox{\boldmath$r$}^T \bar{\mbox{\boldmath$a$}}
\end{displaymath} (305)

となる。

こうして得られた複素自己回帰係数の性質について考えよう。輪郭点列が複素表現され ている場合、原点まわりの回転は各輪郭点列を $e^{i\theta}$ 倍することに対応する。 回転された輪郭点列 $\{e^{i\theta}z_j\}_{j=0}^{N-1}$ に対する $r_k$ は、式( 7.8)から求められる。そこでは $e^{i\theta}$ 同士は打ち消しあい、 $r_k$ は回転によって変わらない。従って、式(7.9)も不変となり、結局、 複素自己回帰係数も回転不変となる。また、式(7.8)から、$r_k$ は輪郭 を追跡する際の始点位置の選び方に依存しない量である。従って、それに基づいて計 算される複素自己回帰係数も始点位置の選び方に依存しない量となる。これらの性質 は、形の記述、表現という意味で好ましいものである。さらに、興味深い性質として 輪郭線の追跡方向(時計回りか反時計回り)を変えると、得られる複素自己回帰係数 が複素共役となる。これは、輪郭点列を逆方向に追跡した場合には、式 (7.3)の右辺は $z$ の添字 $-k$$+k$ としたものとなり、$r_k$ が 複素共役となることから明らかである。



Takio Kurita 平成14年7月3日