R/getStreakBacksweep.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getStreakBacksweep" <- function(streakT, k, mu, x) {
  Y <- matrix(0, length(x), length(k))
  for (i in 1:length(k)) {
    x1 <- exp(-k[i] * (x - mu + streakT))
    x2 <- exp(-k[i] * ((streakT / 2) - (x - mu)))
    x3 <- exp(-k[i] * streakT)
    Y[, i] <- (x1 + x2) / (1 - x3)
  }
  Y
}