# Set up the harmonic acceleration operator
dayrange <- c(0,365)
Lbasis <- create.constant.basis(dayrange,
axes=list("axesIntervals"))
Lcoef <- matrix(c(0,(2*pi/365)^2,0),1,3)
bfdobj <- fd(Lcoef,Lbasis)
bwtlist <- fd2list(bfdobj)
harmaccelLfd <- Lfd(3, bwtlist)
oldpar <- par(no.readonly=TRUE)
plot(harmaccelLfd)
par(oldpar)
Run the code above in your browser using DataLab