Write a script for Matlab that calculates the position of a ball being dropped f

Question

Write a script for Matlab that calculates the position of a ball being dropped from a height of 10 meters with an initial velocity of 0 as a function of time. Have the script accept the value of time from the user using the input function. Then display the result in the command window.

Now the position of the ball as a function of time is known calculate the position of the ball for a set of time points using a ‘for loop.’ Start your time vector at time 0 (the initial time). Continue your time points until the ball hits the ground (at height = 0).

Lastly add a plot and figure inside the ‘for loop’ and plot the position of the ball as the for loop runs.

Explanation / Answer

If u want i can send this file to ur email id.Message me the id. In this program u can input many variables.Choice is yours u can make some variables constant function projectile_motion % This program is used to view the path & distance covered by a particle in % animated form. % % Execute function in command window it will prompt yout to enter values % required by program, after getting input parameters it will calculate the % equations for projectile motion then it will show the animated path and % distance covered by particle. % %% Input Parameters for projrctile motion theta = input('Enter launch angle= '); %Launch angle, theta (degrees) m = input('Enter mass of the body= '); %Mass, m (grams) Fs = input('Enter constant force spring force= '); %Constant force spring force, Fs (N) Ns = input('Enter number of springs= '); %Number of springs, Ns LL = input('Enter distance over which force is applied= '); %Distance over which force is applied, LL (mm) eta = input('Enter efficiency= '); %efficiency, eta g = input('Enter acceleration due to gravity= '); %Acceleration due to gravity, g (m/s^2) N = input('Enter no. of steps= '); %% Calcultion for projectile motion V = sqrt(2.*Ns.*Fs.*LL.*eta./m); %Launch velocity, V (m/s) Vy = V.*sind(theta); %Vertical velocity component, Vy (m/s) Vx = V.*cosd(theta); %Horizontal velocity component, Vx (m/s) tmh = V.*sind(theta)./g; % Time to maximum height, tmh (seconds) ttotal = 2.*tmh; % Total travel time, ttotal (seconds) ymax = (V.^2).*(sind(theta).^2)./(2.*g); %Maximum height, ymax (m) xmax = (V.^2).*(sind(2.*theta))./g; %Maximum range, xmax (m) tinc = ttotal./N; %Time increment, tinc %% Preallocation of answer spaces time = zeros(N,1); %Preallocation for time xpos = zeros(N,1); %Preallocation for X-Position ypos = zeros(N,1); %Preallocation for Y-Position %% Loop for generation of time vector for l = 1:N-1 time(l+1) = time(l)+tinc; end %% Loop for generation of x-position & y-position vector for j = 1:N xpos(j) = time(j).*Vx; ypos(j) = V.*time(j).*sind(theta)-0.5.*g.*time(j).^2; end %% Loop for visualization of plot for k = 1:N hold all h = plot(time(k),xpos(k),'p',time(k),ypos(k),'o',... 'LineWidth',2,'MarkerEdgeColor','g',... 'MarkerFaceColor','y','MarkerSize',2+k); pause(0.6) drawnow end %% Annotation of graph legend('Distance','Path','Location','NorthWest') title('Projetile Trajectory','FontSize',14,'FontWeight','Bold','FontName','Calibri') ylabel('Distance(m)','FontSize',12,'FontWeight','Bold','FontName','Calibri') xlabel('Time(s)','FontSize',12,'FontWeight','Bold','FontName','Calibri') grid on

Related Questions

Navigate

• Browse (All)
• Browse W
• Subjects

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.