这个的例子函数是y=exp(-Ts*n)的自身卷积，源程序只需要修改为相应函数。矩阵的算法相同。