ls()
q()
?dget
dget(yout184-dget.Data)
dget("yout184-dget.Data")
ls()
yout184<-dget("yout184-dget.Data")
source("yout184.s")
history()
t<-seq(0,16384)*2
plot(yout184,t,type="l")
length(t)
length(yout184)
t<-seq(0,16383)*2
t<-seq(1,16384)*2
plot(yout184,t,type="l")
plot(t,yout184,type="l")
#time<-seq(1,16384)*2
time<-seq(1,16384)*2
rm(t)
ls()
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
plot(time,yout184,type="l",ylab="chi-2 (degrees)",xlab="time (fs)")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
graphics.off()
source("scr.s")
source("scr.s")
q()
ls()
aa2.nls<-nls(Ct ~ Ct[1]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1, w=20))
summary (aa2.nls)
summary (aa.nls)
Ct[1]
ls()
aa2.nls<-nls(Ct ~ aa[1]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1, w=20))
summary (aa2.nls)
a.df[1,"Ct"]
aa.df[1,"Ct"]
aa[1]
aa2.nls<-nls(Ct ~ aa.df[1,"Ct"]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1, w=20))
summary (aa2.nls)
plot(time[1:500], aa.df[1:500,2],type="l",ylab="C(t)",xlab="time (fs)")
lines(time[1:500], fitted.values(aa2.nls),type="l",lty=2)
lines(time[1:500], fitted.values(aa.nls),type="l",lty=6)
history()
q()
19.48929*(pi/180)^2
sqrt(19.48929*(pi/180)^2)
.084*180/pi
1.98*100/1000/0.005936777
1.98*100/1000/(.084)^2
source("yout184.s")
plot(Mod(a[1:750]),type="l",ylab="Mod(fft(chi-2))",xlab="frequency bin")
summary (aa2.nls)
summary (aa3.nls)
summary (aa2.nls)
summary (aa.nls)
17.3*32.7
525/32.7
225/32.7
407/17.1
2^10
2^14
2^15
plot(seq(0,749)*1/32768,Mod(a[1:750]),type="l",ylab="Mod(fft(chi-2))",xlab="frequency bin")
plot(seq(0,749)*1/32.768,Mod(a[1:750]),type="l",ylab="Mod(fft(chi-2))",xlab="frequency bin")
plot(seq(0,749)*1/32.768,Mod(a[1:750]),type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
source("yout184.s")
graphics.off()
source("yout184.s")
source("scr.s")
f<-c(seq(0,8192)*1/32.768,seq(8191,1)*1/32.768)
length(f)
Mod(aa[1])
Mod(aa[16384])
Mod(aa[2])
plot(Mod(a),type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
plot(f,Mod(a),type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
plot(f[1:8192],Mod(a)[1:8192],type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
plot(f,Mod(a),type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
plot(f[1:8192],Mod(a)[1:8192],type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
33*c(7,16,30,75,130,25,40,60)
f<-c(seq(0,8192)*1/32.768)
plot(f[1:8193],Mod(a)[1:8193],type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
plot(f[1:750],Mod(a)[1:750],type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
Mod(a[1])
Arg(a[1])
Re(a[1])
Im(a[1])
plot(f[1:8193],Mod(a)[1:8193],type="l",ylab="Mod(fft(chi-2))",xlab="frequency (1/ps)")
graphis.off()
graphics.off()
source("yout184.s")
source("scr.s")
1.98*100/1000*(180/pi)^2*1/14.8^2
(180/pi)^2*1/14.8^2
1/((180/pi)^2*1/14.8^2)
14.8*pi/180
14.8*(pi/180)^2
1/((180/pi)^2*1/14.8)
1.98*100/1000*(180/pi)^2*1/14.8
1.98*100/1000*(180/pi)^2*1/19.5
33*1000*4.184*10^7/(7.7*10^-15)
sqrt(33*1000*4.184*10^7/(7.7*10^-15))
sqrt(33.3*1000*4.184*10^7/(7.7*10^-15))
sqrt(44.3*1000*4.184*10^7/(7.7*10^-15))
w=8;aa2b.nls<-nls(Ct ~ aa.df[1,"Ct"]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1))
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
plot(time[1:500], aa.df[1:500,2],type="l",ylab="C(t)",xlab="time (fs)")
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
w=20;aa2b.nls<-nls(Ct ~ aa.df[1,"Ct"]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1))
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
w=12;aa2b.nls<-nls(Ct ~ aa.df[1,"Ct"]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1))
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
w=14;aa2b.nls<-nls(Ct ~ aa.df[1,"Ct"]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1))
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
plot(time[1:500], aa.df[1:500,2],type="l",ylab="C(t)",xlab="time (fs)")
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
w=16;aa2b.nls<-nls(Ct ~ aa.df[1,"Ct"]*exp(-time/tau)*(cos(w*time) + 1/(tau*w)*sin(w*time)),aa.df[1:500,],start=c(tau=.1))
lines(time[1:500], fitted.values(aa2b.nls),type="l",lty=6)
ls()
rm(ls())
ls()
xx<-ls()
rm(xx)
ls()
?remove
remove(xx)
ls
ls()
xx<ls()
xx<-ls()
xx
remove(xx)
ls()
?rm
remove(list=xx)
xx
xx<-ls()
remove(list=xx)
ls()
rm(xx)
source("yout184.s")
source("yout184.s")
source("scr1.s")
source("yout184.s")
source("yout184.s")
source("scr1.s")
vi *.txt
source("scr1.s")
aa3.nls[3]
aa3.nls
aa3.nls[,"C"]
list(aa3.nls)
nls
?nls
?nlsModel
aa3.nls$getpars()
ls()
aa3.nls$getPars()
aa3.nls$getPars
aa3.nls$
$
list(aa3.nls)
a3.nls[[1]]
aa3.nls[[1]]
aa3.nls[[1]]$getPars()
aa3.nls$getPars()
aa3.nls$$getPars()
aa3.nls[[1]]$getPars()
aa3.nls[[1]]$getPars()[3]
fitted.values(aa3.nls)[1]
source("scr1.s")
source("scr1.s")
source("yout184.s")
source("yout184.s")
graphics.off()
33*c(7,16,30,75,130,25,40,60)
I=7.7*10^15;XX=1.98*100*(180/pi)^2*4.184*10^7/(7.7*10^-15);aa2c.nls<-nls(Ct ~ C*exp(-time/tau)*(cos(sqrt(XX/C - 1/tau^2)*time) + 1/(tau*sqrt(XX/C - 1/tau^2))*sin(sqrt(XX/C - 1/tau^2)*time)),aa.df[1:500,],start=c(tau=.1, C=20))
I=7.7*10^15;XX=1.98*100*(180/pi)^2*4.184*10^7/(7.7*10^-15)
XX
XX/15
sqrt(XX/15)
sqrt(XX/14.85)
?nls
aa2c.nls<-nls(Ct ~ C*exp(-time/tau)*(cos(sqrt(XX/C - 1/tau^2)*time) + 1/(tau*sqrt(XX/C - 1/tau^2))*sin(sqrt(XX/C - 1/tau^2)*time)),aa.df[1:500,],start=c(tau=.1, C=15), trace=TRUE)
aa2c.nls<-nls(Ct ~ C*exp(-time/tau)*(cos(sqrt(XX/C - 1/tau^2)*time) + 1/(tau*sqrt(XX/C - 1/tau^2))*sin(sqrt(XX/C - 1/tau^2)*time)),aa.df[1:500,],start=c(tau=.1, C=12), trace=TRUE)
aa2c.nls<-nls(Ct ~ C*exp(-time/tau)*(cos(sqrt(XX/C - 1/tau^2)*time) + 1/(tau*sqrt(XX/C - 1/tau^2))*sin(sqrt(XX/C - 1/tau^2)*time)),aa.df[1:500,],start=c(tau=.1, C=20), trace=TRUE)
aa2c.nls<-nls(Ct ~ C*exp(-time/tau)*(cos(sqrt(XX/C - 1/tau^2)*time) + 1/(tau*sqrt(XX/C - 1/tau^2))*sin(sqrt(XX/C - 1/tau^2)*time)),aa.df[1:500,],start=c(tau=.15, C=20), trace=TRUE)
XX/15
XX<-XX/10^24
XX
XX/15
sqrt(XX/15)
sqrt(XX/14.85)
XX=1.98*100*(180/pi)^2*4.184*10^7/(7.7*10^-15)*10^-24
sqrt(XX/14.85)
aa2c.nls<-nls(Ct ~ C*exp(-time/tau)*(cos(sqrt(XX/C - 1/tau^2)*time) + 1/(tau*sqrt(XX/C - 1/tau^2))*sin(sqrt(XX/C - 1/tau^2)*time)),aa.df[1:500,],start=c(tau=.15, C=20), trace=TRUE)
summary(aa2c.nls
)
sqrt(XX/12.3)
sqrt(XX/12.3 - 1/(.186)^2)
plot(time[1:500], aa.df[1:500,2],type="l",ylab="C(t)",xlab="time (fs)")
history()
lines(time[1:500], fitted.values(aa2c.nls),type="l",lty=6)
history(-200)
history(100)
graphics.off()
source("yout184.s")
source("yout184.s")
source("yout184.s")
aafit.nls
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
aa2c.nls[1]
fitted.values(aa2c.nls)["tau"]
fitted.values(aa2c.nls)[1]
coef(aa2c.nls)[1]
coef(aa2c.nls)
coef(aa2c.nls)["tau"]
coef(aa2c.nls)["C"]
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
source("yout184.s")
1/tau^2
XX/C
XX/C-1/tau^2
sqrt(XX/C-1/tau^2)
sqrt(XX/C)
1.98*100/1000*(180/pi)^2*1/C
ls()
xx<-ls()
rm(xx)
ls()
xx<-ls()
rm(list=xx)
ls()
rm(xx)
ls()
source("yout184-out.s")
source("yout184.s")
7/19
q()
