rm(list=ls())

library(tuneR)
library(fastICA)

dataMatrix<-readWave("ICATest.wav")

#summary(dataMatrix)
plot(dataMatrix)
#play(dataMatrix)
W<-mono(dataMatrix,"left")
WdownSampled<-downsample(W,11025)  #Wdown is of class "Wave"

a <- fastICA(matrix(c(seq(1,140288),WdownSampled@left),ncol=2), 2, alg.typ = "parallel", fun = "logcosh", alpha = 1,method = "R", row.norm = FALSE, maxit = 200, tol = 0.001, verbose = TRUE)

aWave1<-Wave(a$X,samp.rate=11025,bit=16)
aWave2<-Wave(a$S,samp.rate=11025,bit=16)








WspecObject <- periodogram(Wdown, normalize = TRUE, width = 1024, overlap = 512)
plot(WspecObject)
summary(WspecObject)

WSpecObjectFF<-FF(WspecObject)
#plot(WSpecObjectFF,type="l")

notes <- noteFromFF(WSpecObjectFF, 440)
snotes <- smoother(notes)
melodyplot(WspecObject, snotes)
qnotes <- quantize(snotes, WspecObject@energy, parts = 8)
quantplot(qnotes, expected = rep(c(0, -12), each = 4), bars = 2)

a <- fastICA(matrix(c(seq(1,140288),Wdown@left),ncol=2), 2, alg.typ = "parallel", fun = "logcosh", alpha = 1,method = "R", row.norm = FALSE, maxit = 200, tol = 0.001, verbose = TRUE)
plot(a$S,pch=".",col=2)





Extra:
Wwave<-extractWave(Wdown)
WwaveCC <- melfcc(Wwave, frames_in_rows=FALSE)
WwaveCClifter<-lifter(WwaveCC,lift=0.6,inv=FALSE,htk=FALSE)
WwaveCCHz<-mel2hz(WwaveCClifter)
WwaveCCcep<-lpc2cep(WwaveCC)