phase.partnered <-
function (n=2000, rho=1, gamma=1, sigma=0.10, mu=0) {
ts.part=matrix(nrow=n, ncol=2, NA)
minVal=0
maxVal=2*pi
fs=seq(from=1, to=n/2, by=1)
epsilon=runif(n=n/2)
epsilon [epsilon > 0.5]=1
epsilon [epsilon < 0.5]=-1
delta=acos(rho)*epsilon
phi1=minVal+(maxVal-minVal)*runif(n=n/2)
phi2=phi1+delta
const.num=1/(fs^(gamma/2))
for (t in 1:n) {
ts.part[t,1]=sum(const.num*(sin(2*pi*fs*t/n + phi1)))
ts.part[t,2]=sum(const.num*(sin(2*pi*fs*t/n + phi2)))
}
ts.part=mu+(ts.part/apply(ts.part, 2, sd))*sigma
ts.part=as.data.frame(ts.part)
colnames(ts.part) <- c("timeseries1", "timeseries2")
results=list(rho=rho, gamma=gamma, sigma=sigma, mu=mu, timeseries=ts.part)
class(results)="phasepart"
return (results)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.