Pendulum Experiment , R

Language: Layout:

data <- data.frame(
                    l = c(1.5,1.25,1,0.75,0.5,0.25),  
                   T1 = c(24.97,22.88,20.15,18.15,14.85,11),  
                   T2 = c(25.03,22.47,20.54,18.25,15.25,11.15),
                   T3 = c(25.13,22.63,20.41,18,15.16,11.03)
                  )
data
data1 <- data.frame(
                    l = c(1.5,1.5,1.5, 1.25,1.25,1.25, 1,1,1,  0.75,0.75,0.75,   0.5,0.5,0.5,  0.25,0.25,0.25),  
                   T10 = c(24.97,25.03,25.13,  22.88,22.47,22.63,      20.15,20.54,20.41,   18.15,18.25,18,    14.85,15.25,15.16,     11,11.15,11.03)
                    )
data1

data2 <- data.frame(
                    l = c(1.5,1.5,1.5, 1.25,1.25,1.25, 1,1,1,  0.75,0.75,0.75,   0.5,0.5,0.5,  0.25,0.25,0.25),  
                   T10 = c(24.97,25.03,25.13,  22.88,22.47,22.63,      20.15,20.54,20.41,   18.15,18.25,18,    14.85,15.25,15.16,     11,11.15,11.03)
                    )
data2 <- data.frame(
                    a = c(15,15,15, 25,25,25, 35,35,35, 45,45,45, 55,55,55, 65,65,65),  
                   T5 = c(10.18,10.25,10, 10.37,10.57,10.34, 10.34,10.32,10.63, 10.75,10.5,10.34,10.88,10.56,10.75, 10.75,10.69,10.69)
                    )

str(data2)
data3 <- data.frame(
                    a = data2$a ,
                   T = (data2$T5)/5
                    )
data3

str(data3)

data4 <- data.frame(
                    a = data2$a ,
                   g = ((2*pi/(data3$T))^2)*1
                    )
data4

str(data3)

M = tapply(data4$g, 
           INDEX = data4$a, 
           FUN   = mean)

boxplot(g~a,
data=data4,
        ylab="Aceleração da gravidade (m/s^2)",
        xlab="Ângulo (°)")
abline(h=9.80665, col ="red", lty = 2)
points(M, 
       col="red",
       pch="+",
       cex=2)
lines(M, 
       col="blue",
       cex=2)
legend("bottomleft", legend=c("g= 9,80665 m/s^2", "Média"),
       lty=c(2,NA), col=c("red","red"), lwd=2:2, bty="o", pch=c(NA,"+"))

data5 <- data.frame(
                    l = data1$l ,
                    T= data1$T10/10 
                    )
data5

data7 <- data.frame(
                    l = data1$l ,
                   g = ((2*pi/(data5$T))^2)*data1$l
                    )
data7
M1 = tapply(data7$g, 
           INDEX = data7$l, 
           FUN   = mean)
boxplot(g~l,
data=data7,
        ylab="Aceleração da gravidade (m/s^2)",
        xlab="Comprimento (m)")
abline(h=9.80665, col ="red", lty = 2)
points(M1, 
       col="red",
       pch="+",
       cex=2)
lines(M1, 
       col="blue",
       cex=2)
legend("bottomright", legend=c("g= 9,80665 m/s^2", "Média"),
       lty=c(2,NA), col=c("red","red"), lwd=2:2, bty="o", pch=c(NA,"+"))

data6 <- aggregate(data5[, 2], list(data5$l), mean)

plot(data5$l, data5$T, ylab="Período (s)",
        xlab="Comprimento (m)")
points(data6$Group.1, data6$x,col="red",
       pch="+",
       cex=2)
curve(  2*pi*sqrt(x/9.8), from=0.25, to=2, n=300, add=TRUE, xlab="xvalue", ylab="yvalue", 
             col="blue", lwd=2, main="Plot of (3x^2 + x)"  )
lines(data6$Group.1, data6$x, 
       col="red",
       cex=2)
legend("topleft", legend=c("eq. 8.2", "Média"),
       lty=c(1,NA), col=c("blue","red"), lwd=2:2, bty="n", pch=c(NA,"+"))

[ + ] Show input
Absolute running time: 1.14 sec, cpu time: 0.95 sec, memory peak: 34 Mb, absolute service time: 1,18 sec