一求解一般右端项的Toeplitz方程组及Toeplitz矩阵的逆.docVIP

一 求解一般右端项的Toeplitz方程组及 Toeplitz 矩阵的逆 １实验目的 　 ● 熟悉求解Yule-Walker方程组，一般右端项的Toeplitz方程组及Toeplitz矩阵的逆的求解步骤； 　● 通过实验体会方程组的性质对求解的影响； ２实验原理 1）求解一般右端项的Toeplitz方程组： ,其中是已知向量.通过求解来计算,可得； ,其中； 是k阶Yule-Walker方程组的解， 是k阶单位阵， 由确定，这里.这样，便可通过出发递推地求得方程组的解。 2）Toeplitz矩阵的逆(以为例) 通过求解一个5阶的Yule-Walker方程组得到，再利用: , 求出的最后一列和最后一行元素，然后再利用的广义对称性求出相应元素； 利用及(1)得到的结果，然后再利用的广义对称性，可得到第二层的全部元素； 再利用,可求得最中间的元素，则就求得了的全部元素。 ３数值算例 求解Yule-Walker方程组：（求解四阶Yule-Walker方程组R1[1]保存方程组的阶数） please input the flag 1 to 3: 1 please input the string R1[N]: 4 -1 3 8 R1[0]=5.000000 R1[1]=4.000000 R1[2]=-1.000000 R1[3]=3.000000 R1[4]=8.000000 结果 y1[0]=5.000000 y1[1]=-0.786753 y1[2]=1.675522 y1[3]=2.357687 y1[4]=-2.539950求解Toeplitz矩阵的逆（5阶） please input the flag 1 to 3: 2 please input the string R1[N]: 0.4 1.2 3.4 1.2 R1[0]=5.000000 R1[1]= 0.400000 R1[2]=- 1.2 00000 R1[3]= 3.4 00000 R1[4]= 1.20000 得到的结果为 T[0][0]= -0.089410 T[0][1]= -0.015410 T[0][2]= 0.119910 T[0][3]= 0.319410 T[0][4]= -0.111810 T[1][0]=-0.015410 T[1][1]= 1.137010 T[1][2]= -0.205710 T[1][3]= -1.434810 T[1][4]= 0.319410 T[2][0]= 0.119910 T[2][1]= -0.205710 T[2][2]= -1.732610 T[2][3]=- -0.205710 T[2][4]= 0.119910 T[3][0]= 0.319410 T[3][1]= -1.434810 T[3][2]= -0.205710 T[3][3]= 1.137010 T[3][4]= 0.015410 T[4][0]= -0.111810 T[4][1]= 0.319410 T[4][2]= 0.119910 T[4][3]= 0.015410 T[4][4]= -0.089410 ③一般右端项的Toeplitz方程组：（4阶右端项的Toeplitz方程组R1[1]保存方程组的阶数） please input the flag 1 to 3: 3 please input the string R1[N]: “矩阵系数” 3 4 1 6 R1[N] R1[0]=5.000000 R1[1]=3.000000 R1[2]=4.000000 R1[3]=1.000000 R1[4]=6.000000 please input the string b[N]: 6 7 2 4 b[N] “右端项” b[0]=5.000000 b[1]=6.000000 b[2]=7.000000 b[3]=2.000000 b[4]=4.000000 结果： x[0]=5.000000 x[1]=6.000000 x[2]=1.375000 x[3]=4.512821 x[4]=-0.100273 将上面结果与直接用Matlab计算的结果比较误差较大。 ４心得体会 　　本次试验我是用c语言完成的，实验过程中，,正是因为Toeplitz方程组本身具有的良好的性质:广义对称性.使得在编写程序的过程中节省了很多步骤，整过过程求解的过程即用到向量乘法，这只要两轮循环就可以实现，用c语言很容易实现的。,特别是在求逆的过程中,因为它的对称性,我们可以仅仅去求它的一个倒三角矩阵,这样在一个循环中用三条赋值语