6.遗忘因子法递推算法估计参数.docVIP

6.遗忘因子法递推算法估计参数 考虑一个时变系统: y(k)＋a(k)y(k–1)＝b(k)u(k–1)＋e(k) a(k), b(k)有下列各值: 0300; k300 e(k)是零均值，方差为0.1的高斯分布白噪声。经数字仿真，采用遗忘因子法的递推算法估计参数。 遗忘因子法递推算法为： 循环900次，遗忘因子值可以随意调整。参数变化曲线如下图所示。 α=0.99时的情况： α=0.9时的情况： 源程序： nclude math.h #include stdio.h #include brmul.c #include brinv.c int main() { FILE *fp1,*fp2,*fp; int u[900];double v[900]; int i, j, k; double y[910],y2[1][1]; double t1[2][1],t2[1][2],t3[1][1],t5[2][2],t6[2][2]; double o[2][1], o1[900][2], w[2][1], w1[1][2]; double p[2][2], km[2][1]; double alpha;//遗忘因子 t3[0][0] = 0.0; y2[0][0] = 0.0; alpha = 0.99;