source("fit.s")
source("fit.s")
source("fit.s")
source("fit.s")
source("fit.s")
?qqnorm
?qqnorm
ls()
aa
residuals(aa)
rr<-residuals(aa)
rr
rm(rr)
ls()
xxx
l
l
residuals
?residuals
xxx[,2]-fitted.values(aa)
(xxx[,2]-fitted.values(aa)) - residuals(aa)
history()
qqnorm((xxx[,2]-fitted.values(aa)))
qqnorm(residuals(aa))
readline()
graphics.off()
source("fit.s")
source("fit.s")
source("fit.s")
source("fit.s")
graphics.off()
source("fit.s")
graphics.off()
source("fit.s")
?fft
l
ls()
plot(fft(xxx))
plot(Mod(fft(xxx)))
xxx[1:10,]
plot(Mod(fft(xxx[,2])))
length(xxx)
length(xxx[,2])
plot(Arg(fft(xxx[,2])))
plot(Mod(fft(xxx[,2])))
plot(Mod(fft(xxx[,2])),type=l)
?plot
plot(Mod(fft(xxx[,2])),type="b")
plot(Mod(fft(xxx[,2])),type="l")
plot(Mod(fft(xxx[1:10,2])),type="l")
plot(Mod(fft(xxx[1:50,2])),type="l")
plot(Mod(fft(xxx[,2]))[1:50],type="l")
?max
Mod(fft(xxx[,2]))[1:50]
plot(Mod(fft(xxx[,2])),type="l")
length(Mod(fft(xxx[,2])))
aa
system("vi fit.s")
length(xxx[,2])
xxx[1-10,]
xxx[1-10]
xxx[1-10,,]
xxx[1-10,1:2]
xxx[1,1:2]
length(xxx[,2])
size(xxx[,2])
help.search("matrix")
?matrix
is.matrix(xxx)
is.array(xxx)
is.table(xxx)
plot(Arg(fft(xxx[,2])),type="l")
plot(Arg(fft(xxx[,2]))[1:50],type="l")
plot(Mod(fft(xxx[,2]))[1:50],type="l")
plot(Re(fft(xxx[,2]))[1:50],type="l")
plot(Im(fft(xxx[,2]))[1:50],type="l")
graphics.off()
graphics.off()
graphics.off()
q()
ls()
summary (aa)
l
q()
options.ps()
ps.options()
ls()
fit-ps.s
source("fit-ps.s")
?postscript
postscript("scr.ps")
source("fit-ps.s")
source("fit-ps.s")
q()
postscript("scr.ps")
source("fit-ps.s")
graphics.off()
?sink
unlink("scr.txt")
?capture.output
capture(summary(aa),file="scr.txt")
capture.output(summary(aa),file="scr.txt")
close("scr.txt")
?close
close("scr.txt",type="rw")
close()
?postscript
?plot
l
l
ps.options()
source("fit-ps.s")
graphics.off()
source("fit-ps.s")
library()
postscript
l
graphics.off()
postscript("scr.ps",horiz=TRUE)
source("fit-ps.s")
graphics.off()
?nls
qqnorm
?qqnorm
ls()
Arg(fft(xxx[,2]))[1:50],type="l")
Arg(fft(xxx[,2]))[1:50]
-1.6692990/(2*pi)
x11()
plot(Arg(fft(xxx[,2]))[1:50],type="l")
plot(Im(fft(xxx[,2]))[1:50],type="l")
plot(Mod(fft(xxx[,2]))[1:50],type="l")
Mod(fft(xxx[,2]))[1:50]
sqrt(17597.82852)
Re(fft(xxx[,2]))[1:50]
sqrt(1730)
Im(fft(xxx[,2]))[1:50]
2048/1024
2*16
2^17
2^16
l
17598/1024
17598/512
?fft
?Mod
?Arg
 
