Following the Elliptical Path: How Particles Move in Mathematical Physics

Have you ever watched a satellite orbit a planet or marveled at the way planets orbit the sun? Their paths often trace out more than just perfect circles — they move in ellipses. But how do we describe such motion mathematically? And how can we measure how far an object has moved if it's not traveling in a straight line? Let’s take a look at a mathematical model for this type of motion, break it down, and even compute how far the particle travels using some Octave code (a free alternative to MATLAB). ---

The Position Vector: A Mathematical Story

In physics, the position of a moving particle in 2D can be described with a position vector. Here's a special one that traces out an ellipse:

r(t) = < a*cos(ωt), b*sin(ωt) >

This vector has:
  • a: the horizontal radius of the ellipse
  • b: the vertical radius
  • ω (omega): the angular speed of the particle
  • t: time
---

What Shape Is This?

Let’s prove it's an ellipse. We define:
  • x(t) = a*cos(ωt)
  • y(t) = b*sin(ωt)
Now, divide both sides by their respective coefficients and square:
  • (x/a)^2 = cos^2(ωt)
  • (y/b)^2 = sin^2(ωt)
Add them:

(x/a)^2 + (y/b)^2 = cos^2(ωt) + sin^2(ωt) = 1

That’s the standard equation of an ellipse. So the particle moves around a squished circle — an ellipse! ---

Velocity and Acceleration

Motion isn’t just about position — it’s about how fast and how the speed changes. We find velocity by taking the derivative of the position vector:

v(t) = < -a*ω*sin(ωt), b*ω*cos(ωt) >

And acceleration by taking the derivative of velocity:

a(t) = < -a*ω²*cos(ωt), -b*ω²*sin(ωt) >

Notice how acceleration is just the original position vector scaled by -ω²:

a(t) = -ω² * r(t)

This means the particle is always being pulled toward the center — like gravity or a spring. ---

How Far Does It Travel?

We can’t just subtract final minus initial position — the particle is moving along a curved path. We need to integrate the speed over time. The speed is the length (magnitude) of the velocity vector:

|v(t)| = ω * sqrt( a²*sin²(ωt) + b²*cos²(ωt) )

So the total distance traveled from time t₀ to t₁ is:

Distance = ∫ₜ₀ᵗ₁ ω * sqrt( a²*sin²(ωt) + b²*cos²(ωt) ) dt

This kind of integral doesn’t have a nice simple formula — but we can calculate it numerically. ---

Let’s Code It in Octave

Below is the Octave code to compute this integral numerically.
% Parameters
a = 5;                 % Horizontal radius
b = 3;                 % Vertical radius
omega = 2 * pi;        % One full cycle per second
t0 = 0;                % Start time
t1 = 1;                % End time (1 second = one full loop)

% Define the speed function
speed = @(t) omega .* sqrt(a^2 .* sin(omega * t).^2 + b^2 .* cos(omega * t).^2);

% Compute distance using numerical integration
distance = quad(speed, t0, t1);

% Display result
fprintf('Distance traveled from t = %.2f to t = %.2f is approximately %.5f units.\n', t0, t1, distance);
Note: If you're using a newer version of Octave or MATLAB, you can replace quad with integral. ---

What Does This Teach Us?

  1. Curved paths require calculus to measure distance.
  2. Elliptical motion arises naturally from sine and cosine functions.
  3. With the right tools — math, physics, and a little code — we can describe and compute motion precisely.
---

Try This Yourself:

Modify the parameters in the code:
  • Change values of a and b to make the ellipse wider or taller.
  • Try computing the distance for half a loop, or multiple loops.
  • Plot the path like this:
t = linspace(t0, t1, 1000);
x = a * cos(omega * t);
y = b * sin(omega * t);
plot(x, y);
axis equal;
title('Elliptical Path');
xlabel('x');
ylabel('y');
--- Thanks for exploring math and motion with me! If you’re curious how this connects to satellites, springs, or atoms — you’re already thinking like a physicist. Let me know in the comments what you'd like to explore next.

Comments

Popular posts from this blog

What is Mathematical Fluency?

Unlocking the Group: Cosets, Cayley’s Theorem, and Lagrange’s Theorem

Verifying the Binomial Formula with Mathematical Induction