Re: Projectielbaan
Geplaatst: zo 24 okt 2021, 17:02
in deze is een anlytische berekening altijd benaderend maar met een hoek van 40 graden wel in overeenstemming ongeveer met een numerieke.
in deze is een anlytische berekening altijd benaderend maar met een hoek van 40 graden wel in overeenstemming ongeveer met een numerieke.
Code: Selecteer alles
f = @(t,y) [y(2); -9.81*0.0006*y(2)*(y(2)^2+y(4)^2)^0.5;y(4);-9.81*0.0006*y(4)*(y(2)^2+y(4)^2)^0.5 - 9.81 ];
tspan = [0: 0.01: 15];
yinit = [0, 50*cosd(40), 0, 50*sind(40)];
[t1,y1]=ode45(f, tspan, yinit)
plot(y1(:,1),y1(:,3))
legend()
hold on
f = @(t,y) [y(2); -9.81*0.0006*y(2)^2;y(4);-9.81*0.0006*y(4)*y(2) - 9.81 ];
tspan = [0: 0.01: 15];
yinit = [0, 50*cosd(40), 0, 50*sind(40)];
[t2,y2]=ode45(f, tspan, yinit)
plot(y2(:,1),y2(:,3))
legend('wnvl1','Rik Speybrouck')
mocht ik er 25 naast zitten zou dit 164.87 geven en jullie zitten er dan nog eens een heel eind onder
Het scheelt 25 % in wrijving bij een hoek van 40 graden. Dat is niet meteen een verschil van 25 % in afstanden, tijden of hoogtes.Rik Speybrouck schreef: ↑zo 24 okt 2021, 17:21 mocht ik er 25 naast zitten zou dit 164.87 geven en jullie zitten er dan nog eens een heel eind onder
die formules zijn die numeriek berekend of via een dif vergelijking.Xilvo schreef: ↑zo 24 okt 2021, 17:25Het scheelt 25 % in wrijving bij een hoek van 40 graden. Dat is niet meteen een verschil van 25 % in afstanden, tijden of hoogtes.Rik Speybrouck schreef: ↑zo 24 okt 2021, 17:21 mocht ik er 25 naast zitten zou dit 164.87 geven en jullie zitten er dan nog eens een heel eind onder
In de grafiek van wnvl1 zie je dat in het begin een afwijking ontstaat. Maar de hoek met de horizontaal wordt ook kleiner, dan wordt het verschil in wrijving ook kleiner.
Pas als de dalingshoek toeneemt neemt het verschil weer verder toe.
Welke formules?Rik Speybrouck schreef: ↑zo 24 okt 2021, 17:29 die formules zijn die numeriek berekend of via een dif vergelijking.
van wvnlXilvo schreef: ↑zo 24 okt 2021, 17:30Welke formules?Rik Speybrouck schreef: ↑zo 24 okt 2021, 17:29 die formules zijn die numeriek berekend of via een dif vergelijking.
heb je al eens mijn numerieke benadering bekeken waarvoor ik een excel file on line heb gezet, ik heb ook een verwijzing gedaan naar de paper waarop deze is ontworpen en redelijk goed aansluit met mijn analytische berekening. Ik vind het vreem vinden maar ik denk dat jullie met een onderschatting zitten van de afstand.wnvl1 schreef: ↑zo 24 okt 2021, 17:18 het valt nog mee
vergelijking.jpg
Code: Selecteer alles
f = @(t,y) [y(2); -9.81*0.0006*y(2)*(y(2)^2+y(4)^2)^0.5;y(4);-9.81*0.0006*y(4)*(y(2)^2+y(4)^2)^0.5 - 9.81 ]; tspan = [0: 0.01: 15]; yinit = [0, 50*cosd(40), 0, 50*sind(40)]; [t1,y1]=ode45(f, tspan, yinit) plot(y1(:,1),y1(:,3)) legend() hold on f = @(t,y) [y(2); -9.81*0.0006*y(2)^2;y(4);-9.81*0.0006*y(4)*y(2) - 9.81 ]; tspan = [0: 0.01: 15]; yinit = [0, 50*cosd(40), 0, 50*sind(40)]; [t2,y2]=ode45(f, tspan, yinit) plot(y2(:,1),y2(:,3)) legend('wnvl1','Rik Speybrouck')
je moet eens kijken welke k waarde wordt voorgesteld in het paper degene die ik gebruik en die ook gebruikt wordt in het boek waarvan copy reeds on line gezet.Xilvo schreef: ↑zo 24 okt 2021, 18:024559-simulating-projectile-motion-in-the-air-with-spreadsheets.pdf
Daar staat weinig anders in dan we zelf al hebben bedacht.
Ik vermoed dat wanneer Rick in z'n Excel de definitie voor k van Ukster gebruikt hij ook op het goede resultaat komt.
Maar ik vrees dat we in kringetjes draaien. Ik vind het wel mooi zo.