swh:1:snp:81eadaa089e8253d8469bcef66aa332632c6c669
Tip revision: 113add4c7a9bfd37466a9b41cf2dc060b7b87091 authored by HwB on 06 December 2012, 00:00:00 UTC
version 1.3.1
version 1.3.1
Tip revision: 113add4
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)
}