﻿功能:求解对称Toeplitz矩阵方程组

格式:SolveEquationLevtoeplitz(A,B)

说明:A为n*n的对称Toeplitz矩阵,或者A为n*1(或者1*n)的矩阵(此时A里面的元素表示对称Toeplitz矩阵的第一列或第一行).B为n*1的矩阵.此求解的含义是 求A*X=B中的X

例子:

a1 =[ 1.00000000000000  2.00000000000000  3.00000000000000  4.00000000000000  5.00000000000000  6.00000000000000  ]
a2  =
[ 1.00000000000000  2.00000000000000  3.00000000000000  4.00000000000000  5.00000000000000  6.00000000000000
  2.00000000000000  1.00000000000000  2.00000000000000  3.00000000000000  4.00000000000000  5.00000000000000
  3.00000000000000  2.00000000000000  1.00000000000000  2.00000000000000  3.00000000000000  4.00000000000000
  4.00000000000000  3.00000000000000  2.00000000000000  1.00000000000000  2.00000000000000  3.00000000000000
  5.00000000000000  4.00000000000000  3.00000000000000  2.00000000000000  1.00000000000000  2.00000000000000
  6.00000000000000  5.00000000000000  4.00000000000000  3.00000000000000  2.00000000000000  1.00000000000000  ]
b =
[ 42.1644512294626
  9.86536732402880
  27.7128759434972
  13.9948676871112
  41.3691531128106
  83.4597765391040  ]

SolveEquationLevtoeplitz(a2,b)//执行我们的求解命令,其实同样执行SolveEquationLevtoeplitz(a1,b)返回一样的结果
ans  =
[  -7.1763828263907
  25.0732962624511
  -15.782758437927
  20.5461468410428
  7.35816900029693
  -12.072152586820  ]

//By 2012/3/22