13.1 An introduction to the analysis of brain waves (Page 3/8)

<< Chapter < Page

The art of the pfug

Page 3 / 8

Chapter >> Page >

Equations:

$cos (ϕ + 1) = sin (ϕ + p)$
$sin (ϕ) = - sin (ϕ + p)$
$tan (ϕ) = \frac{sin (ϕ)}{sin ϕ + p}$
$cos (ϕ) = - tan (ϕ + p) \cdot cos (ϕ - p)$

( $tan (ϕ)$ is defined in exercise 1)

Trapezoid rule for estimating area (integration)

Basic method

A useful tool for analyzing curves is finding the area underneath them. When we have an unknown combination of waves, we can estimate the area under the curve using the trapezoid rule . We will use $f (t) = sin (t)$ as an example. If we were to estimate part of the area under the curve with one trapezoid, we might do the following:

We have labeled the two heights, $h_{1}$ and $h_{2}$ , and the length of the base $b$ . The area of the square is:

area of square = (base) \times (height) = b \cdot h_{1}

The area of the top triangle is:

area of triangle = \frac{(base) \times (height)}{2} = \frac{b \cdot (h_{2} - h_{1})}{2}

The total area of the trapezoid is then:

area of trapezoid = b \cdot h_{1} + \frac{b \cdot (h_{2} - h_{1})}{2} = \frac{b \cdot (h_{2} + h_{1})}{2}

If we know that the two points on the $x$ -axis are $t_{1}$ and $t_{2}$ , then $b = t_{2} - t_{1}$ . In the figure above, $t_{1} = . 5$ and $t_{2} = 1.5$ . Then the heights follow from the function: $h_{1} = f (t_{1}) = sin (t_{1})$ and $h_{2} = f (t_{2}) = sin (t_{2})$ . Thus in general, the area of a trapezoid approximating the area under $f$ between the points $t_{1}$ and $t_{2}$ is:

area of trapezoid = \frac{b \cdot (h_{2} + h_{1})}{2} = (t_{2} - t_{1}) \frac{f (t_{1}) + f (t_{2})}{2}

In order to get a good estimate, we split up the domain of the function $f (t)$ into several intervals $[t_{i}, t_{i + 1}]$ . For each interval, we calculate the area of the trapezoid that approximates the area under that curve. For example, we could approximate $f (t)$ over $[0, 1]$ using four equal intervals. This would look like:

In this case, our estimate would be

approximate area of curve = (t_{2} - t_{1}) \frac{f (t_{1}) + f (t_{2})}{2} + \dots + (t_{4} - t_{3}) \frac{f (t_{3}) + f (t_{4})}{2}

As we take smaller and smaller intervals, our approximation will get better, because there will be less space between the trapezoids and the curve. We can prove that the trapezoid rule given `order 2' convergence–that is, if we cut our intervals in half, our error gets four times smaller.

In the general case, if we split up the domain of the function at points ${t_{1}, t_{2}, \dots, t_{n}}$ , then the rule for the estimate is

approximate area of curve = (t_{1} - t_{2}) \frac{f (t_{1}) + f (t_{2})}{2} + \dots + (t_{n} - t_{n - 1}) \frac{f (t_{n - 1}) + f (t_{n})}{2}

This formula can be further reduced, which is the subject of Exercise 2.1.

Coding the trapezoid rule

Here we present a code that uses the trapezoid rule to find the area under any function we provide. We need a vector x that holds the values of the domain, for example x = 0:.01:pi . We then need a vector y that holds the function values at those x points, for example y = sin(x) .

function curve_area = mytrapz(x, y, fast)
% function curve_area = mytrapz(x, y, fast)%
% mytrapz.m performs the trapezoid rule on the vector given by x and y.% Input:
%   x - a vector containing the domain of the function%   y - a vector containing values of the function corresponding to the

curve_area = 0;%loop through and add up trapezoids for as many points as we are givenfor n = 2 : numel(x)

We start the code with zero area under the curve, since we haven't counted anything yet. Then we create a for loop to count each triangle individually. As we see above, more trapezoids leads to better answers, so we want to use as many trapezoids as we possibly can. In this situation, that means using every point in x and y . The function numel simply counts the number of elements in x . We then calculate the area of the current triangle (within the loop):

<< Chapter < Page Page > Chapter >>

Read also:

Get Jobilize Job Search Mobile App in your pocket Now!

100% Free Mobile Applications
Receive real-time job alerts and never miss the right job again

Source: OpenStax, The art of the pfug. OpenStax CNX. Jun 05, 2013 Download for free at http://cnx.org/content/col10523/1.34

Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'The art of the pfug' conversation and receive update notifications?

Ask

	6 Neuroanatomy 06 Head Somatic Visceral Sensory By Stephen Voron Start Quiz
	2 Lec:2 Randomized Controlled Trials By Janet Forrester Start Quiz
	Immunology Practice Test By Sandhills MLT Start Test
	NCE Ch 09 Research and Program Evaluation By Anh Dao Start Quiz
	MAPEH TEST By Edgar Delgado Start Test
	20 AP 20 Blood Vessels Circulation MCQ By OpenStax Start Quiz
	6 BOD Respiratory Exam By Brooke Delaney Start Exam
©flickr: Abraham	3 Biology 3 By Sarah Warren Start Test
	NCE Ch 04 Social and Cultural Foundations By Anh Dao Start Quiz
	Vocabulary for "A Rose for Emily" By Bonnie Hurst Start Quiz

13.1 An introduction to the analysis of brain waves (Page 3/8)

Equations:

A sine wave.

A cosine wave.

Trapezoid rule for estimating area (integration)

Basic method

A single trapezoid.

Estimating the area with four trapezoids.

Coding the trapezoid rule

Read also:

Get Jobilize Job Search Mobile App in your pocket Now!