R/nntsrandominitial.R

In CircNNTSR: Statistical Analysis of Circular Data using Nonnegative Trigonometric Sums (NNTS) Models

nntsrandominitial <- function(M=1) 
{
    if (M<0)
       return("M must be nonnegative")
    res <- complex(M + 1)
    aux <- rnorm(2*(M + 1))
    aux <- sqrt(1/(2 * pi)) * (aux/sqrt(sum(aux^2)))
    res[1]<-Mod(aux[1]+1i*aux[M+2])
    if (M>0){
        for (k in 2:(M+1)){
	    res[k] <- aux[k] + 1i*aux[k+M+1]
        }   
    } else
        res<-1/(2*pi)
    return(res)

}