Trajectory Evaluation

A particle is launched at time t = 0 t = 0 from ground level with initial speed v 0 v_0 at an angle θ \theta with respect to ground. Gravitational acceleration g g is downward.

Let T ( t ) T(t) be the particle's kinetic energy as a function of time, and let V ( t ) V(t) be the particle's gravitational potential energy as a function of time, measured with respect to ground. Let t f t_f be the time at which the particle lands.

If θ = 45 \theta = 45 degrees, what is the following ratio?

0 t f T ( t ) d t 0 t f V ( t ) d t = ? \large{\frac{\int_0^{t_f} T(t) \, dt}{\int_0^{t_f} V(t) \, dt} = \,\, ?}


The answer is 2.0.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

4 solutions

Krishna Karthik
Jun 25, 2020

Here's a python solution.

Neeraj, what I have done here in the code is as follows:

First step: Set the condition to when the object hits the ground, when y > = 0 y >= 0 .

Second step: Calculate the Kinetic and Potential energy of the object during trajectory by finding magnitude of velocity and dividing by 2, and potential by finding m g h mgh .

Third: Increment the value of x x every 1 0 7 10^{-7} seconds. This will be done by multiplying the time by the velocity, does adding small distances to the current distance. On a more technical perspective this is called Euler Integration.

4rth: Do the same steps to the y y value, by calculating velocity using the acceleration and then numerically integrating to get the displacement.

5th: using the computations, numerically integrate to find the integrals of the Kinetic and Potential Energy.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
import math


#time values and increment
time = 0
deltaT = 10**-7 

#x and y
x = 0
y = 0

#velocities in both directions (arbitrary)
xDot = 10
yDot = 10

#arbitrary g value
g = 10


#Kinetic Energy and potential arbitrary initialisation
integralKinetic = 0
integralPotential = 0

while y >= 0:
    #kinetic and potential computation
    Kinetic = (xDot**2+yDot**2)/2
    Potential = g*y

    #acceleration in y
    yDotDot = -g

    #updating x value
    x += xDot*deltaT

    #numerically integrating yDot to update y value
    yDot += yDotDot*deltaT
    y += yDot*deltaT

    #Kinetic energy increment and Potential increment
    dK = Kinetic*deltaT
    dP = Potential*deltaT

    #adding dKs and dPs to get the resulting integral
    integralKinetic += dK
    integralPotential += dP

    time += deltaT

print(integralKinetic/integralPotential)

@Krishna Karthik Thank you so much sir for this solution.
I just upvoted

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

You're welcome bro, cheers!

Krishna Karthik - 11 months, 3 weeks ago

Woah I just found a really weird brilliant glitch... You can actually get the "interesting" option into the negatives such as negative 5 and stuff by just repeatedly clicking lol

Krishna Karthik - 11 months, 3 weeks ago

Log in to reply

@Krishna Karthik yes, nice catch.
But i think after some time it will come to it's original values.
I don't know why but Helpful option is not going into negative. Ha ha


A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

@A Former Brilliant Member Yeah; if you refresh the browser it's simply all increment by only one. Kinda funny though lol🤣

Krishna Karthik - 11 months, 3 weeks ago

Log in to reply

@Krishna Karthik @Krishna Karthik think logical if confused option is in negative then it means, that the person agrees with the solution.

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

@A Former Brilliant Member Lol, that's a funny way to interpret a bug. If the interesting option is negative then that's kinda sad🤣 If the brilliant option is negative then someone is trying to call the solution writer dumb LMAO🤣

Krishna Karthik - 11 months, 3 weeks ago

Log in to reply

@Krishna Karthik @Krishna Karthik do you play chess?
If yes, we can play chess together at chess.com

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

@A Former Brilliant Member Yes, I certainly play chess. Let's play at chess.com!

Krishna Karthik - 11 months, 3 weeks ago

Log in to reply

@Krishna Karthik @Krishna Karthik please give me your username

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

@A Former Brilliant Member It's Kk1356

Krishna Karthik - 11 months, 3 weeks ago

Log in to reply

@Krishna Karthik @Krishna Karthik request sent

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

@A Former Brilliant Member Yup, I accepted. Let's play!

