View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
Gravity Measurement Case Study

4.8 | 4 ratings Rate this file 18 Downloads (last 30 days) File Size: 10.8 MB File ID: #3700 Version:
image thumbnail

Gravity Measurement Case Study


Robert Bemis (view profile)


08 Jul 2003 (Updated )

Scientific imaging case study with MATLAB and toolboxes.

| Watch this File

File Information

The goal of this case study is to show how MATLAB and various toolboxes can be used together to solve an imaging problem. The specific problem shown here is a science experiment. Given a pendulum, measure gravity. The mathematics are well defined by classical physics. Gravity depends on the length of the pendulum and the period of oscillation. An example data set is included in which a ball suspended by a string was used as a pendulum. A captured video of the swinging ball and the time stamps for each video frame are stored in a native MATLAB data file. Another such file contains the predetermined pixel resolution associated with this camera setup.
The Image Acquisition Toolbox was originally used to acquire a video of the oscillating pendulum directly into MATLAB. One advantage of this approach over importing video data from an AVI file was having direct control of the desired image resolution and frame rate of the camera at the device level. It was also convenient being able to work with the camera from within the MATLAB environment instead of having to deal with additional applications and intermediate movie files.

The Image Processing Toolbox was used to find and extract the location of the moving ball in each frame. A circle was then fitted to the ball positions to find the center of rotation, average radius and radial variation or spatial measurement error. The positions were transformed from cartesean to polar coordinates, which produced a damped harmonic motion of angular positions versus time.

The System Identifaction Toolbox was used to fit an autoregressive moving average (ARX) model to the damped sinusoid. This resulted in a discrete-time transfer function description of the harmonic motion. The discrete mathematical model was then transformed to a continuous-time domain transfer function expressed in the frequency domain using Laplace "s" notation. From the complex poles of the transfer function, the radian frequency of oscillation was extracted and then inverted to get the period of oscillation. In addition, variation in time stamps for each frame were determined to estimate temporal measurement error.

Finally, the Symbolic Math Toolbox was used to manipulate the governing equations of motion from physics. This simplified the algebraic task of solving for the measured gravity value. It also simplified the calculus to determine total stack-up error in the gravity measurement.

Required Products Image Acquisition Toolbox
Image Processing Toolbox
Symbolic Math Toolbox
System Identification Toolbox
MATLAB release MATLAB 6.5 (R13)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (5)
29 Aug 2007 Colin Do

20 Jan 2007 Ravi Reddy


06 Jun 2006 Raed alhalaika

i need the mfile code for the tracking object as in picture

Comment only
20 Aug 2005 Max Huber


28 Mar 2005 Chai Yu Kiong

13 Jan 2010 1.1


18 Jan 2010 1.2


02 Feb 2010 1.3

Converted ReadMe.TXT to readme.m (for BSD compliance)

01 Sep 2016

Updated license

Contact us