Need help with my code
5 次查看(过去 30 天)
显示 更早的评论
I am trying do an iterative loop that will input an increasing velocity into 3 different trajectories until the final trajectory gets into tolerance around the initial height. To put it simply, a ball is thrown at a wall, bounces off the wall and hits the ground and then bounces towards the initial throwing point. The loop needs to converge on the initial height plus/minus the tolerance by inputting an increasing velocity until the final height is reached, and then output the number of iterations and the plot of the trajectories. Here is the code I have so far. I am still a newby with loops so not sure how to proceed. The help is appreciated.
% All distances and heights measured in meters, time in seconds, velocity in m/s, angles in degrees and
%%acceleration in m/s^2
H0 = 8;
return_height_tolerance = .001;
theta_0=26;
g=9.81;
dAC=16;
eR=1;
%initial velocities in both directions
v0x=v0*cosd(theta_0) %Trig
v0y=v0*sind(theta_0) %Trig
%Point B = Maximmum height along trajectory path 1
H1=H0+v0y^2/(2*g)
%trajectory 1
H2=H0+dAC*tand(theta_0)-g/(2*v0x^2)*dAC^2 %Trajectory Eq.
%velocities in x and y direction at max height
v1x=v0x %Zero Acceleration
v1y=(sqrt(v0y^2-2*g*(H1-H0)))
%velocity just prior to hitting wall
v2x=v1x
v2y=sqrt(v1y^2+2*g*((H0+H1)-(H0+H2)))
v2=sqrt(v2x^2+v2y^2)
%velocity just after hitting the wall
v3x=v2x*eR
v3y=v2y
v3=sqrt(v3x^2+v3y^2)
theta_3_ref=atand(v3y/v3x)
% Trig: Reference Triangle
theta_3_abs=theta_3_ref+180
% Absolute Angle (Quadrant "3")
% velocity just prior to hitting ground
v4x=v3x
v4y=(sqrt(v3y^2-2*g*(H3-H2)))
v4=sqrt(v4x^2+v4y^2)
theta_4_ref=atand(v4x/v4y)
theta_4_abs=theta_4_ref+180
%velocity just after hitting the ground
v5x=v4x
v5y=v4y*eR
v5=sqrt(v5x^2+v5y^2)
theta_5_ref=atand(v5x/v5y)
theta_5_abs=theta_4_ref+90
%Distance between wall and ground
A=(-g/(2*v3x^2));
B=tand(theta_3_abs);
C=(1/(2*g)*(v4^2-v3^2));
%Quadratic formula for distance between wall and ground
%%based on energy Conservation 2-3 022& Trajectory Eq. 2-3 equations
dCD=(-B+sqrt(B^2-4*A*C))/(2*A)
%trajectory 2
H3=H2+dCD*tand(theta_3_abs)-(g/(2*v3^2*(cosd(theta_3_abs))^2))*dCD^2
%distance between ground and return height
dDA=-dAC-dCD
%trajectory 3
H4=H3+dDA*tand(theta_5_abs)-(g/(2*v4^2*(cosd(theta_5_abs)^2)))*dDA^2
v0(n)=0;
while (H0-return_height_tolerance<H4)<H0+return_height_tolerance
n=0;
n=n+.001;
H2=H0+dAC*tand(theta_0)-g/(2*v0x^2)*dAC^2;
H3=H2+dCD*tand(theta_3_abs)-(g/(2*v3^2*(cosd(theta_3_abs))^2))*dCD^2;
H4=H3+dDA*tand(theta_5_abs)-(g/(2*v4^2*(cosd(theta_5_abs)^2)))*dDA^2;
end
4 个评论
Michael Van de Graaff
2021-12-16
Also, just playing with this, you are using variable before you've defined them all over the place. What is v0?
采纳的回答
Sargondjani
2021-12-16
You need to use a equation solver. I assume it is a non-linear problem, so then use fminsearch or fsolve (or lsqnonlin if you want to use a least squares approach).
Basically what you need to do is:
- code a function that takes has input v0, and gives the distance d from the target (assuming the distance is one dimensional):
function dd = my_distance(par,v0)
% compute the distance using initial speed v0, and all other parameters
dd = ...
end
- Make function handle to use it as the objective:
fun_obj = @(v0)my_distance(par,v0);
- Use the solver, with initial guess v0_ini:
v0_solution = fsolve(fun_obj,v0_ini);
0 个评论
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!