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) >
- 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)
(x/a)^2 = cos^2(ωt)(y/b)^2 = sin^2(ωt)
(x/a)^2 + (y/b)^2 = cos^2(ωt) + sin^2(ωt) = 1
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) >
a(t) = < -a*ω²*cos(ωt), -b*ω²*sin(ωt) >
a(t) = -ω² * r(t)
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) )
Distance = ∫ₜ₀ᵗ₁ ω * sqrt( a²*sin²(ωt) + b²*cos²(ωt) ) dt
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?
- Curved paths require calculus to measure distance.
- Elliptical motion arises naturally from sine and cosine functions.
- 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
aandbto 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
Post a Comment