Flugbahn | R Documentation |
Berechnet die Flugbahn eines Geschosses mit berücksichtigung des Luftwiderstandes. Es wird über den Zeitintervall "t" dv, und dx,dy berechnet.
Flugbahn(v0,t,angle, target_vector, m, k,nb_return_values)
v0 |
Geschosss-Geschwindigkeit Vektor |
t |
Zeit (interval der Berechnung) |
angle |
Schusswinkel |
target_vector |
target_vectorvektor |
m |
Masse des Geschosses |
k |
k-Wert = 0.5cwrho*A: cw ==> Strömungswiderstandskoeffizient, Widerstandsbeiwert, rho ==> Dichte (kg/m^3), A ==> Fläche (m^2); |
nb_return_values |
Anzahl Werte die zurückgegeben werden. Wenn Null gewählt wird dann werden alle brechneten Werte zurückgegenen. |
matrix
matrix t,v_x, v_y, s_x, s_y, distance (to target)
# target_vector
target_vector <- c(100, 20) # m
angle <- atan(target_vector[2]/target_vector[1]) #rad
# Zeitintervall
t <- 1e-5
# Geschoss
v0 = 55 # [m/s]
m <- 0.440 # [kg]
cw <- 0.18; # cw Wert des Geschosses (Nach Tabellen der Geschossform)
rho <- 1.293 # Dichte des Mediums (Luft) [kg/m^3]
A <- (0.22)^2*pi/4 #Projektions Flaeche in Flugrichtung [m^2]
k <- 0.5*cw*rho*A
Flugbahn_c <- Flugbahn(v0,t,angle,target_vector,m,k,1000)[,4:5]
colnames(Flugbahn_c) <- c("sx","sy")
Flugbahn_c|>
plot(type = "l")
#tidyvers
library(tidyverse)
Flugbahn_c|>
as_tibble()|>
ggplot(aes(sx,sy))+
geom_line()+
geom_point(data = tibble(x = target_vector[1],y = target_vector[2]),aes(x,y), color = "red")+
geom_text(aes(x = target_vector[1],y = target_vector[2]*1.1, label = "target_vector"))+
coord_fixed()+
labs(title = "Flugbahn Rcpp", subtitle = paste("v0:",v0,"\nt:",t))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.