How can I model and fit the data of the fastest human to compare to the motion

No answer yet for this question. Ask a Tutor

Question:

How can I model and fit the data of the fastest human to compare to the motion of a mammal which I calculated and plotted it?

pls basic level of python and the only permitted libarey is pylab.

here is the task: #photo 3, flow chart #photo 2, sample code #photo1

here is the task and my code

from pylab import *
###################### Functions are defined here ###########################
# Constants
tts_predator = 7 # Time to give up for the predator (in seconds)
time_step = 0.1 # Time step for position calculation (in seconds)

# Power law equation for maximal running speed
def maximal_running_speed(mass):
return 10 ** (1.47832 + 0.25892 * log10(mass) - 0.06237 * (log10(mass)) ** 2)

# Power law equation for time to reach top speed
def time_to_reach_top_speed(mass):
return 10 ** (0.16 + 0.25 * log10(mass))

# function that asks the user to press "0" button to move through the program
def button(x):
x=float(input("Press 0 to continue: "))
return x

###################### main body of the code here ###########################
print("Hi there! Welcome to St Lucia Public Science Museum! This exhibition is called'Mother Nature, racing to survive'. Today we will explore the amazing life of animals in terms of chase and prey in order to survive.")
print("Have you ever wondered about the potential of animals, animals speed or chase and compareing the speed of fastes human on earth with them ?")
print("If you have, then prepare to be amazed. If you haven't, then prepare to be mindblown!")

# prompts user to enter their patron type
patron_type = float(input("Press 2 if you are a science rookie OR Press 1 if you are a scienceenthusiast: "))

# begin program for "science enthusiast"
if patron_type== 1:
# display intro about animals speed
print("The vast plains of the Serengeti comprise 1.5 million ha of savanna. The annual migrationto permanent water holes of vast herds of herbivores (wildebeest, gazelles and zebras),followed by their predators, is one of the most impressive natural events in the world.")
print("Can you guess how both predator and prey accelerate to chase/aviod their adversay?")

x=float(input("Press 0 to continue: "))

#explains maximal_running_speed
print("By using the function of maximum running speed equation we can calculate the speed of animals")
print()

button(x)

# start of loop so user can do multiple animals searches if they wish

# prompt the user to enter the values for mass_predator and mass_prey

mass_predator = float(input("Take a guess and enter choose a mass for the predator: "))
mass_prey = float(input("Take a guess and enter choose a mass for the prey: "))

# Mass of the predator (in kg)

max_speed_predator = maximal_running_speed(mass_predator)
print("This predators max speed is", max_speed_predator)
time_to_top_speed_predator = time_to_reach_top_speed(mass_predator)
print("The time to reach it's maximum speed is ", time_to_top_speed_predator)
max_speed_prey = maximal_running_speed(mass_prey)
print("This prey max speed is", max_speed_prey)
time_to_top_speed_prey = time_to_reach_top_speed(mass_prey)
maximal_running_speed
print("The time to reach it's maximum speed is ", time_to_top_speed_prey)
print()

print("If you want to see a plot ")

x=float(input("Press 0 to continue: "))

# use user inputs to set up variables used in functions that calculate orbital period,transit time and minimum relative intensity
log_mass = array([1.48, 1.78, 2.08, 2.38, 2.68])
log_time_to_top_speed = array([0.16, 0.41, 0.67, 0.93, 1.18])

# Calculate the values of a and b in the power law equation for time to reach top speed
a, b = polyfit(log_mass, log_time_to_top_speed, 1)

# Calculate the maximal running speed and time to reach top speed for the predator and prey
mass_predator = 45 # Mass of the predator (in kg)
mass_prey = 250 # Mass of the prey (in kg)
max_speed_predator = maximal_running_speed(mass_predator)
max_speed_prey = maximal_running_speed(mass_prey)
time_to_top_speed_predator = time_to_reach_top_speed(mass_predator)
time_to_top_speed_prey = time_to_reach_top_speed(mass_prey)

# Calculate the number of time steps for the chase
num_steps = int(tts_predator / time_step)

# Arrays to store positions and times
positions_predator = zeros(num_steps)
positions_prey = zeros(num_steps)
times = zeros(num_steps)

# Calculate positions and times for each time step
for i in range(num_steps):
# Calculate current time
times[i] = i * time_step

# Calculate current speed for the predator and prey
if times[i] < time_to_top_speed_predator:
speed_predator = times[i] * max_speed_predator / time_to_top_speed_predator
else:
speed_predator = max_speed_predator

if times[i] < time_to_top_speed_prey:
speed_prey = times[i] * max_speed_prey / time_to_top_speed_prey
else:
speed_prey = max_speed_prey

# Calculate current positions for the predator and prey
if i > 0:
positions_predator[i] = positions_predator[i-1] + speed_predator * time_step
positions_prey[i] = positions_prey[i-1] + speed_prey * time_step
else:
positions_predator[i] = speed_predator * time_step
positions_prey[i] = speed_prey * time_step

# Print the positions and times
for i in range(num_steps-1):
print("Time:", times[i], "s - Predator Position:", positions_predator[i], "m - Prey Position:", positions_prey[i], "m")

# Plot the positions of the predator and prey
plot(times, positions_predator, 'o-', label='Predator')
plot(times, positions_prey, 'o-', label='Prey')

# Set the axis labels and title
xlabel('Time (s)')
ylabel('Position (m)')
title('Predator-Prey comprison')
print("As you can see the above graph show the position of a preditor and a prey in time period of 7 second which basically the the time that a mammal reaches it's maximum speed which can determine it's position in chase situation")

# Add a legend
legend()

# Show the plot
show()

if patron_type== 2:
print("so your a science rookie")