每-1 y(i)=∑g(j)〔u(i-)+z(i-i)〕 (8) j=0 如果噪声：()为相关噪声，那么就不能得到一致性估计。 1,4考虑过程中存在零漂f()的情况 由于f()是趋势项，它可表示为： f)=∑w(j)i (9) i=0 其中 j=1,2,…,N 通过实验就可按式(9)来拟合趋势项f(k)。 1.5交替矩阵(Rotation Matrix)法 上述参数模型有时可简写成 y=60+e (10) 其中 y=〔y(1)y(2)…y(N)]T多 pT(1) 中= pT(2) T(N)) -y(0) -y(-1)…-y(1-n),4(0)u(-1)…u(1-m) -y(1) -y(0)… -y(2-n),4(1)u(0)…u(2-m) :;: -y(N-1)-y('-2)…-y(N-n),u(N-1)u(N-2)…u(N-M)5 0T=〔a:…0.b1…baJ eT=〔e(1)e(2)e(N)〕。 式(10)的未知参数可用多种参数估计方法求得。如用最小二乘法的话，参数估计值为 8=(中T中)-'中Ty=(φ*)-1Ty (11) 要使上式有解，必须中T阵为非奇异。如果中T中阵为奇异，则称中？中为病态矩阵。这时 可用Householder?变换的方法求解未知参数，但这方法运算较麻烦。这里介绍一种简单的方 法一交替矩阵T法，其主要思想是将φ·阵经过多次选得的交替矩阵左乘后，变为一个三角阵。 而T?T=【(单位阵)，说明参数估计的误差平方和是不变的。 T阵的一般形式为 260伍 一 ， ‘ 习 。 〔 ‘ 一 ‘ 一 〕 如 果噪 声 旬 为 相关噪 声 ， 那 么就 不能 得到 一 致 性估计 。 考虑过 程 中存在零 漂 的情 况 由干厂 是趋势项 ， 它 可表示 为 护 吞 二 习 却 二 其 中 二 ， ， … ， 通过 实验就可 按式 来拟 合 趋势项 。 交替 矩阵 法 上述参数模 型有时可简写 成 价 其中 二 〔夕 一 〕 甲 伊 卯 一 一 一 一 … 一 … 一 一 ” ， 赵 一 又 一 ， “ 一 · ” “ … ‘ 一 优 ‘ 一 用 一 少 刃 一 一 夕 入 尹 一 · 一 少 一 ， “ 一 一 一 一 、‘ ‘… 苗︸一一 口 二 〔 口 一 。 … 二 〕 〔 … 〕 。 式 的未知参数可 用多种 参数估计方法求 得 。 如 用 最小 二乘 法的 话 ， 参数估计值 为 口 二 价 沪 一 ’ 功 夕 二 功 一 ‘ 价 要使 上 式 有解 ， 必须 沪 功阵为非 奇异 。 如果 功 功阵为 奇异 ， 则称 价 必为病 态 矩 阵 。 这 时 可 用 变换 的方法求解未知参数 ， 但 这 方 法运算较麻 烦 。 这 里介绍 一种 简 单 的方 法一交替矩 阵 法 ， 其主 要 思想是 将功 阵经 过多次 选 得的 交替矩 阵左乘后 ， 变为一 个三角阵 。 而 单 位阵 ， 说 明参数 估 计的误 差平 方和 是不变的 。 阵的 一般形 式为