toeplitz.mata

来自「是一个经济学管理应用软件 很难找的 但是经济学学生又必须用到」· MATA 代码 · 共 26 行

*! version 1.0.0  10sep2004
version 9.0

mata:

numeric matrix Toeplitz(numeric colvector fc, numeric rowvector fr)
{
	real scalar	r, c, i, j
	numeric matrix	T

	r = rows(fc)
	c = cols(fr)
	T = J(r, c, iscomplex(fc) | iscomplex(fr) ? 0i : 0)
	if (r==0 | c==0) return(T)

	T[1,.] = fr
	T[.,1] = fc

	for (i=2; i<=r; i++) {
		for (j=2; j<=c; j++) T[i,j] = T[i-1, j-1] 
	}
	return(T)
}

end