Krishna Karthik - 11 months, 3 weeks ago
Karan Chatrath
Jun 24, 2020

Python code based solution:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
import math

# Parameter initialisation:
m     = 1
u     = 1
theta = math.pi/4
g     = 10

# Time step:
dt    = 10**-6

# Variables initialisations:
t     = 0
x     = 0
y     = 0
dx    = u*math.cos(theta)
dy    = u*math.sin(theta)

# Integrals of KE and PE wrt time initialisation:
TS    = 0
VS    = 0

# Loop for numerical integration:
while y >= 0:

    # Acceleration components:
    ddx = 0
    ddy = -g

    # Kinetic and potential energies:
    T   = 0.5*m*(dx**2 + dy**2)
    V   = m*g*y

    # Numerical integration to obtain velocity components:
    dx  = dx + ddx*dt
    dy  = dy + ddy*dt

    # Numerical integration to obtain position coordinates:
    x   = x + dx*dt
    y   = y + dy*dt

    # Numerical integration to obtain Integrals of KE and PE respectively:
    TS  = TS + dt*T
    VS  = VS + dt*V

    t   = t + dt

# Required ratio:
ANSWER = TS/VS
print(ANSWER)
# 2.0000348700607242

@Karan Chatrath sir thank you so much ,i was struggling for this (i just upvoted)
btw how can i plot α = 3 2 sin 2 θ 1 \large \alpha = \frac {3 }{2\sin^{2} \theta}-1 this in excel.
if i plot that curve in my solution it will be awesome .
thanks in advance

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

This is a pretty straightforward operation in Excel. Refer to the internet to see how to plot a any function y = f ( x ) y=f(x) using Excel. The explanation would be detailed but the steps are simple.

Karan Chatrath - 11 months, 3 weeks ago

Nice one. One could also run while y = 0 y \geq = 0 , eliminating the need to calculate the final time.

Steven Chase - 11 months, 3 weeks ago

Log in to reply

@Steven Chase sir for this problem only?
if we do y > = 0 y>=0
we don't have to calculate time???
is this you want to say

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

If you use the "while" loop, you don't need to pre-calculate the final time

Steven Chase - 11 months, 3 weeks ago

Log in to reply

@Steven Chase @Steven Chase sir i did not understand much,i will be happy if you post the solution of this problem only by the method of python without calculating time.
thanks in advance.

A Former Brilliant Member - 11 months, 3 weeks ago

Yes, in hindsight, choosing y y as the loop variable is a better way of solving the problem. I have edited my solution. Thanks for the feedback.

Karan Chatrath - 11 months, 3 weeks ago

Log in to reply

Yes, very true. Calculating the time required is an unnecessary condition if you know it'll land at y = 0 y = 0 .

Krishna Karthik - 11 months, 3 weeks ago

@Steven Chase ohhh after seeing the edited solution i can feel what Steven sir wants to say 1 hour ago.

A Former Brilliant Member - 11 months, 3 weeks ago

@Karan Chatrath sir how can we make 2D images in python representing 3D things.
Suggest me some website or video to learn.
Thanks in advance.

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

By images, do you mean plots?

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath Yes sir like this problem have https://brilliant.org/problems/more-paraboloid-surface-dynamics/
I also want to post question like this images I like this types of images very much.
How can I learn to make this images?
Thanks in advance.

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

@A Former Brilliant Member You can write a code on python generating the coordinates of that surface and then you can export it to excel and generate a 3D plot. I am not sure how to export data from python to Excel but I am sure it is possible.

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath @Karan Chatrath can you give me 1 code , no matter whatever the plot in 3D you are giving me.
I know how to export it to excel.
But I don't know how to write this code in python.
If you give me one code it will be a good example for me , to how to draw
We can check.
Thanks in advance.

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

@A Former Brilliant Member Okay, I will do that. I do not know how to export data from Python to excel, so I am hoping you will help me with that. That way we both learn something new. Sounds good?

I will type and code and share it after some time.

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

Log in to reply

@A Former Brilliant Member Actually, can you tell me how to export data from python to excel? I will attach the code shortly

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath That opaque helix is one example. But there are much simpler examples

