Заглавная страница Избранные статьи Случайная статья Познавательные статьи Новые добавления Обратная связь FAQ Написать работу КАТЕГОРИИ: АрхеологияБиология Генетика География Информатика История Логика Маркетинг Математика Менеджмент Механика Педагогика Религия Социология Технологии Физика Философия Финансы Химия Экология ТОП 10 на сайте Приготовление дезинфицирующих растворов различной концентрацииТехника нижней прямой подачи мяча. Франко-прусская война (причины и последствия) Организация работы процедурного кабинета Смысловое и механическое запоминание, их место и роль в усвоении знаний Коммуникативные барьеры и пути их преодоления Обработка изделий медицинского назначения многократного применения Образцы текста публицистического стиля Четыре типа изменения баланса Задачи с ответами для Всероссийской олимпиады по праву Мы поможем в написании ваших работ! ЗНАЕТЕ ЛИ ВЫ?
Влияние общества на человека
Приготовление дезинфицирующих растворов различной концентрации Практические работы по географии для 6 класса Организация работы процедурного кабинета Изменения в неживой природе осенью Уборка процедурного кабинета Сольфеджио. Все правила по сольфеджио Балочные системы. Определение реакций опор и моментов защемления |
Метод кратного перерахунку за допомогою ExcelСодержание книги
Поиск на нашем сайте
Задача 1. 1. Обчислити наближене значення інтеграла функції f (x) = esin x cos2x на відрізку [0;1] за узагальненою формулою трапецій з кроками інтегрування h = 0,2 0,1 0,05. 2. Уточнити значення інтеграла і оцінити його похибку методом кратного перерахунку. Розв’язання. 1.Спочатку побудуємо електронну таблицю значень функції f (x) у вузлах інтегрування. Надамо чарункам електронної таблиці таких значень:
Як і раніше, символ ↓ означає копіювання попередніх чарунок. Тут спочатку h = 0,2 і треба копіювати у стовпці А до значення 1, тобто до чарунки А7. В результаті отримаємо таку таблицю:
У стовпці С до отриманої таблиці додамо відповідні коефіцієнти Котеса узагальненої формули трапецій:
Аналогічно дістанемо такі таблиці для h = 0,1 і h = 0,05:
2. На основі цих обчислень побудуємо таблицю кратного перерахунку. Надамо чарункам таких значень:
Тут у стовпці А n – кількість відрізків, на які вузли інтегрування ділять [0;1] (1/n = h), яка подвоюється згідно з методом кратного перерахунку. У стовпці В інтегральна сума, а у стовпці С – значення інтеграла згідно з узагальненою формулою трапецій (10). У стовпці D підрахунок оцінки похибки отриманого значення інтеграла згідно з правилом Рунге (14). Оскільки в узагальненої формули трапецій порядок р дорівнює 2, то тут ділимо на 3 = 2р – 1. В результаті маємо таку таблицю:
Поклавши n = 5, знаходимо з таблиці: In ≈ 0,549344, I2n ≈ 0,563787, I4n ≈ 0,567381; згідно з правилом Рунге R2n (f) ≈ , отримуємо з стовпця D R2n (f) ≈ 0,004814, R4n (f) ≈ 0,001198. Як бачимо, із зростанням n, тобто зменшенням кроку h = 1/ n R(f) зменшується. Тепер розширимо попередню таблицю направо і проведемо в ній наступні перерахунки. Надамо чарункам таких значень:
Символ * у чарунці цієї таблиці означає, що в ній залишилась та ж сама формула, що була до розширення. Далі у стовпці Е знаходимо уточнені значення інтеграла по формулі екстраполяції за Річардсоном (15), тобто значення In,2n та I2n,4n уже порядку 3. Цим започаткований другий перерахунок, а потому в чарунці F28 обчислюємо оцінку похибки R2n,4n (f) знову за правилом Рунге (14): із зростанням кратності перерахунку на одиницю порядок р теж зростає на одиницю і отже при другому перерахунку р = 3, 2р – 1 = 7. Нарешті у чарунці G28 обчислюємо наближене значення інтегралу з залишковим членом порядку 4 за формулою Річардсона. В результаті обчислення в Excel дістаємо таблицю:
Як бачимо з чарунці F28, значення R2n,4n (f) насправді значно краще порядку 3: як було показано у доведенні метод гарантує лише не гірші результати. Отже, наближене значення інтеграла в G28 уже навіть порядку 5 (а не 4), як ми розраховували. Задача 2. Обчислити наближене значення функції f (x) = exsinx на відрізку [0;1] за узагальненою формулою Сімпсона з точністю 10-8, для оцінки похибки використавши метод кратного перерахунку. Розв’язання. Оскільки апріорні оцінки тут вимагають значних обчислень і все одно значно завищені, то спочатку візьмемо найменші можливі n = 2, 4, 8, тобто h = 0,5 0,25 0,125 і оцінимо відповідні похибки за методом кратного перерахунку. (Зауважимо, що за узагальненою формулою Сімпсона n обов’язково парне). Надамо чарункам електронної таблиці таких значень:
В результаті, наприклад, при h = 0,125 отримаємо таку таблицю:
У стовпці С до отриманої таблиці додамо відповідні коефіцієнти Котеса узагальненої формули Сімпсона:
Аналогічно з двома іншими таблицями:
На основі цих обчислень побудуємо таблицю кратного перерахунку, ідентичну таблиці задачі 1:
Тут у стовпці С – значення інтеграла згідно з узагальненою формулою Сімпсона (12). Оскільки в узагальненої формули Сімпсона порядок р дорівнює 4, то тут ділимо на 15 = 2р – 1. В результаті маємо таку таблицю:
Поклавши n = 2, знаходимо з таблиці: In ≈ 0,908185, I2n ≈ 0,909254, I4n ≈ 0,909326; згідно з правилом Рунге R2n (f) ≈ , отримуємо з стовпця D R2n (f) ≈ 7,12176E-05, R4n (f) ≈ 4,81561E-06. Як бачимо, отримані наближені значення інтегралу є недостатньо точними. Тому розширимо попередню таблицю направо і проведемо у ній другий перерахунок:
Символ * у чарунці цієї таблиці означає, як і раніше, що в ній залишилась та ж сама формула, що була до розширення. Далі у стовпці Е знаходимо уточнені значення інтеграла по формулі екстраполяції за Річардсоном (15), тобто значення In,2n та I2n,4n наступного порядку, в чарунці F18 оцінку похибки R2n,4n (f) знову за правилом Рунге (14). Отже, дістаємо таблицю:
Отже, R2n,4n (f) ≈ 1,88137E-07. Порядок похибки зріс на одиницю, проте отримані наближені значення інтегралу все ще є недостатньо точними. Порядок зростає як із зростанням кількості перерахунків, так і з зростанням n. Оскільки всі можливості збільшення кількості перерахунків вичерпані при даних n, то треба покласти n = 16 і провести відповідні додаткові обчислення. При h = 1/16 = 0,0625 отримуємо:
Додамо отримані дані у попередню таблицю кратного перерахунку. Маємо:
Отже, нарешті ми отримали оцінку похибки належного порядку у чарунці F19, задача розв’язана. Це оцінка похибки R4n,8n (f) наближеного значення інтегралу І4n,8n (f), що знаходиться у чарунці Е19: І4n,8n (f) ≈ 0,909330672 (таке значення було отримане після розширення стовпця цієї чарунки в Excel). Насправді отриманий порядок знову більший на одиницю гарантованого відповідною теоремою. Ми можемо тепер провести ще третій перерахунок і подивитись на порядок третього уточненого значення:
Тут р = 6, 2р – 1 = 63. В результаті дістаємо:
Порядок похибки третього уточненого значення не зріс: насправді можна довести, що при зростанні порядку у деякому перерахунку вище гарантованого у наступному перерахунку зростання, як правило, не відбудеться. Питання, тести
1. Однобічні формули чисельного диференціювання – це
2. Нехай значення скінчених різниць Δkу0 знаходяться у діапазоні А1:F1 (Δ1у0 у чарунці В1, Δ2у0 у чарунці С1, …), крок інтерполяції h = 0,1. Тоді наближене значення похідної f ′(х0) згідно з однобічними формулами чисельного диференціювання f ′(х0) ≈ L′n(х0) = (Δу0 – Δ2у0 + Δ3у0 + … + Δnу0) обчислюється такою таблицею: А:
Б:
В:
Г:
3. Нехай значення скінчених різниць Δkу0 знаходяться у діапазоні А1:F1 (Δ1у0 у чарунці В1, Δ2у0 у чарунці С1, …), крок інтерполяції h = 0,1. Наближене значення похідної f ′(х0) обчислене згідно з однобічними формулами чисельного диференціювання f ′(х0) ≈ L′n(х0) = (Δу0 – Δ2у0 + Δ3у0 + … + Δnу0) при n = 4 такою таблицею:
знаходиться у чарунці
4. Нехай значення скінчених різниць Δkу0 знаходяться у діапазоні А1:F1 (Δ1у0 у чарунці В1, Δ2у0 у чарунці С1, …), крок інтерполяції h = 0,1. Наближене значення похідної f ′(х0) обчислене згідно з однобічними формулами чисельного диференціювання f ′(х0) ≈ L′n(х0) = (Δу0 – Δ2у0 + Δ3у0 + … + Δnу0) при n = 4 такою таблицею:
В результаті дістали значення:
Тоді шукане значення дорівнює
5. При застосуванні однобічних формул чисельного диференціювання похибка із зменшенням кроку методу h
6. Нехай функція f (x) лінійна на відрізку [a; b]. Тоді похибка квадратурної формули буде найменшою можливою, якщо для наближеного обчислення застосувати
7. Нехай функція f (x) лінійна на відрізку [a; b]. Тоді найменша можлива похибка квадратурної формули дорівнює
8. Формулою Ньютона – Котеса називають квадратурну формулу, якщо
9. Серед наступних квадратурних формул найкраща (тобто найменша) оцінка для залишкового члена у
10. З алишковий член R(f) квадратурної формули Ньютона – Котеса має найкращий (тобто найбільший) порядок відносно кроку інтегрування для
11. Серед наступних узагальнених квадратурних формул при достатньо малих кроках інтегрування найкраща (тобто найменша) оцінка для залишкового члена у
12. З алишковий член квадратурної формули Ньютона – Котеса має найбільший порядок відносно кроку інтегрування для
13. Отримані оцінки наступного методу є апостеріорними
14 Метод подвійного перерахунку ґрунтується на таких формулах:
15. Формула правила Рунге – це
16. Формула екстраполяції за Річардсоном – це
18. Точний варіант формули екстраполяції за Річардсоном – це
19. У стовпці С наступної таблиці
знаходяться коефіцієнти Котеса
20. Коефіцієнти Котеса узагальненої формули Сімпсона знаходяться у стовпці С наступної таблиці: А:
Б:
В:
Г:
21. У стовпці В наступної таблиці знаходяться наближені значення інтеграла функції, обчислені за узагальненою формулою трапецій з кроками h = 0,2 0,1 0,05 у чарунках В2, В3 і В4 відповідно. Така таблиця здійснює подвійний перерахунок цих даних:
A:
Б:
В:
Г:
21. У стовпці В наступної таблиці знаходяться наближені значення інтеграла функції, обчислені за узагальненою формулою трапецій з кроками h = 0,2 0,1 0,05 у чарунках В2, В3 і В4 відповідно. Ця таблиця здійснює подвійний перерахунок цих даних:
В результаті обчислень отримали:
Оцінка похибки для наближеного значення 0,563787 знаходиться у чарунці
Оцінка похибки для наближеного значення 0,568578 знаходиться у чарунці
22. У стовпці В наступної таблиці знаходяться наближені значення інтеграла функції, обчислені за узагальненою формулою Сімпсона з кроками h = 0,5 0,25 0,125 у чарунках В2, В3 і В4 відповідно. Така таблиця здійснює кратний перерахунок цих даних:
A:
Б:
В:
Г:
Завдання 1. У точках 1) х = 1 та 2) х = 1,1 знайти похідну від функції |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Поделиться: |
Познавательные статьи:
Последнее изменение этой страницы: 2016-04-23; просмотров: 370; Нарушение авторского права страницы; Мы поможем в написании вашей работы!
infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 3.129.73.6 (0.012 с.)