h
history()
57.29578/pi
57.29578*pi
180/pi
Mod(fft(xxx[,2]))[1:20]
Arg(fft(xxx[,2]))[1:20]
-1.6692990/(2*pi)
(2 * pi)/128
cat("A = ", Mod(fft(xxx[,2]))[17]/512, " ;   B = ", Arg(fft(xxx[,2]))[17]/(2*pi))
source("fit.s")
source("fit.s")
?exit
help.search("exit"
)
graphics.off()
source("fit.s")
li -rt
summary(aa)
rm(aa,xx)
rm(xxx)
ls()
source("scr.s)
source("scr.s")
graphics.off()
source("scr.s")
source("scr.s")
source("scr.s")
q()
source("scr.s")
source("scra.s")
ls()
summary(aa)
options()
source("scra.s")
source("scra.s")
?summary
library()
library(base)
source("scra.s")
source("scra.s")
?library
.packages
.packages()
.packages(T)
library(akima)
.packages()
library()
?library
?cat
?summary
source("scra.s")
 aa
summary.lm(aa)
summary.nls(aa)
help.search(summary)
help.search("summary")
summary.nlsList9aa)
summary.nlsList(aa)
?summary.nlsList
help.search(summary.nlsList)
help.search(nls)
help.search("summary.nlsList")
summary.nlsList
aa$
$
aa
list(aa)
nls
?nls
nlsModel
?nlsModel
aa$resid()
aa$resid
aa
aa$dev
aa$model
aa$
$
aa
aa$$data
help("list")
list(aa)
alist(aa)
is.list(aa)
as.list.environment(aa)
as.list.environment(aa,all.names=F)
summary(aa)
l
ls()
last.warning
xxx
rm(xxx,aa)
aa
source("scra.s")
sink
sink()
source("scr.s")
ls()
source("scra.s")
source("scra.s")
source("scra.s")
?print
source("scra.s")
source("scra.s")
source("scr.s")
source("scr.s")
source("scr.s")
Arg(fft(xxx[,2]))[1:20]
plot(Arg(fft(xxx[,2])))
plot(Arg(fft(xxx[,2])),type="l")
mean(Arg(fft(xxx[,2])))
plot((fft(xxx[,2])),type="l")
plot(Mod(fft(xxx[,2])),type="l")
plot(Im(fft(xxx[,2])),type="l")
plot(Arg(fft(xxx[,2])),type="l")
plot(Re(fft(xxx[,2])),type="l")
plot(Re(fft(xxx[,2]))[1:20],type="l")
plot(Im(fft(xxx[,2]))[1:20],type="l")
?tan
Im(fft(xxx[,2]))[17]
Re(fft(xxx[,2]))[17]
atan(-17512.52/-1730.632)
pi/2
atan(-1730.632/-17512.52)
source("scr.s")
source("scr.s")
source("scr.s")
source("scr.s")
source("scr.s")
plot(Arg(fft(xxx[,2])),type="l")
mean(Arg(fft(xxx[,2])))
mean(Arg(fft(xxx[,2]))[1-20])
mean(Arg(fft(xxx[,2]))[15-20])
mean(Arg(fft(xxx[,2]))[15-19])
mean(Arg(fft(xxx[,2]))[17])
m11ean(Arg(fft(xxx[,2]))[17])
plot(Arg(fft(xxx[,2]))[1-30],type="l")
graphics.off()
plot(Arg(fft(xxx[,2]))[1:30],type="l")
lines(Mod(fft(xxx[,2]))[1:30],type="l")
plot(Mod(fft(xxx[,2]))[1:30],type="l")
xxx
xxx[(seq(1,1024)%128)]
?%
?mod
help.search("modulo")
help.search("mod")
?Syntax
?Arithmetic
xxx[(seq(1,1024)%%128)]
xxx[(seq(1,1024)%%128),2]
seq(1,1024)%%128
?seq
seq(1,1024)%%128
seq(1,1024)%%128==0
xxx[seq(1,1024)%%128==0,2]
xxx[seq(1,1024)%%128==0,]
xxx[seq(1,1024)%%64==0,]
xxx[seq(1,1024)%%64 +1==0,]
xxx[seq(1,1024)%%64==0 + 1,]
xxx[seq(1,1024)%%64==0 + 2,]
xxx[seq(1,1024)%%64==0 + 3,]
aa
?Mod
Im(fft(xxx[,2]))[17]/Re(fft(xxx[,2]))[17]
?atan
atan(Im(fft(xxx[,2]))[17]/Re(fft(xxx[,2]))[17])
atan(Im(fft(xxx[,2]))[17]/Re(fft(xxx[,2]))[17])rg
(Arg(fft(xxx[,2]))[17])
tan(Arg(fft(xxx[,2]))[17])
atan(Arg(fft(xxx[,2]))[17])
atan(tan(Arg(fft(xxx[,2]))[17]))
atan(-tan(Arg(fft(xxx[,2]))[17]))
tan(-Arg(fft(xxx[,2]))[17])
acos(Im(fft(xxx[,2]))[17]/Mod(fft(xxx[,2]))[17])rg
?acos
acos(Im(fft(xxx[,2]))[17]/Mod(fft(xxx[,2]))[17])
acos(Re(fft(xxx[,2]))[17]/Mod(fft(xxx[,2]))[17])
asin(Im(fft(xxx[,2]))[17]/Mod(fft(xxx[,2]))[17])
tan(asin(Im(fft(xxx[,2]))[17]/Mod(fft(xxx[,2]))[17]))
tan(-Arg(fft(xxx[,2]))[17])
Arg(fft(xxx[,2]))[17]
atan2(Im(fft(xxx[,2]))[17],Re(fft(xxx[,2]))[17])
atan2(Im(fft(xxx[,2]))[17],Re(fft(xxx[,2]))[17])/(2*pi)
pi/2
atan2(Im(fft(xxx[,2]))[17],Re(fft(xxx[,2]))[17]) - pi/2
atan2(Im(fft(xxx[,2]))[17],Re(fft(xxx[,2]))[17]) + pi/2
source("scr.s")
source("scr.s")
source("scr.s")
source("fit.s")
source("fit.s")
graphics.off()
source("fit-output.s")
source("fit-output.s")
source("fit-output.s")
q()
z=(sin(s)+ 1i*cos(x))
z=(sin(x)+ 1i*cos(x))
z=function(x){aa=(sin(x)+ 1i*cos(x));aa}
z(1)
vi(z)
z=function(x=1,d=0){aa=(sin(x+d)+ 1i*cos(x+d));aa}
z()
z(2)
z(2,2)
z(2,-2)
z(0)
z(seq(0,1024))
fft(z(seq(0,1024)),inverse=T)
plot(Re(fft(z(seq(0,1024)),inverse=T)))
plot(Mod(fft(z(seq(0,1024)),inverse=T)))
plot(Mod(fft(z(seq(0,1024)))))
plot((z(seq(0,1024))))
plot((z(seq(0,1024),1)))
plot(Mod(fft(z(seq(0,1024)))))
plot(Mod(fft(z(seq(0,1024)))),type="l")
plot((fft(z(seq(0,1024)))),type="l")
plot((fft(z(seq(0,1024),d=1))),type="l")
plot((fft(z(seq(0,1024)))),type="l")
plot((fft(z(seq(0,1024),d=1))),type="l")
plot(Mod(fft(z(seq(0,1024),d=1))),type="l")
plot(Mod(fft(z(seq(0,1024),d=0))),type="l")
plot(Arg(fft(z(seq(0,1024),d=0))),type="l")
plot(Arg(fft(z(seq(0,1024),d=1))),type="l")
plot(Arg(fft(z(seq(0,1024),d=pi))),type="l")
plot(Arg(fft(z(seq(0,1024),d=2*pi))),type="l")
plot(Arg(fft(z(seq(0,1024),d=pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=3*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=6*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=5*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=2*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=2*pi))),type="l")
plot(Arg(fft(z(seq(0,1024),d=4*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=0*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=0*pi/2))),type="l")
plot(Arg(fft(z(seq(0,1024),d=0*pi/2),inverse=T)),type="l")
plot(Mod(fft(z(seq(0,1024),d=0*pi/2),inverse=T)),type="l")
seq(0,1024)
z
z1=function(x=1,d=0){aa=(sin(x+d)+ 1i*cos(x));aa}
plot(Mod(fft(z1(seq(0,1024),d=0*pi/2),inverse=T)),type="l")
plot(Mod(fft(z1(seq(0,1024),d=1*pi/2),inverse=T)),type="l")
plot(Mod(fft(z1(seq(0,1024),d=2*pi/2),inverse=T)),type="l")
plot(Mod(fft(z1(seq(0,1024),d=3*pi/2),inverse=T)),type="l")
plot(Mod(fft(z1(seq(0,1024),d=.1*pi/2),inverse=T)),type="l")
plot(Im(fft(z1(seq(0,1024),d=.1*pi/2),inverse=T)),type="l")
plot(Im(fft(z1(seq(0,1024),d=.1*pi/2))),type="l")
z1=function(x=1,d=0){aa=(sin(x+d)+ cos(x));aa}
plot(Mod(fft(z1(seq(0,1024),d=1*pi/2))),type="l")
plot(Mod((z1(seq(0,1024),d=1*pi/2))),type="l")
plot(Mod((z1(seq(0,1024),d=1*pi/2)))[1:10],type="l")
plot(Mod((z1(seq(0,1024),d=1*pi/2)))[1:20],type="l")
plot(Mod((z1(seq(0,1024),d=2*pi/2)))[1:20],type="l")
plot(Mod((z1(seq(0,1024),d=0*pi/2)))[1:20],type="l")
plot(((z1(seq(0,1024),d=0*pi/2)))[1:20],type="l")
plot(((z1(seq(0,1024),d=1*pi/2)))[1:20],type="l")
plot(((z1(seq(0,1024),d=2*pi/2)))[1:20],type="l")
plot(((z1(seq(0,1024),d=0*pi/2)))[1:20],type="l")
plot(((z1(seq(0,1024),d=1*pi/2)))[1:20],type="l")
plot(((z1(seq(0,1024),d=0*pi/2)))[1:20],type="l")
plot(((2*pi/16*z1(seq(0,1024),d=0*pi/2)))[1:20],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=0*pi/2)))[1:20],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=0*pi/2)))[1:256],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=0*pi/2)))[0:256],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=0*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=1*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=2*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-2*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-1*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=1*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=.1*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.1*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.5*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.9*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.99*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-1*pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-pi/2)))[1:64],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:64],type="l")
z1
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:64],type="l")
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:1024],type="l")
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:100],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:100],type="l")
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:100],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:100],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[60:70],type="l")
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[60:70],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[60:70],type="l")
z1
z1=function(x=1,d=0){aa=(sin(x+d));aa}
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[60:70],type="l")
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:1024],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:1024],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.9999*pi/2)))[1:128],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-0*pi/2)))[1:128],type="l")
plot(Mod(fft(z1(2*pi/16*seq(0,1024),d=-0*pi/2)))[1:128],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-0*pi/2)))[1:128],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-0*pi/2))),type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.999*pi/2))),type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.999*pi/2))),type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.999*pi/2)))[1:100],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-0*pi/2)))[1:128],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-0*pi/2)))[1:128],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-0.999*pi/2)))[1:128],type="l")
z1
z
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.1*pi/2)))[1:128],type="l")
plot(((z1(2*pi/16*seq(0,1024),d=-.1*pi/2)))[1:128],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.1*pi/2)))[1:128],type="l")
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.15*pi/2)))[1:128],type="l")
Arg(fft(z1(2*pi/16*seq(0,1024),d=-.15*pi/2)))[1:128]
pi/2
plot(Arg(fft(z1(2*pi/16*seq(0,1024),d=-.12*pi/2)))[1:128],type="l")
Arg(fft(z1(2*pi/16*seq(0,1024),d=-.12*pi/2)))[1:128]
Arg(fft(z1(2*pi/64*seq(0,1024),d=-.12*pi/2)))[1:128]
plot(Arg(fft(z1(2*pi/64*seq(0,1024),d=-.12*pi/2)))[1:128],type="l")
Arg(fft(z1(2*pi/64*seq(0,1024),d=-.12*pi/2)))[1:128]
pi/2-1.710432
plot(Arg(fft(z1(2*pi/64*seq(0,1024),d=-.12*pi/2)))[1:1024],type="l")
q()
source(fit.s)
source("fit.s")
q()
