Hello,I know solution method for solving differentialequations.

  • Thread starter Thread starter HWGXX7
  • Start date Start date
  • Tags Tags
    Method
HWGXX7
Messages
46
Reaction score
0
Hello,

I know solution method for solving differentialequations. But last I was asked to solve the equation that describes a mathematical pendulum.

\frac{d^{2}\theta}{dt^{2}}+\frac{g}{l}.sin(\theta)=0

How can I solve this the with the most common method?

thank you
 
Physics news on Phys.org


Usually you would assume a small angle to get rid of the sin() term so the problem becomes linear.
Formally, you could express the solution in terms of the Jacobi elliptic function.
 


I'am familiar with the fact that I can assume the sin component equal to the angle, only for small angles..But the questions was related to angles great dan 45°, so I may not neglect the sin component.

Is it also possible to express the sinus as a power series with n = 5 terms, to get a good result?
I don't know much about Jacobi elliptic function...

thanks
 


I an concerned about you saying "I know solution method for solving differentialequations". One of the first things you should have learned is that there is NO one method of solving differential equations
In particularly, this is a badly non-linear equation and there is no "method" for solving it. There are a number of ways of either getting an approximate solution or getting qualitative information about the solution.

1) Linearization. The set up of your equation implies that you are taking "0" as the bottom of the swing. As long as you do not swing the pendulum "too hard" you have the angle staying small and you can approximate sin(x) by x. Solve \theta''+ \theta=0.

2) "Quadrature". By the chain rule, if we let \omega= d\theta/dt,
\frac{d^2\theta}{dt^2}= \frac{d\omega}{dt}= \frac{d\omega}{d\theta}\frac{d\theta}{dt}= \omega\frac{d\omega}{dt}
The equation becomes
\omega\frac{d\omega}{d\theta}+ (g/l)sin(\theta)= 0
which is a separable first order equation for \omega. Of course, after solving for \omega you have another first order equation for \theta. That equation can be reduced to an integral but, an "elliptic integral, one that has no elementary anti-derivative. That means that it must be done numerically which is why this is an approximation method.

3) "perturbation". sin(\theta)= \theta- (1/3!)\theta^3+ (1/5!)\theta^5- .... Dropping all except the \theta gives the linear equation of (1). Then subtract that solution from \theta and put the difference into the equation including the \theta^2 term. Continue as long as you wish.

4) "phase plane". Go ahead and integrate \int d\omega= -(g/l)sin(\theta)d\theta to get \omega= (g/l)cos(\theta)+ C and graph \omega against \theta for different values of C. That gives information about the behavior of solutions.
 


I an concerned about you saying "I know solution method for solving differentialequations". One of the first things you should have learned is that there is NO one method of solving differential equations

Well, I meant the methods I learned for solving (high) second order lineair (non) homogenous differential equations.

I know methode of annihilation: by which the differential eqaution is transformed to a differentialpolynoom. For example: (D-a)y(t)=0 has solution: Ce^{-at}

Also possible for non homogenous equations..

I assumed that this a method for solving ODE..

Back to topic: you proposed a good list of possible way's to solve the equation. I migth be able tot handle the problem now.

Thank you!
 
Back
Top