User talk:Reedbeta

From Wikipedia, the free encyclopedia
Latest comment: 23 November 2015 by MediaWiki message delivery in topic ArbCom elections are now open!
Jump to navigation Jump to search

Derivation of the equation for simple harmonic motion

There must be an easier way to derive this.

The desired closed-form solution is: x(t)=Acos(ωt+ϕ) for some constants A, ϕ.

A differential equation for a 1-dimensional harmonic oscillator is:

 F=kx

Since force is mass times acceleration,

d2xdt2=kmx

To solve this second-order differential equation, we introduce a new variable v, reducing it to a system of linear first-order differential equations. Also writing ω2=k/m for convenience:

dxdt=v
dvdt=ω2x

Rewrite this system using matrices:

ddt[xv]=[01ω20][xv]

The coefficient matrix has characteristic polynomial

 det(AλI)=λ2+ω2

so its complex eigenvalues are ±ωi. A complex eigenbasis is then:

={[iω1],[iω1]},

which leads to the general solution of our system,

x(t)=c1eωitiω+c2eωitiω

Here, c1 and c2 are the (complex) coordinates of [x0,v0]T (that is, the initial state vector of the system) with respect to the complex eigenbasis .

Applying Euler's formula:

x(t)=ic1ω(cosωt+isinωt)+ic2ω(cosωt+isinωt)

Using the trigonometric identities  cosa=cosa and  sina=sina, this simplifies to:

x(t)=i(c2c1)ωcosωt+c1+c2ωsinωt

Now, let us return to the coordinates c1 and c2. We know that

[x0v0]=[iωiω11][c1c2]

Inverting the matrix, we have

[c1c2]=12[ωi1ωi1][x0v0] where x0 and v0 are both real.

Substituting back into the equation for x(t), we have

x(t)=x0cosωt+v0ωsinωt

Another trig identity gives:

x(t)=x02+v02ω2sin(ωt+θ), where θ=arctanx0ωv0.

If we then let

A=x02+v02ω2=x02+mv02k and
ϕ=θπ2=arctanx0ωv0π2,

then we finally end up at the desired solution

 x(t)=Acos(ωt+ϕ).

Gradient of an ellipsoid

The basic ellipsoid equation is

f(x,y,z)=x2a2+y2b2+z2c21

So the gradient is

f(x,y,z)=2xa2,2yb2,2zc2

ArbCom elections are now open!

{{Wikipedia:Arbitration Committee Elections December 2015/MassMessage}} MediaWiki message delivery (talk) 13:01, 23 November 2015 (UTC)Reply

ArbCom Elections 2016: Voting now open!

Template:Ivmbox

ArbCom 2018 election voter message

Template:Ivmbox