Code covered by the BSD License  

Highlights from
Gravity Measurement Case Study

4.75

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

Gravity Measurement Case Study

by

 

08 Jul 2003 (Updated )

Scientific imaging case study with MATLAB and toolboxes.

| Watch this File

File Information
Description

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 Processing Toolbox
System Identification Toolbox
Image Acquisition Toolbox
Symbolic Math 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

Good

06 Jun 2006 Raed alhalaika

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

20 Aug 2005 Max Huber

Excellent.

28 Mar 2005 Chai Yu Kiong  
Updates
13 Jan 2010

copyright

18 Jan 2010

copyright

02 Feb 2010

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

Contact us