The Basics of MRI

Chapter 2

THE MATHEMATICS OF NMR



Exponential Functions

The number 2.71828183 occurs so often in calculations that it is given the symbol e. When e is raised to the power x, it is often written exp(x).

ex = exp(x) = 2.71828183x

Logarithms based on powers of e are called natural logarithms. If

x = ey

then

ln(x) = y,

Many of the dynamic MRI processes are exponential in nature. For example, signals decay exponentially as a function of time. It is therefore essential to understand the nature of exponential curves. Three common exponential functions are

y = e-x/t 
y = (1 - e-x/t
y = (1 - 2e-x/t

where t is a constant. 

Trigonometric Functions

The basic trigonometric functions sine  and cosine  describe sinusoidal functions which are 90o out of phase.

The trigonometric identities are used in geometric calculations. 

Sin(q) = Opposite / Hypotenuse
Cos(q) = Adjacent / Hypotenuse
Tan(q) = Opposite / Adjacent

Three additional identities are useful in understanding how the detector on a magnetic resonance imager operates.

Cos(q1) Cos(q2) = 1/2 Cos(q1 - q2) + 1/2 Cos(q1 + q2)

Sin(q1) Cos(q2) = 1/2 Sin(q1 + q2) + 1/2 Sin(q1 - q2)

Sin(q1) Sin(q2) = 1/2 Cos(q1 - q2) - 1/2 Cos(q1 + q2)

The function sin(x) / x occurs often and is called sinc(x). 

Differentials and Integrals

A differential can be thought of as the slope of a function at any point. For the function

the differential of y with respect to x is

An integral is the area under a function between the limits of the integral.

An integral can also be considered a summation; in fact most integration is performed by computers by adding up values of the function between the integral limits. 

 

Vectors

A vector is a quantity having both a magnitude and a direction.  The magnetization from nuclear spins is represented as a vector emanating from the origin of the coordinate system. Here it is along the +Z axis.  In this picture  the vector is in the XY plane between the +X and +Y axes. The vector has X and Y components and a magnitude equal to

( X2 + Y2 )1/2

Matrices

A matrix is a set of numbers arranged in a rectangular array.  This matrix has 3 rows and 4 columns and is said to be a 3 by 4 matrix.

To multiply matrices the number of columns in the first must equal the number of rows in the second.  Click sequentially on the next start buttons to see the individual steps associated with the multiplication.     

Coordinate Transformations

A coordinate transformation is used to convert the coordinates of a vector in one coordinate system (XY) to that in another coordinate system (X"Y"). 

Convolution

The convolution of two functions is the overlap of the two functions as one function is passed over the second.  The convolution symbol is . The convolution of h(t) and g(t) is defined mathematically as

The above equation is depicted for rectangular shaped h(t) and g(t)functions in this animation. 

Imaginary Numbers

Imaginary numbers are those which result from calculations involving the square root of -1. Imaginary numbers are symbolized by i.

A complex number is one which has a real (RE) and an imaginary (IM) part. The real and imaginary parts of a complex number are orthogonal.

Two useful relations between complex numbers and exponentials are

e+ix = cos(x) +isin(x)
and
e-ix = cos(x) -isin(x).

The quantity e+ix is said to be the complex conjugate of e-ix. In other words, the complex conjugate of a complex number is the number with the sign of the imaginary component changed.

Fourier Transforms

The Fourier transform (FT) is a mathematical technique for converting time domain data to frequency domain data, and vice versa. 

The Fourier transform will be explained in detail in Chapter 5.


Problems

  1. In a right triangle the hypotenuse is 5 cm, and the remaining two sides are 3 cm and 4 cm. What is the size of an angle opposite the 3 cm long side? 
  2. What is the integral of y between 0 and 5 where y = 3x2 + 3 ? 
  3. What is the slope of the function y = 3x2 + 3 at x = 2 ? 
  4. What is the product of eiA with eiB
  5. You have some laboratory data which has the functional form y = e-x/t. What could you plot as a function of x to make the data linear? 
  6. What is the product of these two matrices?  

Go to the: [next chapter | chapter top | previous chapter | cover ]

Copyright © 1996-99 J.P. Hornak.
All Rights Reserved.