10992

Extended Formulas (Closed)

Лекция

Математика и математический анализ

Extended Formulas Closed If we use equation 2.5 N 1 times to do the integration in the intervals x1; x2; x2; x3; xN 1; xN and then add the results we obtain an extendedr or compositer formula for the integral from x1 to xN. Extended trapezoidal rule: In this method the area under the curve is approximated by sums of trapezoids areas under the curve see Fig. 2.3.. Figure 2.3. Extended trapezoidal rule. Trapezoid formul...

Английский

2013-04-03

145 KB

0 чел.

Extended Formulas (Closed)

If we use equation (2.5) N – 1 times to do the integration in the intervals (x1; x2); (x2; x3); … (xN - 1; xN), and then add the results, we obtain an “extended” or “composite” formula for the integral from x1 to xN.

Extended trapezoidal rule:

In this method the area under the curve is approximated by sums of trapezoids areas under the curve (see Fig. 2.3.).

Figure 2.3. Extended trapezoidal rule.

Trapezoid formula is

  (2.9)

Here we have written the error estimate in terms of the interval b – a and the number of points N instead of in terms of h.  This is clearer, since one is usually holding a and b fixed and wanting to know (e.g.) how much the error will be decreased by taking twice as many steps (in this case, it is by a factor of 4). In subsequent equations we will show only the scaling of the error term with the number of steps.

If we apply equation (2.5) to successive, nonoverlapping pairs of intervals, we get the extended Simpson’s rule:

 (2.10)

Notice that the 2/3, 4/3 alternations continues throughout the interior of the evaluation. Many people believe that the wobbling alternation somehow contains deep information about the integral of their function that is not apparent to mortal eyes. In fact, the alternation is an artifact of using the building block (2.6).

The trapezoidal rule is two orders lower than Simpson’s rule, however, its contribution to the integral goes down as an additional power of N (since it is used only twice, not N times). This makes the resulting formula of degree one less than Simpson.

We can construct open and semi-open extended formulas by adding the closed formulas, evaluated for the second and subsequent steps, to the extrapolative open formulas for the first step. The resulting formula for an interval open at both ends is as follows. This one is known as the extended midpoint rule

(2.11)

This rule is shown in Figure 2.4.

Figure 2.4. Midpoint extended formula with equally spaced abscissas compute the integral of a function y = x3.

Example

Suppose we want to evaluate the integral  The definite integral is the area under the curve y = ex, above the x-axis and between the lines x = 1 and
x = 2.

It is obvious the above approximation is an underestimate for the actual value. The actual value can be easily obtained as follows

We can find an approximation by Left Sum Approximation formula:

(2.12)

Figure 2.5. Left Sum Approximation.

In this rule we divide the interval [a, b] into subintervals with same length (h). Construct rectangle over each subinterval, the height of the rectangle being the function value at the left end point of the subintervals. The sum of the areas of these rectangles will provide an approximation to the area under the curve on the interval (see Fig 2.5).

Right Sum Approximation

We can find another approximation to the area under the graph. Divide the interval [a, b] into subintervals and construct rectangle over each subinterval, the height of the rectangle being the function value at the right end point of the subinterval.

The sum of the areas of these rectangles will provide an approximation to the area under the curve and is shown in Fig. 2.6.


Figure 2.6. Right Sum Approximation.

Right Sum Approximation formula is:

(2.13)

The following table gives the number of subintervals and the corresponding approximations.

Number of subintervals

Approximate value of integral

Left Sum Approximation

Midpoint Approximation

Right Sum Approximation

Trapezoidal

Approximation

Simpson’s Approximation

10

4.44112722

4.66882868

4.90820465

4.25683746

4.67077443

100

4.64745932

4.67075481

4.69416706

4.62903035

4.67077427

In this table we can compare accuracy of different approximation formulas and see how the increase in the number of subintervals affects the approximation.

As the number of subintervals increases, the approximations gets better.

2.6. Simple Monte Carlo Integration

Suppose that we pick N random points, uniformly distributed in an interval [a, b]. Call them x1, x2, ... xN. Then the basic theorem of Monte Carlo integration estimates the integral of a function f over the interval
[a, b],

. (2.14)

In formula (2.14) we note that each random point xi is actually a random variable which is uniformly distributed on [a, b]

. (2.15)

Uniform deviates  are just random numbers that lie within a specified range [0, 1], with any one number in the range just as likely as any other. Then f(xi) is also a random variable as a function of random point xi.

There is an another advanced technique of Monte Carlo integration (see Fig. 2.7).

Figure 2.7. Monte Carlo integration.

Random points are chosen within the area A. The integral of the function f is estimated as the area of A multiplied by the fraction of random points that fall below the curve f:

. (2.16)

Where S(A) is the area of A;  n is the number of random points that lies below the curve f; N is the total number of random points.

The Monte Carlo method thus becomes a deterministic quadrature scheme — albeit a simple one — whose fractional error O(.) decreases at least as fast as .

2.7. Questions

1. What is a problem of a numerical integration?

