5.递推最小二乘法程序.docVIP

最小二乘参数的递推算法 根据递推公式用循环实现递增。如果满足所要求的精度，则跳出循环。 递推公式如下： 源程序为：“ls.dsw ” 。输入信号M序列存在数据文件“data1.txt”中 ，均值为0,方差为0.1的高斯白噪声存在数据文件“data2.txt”中 。估计结果存在数据文件“data3.txt”中。 源程序： math.h #include stdio.h #include brmul.c #include brinv.c int main() { FILE *fp1,*fp2,*fp; int u[600]; double v[600]; int i, j, k, count; double y1[4][1],y[510],y2[1][1]; double t1[4][1],t2[1][4],t3[1][1],t5[4][4],t6[4][4]; double o[4][1], o1[500][4], w[4][1], w1[1][4]; double p[4][4], km[4][1], er[4][1]; double maxe; t3[0][0] = 0.0; y2[0][0] = 0.0; fp1 = fopen(data1.txt,r); fp2 = fopen(data2.txt,r); for( i = 0;