Steven Chase - 11 months, 2 weeks ago

Unit hemisphere centered at the origin. The data is exported to 'HEM.txt'. Copy that to Excel and try 3D plotting the vectors.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
import math

dtheta = math.pi/100
dphi   = math.pi/100

theta  = 0

f = open('HEM.txt', 'w')

while theta <= math.pi/2: 

   phi = 0
   while phi <= 2*math.pi:

      X = math.sin(theta)*math.cos(phi)
      Y = math.sin(theta)*math.sin(phi)
      Z = math.cos(theta)

      f.write(str(X)+' '+str(Y)+' '+str(Z)+'\n')

      phi = phi + dphi


   theta = theta + dtheta

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Neeraj Anand Badgujar Hope this is helpful.

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath I am getting this

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

@A Former Brilliant Member Yes, that is not a 3D plot. How to generate 3D plots in Excel is something I am unsure of at this moment. But the code I provided does generate the required data.

Maybe @Steven Chase can provide some insight into this.

Some context: This is an unsuccessful attempt at plotting a unit hemisphere centred at the origin, as a 3D surface.

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath @Karan Chatrath sir i think Steven sir is sleeping now .
Hope he will give us a reply in 3-4 hours .

A Former Brilliant Member - 11 months, 2 weeks ago

@Karan Chatrath @Karan Chatrath sorry sir . You have done your work from your side is perfect but i failed for my contribution.
but i am trying that how to make 3-D plots from web

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

@A Former Brilliant Member Oh, no worries. 3D plots on Excel are tricky. I couldn't figure this out myself and am curious to learn about it. I will continue to try searching for an answer later on.

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath @Karan Chatrath can we make in MATLAB??

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

@A Former Brilliant Member Yes

Karan Chatrath - 11 months, 2 weeks ago

Log in to reply

@Karan Chatrath @Karan Chatrath wow it's very Beautifull

A Former Brilliant Member - 11 months, 2 weeks ago

@Karan Chatrath Regarding 3D plots in Excel, as I do them, Excel is very dumb. It just plots points. The code script is where all the brains are.

Steven Chase - 11 months, 2 weeks ago

Log in to reply

@Steven Chase @Steven Chase i laughed so hard when i read that "Excel is very dumb"
sir according to you which tool is best for making 3-D plots except tools like MATLAB which are very expensive??

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

@A Former Brilliant Member There is an add-on pack for Python called sci-py that has all kinds of extra features like plotting. That might be something to look at.

Steven Chase - 11 months, 2 weeks ago

Log in to reply

@Steven Chase @Steven Chase How can i install it in my computer ??
Did u have installed??

A Former Brilliant Member - 11 months, 2 weeks ago

From this Problem we know that the generalise value of 0 t f T ( t ) d t 0 t f V ( t ) d t = α \frac{\int_{0}^{t_{f} } T(t) dt} {\int_{0}^{t_{f} } V(t) dt } = \alpha
comes out to be sin θ + 4 sin 3 θ 3 2 sin 3 θ 2 sin 3 θ 3 = α \Large \frac{\sin \theta+ \frac{4 \sin^{3} \theta }{3} -2\sin^{3} \theta}{ \frac{2 \sin^{3} \theta }{3}} = \alpha
simpifying the α \alpha leads to α = 3 2 sin 2 θ 1 \large \alpha = \frac {3 }{2\sin^{2} \theta}-1
Now ,lets put θ = 45 \theta =45 degrees α = 2 \boxed{ \alpha = 2}

I will be happy if anyone will post the solution of this problem by solving it with python programming.
Thanks in advance

A Former Brilliant Member - 11 months, 3 weeks ago

Log in to reply

Just solved this problem. I have posted a python based solution.

Karan Chatrath - 11 months, 3 weeks ago

We can show (see this ) that the ratio is

η = 3 2 sin 2 θ 1 \eta=\dfrac {3}{2\sin^2 \theta}-1

Here θ = 45 ° \theta=45\degree

So η = 3 2 × 1 2 1 = 3 1 = 2 \eta=\dfrac {3}{2\times \frac 12}-1=3-1=\boxed 2 .

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...