2. What is possible to say about error of left sum, right sum and midpoint approximation formulas?

3. How to increase an accuracy of calculations in the problem of numerical integration?

4. Compare both advantages and defects of trapezoidal and Monte Carlo integration formulas.

PAGE  1


EMBED PBrush  

EMBED Word.Picture.8  

EMBED Word.Picture.8  

EMBED Word.Picture.8  


 

А также другие работы, которые могут Вас заинтересовать

40766. Стилістика як мовознавча дисципліна Стилістична система української мови. Стилістична норма 670.91 KB
  Стилістика як мовознавча дисципліна Стилістична система української мови. Стилістична норма Мета: поглибити знання студентів про систему стилів ознайомити студентів із стилістичними нормами набуття ними знань щодо стилістичних особливостей сучасної української літературної мови; вивчення особливостей функціонування офіційноділового наукового публіцистичного художнього та розмовного стилів розкрити усі параметри їх характеристик їх різновиди; розвивати вміння будувати висловлювання за даним стилем; виховувати мовне чуття яке...
40767. Ординалістська теорія поведінки споживача. Гранична норма заміщення благ: суть і методика обчислення 61.08 KB
  Криві байдужості як спеціальний інструментарій мікроекономічного аналізу. Карта байдужості. Із цієї гіпотези можна зробити висновок що криві байдужості випуклі до початку координат. Крива байдужості Залежність сукупної корисності від кількості двох товарів TU крива споживчої байдужості однакової сукупної корисності.
40768. МОДЕЛЮВАННЯ ПОВЕДІНКИ СПОЖИВАЧА НА РИНКУ ТОВАРІВ 138.97 KB
  Тоді формула бюджетного обмеження матиме вигляд: Гранична норма заміщення благом Х витрат споживача на інші блага виражається в грошовій формі й означає міру готовності пожертвувати можливістю витратити бюджет на інші блага задля споживання додаткової одиниці Х залишаючись на даній кривій байдужості.Побудова лінії Енгеля на основі лінії дохід ...
40769. Попит, пропозиція їх взаємовідносини 78.95 KB
  Такими є основні контури моделі в теорії попиту і пропозиції що посідає одне з найважливіших місць у поясненні поведінки ринкових субєктів як на мікро так і на макрорівні. Функція на якій базується модель попиту має вигляд 5.1 де QD обсяг величина попиту на товар одиниць за проміжок часу; P ціна на товар грош. Обсяг величина попиту є саме такою кількістю товару яку бажає і може придбати споживач споживачі в певний проміжок часу за деякою ціною з діапазону можливих цін за інших незмінних умов.
40770. Економічна природа лізингового кредиту 32.86 KB
  Об'єктом предметом лізингу є будьяке нерухоме і рухоме майно яке може бути віднесене до основних фондів: машини устаткування транспортні засоби обчислювальна та інша техніка будинки споруди системи телекомунікацій тощо не заборонене законодавством до вільного обігу на ринку. Суб'єктами лізингу як правило є три особи: постачальник або продавець предмету лізингу підприємство організація та інші суб'єкти господарювання які здійснюють виробництво і реалізацію машин і обладнання; лізингодавець суб'єкт господарювання який...
40771. Валютні відносини та валютні системи 444.55 KB
  Поняття та види валюти. У цей момент національні гроші набувають форми валюти. Використання валюти для обслуговування зовнішньоекономічних відносин надає останнім специфічної форми валютних відносин. Залежно від економічного призначення валютних платежів у валютних відносинах можна виділити такі складові: міжнародні розрахунки за зовнішньоторговельними операціями; одержання і погашення різних видів міжнародних позичок; переміщення валюти при здійсненні зовнішніх інвестицій; переміщення валюти при наданні економічної і технічної...
40772. Логос античної і середньовічної філософії 40.42 KB
  Античні філософи Анаксагор Анаксімандр Анаксімен Аристокл Аристотель Геракліт Гіппій Горгій Демокріт Демокріт Емпедокл Епіктет Епікур Зенон Ксенофан Левкіп Лукрецій Кар Марк Аврелій Марк Терренцій Варрон Меліс Панецій Парменід Піфагор Платон Плотін Порфирій Посідоній Прарменід Продік Прокл Протагор Сенека Сократ Фалес Ціцерон Антична філософія. Давньогрецька філософія класичного періоду. Філософія Арістотеля: логіка метафізика етика. Елліністична і римська філософія...
40773. ТЕХНОГЕННІ НЕБЕЗПЕКИ ТА ЇХ НАСЛІДКИ 48.43 KB
  Техногенні небезпеки та їх вражаючі фактори Техногенні небезпеки це небезпеки повязані з використанням транспортних засобів експлуатацією транспортних комунікацій використанням горючих легкозаймистих і вибухонебезпечних речовин та матеріалів хімічних речовин небезпечних гідродинамічних обєктів та різними видами випромінювання. Станом на 2005 рік Державний реєстр потенційно небезпечних обєктівПНО містить відомості про 9000 обєктів до переліку яких входять промислові підприємства шахти карєри магістральні газо нафто і...