swh:1:snp:81eadaa089e8253d8469bcef66aa332632c6c669
Tip revision: fdf16693b000f3e309c56091892c61f0ec9fd670 authored by Hans W. Borchers on 21 November 2017, 16:15:00 UTC
version 2.1.1
version 2.1.1
Tip revision: fdf1669
hankel.R
##
## h a n k e l . R
##
hankel <- function(a, b) {
if (!is.vector(a))
stop("Argument 'a' must be a numeric or complex vector.")
n <- length(a)
if (missing(b)) b <- c(a[n], rep(0, n-1))
if (!is.vector(b))
stop("Argument 'b' must be a numeric or complex vector.")
m <- length(b)
if (a[n] != b[1])
warning("a[n] not equal to b[1], b[1] set to a[n].")
H <- matrix(0, n, m)
for (i in 2:(n+m))
H[row(H)+col(H) == i] <- if (i <= n+1) a[i-1] else b[i-n]
return(H)
}